r - glmmPQL crashes on inclusion of corSpatial object -


link data (1170 obs, 9 variables, .rd file)

simply read in using readrds(file).

i´m trying setup glmm using glmmpql function mass package including random effects part , accounting spatial autocorrelation. however, r (version: 3.3.1) crashes upon execution.

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) 

what tried far:

  • reduce number of observation
  • play corspatial parameters (range , nugget)
  • reduce number of fixed predictors
  • execute code on windows, linux (debian) , mac r installations

while no error message on local pc (rstudio crashes), running script on server returns following error message:

r: malloc.c:3540: _int_malloc: assertion (fwd->size & 0x4) == 0' failed. aborted

i'd use inla package model this, allows use spatially correlated random effects. required code bit long place here. therefore i've place in document on http://rpubs.com/inbostats/spde


Comments