## 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))