
## handle the model comparison results from mplus

v500<-scan('all_results500v.txt')

v500m<-matrix(v500, nrow=1500, byrow=T)

inv500<-scan('all_results500inv.txt')

inv500m<-matrix(inv500, nrow=1497, byrow=T)

v500f<-v500m[c(1:523, 525:1034, 1036:1279, 1281:1500), c(40:42, 169)]
inv500f<-inv500m[, c(40:42, 163)]

mean(v500f[,4])
var(v500f[,4])
ks.test(v500f[,4],'pchisq',23)

mean(inv500f[,4])
var(inv500f[,4])
ks.test(inv500f[,4],'pchisq',27)

chi.diff<-inv500f[,4]-v500f[,4]
mean(chi.diff)
var(chi.diff)

ks.test(chi.diff, 'pchisq', 3)


## for the sample size 100
## function to delete an element
delOne<-function(x, a){
  nn<-NULL
  for (i in 1:length(a)){
    nn<-c(nn, which(x==a[i]))
  }
  x<-x[-nn]
  x
}

ind<-scan("ind.txt")

inv<-scan("all_results100inv.txt")
invm<-matrix(inv, ncol=179, byrow=T)
err.inv<-scan("error-100-inv.log")

ind.inv<-delOne(ind,err.inv)
chi.inv<-cbind(ind.inv,invm[,163])

v<-scan("all_results100v.txt")
vm<-matrix(v, ncol=185, byrow=T)
err.v<-scan("error-100-v.log")
ind.v<-delOne(ind,err.v)

chi.v<-cbind(ind.v, vm[,169])

## create two NA matrix
chi.inv.temp<-array(NA, dim=c(1500,2))
chi.v.temp<-array(NA, dim=c(1500,2))
i<-1
for (i in 1:1500){
  temp<-chi.inv[chi.inv[,1]==i ,]
  if (length(temp)!=0){
    chi.inv.temp[i,]<-temp
  }

  temp<-chi.v[chi.v[,1]==i ,]
  if (length(temp)!=0){
    chi.v.temp[i,]<-temp
  }
}

chi.all<-cbind(chi.inv.temp[,2], chi.v.temp[,2], chi.inv.temp[,2]-chi.v.temp[,2])

chi.all<-chi.all[chi.all[,3]>0, ]

apply(chi.all, 2, mean, na.rm=T)
apply(chi.all, 2, var, na.rm=T)

ks.test(chi.all[,1],'pchisq',27)
ks.test(chi.all[,2],'pchisq',24)
ks.test(chi.all[,3],'pchisq',3)