Hello everybody,

consider this dataset:

*arti <- data.frame(
                    fieldID = as.factor(rep(1:25,each=4)),       # Object
ID
                    wave = as.factor(rep(c(1990,2000,2010,2020),25)), #
Year of survey
                    event = sample(c(0,0,0,0,1),100,r=T), # Event has
happened 1 and or not 0
                    height = rep(runif(25)*10,each=4)+rep(0:3,25), # some
explanatory variable
                    xcor = rep(1:25,each=4), # x-coordinates of object
                    ycor = rep(1:25,each=4) # y-coordinates of object
                  )*

with this model:

library(MASS)

*glmmfit <- glmmPQL (event ~ height, random=~1|fieldID, data=arti,
family="binomial",
correlation=corExp(form=~ycor+xcor|fieldID/wave),niter=50)*

and the resulting error message

*Iteration 1
Error in corFactor.corSpatial(object) :
 NA/NaN/Inf in foreign function call (arg 1)
In addition: Warning messages:
1: In min(unlist(attr(object, "covariate"))) :
 no non-missing arguments to min; returning Inf
2: In min(unlist(attr(object, "covariate"))) :
 no non-missing arguments to min; returning Inf*


I would be very grateful, if anyone could give me some pointers on what the
hell is going on here or what am I doing wrong? I already tried everything
imaginable, at least what is imaginable to me.

Greetings
Gregor
(Mastersstudent Forest Science - Freiburg (Germany))

	[[alternative HTML version deleted]]


