Valores ajustados del resultado de la función "lme

Question

He ajustado un modelo lineal mixto con R de la siguiente manera.

lmme3cs = lme(LAZ ~ group_k + x1 + x2 + x4 + x6 + x1_k + x2_k + x4_k + x6_k,
              random = ~1|SECTORID/CHILDUID,
              correlation = corCompSymm(form = ~1|SECTORID/CHILDUID), data=dat2c)

Tiene 50043 observaciones, pero el objeto 'lmme3cs$fitted' es una matriz de 50043 por 3 que tiene tres columnas fijas, SECTORID y CHILDUID. ¿Qué son estas tres columnas?

Answer 1

Son los valores ajustados a diferentes niveles de agrupación. Por ejemplo, si se ajusta

> library(nlme)
> fmOxide <- lme(Thickness ~ 1, Oxide, ~1|Lot/Wafer)

se puede obtener el resultado global previsto mediante

> fitted(fmOxide, level = 0)

y es 2000.153. Los resultados previstos para cada lote son:

> fitted(fmOxide, level = 1)

y hay ocho salidas como hay ocho Lotes. Los resultados previstos para cada Oblea, que están anidados dentro de los Lotes, son:

> fitted(fmOxide, level = 2)

Si se ejecuta

> fitted(fmOxide, level = 0:2)

se obtienen los resultados previstos para cada nivel de agrupación, es decir, fmOxide$fixed.