這篇“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語言怎么實現模擬表型數據”這篇文章的內容，相信大家都有了一定的了解，希望小編分享的內容對大家有幫助，若想了解更多相關的知識內容，請關注億速云行業資訊頻道。
免責聲明：本站發布的內容（圖片、視頻和文字）以原創、轉載和分享為主，文章觀點不代表本網站立場，如果涉及侵權請聯系站長郵箱：is@yisu.com進行舉報，并提供相關證據，一經查實，將立刻刪除涉嫌侵權內容。
