# 遗传增益(genetic gain)如何计算标准误？

Δ G = h 2 S = h 2 ∗ ( μ s − μ p )

Δ G 2 = 100 ∗ h 2 S / μ p = 100 ∗ h 2 ∗ ( μ s − μ p ) / μ p

// An example

pm.esr<-echidna(fixed=height~1+Prov,

random=~Female+Block+Female:Block,

residual=~units,

es0.file='Provenance.es0')

// heritability

pin(pm.esr, mulp=c(h2~4*V3/(V1+V3+V4)))

#   Term Estimate    SE

# 1   h2    0.260 0.110

// count mean

Fem.mv<-aggregate(height~Female,data=df,mean)

mu.p<-mean(Fem.mv$height) raneff<-coef(pm.esr)$random

gca<-raneff %>% dplyr::filter(Term=='Female') %>%

arrange(-Effect)

sfem<-gca$Level[1:3] Fem.s<-dplyr::filter(Fem.mv,Female %in% sfem) mu.s<-mean(Fem.s$height)

pin(pm.esr, mulp=c(h2~4*V3/(V1+V3+V4),

delG.m~(mu.s-mu.p)*4*V3/(V1+V3+V4),

delG2.p~100*(mu.s-mu.p)*4*V3/((V1+V3+V4)*mu.p)

))

#      Term Estimate     SE

# 1       h2    0.260 0.1100

# 11  delG.m    0.181 0.0767

# 12 delG2.p    1.733 0.7320

AFEchidna包免费对外开放，仅用于学术研究，不可用于商业研究，否则造成的后果自负。

Zhang WH, Wei RY, Liu Y, Lin YZ. AFEchidna is a R package for genetic evaluation of plant and animal breeding datasets. BioRxiv. DOI: 10.1101/2021.06.24.449740.

https://wap.sciencenet.cn/blog-1114360-1305931.html

## 全部精选博文导读

GMT+8, 2022-1-27 22:47