Резюмировать статистику после lapply

У меня есть фрейм данных с несколькими столбцами и двумя разными группами — см. ниже.

set.seed(123) 
d <- data.frame(
  q1 = rnorm(20),
  q2 = rnorm(20),
  q3 = rnorm(20),
  group = sample(c("A", "B"), size = 20, replace = TRUE))

Я использую lapply для вычисления ttest для каждого столбца между двумя группами, как указано ниже:

lapply(d[,-4], function(i) t.test(i ~ d$group))

lapply возвращает для каждого столбца результаты, содержащие несколько статистических данных (я только что сообщил столбец q1)

$q1

    Welch Two Sample t-test

data:  i by d$group
t = -0.76262, df = 17.323, p-value = 0.4559
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
 -1.2294678  0.5759458
sample estimates:
mean in group A mean in group B 
    -0.05443279      0.27232820 

Я хочу резюмировать основную информацию о статистике (t, df, pvalue) в виде одной таблицы для каждого столбца (q1, q2, q3...)

Стоит ли изучать PHP в 2023-2024 годах?
Стоит ли изучать PHP в 2023-2024 годах?
Привет всем, сегодня я хочу высказать свои соображения по поводу вопроса, который я уже много раз получал в своем сообществе: "Стоит ли изучать PHP в...
Поведение ключевого слова "this" в стрелочной функции в сравнении с нормальной функцией
Поведение ключевого слова "this" в стрелочной функции в сравнении с нормальной функцией
В JavaScript одним из самых запутанных понятий является поведение ключевого слова "this" в стрелочной и обычной функциях.
Приемы CSS-макетирования - floats и Flexbox
Приемы CSS-макетирования - floats и Flexbox
Здравствуйте, друзья-студенты! Готовы совершенствовать свои навыки веб-дизайна? Сегодня в нашем путешествии мы рассмотрим приемы CSS-верстки - в...
Тестирование функциональных ngrx-эффектов в Angular 16 с помощью Jest
В системе управления состояниями ngrx, совместимой с Angular 16, появились функциональные эффекты. Это здорово и делает код определенно легче для...
Концепция локализации и ее применение в приложениях React ⚡️
Концепция локализации и ее применение в приложениях React ⚡️
Локализация - это процесс адаптации приложения к различным языкам и культурным требованиям. Это позволяет пользователям получить опыт, соответствующий...
Пользовательский скаляр GraphQL
Пользовательский скаляр GraphQL
Листовые узлы системы типов GraphQL называются скалярами. Достигнув скалярного типа, невозможно спуститься дальше по иерархии типов. Скалярный тип...
1
0
56
1
Перейти к ответу Данный вопрос помечен как решенный

Ответы 1

Ответ принят как подходящий

Вы можете снова использовать lapply() для извлечения каждого параметра и bind_rows():

library(dplyr)
lapply(l, function(x) {
  data.frame(t = x$statistic,
             df = x$parameter,
             pv = x$p.value) # returns a dataframe for each element in l
}) %>% bind_rows()

#           t        df         pv
# 1 -1.031983 13.533116 0.32017136
# 2 -2.458574  9.771018 0.03427922
# 3  1.421821 11.416813 0.18181697

Вы можете сделать это одним выстрелом:

lapply(d[,-4], function(i) {
  res <- t.test(i ~ d$group)
  data.frame(t = res$statistic,
             df = res$parameter,
             pv = res$p.value)
  }) %>% bind_rows()

Если вы хотите сохранить ссылку на имена столбцов, передайте .id на bind_rows():

lapply(d[,-4], function(i) {
  res <- t.test(i ~ d$group)
  data.frame(t = res$statistic,
             df = res$parameter,
             pv = res$p.value)
}) %>% bind_rows(.id='id')
#   id          t       df        pv
# 1 q1 -0.7626249 17.32329 0.4559469
# 2 q2 -1.6467070 17.73117 0.1172263
# 3 q3  0.5288851 13.01589 0.6057874

Пример:

set.seed(123) 
d <- data.frame(
  q1 = rnorm(20),
  q2 = rnorm(20),
  q3 = rnorm(20),
  group = sample(c("A", "B"), size = 20, replace = TRUE))
l <- lapply(d[,-4], function(i) {
  t.test(i ~ d$group)

  })

вы можете избежать bind_rows, используя sapply ? data.frame(sapply(d[,-4], function(i) { res <- t.test(i ~ d$group); c(t = res$statistic, df = res$parameter, pv = res$p.value) }))

Ronak Shah 22.05.2019 13:43

Знаете ли вы, как сообщать о colnames?

Al14 22.05.2019 15:15

Забудьте об этом, смотрите мое обновление, для .id есть параметр bind_rows(), который дает то, что вы хотите.

RLave 22.05.2019 15:37

Другие вопросы по теме

Передача более одного аргумента в FUN of lapply
Линейная интерполяция по списку фреймов данных
Удаление элемента из всех встроенных списков в списке в R
Есть ли функция, подобная lapply, для объекта S4?
Как запустить несколько линейных регрессий с разными независимыми переменными и зависимыми переменными, добавляя стандартизированные коэффициенты в R?
Поиск элементов вектора на основе числовых позиций, определенных в другом векторе в R
Как преобразовать список выходных данных регрессии в фреймы данных с помощью broom::tidy с помощью функции lapply?
Выберите/получите имена всех столбцов, которые имеют отрицательное значение от 0 до 10
Как использовать цикл for или Laaply для усечения случайной переменной t на биржевых данных для имитации цен опционов с использованием библиотеки foptions
Создайте несколько фреймов данных, отфильтровав подмножества, которые не равны значению в R