前段時間，有用戶咨詢過配合力試驗數(shù)據(jù)缺區(qū)怎么處理的問題，在此簡單梳理一下思路，備忘。

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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