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

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

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

#Read SoybeanGenotype data fileD<-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 QTLu <-rep(0,31620) #marker effectsu[QTL] <- 1g <-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语言怎么实现模拟表型数据"这篇文章的内容，相信大家都有了一定的了解，希望小编分享的内容对大家有帮助，若想了解更多相关的知识内容，请关注行业资讯频道。