博文
SAS 处理裂区试验代码
||
建立了一个数据分析500人群，感兴趣的老师、同学
添加微信号： joinVSNC 为好友，拉你进入数据分析群
这里的数据是GenStat里面的Oats.ghs
data split;
input block$ w$ s$ v$ n$ y;
cards;
1 1 1 Marvellous 0.6cwt 156
1 1 2 Marvellous 0.4cwt 118
1 1 3 Marvellous 0.2cwt 140
1 1 4 Marvellous 0cwt 105
1 2 1 Victory 0cwt 111
1 2 2 Victory 0.2cwt 130
1 2 3 Victory 0.6cwt 174
1 2 4 Victory 0.4cwt 157
1 3 1 Goldenrain 0cwt 117
1 3 2 Goldenrain 0.2cwt 114
1 3 3 Goldenrain 0.4cwt 161
1 3 4 Goldenrain 0.6cwt 141
2 1 1 Marvellous 0.4cwt 104
2 1 2 Marvellous 0cwt 70
2 1 3 Marvellous 0.2cwt 89
2 1 4 Marvellous 0.6cwt 117
2 2 1 Victory 0.6cwt 122
2 2 2 Victory 0cwt 74
2 2 3 Victory 0.2cwt 89
2 2 4 Victory 0.4cwt 81
2 3 1 Goldenrain 0.2cwt 103
2 3 2 Goldenrain 0cwt 64
2 3 3 Goldenrain 0.4cwt 132
2 3 4 Goldenrain 0.6cwt 133
3 1 1 Goldenrain 0.2cwt 108
3 1 2 Goldenrain 0.4cwt 126
3 1 3 Goldenrain 0.6cwt 149
3 1 4 Goldenrain 0cwt 70
3 2 1 Marvellous 0.6cwt 144
3 2 2 Marvellous 0.2cwt 124
3 2 3 Marvellous 0.4cwt 121
3 2 4 Marvellous 0cwt 96
3 3 1 Victory 0cwt 61
3 3 2 Victory 0.6cwt 100
3 3 3 Victory 0.2cwt 91
3 3 4 Victory 0.4cwt 97
4 1 1 Marvellous 0.4cwt 109
4 1 2 Marvellous 0.6cwt 99
4 1 3 Marvellous 0cwt 63
4 1 4 Marvellous 0.2cwt 70
4 2 1 Goldenrain 0cwt 80
4 2 2 Goldenrain 0.4cwt 94
4 2 3 Goldenrain 0.6cwt 126
4 2 4 Goldenrain 0.2cwt 82
4 3 1 Victory 0.2cwt 90
4 3 2 Victory 0.4cwt 100
4 3 3 Victory 0.6cwt 116
4 3 4 Victory 0cwt 62
5 1 1 Goldenrain 0.6cwt 96
5 1 2 Goldenrain 0cwt 60
5 1 3 Goldenrain 0.4cwt 89
5 1 4 Goldenrain 0.2cwt 102
5 2 1 Victory 0.4cwt 112
5 2 2 Victory 0.6cwt 86
5 2 3 Victory 0cwt 68
5 2 4 Victory 0.2cwt 64
5 3 1 Marvellous 0.4cwt 132
5 3 2 Marvellous 0.6cwt 124
5 3 3 Marvellous 0.2cwt 129
5 3 4 Marvellous 0cwt 89
6 1 1 Victory 0.4cwt 118
6 1 2 Victory 0cwt 53
6 1 3 Victory 0.6cwt 113
6 1 4 Victory 0.2cwt 74
6 2 1 Goldenrain 0.6cwt 104
6 2 2 Goldenrain 0.4cwt 86
6 2 3 Goldenrain 0cwt 89
6 2 4 Goldenrain 0.2cwt 82
6 3 1 Marvellous 0cwt 97
6 3 2 Marvellous 0.2cwt 99
6 3 3 Marvellous 0.4cwt 119
6 3 4 Marvellous 0.6cwt 121
;
proc print data=split;
run;
proc anova;
class block v n;
model y= block v n block*v v*n;
test H=v e=block*v;
run;
https://blog.sciencenet.cn/blog-2577109-985634.html
上一篇：不同试验设计广义遗传率的分析
下一篇：ASReml-R多性状分析，怎么设置初始值？
扫一扫，分享此博文
全部作者的其他最新博文
全部精选博文导读
相关博文
- • 聚英才 建高地 | 北京理工大学“特立青年学者”全球招聘开启
- • 700年后日本或濒临灭绝？日本学者推算预测：届时或仅剩1名15岁以下孩子
- • [转载]【同位素视角】非英语母语学者如何区分’e.g.’, ‘i.e.’, ‘namely’与‘such as’等混淆难题
- • 美国佐治亚大学等机构学者：刈割策略对Bulldog 805紫花苜蓿+Tifton 85狗牙根混播草地产量及品质的影响
- • 美国堪萨斯州立大学、密苏里大学等机构学者研究成果：土壤水分管理策略和品种多样性对紫花苜蓿产量、营养品质和农场盈利能力的影
- • 德国、捷克草业科学学者长期放牧实验：异质草地斑块中的土壤有机碳储量和地下生物量