作者anovachen ( )
看板Statistics
标题Re: [问题] 各位大大 救救我 简单的生统 我却卡关><
时间Sun Oct 20 01:35:30 2013
应该是要以两因子两水准的方式分析,
这样交互作用(假设AB交互作用是正的)可以图解成:
^ .B+
| .
| .
| .
| B+.
| .B-
| .
| B-.
---------->
A- A+
A因子(A细菌) B因子(B细菌)
对照组 - -
A细菌组 + -
B细菌组 - +
AB细菌组 + +
可是two-way又repeated measure,有对应的无母数方法能用吗?
如果资料不服从常态性和同质变异数假设,
通常可以用Box-Cox转换,
但万一转换後还是不像常态分配就很麻烦了Orz
(但是,很多生命科学相关领域的人不作常态性检验,
就直接跑t-test和ANOVA,硕士论文还是照样写...)
--
以下是我分析的结果,使用SAS写的程式码...
data bac;
input id group $ A $ B $ d1 d2 d3 d4;
cards;
1 Con N N 1 3 4 5
...略
11 A Y N 2 3 4 5
...略
21 B N Y 2 2 3 3
...略
31 AB Y Y 1 2 3 4
...略
;
(考虑到原PO可能不愿将资料做为示范用途,
故在此不公开原PO的资料)
proc glm data=bac;
class A B;
model d1-d4=A B A*B / nouni;
repeated day 4 /summary printe ;
run;
可参考此文:
http://r97846001.blog.ntu.edu.tw/2010/08/23/sasrepeated-measures-anova_demo02/
http://www.pt.ntu.edu.tw/hmchai/SAS/SAScontinuous/SASrepeated.htm
PS:本程式码有点问题,运算时间过长无法停止,
可能要请高手指点一下...
※ 编辑: anovachen 来自: 111.255.12.110 (10/20 02:36)
※ 编辑: anovachen 来自: 111.255.12.110 (10/20 03:41)
1F:→ andrew43:不只硕士论文照样写, 连paper都照样发表. XD 10/23 20:38