[R-es] Efectos fijos y aleatorios en un modelo lineal

[Fernando Piston]

Muchas Gracias Luciano por tu pronta respuesta.

Como bien decis cuesta imaginar como se puede considerar un factor como 
aleatorio cuando se tienen solo 2 niveles en el caso de Year y 3 en el 
caso de Rep. Sin embargo, aunque seria deseable tener mas niveles para 
el factor Year (más años de ensayo) su efecto es, en cualquier caso, 
aleatorio.

Saludos
Fernando

On 19/03/12 15:22, Luciano Selzer wrote:
> Fernando:
> Más allá de la solución que alguien de la lista te pueda dar me
> gustaría darte una advertencia. Según varios autores, los efectos
> aletarorios se pueden determinar correctamente a partir de 6 niveles.
> Así que con 3 y 2 niveles es muy dificil que se pueda determinar
> correctamente la varianza de cada uno.
>
> Saludos
> Luciano
>
>
>
> El día 19 de marzo de 2012 10:21, Fernando Piston
> <fpiston en ias.csic.es>  escribió:
>> Hola a todos,
>>
>> Tengo algunas dudas sobre como introducir en un modelo lineal factores con
>> efectos fijos y aleatorios. Mi diseño es el siguiente:
>>
>> Factores
>>
>> Line: 40 líneas de trigo (Fijo)
>>
>> Rep: 3 Bloques (Aleatorio)
>>
>> Year: 2 Años (Aleatorio)
>>
>>
>> Variable dependiente
>>
>> alpha.ug.mg
>>
>>
>> Nota: Adjunto tabla de datos
>>
>>
>> Dentro del diseño, 'Rep' y 'Year' se considera de efectos aleatorios y
>> 'Line' es de efecto fijo.
>>
>>
>> Lo he intento llevar a cabo con la función 'lmer' del paquete 'lme4',
>>
>> lmer(alpha.ug.mg ~ Line + (1|Year) + (1|Rep%in%Year) + (1|Line:Year), data =
>> Dataset)
>>
>>
>> pero el modelo anterior me considera como error del factor 'line' el error
>> de los residuales en lugar del error de Year*Line.
>>
>> ¿como puede decirle al modelo que para 'line' me considere el error de
>> 'Year*Line'?
>>
>>
>> Gracias
>>
>>
>>
>> _______________________________________________
>> R-help-es mailing list
>> R-help-es en r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-help-es
>>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es