¿Cómo ajustar el modelo para diseño cruzado y anidado usando la función lme en R?

8

Supongamos que A es un factor fijo y B es un factor aleatorio. A y B están cruzados. Otro factor aleatorio C está anidado en B. Teniendo en cuenta las interacciones AB y AC, ¿cómo debo ajustar el modelo usando la función lme en el paquete nlme? Gracias.

r mixed-model lme4-nlme Knightgu
fuente
Usted menciona el paquete lme en el título y nlme en su pregunta. ¿Cuál es el paquete correcto?
João Daniel
1
Y .... ¿dónde está el código?
DWin
Dije la función lme en el paquete nlme. lme no es un paquete. Y me pregunto cómo especificar los argumentos en esta función para mi problema.
¿Has probado lmer()en el paquete lme4? Ver, por ejemplo, este capítulo del libro
invitado
Muchas gracias. Sí, la función lmer () puede manejar fácilmente efectos aleatorios cruzados. Es por eso que el paquete lme4 es una gran mejora para el paquete nlme. Más tarde me di cuenta de que era muy difícil usar la función lme () para ajustar los efectos aleatorios cruzados. Ya que tenemos el paquete más nuevo, ¿por qué no usarlo?？
Knightgu

Respuestas:

5

La parte difícil aquí es especificar la fórmula correcta.

library(nlme)
# Assume y is your response variable.
lme(fixed = y ~ A, random = ~ A | B/C)

Aquí B/Csignifica C anidado en B y es la abreviatura de B + B:C. La interacción con A está presente porque A se repite en la parte aleatoria.

Para una introducción básica a nlme, recomiendo este documento .

gui11aume
fuente
Muchas gracias por tu solución. Creo que tienes razón. Usaré un ejemplo práctico para verificarlo. Por cierto, muy buena referencia. Gracias.
Knightgu
1
He probado la fórmula que mencionaste usando R. Desafortunadamente, el conjunto de datos que usé no justificó la fórmula. La función lme () aplicada simplemente no convergió a una solución numérica. Parece que la función lmer () en el paquete lme4 hace lo correcto.
Knightgu