glmmPQL在包含corSpatial对象时崩溃

Question

Link to data（1170个障碍，9个变量，.Rd文件）

只需使用readRDS(file)。

即可阅读

我正在尝试使用glmmPQL包中的MASS函数设置GLMM，包括随机效果部分并考虑空间自相关。但是，R（版本：3.3.1）在执行时崩溃。

library(nlme)

# setup model formula
fo <- hail ~ prec_nov_apr + t_min_nov_apr + srad_nov_apr + age

# setup corSpatial object
correl = corSpatial(value = c(10000, 0.1), form = ~ry + rx, nugget = TRUE,
                    fixed = FALSE, type = "exponential")
correl = Initialize(correl, data = d)

# fit model
fit5 <- glmmPQL(fo, random = ~1 | date, data = d, 
                correl = correl, family = binomial)

到目前为止我尝试了什么：

• 减少观察次数
• 使用corSpatial参数（范围和金块）
• 减少固定预测变量的数量
• 在Windows，Linux（Debian）和Mac R安装上执行代码

虽然我的本地电脑上没有收到任何错误消息（RStudio只是崩溃），但在服务器上运行脚本会返回以下错误消息：

R: malloc.c:3540: _int_malloc: Assertion (fwd->size & 0x4) == 0' failed. Aborted

Answer 1

我使用INLA包对此进行建模，因为它允许使用空间相关的随机效果。所需的代码有点太长了，无法放在这里。因此，我已将其放在http://rpubs.com/INBOstats/spde

上的文档中