[R] fitting generalized linear models using glmmPQL

[Eduardo.Garcia@uv.es]

Hi, I have tried to run the following (I know it's a huge data set but 
I tried to perform it with a 1 GB RAM computer):

library(foreign)
library(MASS)
library(nlme)
datos<-read.spss(file="c:\\Documents and 
Settings\\Administrador\\Escritorio\\datosfin.sav",to.data.frame=TRUE)
str(datos)

`data.frame':   1414 obs. of  5 variables:
 $ POB     : Factor w/ 6 levels "CHI","HOS","HYR",..: 1 1 1 1 1 1 1 1 
1 1 ...
 $ CLON    : num  1 1 1 1 1 1 1 1 2 2 ...
 $ TEMP    : Factor w/ 2 levels "20 C","25 C": 1 1 1 1 2 2 2 2 1 1 ...
 $ SALINITA: Factor w/ 2 levels "15 g/l","30 g/l": 1 1 2 2 1 1 2 2 1 
1 ...
 $ DE      : num  17 0 7 1 15 28 4 14 13 16 ...
 - attr(*, "variable.labels")= Named chr  "" "" "" "" ...
  ..- attr(*, "names")= chr  "POB" "CLON" "TEMP" "SALINITA" ...

datos$CLON<-as.factor(datos$CLON)

modelo1<-glmmPQL(DE ~ (POB/CLON)*TEMP*SALINITA, data = datos, random = 
~ 1|CLON, family = poisson)

And I have obtained the following:

Error: NA/NaN/Inf in foreign function call (arg 1)
In addition: Warning message:
step size truncated due to divergence 

This is the first time I've observed such a message and I have no idea 
about what does it mean. Is it possible that the process failed 
because of the size of the data set (180 levels of the CLON factor)? 
Or maybe is it a syntax problem?

Thank you in advance.


********************************
Eduardo Moisés García Roger

Institut Cavanilles de Biodiversitat i Biologia
Evolutiva - ICBIBE.
Tel. +34963543664
Fax  +34963543670