R：“for”循环，混合效果模型崩溃

Question

我在分析期间R崩溃时出现“for”循环功能有问题。我使用“for”循环随机地对较大的数据集进行子样本处理（维护嵌套设计），将嵌套的混合模型应用于子样本，并提取方差分量。该循环旨在运行999次迭代，但R将在循环完成之前崩溃。我从Windows获得此错误消息“R for Windows GUI前端遇到问题，需要关闭”。我已经尝试使用32位和64位Windows系统的多台计算机，所以我开始相信它可能是我编写的代码。

R=999
var10<-numeric(R)
for (i in 1:R){ 
not10=subset(not, not$pat==sample(not$pat,replace=TRUE,10))
mod10=lme(fixed=sl~1,random=~1|pat/mat,data=not10,method="REML")
VarCorr(mod10)
sire10 <- as.numeric(VarCorr(mod10)[2,1])    
sire10
dam10 <- as.numeric(VarCorr(mod10)[4,1])      
dam10
wg10 <- as.numeric(VarCorr(mod10)[5,1])     
wg10
var10[i]<- sire10/(dam10+sire10+wg10)
}

有没有人对修复此问题有任何见解？有没有更好的方法来编码我的分析以获得更有效的处理？

谢谢！