r - 如何从半同胞设计中计算遗传力

Question

我正在尝试测量从半同胞设计收集的一组数据的性状，开花时间（FT）的遗传力。数据包括每株母株的 FT 和来自该母株的 2 个半兄弟姐妹，用于约 150 个不同的母系 (ML)。亲子关系不明。

我试过了：

1. 用母体 FT 和平均同胞 FT 的回归估计遗传力，并将斜率加倍。这工作得很好，并产生了 0.14 的估计值。

2. 运行 ANOVA 并使用 ML 之间的变异来估计加性遗传方差。从这个 powerpoint的幻灯片 25和这个线程中得到了这个想法，在方差计算之内和之间

   fit = lm(FT ~ ML, her)

   anova(fit)

她是数据集，在这种情况下，它只包括半同胞的 FT 值（我排除了母 FT 值以尝试遗传性）

从 ANOVA 输出中，我使用“ML”术语均方作为 ML 之间的变异，这也等于加性遗传方差的 1/4，因为半同胞之间的相关系数为 0.25。这个值原来是 0.098。此外，通过将其乘以 4，我可以获得加性遗传方差。我使用“残差”均方作为除“ML”术语之外的所有可变性。因此，所有方差减去加性遗传方差的 1/4。结果是 1.342。

然后尝试将可遗传性计算为Va/Vp = (4*0.098)/(1.342 + 0.098) = 0.39

这与我的斜率估计完全不同，我不确定我的推理是否正确。

我已经尝试过使用 R 的 sommer 和 heitability 包，但在使用半同胞设计时都没有成功，也没有找到任何一个包的半同胞设计示例。

有什么建议么？

标签： r