0

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!

1
  • Bienvenido a Stack Overflow en español, te sugiero que hagas el recorrido de bienvenida y de paso ganes tu primer medalla, también es muy importante que leas Cómo preguntar para poder mejorar tu pregunta y que sea bien recibida por la comunidad mejorando así, tus chances de obtener buenas respuestas. el 30 dic. 2022 a las 20:13

0

Tu Respuesta

Al pulsar en “Publica tu respuesta”, muestras tu consentimiento a nuestros términos de servicio, política de privacidad y política de cookies

Examina otras preguntas con la etiqueta o formula tu propia pregunta.