Question

我已经看到了非常类似的问题，但仍然无法解决这个简单的问题。

我想在分组变量上运行配对t检验，并将结果存储在数据帧中。包tidy中的命令broom执行此操作，但是当我运行下面的代码时，输出对于每个时间点都是相同的。我知道我可以使用循环来完成它，但我想在这里使用tidy并了解正在发生的事情。

这是一个可重复的例子：

library(dplyr)
library(broom)
df <- data.frame (time1=rep(1:4,30),
              sbp=runif(120, 100, 150),
              sbp1=runif(120, 120, 170))

R>head(df)
  time1      sbp     sbp1
1     1 146.9411 151.4842
2     2 102.7216 139.3380
3     3 125.7126 167.1806
4     4 126.0086 146.2177
5     1 149.9213 139.7968
6     2 117.6843 135.2726   



 z<-df %>%
 arrange(time1) %>%
 group_by(time1) %>%
 do(tidy(t.test(df$sbp,df$sbp1,paired=TRUE)))

每个时间点的结果输出相同：

R>head(z)
Source: local data frame [4 x 7]
Groups: time1 [4]

  time1  estimate statistic      p.value parameter  conf.low conf.high
  (int)     (dbl)     (dbl)        (dbl)     (dbl)     (dbl)     (dbl)
1     1 -22.48413 -11.69648 1.660608e-21       119 -26.29047 -18.67779
2     2 -22.48413 -11.69648 1.660608e-21       119 -26.29047 -18.67779
3     3 -22.48413 -11.69648 1.660608e-21       119 -26.29047 -18.67779
4     4 -22.48413 -11.69648 1.660608e-21       119 -26.29047 -18.67779

我确定有一个简单的修复，但不确定它是什么，任何建议都会很棒。感谢

Answer 1

如果您想引用dplyr内部使用的数据集，则需要使用.：

z <- df %>%
       arrange(time1) %>%
       group_by(time1) %>%
       do(tidy(t.test(.$sbp, .$sbp1,paired=TRUE)))     
Source: local data frame [4 x 7]
Groups: time1 [4]

  time1  estimate statistic      p.value parameter  conf.low conf.high
  (int)     (dbl)     (dbl)        (dbl)     (dbl)     (dbl)     (dbl)
1     1 -22.45646 -6.339261 6.307292e-07        29 -29.70157 -15.21135
2     2 -19.85310 -4.346528 1.550448e-04        29 -29.19485 -10.51136
3     3 -21.00503 -5.996117 1.609850e-06        29 -28.16968 -13.84037
4     4 -23.56341 -6.037170 1.438441e-06        29 -31.54606 -15.58077

引用df$sbp您忽略dplyr完成的分组，并且每次都使用它来完成数据集。这解释了为什么拟合t检验每次都相同。

R ttest使用整齐的

1 个答案: