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

博文

关于联合方差分析的讨论-2-联合方差分析的计算方法

已有 10432 次阅读 2016-3-9 14:48 |个人分类:农学统计|系统分类:科研笔记

   注:这里的公式,没法编辑,如果有需要,可以发我邮件:dengfei_2013@163.com。

   这里的计算方法主要是同质性检验(即:方差齐次性检验),联合方差分析中的同质性检验不是方差的检验,而是误差的检验,下面会有解释。

   同一个地点的误差同质性检测,虽然也是用bartlet检测,但不是检测的方差,而是检测的误差方差,是不一样的。在进行多点区试联合方差分析之前,需要首先对每个点的误差方差进行同质性检验(即齐性检验),如果方差同质,就可以将各试点结果列入汇总表。Bartlett的卡方计算公式为:

上个公式中的Vi = ni – 1,其中ni为样本容量,而C为矫正数:

]

 

如果将自然对数换成以10为底的对数(即由ln变成log),在一年多点试验中可以写作:

DF:为各试验点随机误差的自由度

L:为测试点个数

:为误差均方的平均数,即(

:为各试验点误差的均方

为各试验点随机误差均方MSe的平均数的对数

为各试验点误差均方MSe对数的和

如果小于,则说明各试验点方差齐次,否则就说明各试验点误差方差不齐次,不能进行联合方差分析,需要提出个别误差偏大的点,等到误差同质后再进行联合分析。

需要注意的是:

1、  联合方差分析测试的前提不是各个地点方差的齐次与否,而是各个地点误差的方差齐次与否,这里用计算机程序bartlet.test来测试是不符合实际的,因为它所计算的是各点的方差是否一致,而不是各点的误差方差是否一致。比如方差分析表中:


这是一个地点的9个品种、4个区组的方差分析结果,该地点的总方差平衡和是1.5775,自由度是35。误差的平方和是0.30983,自由度是24,这里如果检测9个地点的方差齐次,就是用的平方和1.5775

2、  如果本次试验有10个地点,这里的是误差的均方,比如上个方差表中,误差的均方为0.01291,需要对每一个地点进行单因素随机区组的方差分析,得到每个地点的误差均方,然后带入公式,计算出

如果不满足误差方差的齐次性,就不能进行联合方差分析。


$$



https://wap.sciencenet.cn/blog-2577109-961568.html

上一篇:关于联合方差分析的讨论-1 联合方差分析的定义
下一篇:关于联合方差分析的讨论3-联合方差分析的终结
收藏 IP: 114.250.168.*| 热度|

0

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

数据加载中...

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

GMT+8, 2024-5-14 09:43

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部