Как запустить leveneTest для 5834 генов одновременно

У нас есть 35 мезенхимальных стволовых клеток (МСК) одноклеточных данных RNA-Seq, и мы хотели бы сравнить гетерогенность экспрессии генов между различными условиями культивирования (т.е. гипоксией и нормоксией). Другими словами, мы хотели бы выявить гены, более гомогенные при гипоксии, чем при нормоксии.

Я знаю, как запустить leveneTest с помощью пакета car на R для одного гена (см. два примера ниже). Однако я не знаю, как запустить leveneTest для всех 5834 генов одновременно. Не могли бы вы мне помочь? Посмотрите наши данные с: http://68.181.92.180/~Gary/temporal/Log2_Transpose_CPM_MSC35Sample_CPM5Sample10_5834Gene.txt

Что еще более важно, как преобразовать все 5834 результата leveneTest в таблицу с двумя столбцами? Первый столбец — это символ гена, а второй столбец — значение P. Большое спасибо.

setwd("/Volumes/TOSHIBA EXT/20170305_LeveneTest/car")

require(car)

msc <- read.table("Log2_Transpose_CPM_MSC35Sample_CPM5Sample10_5834Gene.txt", header = T)

Тест:

leveneTest(CEBPB ~ Culture, data = msc, center=median)

Levene's Test for Homogeneity of Variance (center = median)

       Df F value  Pr(>F) 
 group  1  6.5486 0.01527 *
       33      

Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

А также:

leveneTest(PLIN2 ~ Culture, data = msc, center=median)

Levene's Test for Homogeneity of Variance (center = median)

  Df F value Pr(>F)
group  1  0.1804 0.6738
      33