1

使用 R 包 nlme 中的函数 gls() 进行模拟时遇到问题。我想对装有 gls() 的线性模型强加一个 AR(1) 相关结构。当我在模拟框架中执行此操作时（例如，将整个事情重复 10,000 次），R 通常（但并非总是）在某个时候突然崩溃。我在 R 中没有收到错误消息，但 R 突然终止并且 Windows 报告“R for Windows GUI 前端已停止工作”。我正在使用 OS Windows 7 和 R 版本 3.0.0 以及最新版本的 nlme (3.1-109)。

这是我的（简化的）代码：

library(nlme)

group <- rep(c("A", "B", "C", "D"), each=5)
person <- as.factor(rep(c(1:5), times=4))
value <- c(8.50, 11.69, 9.07, 8.82, 10.63, 10.87, 10.02, 9.84, 9.82, 9.68,
           7.55, 12.19, 9.16, 8.59, 9.98, 11.16, 9.63, 9.73, 10.15, 9.06)
dat <- data.frame(group, person, value)

xyz <- numeric(1)

runs <- 1E4

for(r in 1:runs){
  model <-  gls(value ~ group - 1, data=dat, corAR1(form=~1|person))
  if(summary(model)$corBeta[1,2] < 0.9){xyz <- xyz + 1}
}

xyz

不同的数据集也会出现问题，而不仅仅是本例中的特定数据集。此外，我在不同（但相同）的工作站和不同的 R 会话中反复运行上述代码。R 在大约 50% 的情况下会崩溃。不幸的是，我不知道是什么导致了这种情况发生（或没有发生），如果你能帮助我，我将不胜感激。

4

0 回答 0