0

所以我知道 R 有一个内置函数 t.test....基本上只是 t.test(y1,y2) 等但是我无法访问我想要比较的数据。

所以我基本上有两个单独的数据输出第一个被称为“数据”，它的输出类似于这个

Time    Kilometres
0   0
1   0.05
2   0.1
3   0.15
4   0.2
5   0.25
6   0.3
7   0.35
8   0.4
9   0.45
10  0.5

另一个输出称为“hunt”并具有此输出

cuts: [20,25)
   Time Kilometres
21   20        7.3
22   21        8.4
23   22        9.5
24   23       10.6
25   24       11.7
------------------------------------------------------------ 
cuts: [25,30)
   Time Kilometres
26   25       12.8
27   26       13.9
28   27       15.0
29   28       16.1
30   29       17.2
------------------------------------------------------------ 
cuts: [30,35)
   Time Kilometres
31   30       18.3
32   31       19.4
33   32       20.5
34   33       21.6
35   34       22.7

我的问题是，是否可以为每次切割进行单独的 T.tests。就像在将每个切割与我的第一个数据“数据”进行比较时，为每个切割获取单独的 p 值。

所以 cut 的 p 值：[0,5] = cut:[5:10] = etc 再次感谢

4

1 回答 1

0

到目前为止，尚不清楚您要在每个组中测试什么。可以用作两组t.test测试或一组测试。每次剪辑只有一组。当用作一组测试时，需要对运行测试的平均值有一个假设值，但尚不清楚哪种参考值是合适的。我想知道您是否真正想做的是在每个切割类别中进行线性回归以测试趋势？这将隐含地针对零斜率进行测试。这是经过轻微测试的代码：

 lapply( split( dat, cut.grp) , 
         function(dgrp) summary(lm( Kilometres ~ Time, data=dgrp))$coefficients [ ,"Pr(>|t|)"] )
于 2013-03-26T06:31:55.370 回答