1. window=100k，step=2k 统计每个window的snp密度，然后用mixtools的normalmixEM（两个组分的混合模型）统计snp的分布模式。
R command：
library(mixtools)
SNPdensity=read.table("snp.density.file")
mixmdl=normalmixEM(SNPdensity)
plot(mixmdl,which=2)
mixmdl$mu mixmdl$sigma mixmdl$lambda
本文共 300 字，大约阅读时间需要 1 分钟。
1. window=100k，step=2k 统计每个window的snp密度，然后用mixtools的normalmixEM（两个组分的混合模型）统计snp的分布模式。
R command：
library(mixtools)
SNPdensity=read.table("snp.density.file")
mixmdl=normalmixEM(SNPdensity)
plot(mixmdl,which=2)
mixmdl$mu mixmdl$sigma mixmdl$lambda
转载于:https://www.cnblogs.com/yangyongzhi/p/6239244.html