статистика по списку фреймов данных

Question

У меня есть список из двух фреймов данных d $ 1 для пациентов ctrl, d $ 2 для больных пациентов. Каждый df содержит микробы Относительная численность от 3 пациентов:

List of 2
 $ CTRL  :'data.frame': 3 obs. of  18107 variables:
  ..$ Azorhizobium caulinodans                                           : num [1:3] 1.48e-07 1.62e-06 1.05e-06
  ..$ Buchnera aphidicola                                                : num [1:3] 9.63e-07 1.01e-06 8.09e-07
  ..$ Cellulomonas gilvus                                                : num [1:3] 1.63e-06 5.39e-07 4.05e-07
  ..$ Dictyoglomus thermophilum                                          : num [1:3] 2.30e-06 3.17e-06 1.34e-06
  ..$ Pelobacter carbinolicus                                            : num [1:3] 9.63e-07 3.70e-06 1.38e-06
  ..$ Shewanella colwelliana                                             : num [1:3] 9.63e-07 1.89e-06 1.62e-07
  ..$ Myxococcus fulvus                                                  : num [1:3] 1.78e-06 4.65e-06 1.50e-06
$ SICK:'data.frame':    3 obs. of  18107 variables:
  ..$ Azorhizobium caulinodans                                           : num [1:3] 4.24e-07 0.00 1.28e-06
  ..$ Buchnera aphidicola                                                : num [1:3] 5.45e-07 6.02e-07 4.47e-07
  ..$ Cellulomonas gilvus                                                : num [1:3] 3.03e-07 0.00 2.23e-07
  ..$ Dictyoglomus thermophilum                                          : num [1:3] 6.66e-07 2.75e-06 1.96e-06
  ..$ Pelobacter carbinolicus                                            : num [1:3] 9.69e-07 1.72e-07 1.62e-06
  ..$ Shewanella colwelliana                                             : num [1:3] 1.76e-06 6.02e-07 3.91e-07
  ..$ Myxococcus fulvus                                                  : num [1:3] 6.66e-07 8.60e-07 1.56e-06

Я хотел бы рассчитать некоторую статистику для каждого таксона (CTRL против SICK) и сохранить результаты для каждой ошибки в виде отдельного df (results.mw). Я пытался:

results.mw = lapply(mylist, function(d, l)
  {
  # Run wilcoxon by column
    as.data.frame(wilcox.test(d, l, exact = F)$p.value)
  }, d$"CTRL", l$"SICK")

, но я получаю сообщение об ошибке

Error in FUN(X[[i]], ...) : unused argument (l$SICK)

Answer 1

Вам нужно перебрать таксоны вместо исходного списка, который содержит два фрейма данных. Ниже я немного отредактировал код, он должен выполнить попарный тест. Я смоделировал данные, чтобы иметь что-то похожее на то, что у вас ..

# create data function
makeData = function(){
df = data.frame(matrix(rnorm(1000*3),3,1000))
colnames(df) = paste("S",1:1000,sep="_")
rownames(df) = letters[1:3]
return(df)
}
# create two data.frames
mylist = list(
       CTRL=makeData(),SICK=makeData()
)
# check 
str(mylist)
# although you said species are the same
# just to be sure
# we take intersection of species names
SPECIES = intersect(names(mylist$CTRL),names(mylist$CTRL))
# loop through species
p = sapply(SPECIES, function(i)
  {
  # Run wilcoxon by species
    wilcox.test(mylist$CTRL[,i],mylist$SICK[,i],exact=F)$p.value
  })
# gives you p-value by species
head(as.data.frame(p))