等待一周后,三个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)的采样。
Chain 1,Iteration:1/5000 [0%](Warmup)Chain 1,Iteration: 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'NOW(CHAIN 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'NOW(CHAIN 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'