2013-03-26 47 views
0

所以我知道R有內置的函數t.test ....這基本上只是t.test(y1,y2)等 但我有訪問我想要比較的數據時遇到問題。T.test in R program for multiple data sets

所以我基本上有兩種獨立的數據輸出 第一個是隻稱爲「數據」和它具有類似於此

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 

的輸出的另一輸出被稱爲「尋線」,並具有該輸出

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值,同時將每個剪輯與我的第一個數據(稱爲「數據」)進行比較。

所以削減的p值:[0,5] = 削減:[5:10] = 等 再次感謝

+2

我知道這不是你問的問題,但你有沒有考慮過使用方差分析?你可以使用Tukey的測試來比較兩種治療方法。 – 2013-03-26 06:16:04

+0

多數民衆贊成在酷...很抱歉,我對統計數據空白....但會給我單獨的p值以及...? – 2013-03-26 06:32:21

+1

如果你不熟悉統計學,也許你可以從統計學家那裏尋求當地的幫助。 – 2013-03-26 08:09:43

回答

0

到目前爲止,尚不清楚要各組內測試的內容。 t.test可以用作兩組測試或一組測試。每個剪輯只有一個組。當作爲單組測試使用時,測試運行的平均值需要有一個假設值,但不清楚什麼樣的參考值是合適的。我想知道你真正想做的是在每個切分類別內進行線性迴歸以測試趨勢嗎?這將隱含地測試零斜率。這是測試代碼:

lapply(split(dat, cut.grp) , 
     function(dgrp) summary(lm(Kilometres ~ Time, data=dgrp))$coefficients [ ,"Pr(>|t|)"])