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
>
[[alternative HTML version deleted]]