Estoy realizando el análisis de un experimento time course con dos factores y una variable de respuesta.
Mi tabla tiene esta estructura esquemática
Individuos | Tiempo |Sexo | Gen | Respuesta
1,2,3,4...44 | 0,5,10...55 | M-F | WT-Mut | Variable numérica a analizar
Cada uno de los 44 individuos se mide en cada punto temporal y se clasifica según el Sexo y el Gen, que son variables dicotómicas. Adjunto datos de ejemplo:
Individuals <- rep(c(1:44), each=12)
Time <- rep(c(0,5,10,15,20,25,30,35,40,45,50,55), 44)
Sex <- rep(c("M","F"), each=2, length.out=528)
Gen <- rep(c("WT","Mut"), length.out=528)
Response <- sample(c(0:50), 528)
WORK <- data.frame(Individuals, Time, Sex, Gen, Response)
Quiero analizar si hay diferencias estadísticamente significativas debidas al sexo y al gen en cada momento. Estoy utilizando las funciones shapiro_test, levene_test, anova_test, pairwise_t_test y cohens_d del paquete rstatix.
Creo que mi problema está en cómo interpretan estas funciones la variable "Tiempo", porque tengo otras pruebas con la misma estructura (pero sin tiempo) y funciona bien. Ya he probado a convertir la variable Time en factor, pero no funciona.
He probado muchas cosas diferentes y con todas obtengo un error. Aquí tenéis mi última y sencilla prueba:
TRABAJO %>% group_by(Tiempo, Sexo, Gen) %>% shapiro_test(Respuesta)
WORK %>% levene_test(Response ~ Time*Sex*Gen, center = 'mean')
WORK %>% anova_test(Respuesta ~ Tiempo*Sexo*Gen)
WORK %>% group_by(Time,Sex) %>% pairwise_t_test(Response~Gen, p.adjust.method = "holm") #En realidad este es el único que funciona
WORK %>% group_by(Time, Sex) %>% cohens_d(Response~Gen, hedges.correction = TRUE)
Por favor, ¿alguien puede aportar alguna idea de cómo resolverlo de forma sencilla?
¡Gracias de antemano!