Buenas tardes


Quiero usar la función *corSymm* de la librería *nlme*  para hacer las
estimaciones de los parámetros  de  un *Linear mixed-effects model. *Esta
metodología es propuesta por  por Pinheiro, J.C., and Bates, D.M. (2000)
"Mixed-Effects Models in S and S-PLUS", Springer.


Quiero obtener una estimación similar a la que genera el siguiente  script:

*library(nlme)*

*spatDat <- data.frame(res=c(2,3,4,2,4), x = c(0,0.25,0.5,0.75,1), y =
c(0,0.25,0.5,0.50,0.75), pob=c(12,23,34,45,32), otro3=c(1,1,1,1,1))*

*cs1Exp <- corExp( 1, form = ~ x + y )*

*cs1Exp <- Initialize( cs1Exp, spatDat )*

*pruebaq <- glmmPQL(res~ x + y + offset(log(pob)), data = spatDat, family =
poisson, random = ~1 |otro3,  correlation = cs1Exp)*

*pruebaq*



*Mi interés recae sobre el siguiente script:
*

*
*

mat<-corMatrix( cs1Exp ) # Estoy partiendo de una matriz de covarianzas
cualquiera. Por comodidad parto de la matriz definida en el anterior
ejemplo.

superior<-t(mat)[lower.tri(t(mat))]

super <- corSymm( value = superior, form = ~1|otro3)

super<- Initialize(super, data = spatDat)

pruebaq2 <- glmmPQL(res~ x + y + offset(log(pob)), data = spatDat, family =
poisson, random = ~1 |otro3,  correlation = super)


*El error que me genera es el siguiente: *

*iteration 1*

*Error en lme.formula(fixed = zz ~ x + y + offset(log(pob)), random = ~1 |  :
*

*  nlminb problem, convergence error code = 1*

*  message = false convergence (8)*

*
*

*Quiero saber cómo corregir este error o algún guía sobre este. *

*
*

*Agradezco su atención *

	[[alternative HTML version deleted]]

