【发布时间】:2013-06-14 08:37:17
【问题描述】:
我在 R 中运行模拟,在一个 for 循环中。在循环的每一步中,我都在尝试使用 VGAM 包的 vglm 函数来拟合 beta 二项分布:
vglm(data ~ 1, betabinomial)
但是,对于我生成的某些数据,beta 二项式模型将不适合,并且函数会崩溃。因为它崩溃了,所以 R 退出循环并停止。我想简单地转到循环中的下一步。我知道功能 tryCatch 但在我的情况下未能使用它。 目前,我的代码如下(注意for循环嵌套在另一个for循环中):
for (i in 1:nsimu) {
## Some code to generate Pos and Neg
data <- cbind(Pos,Neg)
if (sum(Pos)==0) {
##Go to next step
}
else {
vglm(data ~ 1, betabinomial)
}
}
非常感谢您的任何帮助或建议。
史蒂芬
【问题讨论】: