Как я могу получить значение p.value всех регрессий, используя функцию lapply - PullRequest
Новая
       15

Как я могу получить значение p.value всех регрессий, используя функцию lapply

1 голос
/ 23 марта 2020

Я использую пример с mtcars, чтобы позже использовать его при анализе моих данных. Я могу получить список резюме и функцию коэф для каждой регрессии. Мой вопрос: как я могу получить значение p.value всех регрессий в списке, используя функцию lapply?

Вот мой код 

    library(data.table)
    regressions <- 
    data.table(mtcars)[, 
          .(Myregressions = lapply(.SD, function(x) summary(lm(mpg ~ x)))), 
          .SDcols = -1]

Regressions$MyRegressions
Regressions[, lapply(MyRegressions, coef)]

1 Ответ

1 голос
/ 23 марта 2020

Мы можем l oop над 'MyRegressions' и извлечь столбец 

regressions[, lapply(Myregressions, function(x) coef(x)[, "Pr(>|t|)"])]

Если нам нужно извлечь только выбранные значения 

regressions[, list(lapply(Myregressions, function(x) {
      x1 <- coef(x)[, "Pr(>|t|)"]
      x1[x1 < 0.05]}))]$V1
...