配合力試驗(yàn)缺區(qū)數(shù)據(jù)分析思路

前段時(shí)間，有用戶(hù)咨詢(xún)過(guò)配合力試驗(yàn)數(shù)據(jù)缺區(qū)怎么處理的問(wèn)題，在此簡(jiǎn)單梳理一下思路，備忘。

上圖中是一個(gè)NCⅡ配合力試驗(yàn)，父本5個(gè)，母本4個(gè)，區(qū)組3個(gè)，每區(qū)組1個(gè)觀測(cè)值。A1B1組合在區(qū)組3的觀測(cè)值缺失1個(gè)。

缺區(qū)估計(jì)法、平均值填入法，都是“補(bǔ)齊”數(shù)據(jù)的策略，由于對(duì)后續(xù)數(shù)據(jù)分析有影響，因此不建議采用。

線性混合模型用于配合力缺區(qū)分析，本質(zhì)上跟本文給出的方法是一樣的，但是，由于矩陣求逆等計(jì)算的中間誤差或者迭代誤差，計(jì)算結(jié)果準(zhǔn)確性會(huì)略差，甚至方程組無(wú)解。往白里說(shuō)就是，一些平均數(shù)可以直接計(jì)算得出，而不是非要去解線性方程組，結(jié)果可能反而算不出來(lái)了。

下面是針對(duì)方差分析缺區(qū)情況博思公司給出的應(yīng)對(duì)參考建議，對(duì)配合力試驗(yàn)缺區(qū)情況同樣適用（配合力試驗(yàn)也需要進(jìn)行方差分析）。

博思公司提出的方差分析缺區(qū)情況應(yīng)對(duì)參考建議

1、缺區(qū)，并沒(méi)有影響試驗(yàn)材料的客觀試驗(yàn)表現(xiàn)，需要進(jìn)行相應(yīng)的平方和和自由度計(jì)算調(diào)整。

2、缺幾個(gè)自由度就減幾個(gè)，平方和采用累加的思路，缺了就不加，跟未缺區(qū)的時(shí)候一樣分析，試驗(yàn)分析參數(shù)進(jìn)行對(duì)應(yīng)調(diào)整。

3、如果缺區(qū)導(dǎo)致試驗(yàn)中某效應(yīng)的平均數(shù)無(wú)法計(jì)算，那就不能進(jìn)行方差分析，這可以作為缺區(qū)情況下是否能夠進(jìn)行方差分析的一個(gè)判斷依據(jù)。

4、F分布要求遵循正態(tài)分布，個(gè)別數(shù)據(jù)缺失，可以理解為對(duì)分布沒(méi)有影響。如果數(shù)據(jù)缺失嚴(yán)重，就不要進(jìn)行方差分析，F測(cè)驗(yàn)了，直接進(jìn)行品種平均數(shù)大小排序即可。

5、如果考慮采用缺區(qū)估計(jì)方式填值處理，應(yīng)在試驗(yàn)數(shù)據(jù)缺區(qū)未影響到效應(yīng)的平均數(shù)計(jì)算的情況下進(jìn)行，否則不宜采用估計(jì)方式。

6、后續(xù)進(jìn)行的多重比較，采用試驗(yàn)中真實(shí)獲得的品種平均數(shù)。

數(shù)據(jù)具體分析：

缺失的數(shù)據(jù)，對(duì)如下平均數(shù)的計(jì)算有影響。

• 總平均
• A1B1的平均數(shù)
• A1的平均數(shù)
• B1的平均數(shù)
• 區(qū)組3的平均數(shù)

所以在計(jì)算這5個(gè)平均數(shù)的時(shí)候進(jìn)行調(diào)整即可。就是分母多減去1.

線性混合模型構(gòu)建計(jì)算矩陣，采用0和1來(lái)標(biāo)注是否有觀測(cè)值。所以其對(duì)缺區(qū)的處理道理是一樣的。

圖中給出的計(jì)算結(jié)果，

A1B1的平均數(shù)是30=（31+29）÷2

A1B2的平均數(shù)是26=（27+25+26）÷3

專(zhuān)業(yè)的育種信息化研發(fā)團(tuán)隊(duì)，提供給客戶(hù)的不僅是產(chǎn)品，更有強(qiáng)大而專(zhuān)業(yè)的技術(shù)服務(wù)保障。如你在育種數(shù)據(jù)分析、育種標(biāo)簽打印方面有困惑，或有育種數(shù)據(jù)管理、育種標(biāo)簽需求，歡迎聯(lián)系我們。