如何解决重复测量方差分析
我有两种类型的数据集。
- 数据集包含 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)
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。