popgen-notes开发者指南：R代码实现遗传漂变模拟与选择系数计算

popgen-notes开发者指南：R代码实现遗传漂变模拟与选择系数计算

【免费下载链接】popgen-notesPopulation genetics notes项目地址: https://gitcode.com/gh_mirrors/po/popgen-notes

popgen-notes是一个专注于群体遗传学研究的开源项目，提供了丰富的R代码工具，帮助开发者和研究者轻松实现遗传漂变模拟与选择系数计算。本指南将带你快速掌握这两项核心功能的使用方法，无需深入复杂的理论细节。

遗传漂变模拟：从理论到实践

遗传漂变是群体遗传学中的关键概念，指由于随机抽样误差导致基因频率在世代间的随机波动。popgen-notes项目中的R代码为我们提供了直观的模拟工具。

核心模拟函数解析

在项目的Rcode/Genetic_drift目录下，WF_loss_of_het.R文件实现了基于 Wright-Fisher 模型的遗传漂变模拟。核心函数wf()接受种群大小(N)、模拟世代数(ngens)和初始等位基因频率(p0)等参数，返回各世代的等位基因分布矩阵：

wf <- function(N, ngens, p0=1/3, mu=0) { N <- 2*N # 二倍体调整 gns <- matrix(NA, nrow=ngens, ncol=N) gns[1, ] <- sample(1:2, N, replace=TRUE, prob=c(p0, 1-p0)) for (i in 2:ngens) { gns[i, ] <- gns[i-1, sample(1:N, N, replace=TRUE)] # 突变处理代码省略 } gns }

可视化遗传漂变效应

通过het()函数计算杂合度，并使用plot()函数可视化随世代变化的杂合度曲线：

het <- function(x) { tbl <- table(x) 1 - sum((tbl/sum(tbl))^2) }

下图展示了不同种群大小下杂合度的衰减情况，种群越小，遗传漂变效应越明显，杂合度下降越快：

快速上手模拟

只需几行代码即可完成一次遗传漂变模拟：

# 模拟种群大小为500，150代的遗传漂变 my.sims <- replicate(40, wf(N=500, ngens=150, p0=0.3)) # 计算并绘制杂合度变化 h <- apply(my.sims, 1, het) plot(h, type='l', xlab='世代', ylab='杂合度')

选择系数计算：揭示自然选择的力量

选择系数是衡量自然选择强度的重要参数，popgen-notes提供了多种计算和可视化选择系数的工具。

固定概率与选择系数关系

在Rcode/prob_fix_diffusion.R文件中，实现了基于扩散近似的选择系数计算。核心公式为：

prob.fix <- (1-exp(-s))/(1-exp(-2*N*s))

其中s为选择系数，N为有效种群大小。该公式描述了一个新突变等位基因的固定概率与选择系数的关系。

选择系数可视化

通过plot()函数可以直观展示不同选择系数下的固定概率：

s <- seq(-0.5e-3, 1e-3, length=1000) plot(s, prob.fix, type="l", lwd=3, xlab="选择系数, s", ylab="固定概率")

下图展示了不同有效种群大小下，选择系数与固定概率的关系曲线：

实际应用案例

在Rcode/cline.R文件中，选择系数被用于模拟等位基因频率的地理梯度变化：

sel = c(0.1, 0.01, 0.001) # 对非 resident 等位基因的选择系数

通过调整选择系数，可以模拟不同强度的自然选择对种群遗传结构的影响。

项目资源与扩展

主要代码目录

• 遗传漂变模拟：Rcode/Genetic_drift/
• 选择系数计算：Rcode/prob_fix_diffusion.R
• 群体遗传结构分析：Rcode/cline.R

进一步学习

要深入了解这些工具的理论背景，可以参考项目中的文档：

• 群体遗传学核心概念：define_terms.md
• 数学背景知识：math_background/Math-background.tex

开始使用popgen-notes

要开始使用popgen-notes项目，首先克隆仓库：

git clone https://gitcode.com/gh_mirrors/po/popgen-notes

进入项目目录后，可以直接运行R代码文件进行遗传漂变模拟和选择系数计算。项目中的代码经过精心优化，既适合教学演示，也可用于实际研究分析。

通过popgen-notes提供的R代码工具，你可以轻松探索群体遗传学的奥秘，从遗传漂变到自然选择，从理论模型到实际数据，全方位掌握群体遗传学的核心概念和分析方法。

【免费下载链接】popgen-notesPopulation genetics notes项目地址: https://gitcode.com/gh_mirrors/po/popgen-notes