Estoy intentando ejecutar un loop para ejecutar múltiples veces un test de shapiro. Para ello me he basado del siguiente link (https://stackoverflow.com/questions/58253937/loop-for-shapiro-wilk-normality-test-for-multiple-variables-in-r).
Siguiendo los pasos, realizo el subset de mi variable de interés:
Partiendo de mi database (checking 3835x9 variables):
Estadia Edad TALLA_D PESO_D IMC `Linfocitos totales` `Albumina total` `Colesterol total` Conut_stage
18 80 NA 67.5 NA 0.89 3.3 123 Light (2-4)
6 87 NA NA NA 6.4 3.9 120 Light (2-4)
2 66 NA NA NA 1.37 4.1 153 Low (0-1)
14 63 NA NA NA 1.04 3.8 239 Low (0-1)
2 58 180 129 40 2.44 NA NA NA
3 82 NA NA NA 1.18 4.2 130 Light (2-4)
5 18 NA NA NA 0.9 4.1 130 Light (2-4)
9 34 151 80 35 0.7 4 220 Light (2-4)
5 76 NA NA NA 1.02 NA NA NA
44 75 NA NA NA 1.42 3.6 165 Low (0-1)
My_list <- split(checking, f = list(checking$Conut_stage))
Preparo la lista: loop_Shapiro_1 <- list()
Y realizo el loop:
for (name in names(My_list)){
My_sub <- My_list[[name]]
loop_Shapiro_1[[name]] <- lillie.test(My_sub$Edad)
}
El problema es que me gustaría realizaro para todas las columnas que tengo en el dataframe (que son 9), en vez de ir realizando el loop por cada variable Edad, IMC, Estadia, etc.
Para ello he intentado cambiar el código del loop sin especificar la columna respecto a la que quiero realizar el loop:
for (name in names(My_list)){
My_sub <- My_list[[name]]
loop_Shapiro_1[[name]] <- lillie.test(My_sub)
}
Sin embargo obtengo 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 9 but subscript `complete.cases(x)` has size 960.
Run `rlang::last_error()` to see where the error occurred.
¿Qué deberia cambiar para que se hiciera el loop con todas las variables del dataframe?
Por otro lado estoy intentando realizar el mismo loop con el test de levene, el cual ni me funciona:
Recojo el nombre de las variables que quiero realizar el test de levene
name_col<-colnames(checking[1:8])
Creo la lista para recoger el output: loop_levene <- list()
Ejecuto el loop:
for (val in name_col) {
loop_levene<-leveneTest(val ~ Conut_stage, BD)
}
Y obtengo el siguiente error:
Error in model.frame.default(form, data) :
variable lengths differ (found for 'Conut_stage')
¿Qué estoy realizando mal? Gracias de antemano
EDIT1: siguiendo las sugerencias en los comentarios de este post he cambiado ambos códigos:
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 también he cambiado el código (siguiendo los comentarios del post):
name_col<-colnames(checking[1:8])
loop_levene <- list()
for (val in name_col) {
f <- formula(paste(val, '~ Conut_stage'))
loop_levene <- leveneTest(f, BD)
}
El cual obtengo el resultado deseado pero para la última variable a medir, es decir Colesterol total. Entiendo que se está sobreescribiendo por encima ¿Cómo puedo solucioanr este problema? Gracias de antemano.