Hola.

Leyendo tu primera pregunta, yo diría que los valores predichos para el 
nivel superior si son las medias que obtienes cuando evaluas los valores 
predichos del nivel inferior en las categorías del efecto aleatorio.

Con el ejemplo que viene en  lmer sería

library(lme4)
(fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy))

predichos <- tapply(fitted(fm1), sleepstudy$Subject, mean)
observados <- tapply(sleepstudy$Reaction,sleepstudy$Subject,mean)
residuos <- observados-predichos

# predichos vs residuos^2
plot(predichos,residuos^2)

# observados vs predichos
plot(observados,predichos)
abline(a=0,b=1,col="red")

Saludos



El 09/05/13 22:09, Joan Molibo escribió:
> Muchas gracias, he encontrado un curso online gratis que tiene buena pinta.
>
>   
>
> http://www.bristol.ac.uk/cmm/research/lemma/
>
>   
>
> No obstante sigo con mi duda después de haber estudiado más sobre el tema;
> cuando termine si continua la reformularé.
>
>   
>
> Gracias.
>
>   
>
> Joan.
>
>   
>
> De: marticasals@gmail.com [mailto:marticasals@gmail.com] En nombre de Martí
> Casals
> Enviado el: Miércoles, 08 de Mayo de 2013 12:54
> Para: JOANMOLIBO@terra.es
> CC: r-help-es@r-project.org
> Asunto: Re: [R-es] lme4 y residuales
>
>   
>
> Hola Joan,
>
> te sugiero mirar este link, creo que te puede ir bien.
>
> http://modeloslineales.wikidot.com/
>
>   
>
> Saludos,
>
>   
>
> Martí
>
>   
>
> 2013/5/8 JOANMOLIBO@terra.es <JOANMOLIBO@terra.es>
>
> Buenos días;
>
> Estoy intentando aprender algo de sobre  los modelos lineales mixtos. Estoy
> siguiendo el libro "Multilevlel analysis: techniques and applications" de
> Joop Hox. Estoy intentado reproducir las tablas y gráficos que trae el libro
> con el paquete lmer.
>
> En este libro se indica que en la comprobación de los supuestos del modelo
> hay que estudiar los gráficos de residuales tanto del nivel 1 como de los
> superiores.
>
> En el primer conjunto de datos que se analiza hay dos niveles. Obtener el
> los residuales y los valores ajustados  por al nivel inferior no ha sido
> problema con las funciones residuals y fitted.
> Para el segundo nivel, la función ranef me extrae los residuales. El
> problema lo tengo con los valores predichos,  y no puedo realizar el gráfico
> de los residuales frente a los valores predichos.
>
> Mi pregunta es: ¿cómo puedo realizar tal gráfico para el nivel 2? ¿Cómo
> obtengo los valores predichos?¿Serían estos las medias predichas para cada
> grupo del nivel 2?
>
> Gracias.
>
> Joan.
>
>
>          [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
>
>
>
>
>
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es


	[[alternative HTML version deleted]]

