in0 <- function(file){ estimates0 <- as.data.frame(matrix(scan(file=file), ncol=66, byrow=TRUE)) names(estimates0) <- c("rep",paste("ex",1:4,".mu", sep=""), paste("in",1:4,".mu", sep=""), paste("b.ex",2:4, sep=""), paste("b.in",2:4, sep=""), paste("ex",1:4,".res", sep=""), paste("in",1:4,".res", sep=""), "y.mu","b.f1","b.f2","f1.var","f1f2.cov","f2.var","y.var", paste("ex",1:4,".mu.se", sep=""), paste("in",1:4,".mu.se", sep=""), paste("b.ex",2:4,".se", sep=""), paste("b.in",2:4,".se", sep=""), paste("ex",1:4,".res.se", sep=""), paste("in",1:4,".res.se", sep=""), "y.mu.se","b.f1.se","b.f2.se","f1.var.se","f1f2.cov.se","f2.var.se","y.var.se", "H0LogLik", "df", "AIC", "BIC", "AdjBIC", "Entropy", "Condition") estimates0 } in1 <- function(file){ estimates1 <- as.data.frame(matrix(scan(file=file), ncol=68, byrow=TRUE)) names(estimates1) <- c("rep",paste("ex",1:4,".mu", sep=""), paste("in",1:4,".mu", sep=""), paste("b.ex",2:4, sep=""), paste("b.in",2:4, sep=""), paste("ex",1:4,".res", sep=""), paste("in",1:4,".res", sep=""), "y.mu","b.f1","b.f2","b.s","f1.var","f1f2.cov","f2.var","y.var", paste("ex",1:4,".mu.se", sep=""), paste("in",1:4,".mu.se", sep=""), paste("b.ex",2:4,".se", sep=""), paste("b.in",2:4,".se", sep=""), paste("ex",1:4,".res.se", sep=""), paste("in",1:4,".res.se", sep=""), "y.mu.se","b.f1.se","b.f2.se","b.s.se","f1.var.se","f1f2.cov.se","f2.var.se","y.var.se", "H0LogLik", "df", "AIC", "BIC", "AdjBIC", "Entropy", "Condition") estimates1 }