Сбой glmmPQL при включении объекта corSpatial

Question

Ссылка на данные (1170 объектов, 9 переменных, файл.Rd)

Просто прочитайте это, используя readRDS(file),

Я пытаюсь настроить GLMM, используя glmmPQL функция от MASS пакет, включающий часть случайных эффектов и учет пространственной автокорреляции. Тем не менее, 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