¿Puedo "descartar" los valores p de mi modelo lineal mixto?

Question

Tengo un modelo de efectos mixtos lineales como el siguiente en R:

modelo.int = lmer(valor ~ grupo * categoría * tiempo + (1|paciente), data=modeldata)

anova(modelo.int) luego me da la siguiente lista de valores p:

1. grupo
2. categoría
3. tiempo
4. grupo * categoría
5. grupo * tiempo
6. categoría * tiempo
7. grupo * categoría * tiempo

Dado que utilizo el modelo mixto en 5 conjuntos de datos, esto resulta en un total de 5*7 = 35 valores p. Ahora sabemos que un mayor número de valores p resulta en un mayor porcentaje de errores de tipo I, por lo que deberíamos corregir para comparaciones múltiples.

Sin embargo, entre los 7 valores p, los únicos que conciernen a mi hipótesis son los dos marcados en negrita. Los otros valores p no me interesan. Por ejemplo, ya sé que tiempo y grupo tienen un efecto significativo en valor por diseño del estudio, ¡esa es específicamente la razón por la cual los incluí como efectos fijos en el modelo de efectos mixtos lineales!

Mi pregunta es, ¿puedo simplemente reportar los 2 valores p concernientes a mi hipótesis en mi artículo y descartar el resto? Y en caso afirmativo, ¿debo corregir para comparaciones múltiples solo entre los 2*5 = 10 valores p concernientes a mi hipótesis o para los 35?

Answer 1

Mi pregunta es, ¿puedo simplemente reportar los 2 valores p relacionados con mi hipótesis en mi documento y descartar el resto? Y si es así, ¿debo corregir para comparaciones múltiples solo entre los 2*5 = 10 valores p relacionados con mi hipótesis o para los 35 en total?

Ten cuidado. Podría dar la impresión de que no incluiste los otros términos en tus modelos, y hacer eso es un gran error. Entiendo que puede haber mucha salida que no sea relevante para tu punto, por lo que tal vez incluye solo los términos de interés en el cuerpo de tu documento, pero asegúrate de incluir una nota al pie que refiera a un apéndice donde muestres la salida completa.