我想在嵌套数据框中执行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
这是否提供您想要的输出?
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