"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.
5 天 后
不好意思,那个地方我朋友写错了。如图,应该是对于每个时间点来说,有两种老鼠,每一种老鼠有六只,所以一共有12个measurements.每一次测量需要杀了老鼠,所以每个时间点的老鼠都是不同的,因此不需要ID.请问大家对于这种问题有什么好的方法吗?
ANOVA需要很多条件。可以用 kruskal- wallis。
谢谢楼上的答复,但无论是ANOVA还是kruskal- wallis,都是把时间当成一个因素,把对老鼠的处理当成另一个因素,有没有什么方法能检验mRNA的含量随时间的变化呢?
一个简单的办法是对老鼠的处理分组讨论。
有时间因素吗? 也许要用纵向数据分析  分层模型

还不太懂背景  只是随便说说 不一定正确
试一下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 表示此观测值含有大的标准化残差
请问纵向数据分析,分层模型的英文名是什么?我可以自己去查查
不太知道minitab里的GLM是怎样运行的,选择了什么link呢?
从大家的回复中得到点启发,可以用ANCOVA来做吗?把老鼠类型当成qualitative factor,把时间当成quantitative factor,这样有道理吗?
想到一种比较复杂的方法

方程一: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

我比较担心的是,变量太多,数据较少,用此种方法可行吗?可信度如何?
10 个月 后