Question

我的数据如下：

> Group Time Result
>     A 1   933
>     A 2   992
>     A 3   1007
>     A 4   1428
>     A 5   1068
>     A 6   721
>     A 1   1175
>     A 2   900
>     A 3   875
>     A 4   1775
>     A 5   986
>     A 6   963
>     A 1   1394
>     A 2   958
>     A 3   919
>     A 4   1103
>     A 5   940
>     A 6   919
>     C 2   1127
>     C 3   990
>     C 4   1033
>     C 5   1073
>     C 6   817
>       ....

我总共有50个小组，每个grous有6个时间点，所以我希望将P值结果作为一个表，使用t.test：

>     Time：1                    
>       A   B   C   D   E
>     A /   N   N   N   N
>     B N   /   N   N   N
>     C N   N   /   N   N
>     D N   N   N   /   N
>     E N   N   N   N   /
>                                                                   
>     Time：2                    
>       A   B   C   D   E
>     A /   N   N   N   N
>     B N   /   N   N   N
>     C N   N   /   N   N
>     D N   N   N   /   N
>     E N   N   N   N   /

这可能吗？

感谢您的帮助。

Answer 1

您可以尝试：

library(tidyverse)
d %>% 
  split(.$Time) %>% 
  map(~pairwise.t.test(.$Result, .$Group, p.adjust.method = "none"))

R应用t.test将多个组和exprot结果作为表进行比较

1 个答案: