育种数据分析之放飞自我分享 http://blog.sciencenet.cn/u/yijiaobai 关注:生物统计,数量遗传,混合线性模型,生物信息,R,Perl,Python,GWAS,GS相关方法,文章及代码

博文

方差分析表 组合 区组 合并转化(写给金金)

已有 1709 次阅读 2018-9-18 17:46 |个人分类:农学统计|系统分类:观点评述

> library(agridat)
> data("oats")
> head(oats) 
 Blocks Nitrogen Subplots    Variety Wplots yield
1      1  0.6_cwt        1 Marvellous      1   156
2      1  0.4_cwt        2 Marvellous      1   118
3      1  0.2_cwt        3 Marvellous      1   140
4      1    0_cwt        4 Marvellous      1   105
5      1    0_cwt        1    Victory      2   111
6      1  0.2_cwt        2    Victory      2   130

> mod = aov(yield ~ Blocks + Variety*Nitrogen,data=oats)
> summary(mod)               
  Df Sum Sq Mean Sq F value   Pr(>F)    
Blocks            5  15875    3175  12.489 4.09e-08 ***
Variety           2   1786     893   3.513   0.0366 *  
Nitrogen          3  20021    6674  26.251 1.13e-10 ***
Variety:Nitrogen  6    322      54   0.211   0.9719    
Residuals        55  13982     254                     
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

想要结果类型如下:

图片.png

方法: Excel

将方差分析表, 放到Excel里面, 首先看组合, 组合的自由度 = P1 + P2  + P1XP2

即: 92 = 2 + 30 + 60

同理, 平方和SS也是一样的

MS = SS/DF

F = MS/MS(误差)

P值查表即可.

Excel结果:

图片.png

我将模板放到附件里面, 里面有具体的公式.

aov.xlsx



http://blog.sciencenet.cn/blog-2577109-1135556.html

上一篇:R语言求解混合线性方程组(有系谱)
下一篇:chrom markdown-here 插件测试

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2021-11-28 16:47

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部