如何对以下数据做t test 用R或者sas均可

原始数据为txt格式如下， 请问各位大侠，如何用R或者SAS 两者之一均可 做出t test，要求ttest是 day-4 day2 day10 各进行一次，每次都比较 dose1和dose0、 dose5和dose0.
然后再对此数据做ANOVA分析。 各位大侠帮忙出出主意啊! 悬赏30个币！

day

data ex;
input
day        dose        rat_number        foodintake;
cards;
-4        0        1        16.784
-4        0        2        16.902
-4        0        3        15.504
-4        0        4        17.635
-4        0        5        17.086
-4        0        6        14.311
-4        0        7        13.817
-4        1        11        15.97
-4        1        12        14.767
-4        1        13        15.492
-4        1        14        17.261
-4        1        15        16.075
-4        1        16        16.821
-4        1        17        17.543
-4        1        18        10.62
-4        5        21        16.3
-4        5        22        18.369
-4        5        23        17.261
-4        5        24        15.757
-4        5        25        14.81
-4        5        26        15.006
-4        5        27        13.178
-4        5        28        15.352
2        0        1        16.139
2        0        2        16.614
2        0        3        14.29
2        0        4        15.752
2        0        5        16.424
2        0        6        17.301
2        0        7        11.643
2        1        11        14.505
2        1        12        12.899
2        1        13        16.876
2        1        14        14.921
2        1        15        16.32
2        1        16        13.401
2        1        17        12.869
2        1        18        12.108
2        5        21        11.605
2        5        22        14.63
2        5        23        13.899
2        5        24        9.394
2        5        25        13.09
2        5        26        13.285
2        5        27        12.524
2        5        28        15.351
10        0        1        14.805
10        0        2        16.513
10        0        3        12.158
10        0        4        17.465
10        0        5        16.428
10        0        6        14.295
10        0        7        15.993
10        1        11        15.37
10        1        12        15.425
10        1        13        17.87
10        1        14        0
10        1        15        16.734
10        1        16        15.245
10        1        17        16.241
10        1        16        13.825
10        5        21        15.498
10        5        22        18.901
10        5        23        17.425
10        5        24        18.078
10        5        25        18.04
10        5        26        18.903
10        5        27        16.789
10        5        28        18.314
;
run;
/*ttest*/
proc ttest data=ex(where=(dose in(0,5)));
by day;
class dose;
var foodintake;
run;
/*ttest*/
proc ttest data=ex(where=(dose in(0,1)));
by day;
class dose;
var foodintake;
run;


/*anova*/
proc glm data=ex;
class day dose;
model foodintake=day|dose;
run;

day        dose        rat number        foodintake
-4        0        1        16.784
-4        0        2        16.902
-4        0        3        15.504
-4        0        4        17.635
-4        0        5        17.086
-4        0        6        14.311
-4        0        7        13.817
-4        1        11        15.97
-4        1        12        14.767
-4        1        13        15.492
-4        1        14        17.261
-4        1        15        16.075
-4        1        16        16.821
-4        1        17        17.543
-4        1        18        10.62
-4        5        21        16.3
-4        5        22        18.369
-4        5        23        17.261
-4        5        24        15.757
-4        5        25        14.81
-4        5        26        15.006
-4        5        27        13.178
-4        5        28        15.352
2        0        1        16.139
2        0        2        16.614
2        0        3        14.29
2        0        4        15.752
2        0        5        16.424
2        0        6        17.301
2        0        7        11.643
2        1        11        14.505
2        1        12        12.899
2        1        13        16.876
2        1        14        14.921
2        1        15        16.32
2        1        16        13.401
2        1        17        12.869
2        1        18        12.108
2        5        21        11.605
2        5        22        14.63
2        5        23        13.899
2        5        24        9.394
2        5        25        13.09
2        5        26        13.285
2        5        27        12.524
2        5        28        15.351
10        0        1        14.805
10        0        2        16.513
10        0        3        12.158
10        0        4        17.465
10        0        5        16.428
10        0        6        14.295
10        0        7        15.993
10        1        11        15.37
10        1        12        15.425
10        1        13        17.87
10        1        14        0
10        1        15        16.734
10        1        16        15.245
10        1        17        16.241
10        1        16        13.825
10        5        21        15.498
10        5        22        18.901
10        5        23        17.425
10        5        24        18.078
10        5        25        18.04
10        5        26        18.903
10        5        27        16.789
10        5        28        18.314

说实话，没太看懂需求。

作配对分析的话，一般是
（１）　用差值与0比较。
（２）　用paired var1*var2 形式做。

webgu 发表于 2013-1-18 18:29
说实话，没太看懂需求。

作配对分析的话，一般是

谢谢 你好，要求是conduct  t-tests oneach day （day-4 day2, day10)  comparing the control group versus Treatment group 1(Dose1) . Compare the control group to the treatment group 2(Dose 5)
以上就是题目的全部内容，你能帮我看一下思路吗？谢谢！

ziyenano 发表于 2013-1-18 17:25
data ex;
input
day        dose        rat_number        foodintake;

谢谢您！ 我会自己好好研究的！ 非常感谢 不知如何将悬赏30币给您？

zcyijia 发表于 2013-1-18 16:42
day        dose        rat number        foodintake
-4        0        1        16.784
-4        0        2        16.902

刚刚选为您的最佳答案了! 悬赏点数应该就过去了吧！非常感谢您！

ziyenano 发表于 2013-1-18 16:19
data ex;
input
day        dose        rat_number        foodintake;

还有个问题请教您，比如在比较dose1和dose0中，day-4里面的t-test结果是

                 Method           Variances        DF    t Value    Pr > |t|

                                   Pooled           Equal            13       0.44      0.6644
                                   Satterthwaite    Unequal      12.281       0.46      0.6563
请问 如果想看ttest结果，以上就足够了对吗？ p值大于0.05说明差异不显著对吗？
但是为什么一个方差均等的pooled检验 一个是方差不等的检验呢？需要看哪个好呢？谢谢！

zcyijia 发表于 2013-1-19 23:44
还有个问题请教您，比如在比较dose1和dose0中，day-4里面的t-test结果是

                 Method     ...

下面还有一个方差齐性的检验；
如果方差相等，则参照equal的检验结果；
反之参照unequal的结果~