对于田间试验和农业数据分析,请加入下面QQ群:VSNC GenStat 交流群 203276907
快速加入群,可以点击这个链接,里面有大量的资料可以下载和交流:点击链接加入群【VSNC GenStat 交流群】:https://jq.qq.com/?_wv=1027&k=42XYzLq
联合方差分析,是用方差分析的方法。首先要考虑各地点见的误差同质性。
数据:
数据如下:
> head(lianhe)
loc cul block yield
1 L1 1 1 2.205
2 L1 1 2 1.925
3 L1 1 3 2.445
4 L1 1 4 2.195
5 L2 1 1 2.400
6 L2 1 2 1.700
代码如下:
head(lianhe)
a
a
# 设置两个空变量
ms
DF
for(i in 1:length(a)){ #写一个循环,提取每个地点的自由度和误差均方
rr = lianhe[which(lianhe$loc == a[i]),]
result
DF[i] = result$df.residual;
ms[i] = sum(residuals(result)^2/DF[i])
}
msdata
msdata
# > msdata
#