Индексирование переменных с помощью lapply и purrr :: map

Question

У меня есть data.frame:

dataset=structure(list(var1 = c(28.5627505742013, 22.8311421908438, 
95.2216156944633, 
43.9405107684433, 97.11211245507, 48.4108281508088), var2 = c(32.9009465128183, 
54.1136392951012, 69.3181485682726, 70.2100433968008, 44.0986660309136, 
62.8759404085577), var3 = c(89.6971945464611, 67.174579706043, 
37.0924087055027, 87.7977314218879, 29.3221596442163, 37.5143952667713
), var4 = c(41.5336912125349, 98.2095112837851, 80.7970978319645, 
91.1278881691396, 66.4086666144431, 69.2618868127465), var5 = c(33.9312525652349, 
88.1815139763057, 98.4453701227903, 25.0217059068382, 41.1195872165263, 
37.0983888953924), var6 = c(39.813664201647, 80.6405956856906, 
30.0273275375366, 34.6203793399036, 96.5195455029607, 44.5830867439508
), kmeans = structure(c(2L, 1L, 3L, 1L, 3L, 1L), .Label = c("1", 
"2", "3"), class = "factor")), .Names = c("var1", "var2", "var3", 
"var4", "var5", "var6", "kmeans"), row.names = c(NA, 6L), class = "data.frame")

Whit lapply и purrr::map in data.frame результат в порядке. См:

lapply(dataset[c(1:6)],shapiro.test)

purrr::map(dataset[c(1:6)],shapiro.test)

Ok. Теперь я хочу применить это в списке:

создать список (mylist):

set.seed(1234)
for(i in 1:6){
names<-paste0('var',i)
assign(names,runif(30,20,100))
}

dataset<-do.call(
cbind.data.frame,
mget(ls(pattern='*va'))
)

cluster<-kmeans(dataset,3)
dataset$kmeans<-as.factor(cluster[['cluster']])
mylist<-split(dataset,dataset$kmeans)
names(mylist)<-paste0('dataset',seq_along(mylist))

создать функцию (f):

f<-function(x){
  apply(x,2,shapiro.test)
}

после, примените эту функцию в lapply и purrr::map:

lapply(mylist[c(1:6)],f)
#Error: is.numeric(x) is not TRUE

purrr::map(mylist[c(1:6)],f)
#Error: is.numeric(x) is not TRUE

попробуйте это:

lapply(mylist[c(1:6)],function(x){
    lapply(x,shapiro.test)
})
#Error: is.numeric(x) is not TRUE 

lapply(mylist[c(1:6)],function(x){
  lapply(x,f)
})
#Error in apply(x, 2, shapiro.test) : dim(X) must have a positive length 

mylist[c(1:6)]%>%
  map(~map(.,shapiro.test))
#Error: is.numeric(x) is not TRUE

mylist[c(1:6)]%>%
  map(~map(.,f))
#Error in apply(x, 2, shapiro.test) : dim(X) must have a positive length

Что не так?

Answer 1

Документация apply гласит:

Если X не массив, а объект класса с ненулевым dim значение (например, фрейм данных), применить попытки привести его к массиву через as.matrix, если он двумерный (например, фрейм данных) или через as.array.

Поскольку последний столбец каждого фрейма данных является фактором (столбец cluster), вызов as.matrix приведет всю строку к символьному вектору (который не принимается в качестве входных данных для shapiro.test)

Это будет работать, если вы выберете только числовые столбцы в функции apply

f<-function(x){
  apply(x[c(1:6)],2 , shapiro.test)
}

lapply(mylist, f)

Примечание: попробуйте запустить as.matrix(dataset[1,]) и as.matrix(dataset[1,c(1:6)]), чтобы увидеть разницу.