Question

我有两种类型的数据集。

数据集包含 10 个枣品种和在 4 个不同成熟阶段收获的果实。品种效果是固定的，收获阶段的效果是随机的。
来自 2 个枣品种的果实在同一成熟阶段收获，进行 4 次采后处理，并储存在冷库中。每隔一周对水果进行定期分析，最长可达 35 天。品种和处理的效果是固定的，而贮藏间隔的效果是随机的。我想使用重复测量方差分析技术来分析来自这两种情况的数据。请帮助我使用可行的 R 命令和用于上述目的的有用软件包。早日回复将不胜感激。提前致谢。

问候

扎法尔·伊克巴尔

第一种情况的数据集

head(data, 10) 
sub  varieties stages width 
1    1     Anokhi     S1 27.08 
2    1     Anokhi     S1 27.15 
3    1     Anokhi     S1 26.87 
4    1 Desikarela     S1 26.21 
5    1 Desikarela     S1 23.47 
6    1 Desikarela     S1 23.77 
7    1   Dilbahar     S1 30.38 
8    1   Dilbahar     S1 26.80 
9    1   Dilbahar     S1 29.69 
10   1       Gurh     S1 20.57

str(data, 10) 'data.frame': 
120 obs. of 4 variables: 
 sub:int1111111111... 
 varieties: chr "Anokhi" "Anokhi" "Anokhi" "Desikarela" ... 
 stages:chr"S1""S1""S1""S1"... 
 width : num 27.1 27.1 26.9 26.2 23.5 ...

第二种情况的数据集

head(data, 10) 
sub var trt stor TSS 
1 1 Umran NoHWT 0day 10.00 
2 1 Umran NoHWT 0day 9.70 
3 1 Umran NoHWT 0day 9.80 
4 1 Umran HWT1 0day 10.20 
5 1 Umran HWT1 0day 9.75 
6 1 Umran HWT1 0day 9.86 
7 1 Umran HWT2 0day 9.65 
8 1 Umran HWT2 0day 9.72 
9 1 Umran HWT2 0day 9.81 
10 1 Umran HWT3 0day 9.45

str(data, 10) 
'data.frame': 144 obs. of 5 variables:
 sub:int1111111111...
 var : chr "Umran" "Umran" "Umran" "Umran" ...
 trt:chr"NoHWT""NoHWT""NoHWT""HWT1"...
 stor: chr "0day" "0day" "0day" "0day" ...
 $ TSS : num 10 9.7 9.8 10.2 9.75 9.86 9.65 9.72 9.81 9.45 ...

我在 R 软件中使用以下命令进行重复测量方差分析。这些是否可行或建议我最好的替代方案。

library(car)
options(contrasts = c("contr.sum", "contr.poly"))
model<-lm(vitc~varieties+stages+varieties:stages, data=data)
Anova(model, type="III")
library(agricolae)
(HSD.test(model, "varieties")$groups)
(HSD.test(model, "stages")$groups)

重复测量方差分析

0 个答案: