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
Post a Comment