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 <[email protected]> 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 > [email protected] > https://stat.ethz.ch/mailman/listinfo/r-help-es > [[alternative HTML version deleted]] _______________________________________________ R-help-es mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-help-es
