> cordat <- readLowTriMat(file = "MatricesReduced.txt", no.var = 15, na.strings="NA") Read 1200 items > N <- c(5689,8597, 24396, 8870, 36435, 34687, 56308, 34299, 53760, 68427) > stage1fixed <- tssem1(Cov = cordat, n = N, method = "FEM") > summary(stage1fixed) > stage1random <- tssem1(Cov = cordat, n = N, method = "REM", + RE.type = "Diag") > summary(stage1random) > firstRun <- reRun(stage1random) > secondRun <- reRun(firstrun)