jingju11
"For each time point, I have 6 measurement.Two strains and each strain has three measurements from three rats. .."
2 strains*6measures*5 time points=60 measures;
so how to understand the latter sentence 'Two strains and each strain has three measurements from three rats'?
since each rat had multiple measures, it is necessary to give the IDs of those rats measured. we need to differentiate
measures from different mouse at the same time and from the same mouse at different time.
vickie
不好意思,那个地方我朋友写错了。如图,应该是对于每个时间点来说,有两种老鼠,每一种老鼠有六只,所以一共有12个measurements.每一次测量需要杀了老鼠,所以每个时间点的老鼠都是不同的,因此不需要ID.请问大家对于这种问题有什么好的方法吗?
JaneQ
ANOVA需要很多条件。可以用 kruskal- wallis。
vickie
谢谢楼上的答复,但无论是ANOVA还是kruskal- wallis,都是把时间当成一个因素,把对老鼠的处理当成另一个因素,有没有什么方法能检验mRNA的含量随时间的变化呢?
JaneQ
一个简单的办法是对老鼠的处理分组讨论。
ypchen
有时间因素吗? 也许要用纵向数据分析 分层模型
还不太懂背景 只是随便说说 不一定正确
ice024
试一下minitab中GLM,结果见下面:
一般线性模型: mRNA 与 strain
因子 类型 水平数 值
strain 固定 2 gk, wky
mRNA 的方差分析,在检验中使用调整的 SS
来源 自由度 Seq SS Adj SS Adj MS F P
time 1 5.244 5.461 5.461 7.41 0.009
strain 1 15.071 15.071 15.071 20.45 0.000
误差 56 41.278 41.278 0.737
合计 58 61.592
S = 0.858544 R-Sq = 32.98% R-Sq(调整) = 30.59%
项 系数 系数标准误 T P
常量 2.8470 0.2600 10.95 0.000
time -0.05356 0.01968 -2.72 0.009
mRNA 的异常观测值
拟合值 标准化
观测值 mRNA 拟合值 标准误 残差 残差
5 5.69704 3.13829 0.22216 2.55875 3.09 R
6 5.49991 3.13829 0.22216 2.36162 2.85 R
48 3.38747 1.27024 0.22597 2.11724 2.56 R
R 表示此观测值含有大的标准化残差
vickie
请问纵向数据分析,分层模型的英文名是什么?我可以自己去查查
vickie
不太知道minitab里的GLM是怎样运行的,选择了什么link呢?
vickie
从大家的回复中得到点启发,可以用ANCOVA来做吗?把老鼠类型当成qualitative factor,把时间当成quantitative factor,这样有道理吗?
vickie
想到一种比较复杂的方法
方程一:y=a1+a2x1+b1x2+b2x2^2+b3x1x2+b4x1x2^2
x1 代表老鼠的种类,for gk, x1=0, for wky, x1=1
x2代表时间,
when x1=0, y=a1+b1x2+b2x2^2
when x1=1,y=a1+a2+b1x2+b2x2^2+b3x2+b4x2^2
这样要判断这两条曲线是否相同,可以用所有数据fit方程一,检验whether a2=0, b3=0, b4=0
我比较担心的是,变量太多,数据较少,用此种方法可行吗?可信度如何?
bailanhong