計(jì)算家系平均遺傳力時(shí)锥腻，如果數(shù)據(jù)不平衡（小區(qū)株樹不一致）需要用到調(diào)和小區(qū)株樹（nk）假哎。
- 單點(diǎn)
- 多點(diǎn)
nk的計(jì)算在黃少偉的論文中（黃少偉; 鐘偉華 用小區(qū)平均值估算單株遺傳力的方法. 華南農(nóng)業(yè)大學(xué)學(xué)報(bào) 1998, 19, 76–81.）有介紹，計(jì)算公式為：
調(diào)和小區(qū)株樹.png
關(guān)鍵點(diǎn)是得到具體方法每個(gè)小區(qū)的株樹匣砖，加起來然后用地點(diǎn)數(shù)科吭、區(qū)組數(shù)和家系數(shù)的積除一下即可得。
- SPSS
用general linear model中的univariate分析猴鲫，響應(yīng)變量選擇一個(gè)最好沒有缺失值的變量对人，單點(diǎn)時(shí)自變量選擇family和block，多點(diǎn)則把site也選上拂共。當(dāng)然我們并不真的要做回歸或ANOVA牺弄，所以要在options中勾選描述性統(tǒng)計(jì)，這個(gè)才是我們需要的宜狐。出來的結(jié)果是這樣的：
output.png
我們需要最右列N势告，復(fù)制到Excel中，注意把所有Total的行去掉抚恒，按照公式計(jì)算就可以了咱台。
- R
n_s=2; n_r=5; n_f <- length(levels(df$Fam))
df_freq <- table(df$Fam, df$Rep) %>% as.data.frame() %>% filter(Freq!=0)
n_p <- n_r*n_f/sum(1/df_freq$Freq)
harm_plot <- function(Fam, Rep) {
n_f <- length(levels(Fam))
n_r <- length(levels(Rep))
df_freq <- table(Fam, Rep) %>% as.data.frame() %>% filter(Freq!=0)
n_p <- n_r*n_f/sum(1/df_freq$Freq)
return(n_p)
}
