[R-es] Estudio posthoc de la interacción

Lun Feb 19 12:36:20 CET 2018

Estimado Javier Marcuzzi,

Antes que nada muchas gracias por sus sugerencias.

Le comento que soy estudiante de posgrado en Mejoramiento animal y el
ejemplo que planteo es simplemente un análisis apenas preliminar de efectos
fijos, que busca estimar la influencia de recuento de células somáticas
(entre otras variables) sobre producción de leche.

Entiendo que el efecto rodeo-año-estación está de alguna manera contemplado
en el ejemplo, con el tambo, año de parto y mes de parto, aunque por error
no los especifiqué correctamente en el modelo.

Además, para ser más directo en el ejemplo, solo mencioné que el efecto de
la interacción dio significativo, pero también lo fue para el resto de las
variables explicativas consideradas.

Con el afán de plantear un ejemplo sencillo del problema que quería
analizar, éste resultó sin dudas no del todo correcto y complejo de
interpretar.

Volviendo a mi duda inicial y para abstraerme del ejemplo anterior,
supongamos que tengo el siguiente modelo de efectos fijos  y que todos los
efectos fueron significativos:

y = a +  b + c*d

La duda es, ¿cómo analizar dentro de la interacción c*:*d, la combinación
de niveles  que fue diferente estadísticamente para la variable 'y.'?

El planteo sería:

modelo<-Anova(lm(y~a+b+c*d,data=a),type="III")

Para analizar la interacción, utilicé el paquete phia:

interacc_medias<-interactionMeans(modelo)

Pero me produjo el siguiente error:

*Error in terms.default(model) : no terms component nor attribute*

Me podrían ayudar a estudiar la interacción posthoc con este paquete u otro?

De nuevo muchas gracias por su ayuda y disculpas si el ejemplo original no
fue del todo correcto.

Muchas gracias.

Rodrigo.

<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail&utm_term=icon>
Virus-free.
www.avast.com
<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail&utm_term=link>
<#m_230063116732962769_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>

El 18 de febrero de 2018, 12:56, Javier Marcuzzi <
javier.ruben.marcuzzi en gmail.com> escribió:

> Estimado Rodrigo López Correa
>
> Yo realizo ese tipo de análisis en mi vida profesional.
>
> ¿Qué título tiene o que estudia?
>
> Lo que usted escribe tiene sentido pero biológicamente está con errores,
> en su modelo, la producción de leche depende del tambo, mes de parto, año
> de parto, lactancia, recuento de células somáticas y relación entre
> lactancia y células somáticas. Luego dice que tiene relación significativa
> entre la lactancia y células somáticas.
>
> Biológicamente las células somáticas afectan la producción de leche, la
> producción de leche afecta las lactancias, entrando muy fino puede haber
> diferencias en la histología, pero un cambio por ejemplo ambiental afecta
> el recuento de células somáticas. Si mira algún artículo del Journal Dairy
> Science su modelo no tiene algo típico como leche = hys + ..., una simple
> pregunta como si al envejecer la vaca aumentan las células somáticas, es lo
> único a lo que apunta claramente su modelo a responder desde mi punto de
> vista, el resto de su modelo tiene problemas biológicos, le recomiendo
> definir un poco más ese aspecto de acuerdo a lo que buscas, R le puede
> procesar información y dar resultados, pero las vacas se manejan con otras
> reglas.
>
> Usted tiene un error anterior, el análisis de las variancias puede ser
> bueno, pero usted no está reflejando la producción de leche, no le falta
> mucho para escribir el modelo en una forma adecuada, luego lo puede
> complicar todo lo que desea, pero le recomiendo ver ese aspecto antes de
> continuar.
>
> Javier Rubén Marcuzzi
>
>
>
> El 18 de febrero de 2018, 12:13, Rodrigo López Correa <rod99hare en gmail.com
> > escribió:
>
>> Estimados,
>>
>> Estoy analizando el siguiente modelo para producción de leche y he
>> utilizado el paquete car para plantear el siguiente modelo de efectos
>> fijos:
>>
>>
>> modelo<-Anova(lm(leche~tambo_ypart+mpart+lact+rcs_categ3+lact*:*
>> rcs_categ3,data=a),type="III")
>>
>> El resultado dio siginificativa la interacción entre lact y  rcs_categ3.
>>
>> Para analizar la interacción, utilicé el paquete phia:
>>
>> interacc_medias<-interactionMeans(modelo)
>>
>> Pero me produce el siguiente error:
>>
>> *Error in terms.default(model) : no terms component nor attribute*
>>
>>
>> Me podrían ayudar a estudiar la interacción posthoc con este paquete u
>> otro?
>>
>>
>> Muchas gracias!
>>
>> Rodrigo.
>>
>>
>> <https://www.avast.com/sig-email?utm_medium=email&utm_source
>> =link&utm_campaign=sig-email&utm_content=webmail&utm_term=icon>
>> Virus-free.
>> www.avast.com
>> <https://www.avast.com/sig-email?utm_medium=email&utm_source
>> =link&utm_campaign=sig-email&utm_content=webmail&utm_term=link>
>> <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
>>
>>         [[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
>>
>
>

-- 
*Dr. Rodrigo López Correa.*

Montevideo.
Uruguay.
Cel: 099 660 549.

-- 
*Dr. Rodrigo López Correa.*

Montevideo.
Uruguay.
Cel: 099 660 549.

	[[alternative HTML version deleted]]