R语言怎么实现模拟表型数据

这篇“R语言怎么实现模拟表型数据”文章的知识点大部分人都不太理解，所以小编给大家总结了以下内容，内容详细，步骤清晰，具有一定的借鉴价值，希望大家阅读完这篇文章能有所收获，下面我们一起来看看这篇“R语言怎么实现模拟表型数据”文章吧。

GWAS根据遗传率和QTL的数量产出模拟表型数据

#Read Soybean
Genotype data file
D<-read.big.matrix("GN.txt", type="char", sep="\t",head = TRUE)
dim(D)
D=D[,2:31261]
D1=as.data.frame(as.matrix(D))
D2=t(D1)
QTL <-100*(1:20) #pick 20 QTL
u <-rep(0,31620) #marker effects
u[QTL] <- 1
g <-as.vector(crossprod(D2,u))
h3 <- 0.6      #heritability 
y <- g +rnorm(346,mean=0,sd=sqrt((1-h3)/h3*var(g)))
#Saving simulated filewrite.table(y,"H60Q20.txt", sep="\t")

以上就是关于“R语言怎么实现模拟表型数据”这篇文章的内容，相信大家都有了一定的了解，希望小编分享的内容对大家有帮助，若想了解更多相关的知识内容，请关注亿速云行业资讯频道。