## model 2
res<-scan("1g-2-500.txt")
res.mat<-matrix(res, byrow=T, ncol=66)
res.f<-res.mat[1:1000,26:28]
 
cbind(apply(res.f,2,mean),apply(res.f,2,var))
 
## for model 3 to 6
res<-scan("1g-5-500.txt")
 
res.mat<-matrix(res, byrow=T, ncol=96)
 
res.f<-res.mat[1:1000,41:43]
 
cbind(apply(res.f,2,mean),apply(res.f,2,var))
 
 
## for model 7
res<-scan("1g-7-500.txt")
res.mat<-matrix(res, byrow=T, ncol=126)
res.f<-res.mat[1:1000,56:58]
 
cbind(apply(res.f,2,mean),apply(res.f,2,var))