Escribo de nuevo la pregunta ya que en la anterior no expuse la duda correctamente y ya se resolvió la pregunta. No obstante el código me sigue sin funcionar, es por ello por lo que voy a plantear la duda desde el principio:
Tengo un dataframe de 3825 por 9 columnas; este está depositado en este link. Me gustaría saber la normalidad y la homocedasticidad de cada una de las variables en función de la variable Conut_stage
.
Para ello he ejecutado los siguientes códigos, que no funcionan correctamente:
Hago una lista de cada variable en función de Conut_stage
My_list <- split(checking, f = list(checking$Conut_stage))
Realizo el loop:
loop_Shapiro <- list()
variables <- colnames(My_list[[1]])[unlist(lapply(My_list[[1]], is.numeric))]
for (name in names(My_list)){
My_sub <- My_list[[name]]
for (v in variables) {
nombre <- paste0(name, '_', v)
loop_Shapiro[[nombre]] <- lillie.test(My_sub[,v])
}
}
Obteniendo el siguiente error:
Error: Must subset columns with a valid subscript vector.
i Logical subscripts must match the size of the indexed input.
x Input has size 1 but subscript `complete.cases(x)` has size 960.
Run `rlang::last_error()` to see where the error occurred.
Para el segundo test realizo el siguiente código:
Selecciono las variables a las que quiero realizar el test de levene:
name_col<-colnames(checking[1:8])
loop_levene <- list()
for (val in name_col) {
f <- formula(paste0('`', val, '`', ' ~ Conut_stage'))
loop_levene <- leveneTest(f, checking)
}
Y este código funciona bien pero me almacena solamente el resultado de la última variable, es decir LBXTCA
. Me gustaría saber cómo hacer para que no se sobreescriban los resultados.
¿Cómo podría solucionarlo?