R 语言 多项式分布的极大似然估计

本文介绍了在R语言中如何进行多项式分布的极大似然估计,特别是针对数量遗传学中的哈迪-温伯格平衡理论。通过建立似然函数,求解基因频率p的极大似然估计,并使用HardyWeinberg包中的Mourant数据集进行实际计算,验证理论结果。
摘要由CSDN通过智能技术生成

这里介绍一个数量遗传学中的多项式分布,哈迪-温伯格平衡(hardy-weinberg equilibrium)。该定律是说在一个大且随机交配的种群中,基因频率和基因型频率在没有迁移、突变和选择的条件下会保持不变。此时若等位基因M的频率为p,等位基因N的频率为q = 1-p,则该基因位点三种基因型频率为:

在这里插入图片描述
在这里插入图片描述
还是先写出似然函数:
在这里插入图片描述
对该函数取对数,计算极值,得到的就是p的极大似然估计:
在这里插入图片描述
下面用代码说明,这里用到了HardyWeinberg包中的Mourant数据集,该数

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Mrrunsen

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值