[R-es] Duda glmer

[Guillermo.Vinue en uv.es]

Hola,

Tengo una base de datos con estructura jerárquica en la que quiero
clasificar observaciones en distintas categorías.

En el caso más simple, tengo una variable con dos categorías (variable
Y1) y dentro de cada una de ellas hay otras dos categorías (variable
Y2). Además tengo una variable explicativa cuantitativa discreta X.

El banco de datos sería de este tipo:

X Y1 Y2
5  0  1
9  0  0
2  0  1
8  0  0
9  0  0
6  0  0
4  1  0
5  1  1
5  1  1
6  1  1
0  1  1
6  1  0

El enfoque estadístico que quiero plantear es un modelo lineal
generalizado mixto con efectos aletorios para tener en cuenta esta
estructura jerárquica o multinivel. 

Para analizar si este enfoque es adecuado he simulado unos datos. En
primer lugar, trato de clasificar las observaciones en las categorías de
Y1 mediante un glm y obtengo estimadores de los coeficientes que son muy
similares a los valores simulados, como era esperable.

Sin embargo, cuando trato de utilizar un glmer con un intercepto
aleatorio y pendientes fijas para tratar de distinguir entre categorías,
obtengo estimadores muy alejados de los valores simulados. 

He leído mucha de la bibliografía específica (por ejemplo, los capítulos
relacionados con los glmer de los libros de Bates o Goldstein) y he
buscado por internet pero no logro encontrar el error en mi procedimiento.

Adjunto el código de la simulación que llevo a cabo en R por si alguien
me pudiera hacer alguna indicación.

Muchas gracias de antemano.

Un saludo,

Guillermo
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: simulation_glmer.R
Type: application/octet-stream
Size: 2498 bytes
Desc: 
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20150610/4f2e37e6/attachment.obj>