Javier,
creo entender que en tu ejemplo, los factores Casa y Auto son fijos, no?
Si es así, la salida del modelo lmer debería darte los errores
asociados a estos factores.
En todo caso, para ayudarte sería conveniente que nos des tu modelo.
Por ejemplo algo del tipo:
fit=lmer(Impuestos ~ Casa + Auto + (1|hogar), data=datos))
Un saludo. Olivier
--
____________________________________
Olivier G. Nuñez
Email: onunez en iberstat.es
Tel : +34 663 03 69 09
Web: http://www.iberstat.es
____________________________________
El 24/02/2010, a las 12:09, Javier Marcuzzi escribió:
> Tengo una duda, hay unos modelos mixtos que se resuelven con lme4,
> no tengo
> problemas con esto, pero necesito tener el resultado para cada uno
> de los
> factores, lo expreso con este ejemplo aunque no es código r creo
> que se
> entendería
>
>
>
> Factores como ser casa y auto, tengo casa 1, casa 2, casa 3, y auto 1,
> auto2, auto 3.
>
>
>
> Con lme4 y utilizando ranef obtengo por ejemplo
>
>
>
> Casa 1 2
>
> Casa 2 2,1
>
> Casa 3 2
>
> Auto 1 0,2
>
> Auto 2 0,18
>
> Auto 3 0,19
>
>
>
> Supongamos que el valor de Casa 1 es el resultado de información de
> impuestos abonados durante 2 años.
>
> De la misma forma suponemos que el valor de Casa 2 es el resultado de
> información sobre impuestos abonados durante 1 año, porque antes la
> casa no
> existía, fue construida hace poco.
>
>
>
> Por lo cual la exactitud del cálculo es mayor cuándo se tienen más
> registros, en este caso Casa 1 tiene un valor más confiable que
> Casa 2.
>
>
>
> Lo que necesito es tener como resultado final es por ejemplo
>
>
>
> Casa 1 2 1
>
> Casa 2 2,1 0,5
>
> Casa 3 2 0,5
>
> Auto 1 0,2 1
>
> Auto 2 0,18 0,2
>
> Auto 3 0,19 0,4
>
>
>
> Donde la tercer columna es la exactitud.
>
>
>
> ¿Hay alguna forma de extraer esta utilizando lme4, como por ejemplo
> ranef
> para los otros valores?
>
> El término específico es accuracy o reliability, pero no enconté
> esto dentro
> de la documentación lme4.
>
>
>
>
>
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es