等待一周后运行三个MC链后,对象未保存,给出错误“错误:is.atomic(x)不为TRUE” . 我已多次运行此代码或其变体,因此它可能不是我的语法,但可能是关于对象大小或内部stan_glmer错误?日志在下面,感谢您的帮助(关于如何保存中间链的建议也很有用!)

stan.T.mod <- stan_glmer(form2, iter=5000, warmup=5000, thin=1, 
  family = binomial(link = "logit"), prior = student_t(location=0, df = 1, scale = 2.5), 
  prior_intercept = student_t(location=0, df = 1, scale=2.5),
  prior_covariance = decov(shape = 10, scale = 1), 
  chains = 3, seed = 0305991, data=Tbinary.dt)

模型'bernoulli'NOW(CHAIN 1)的采样 . 链1,迭代:1/5000 [0%](热身)链1,迭代:500/5000 [10%](热身)链1,迭代:1000/5000 [20%](热身)链1,迭代: 1500/5000 [30%](热身)链1,迭代次数:2000/5000 [40%](热身)链1,迭代次数:2500/5000 [50%](预热)链1,迭代次数:3000/5000 [60 %](热身)链1,迭代:3500/5000 [70%](热身)链1,迭代次数:4000/5000 [80%](热身)链1,迭代次数:4500/5000 [90%](预热)链1,迭代次数:5000/5000 [100%](预热)经过时间:394305秒(预热)0秒(采样)394305秒(总计)模型'bernoulli'现在采样(链2) . 链2,迭代:1/5000 [0%](预热)链2,迭代:500/5000 [10%](预热)链2,迭代:1000/5000 [20%](预热)链2,迭代: 1500/5000 [30%](热身)链2,迭代次数:2000/5000 [40%](热身)链2,迭代次数:2500/5000 [50%](预热)链2,迭代次数:3000/5000 [60] %](热身)链2,迭代次数:3500/5000 [70%](热身)链2,迭代次数:4000/5000 [80%](热身)链2,迭代次数:4500/5000 [90%](预热)链2,迭代次数:5000/5000 [100%](预热)经过时间:434251秒(预热)0秒(采样)434251秒(总计)模型'bernoulli'现在采样(链3) . 链3,迭代:1/5000 [0%](热身)链3,迭代:500/5000 [10%](热身)链3,迭代:1000/5000 [20%](热身)链3,迭代: 1500/5000 [30%](热身)链3,迭代次数:2000/5000 [40%](热身)链3,迭代次数:2500/5000 [50%](预热)链3,迭代次数:3000/5000 [60] %](热身)链3,迭代次数:3500/5000 [70%](热身)链3,迭代次数:4000/5000 [80%](热身)链3,迭代次数:4500/5000 [90%](预热)链3,迭代次数:5000/5000 [100%](预热)经过时间:345070秒(预热)0秒(采样)345070秒(总计)错误:is.atomic(x)不是TRUE

save(stan.T.mod, file=paste0(outpath, "stan.GCMfixedmods"));

保存错误(stan.T.mod,file = paste0(outpath,“stan.GCMfixedmods”)):找不到对象'stan.T.mod'stan.T.mod错误:对象'stan.T.mod'未找到