各位大神,我为了读取产生数据的models写了一段代码,但是单独拿出来运行没有问题,整体运行系统只运行部分,也并没有报错,想请教是怎么回事?
outputDir <- "D:/data sim/output/"
j <- 1
repNum <- 1
Res_All <- matrix(0, 2, 4)
RepAna <- function(i){
while (j <= repNum){
file <- paste("D:/data sim/output/CFA_", i, "_", j, ".out", sep = "")
ModRes <- readModels(file)
cfi <- 0
tli <- 0
rmsea <-0
a <- 0 #a is the number of fully proper model
CR <- 0
if (is.null(ModRes$warnings[1][[1]]) & is.null(ModRes$errors[1][[1]])){
#exclude non-converge results
a <- a+1
if (ModRes$summaries$CFI > 0.95){
cfi <- cfi+1
}
if (ModRes$summaries$TLI > 0.95){
tli <- tli+1
}
if (ModRes$summaries$RMSEA_Estimate < 0.08){
rmsea <- rmsea+1
}
MRs <- ModRes$parameters$unstandardized
MA <- MRs[,c(2,3)]
MA1 <- subset(MA, param == "CR")
CR <- MA1$est + CR
}
j <- j+1
}
CFI <- cfi/repNum
TLI <- tli/repNum
RMSEA <- rmsea/repNum
CR <- CR/a
Res_All[i,] <- cbind(CFI, TLI, RMSEA, CR)
}
for (i in 1:2){
RepAna(i)
}
系统运行只显示
Reading model: D:/data sim/output/CFA_1_1.out
Reading model: D:/data sim/output/CFA_2_1.out
但是ModRes却没有显示,如果一句一句的运行又没有问题。请问各位大神问题到底在哪里呀?