如何在嵌套数据帧中执行pairwise.t.test

问题描述 投票:1回答:1

我想在嵌套数据框中执行pairwise.t.test。有可能吗?

library(dplyr)
library(broom)
library(tidyr)
library(purrr)
a <- mtcars%>%
  dplyr::select(disp, gear,am)%>%
  nest(-am)%>%
  mutate(t_test = map(data, ~ pairwise.t.test(disp, gear, p.adjust = "bonferroni")))

# Error in factor(g) : object 'gear' not found
r purrr
1个回答
0
投票

这是否提供您想要的输出?

library(dplyr)
library(broom)
library(tidyr)
library(purrr)

mtcars_sub <- mtcars %>% 
  dplyr::select(disp, gear, am) %>% 
  group_by(am) %>% 
  nest()  


t_tests <- function(df){
  pairwise.t.test(df$disp, df$gear, p.adjust = "bonferroni")
}

map(mtcars_sub$data, t_tests)

#> [[1]]
#> 
#>  Pairwise comparisons using t tests with pooled SD 
#> 
#> data:  df$disp and df$gear 
#> 
#>   4    
#> 5 0.048
#> 
#> P value adjustment method: bonferroni 
#> 
#> [[2]]
#> 
#>  Pairwise comparisons using t tests with pooled SD 
#> 
#> data:  df$disp and df$gear 
#> 
#>   3     
#> 4 0.0027
#> 
#> P value adjustment method: bonferroni

让我知道这是否不是您需要的正确输出。

贷给@ ronak-shah

https://stackoverflow.com/a/62408521/8742237

© www.soinside.com 2019 - 2024. All rights reserved.