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

发布时间:2022-03-21 10:11:14 作者:iii
来源:亿速云 阅读:287

这篇“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语言怎么实现模拟表型数据”这篇文章的内容,相信大家都有了一定的了解,希望小编分享的内容对大家有帮助,若想了解更多相关的知识内容,请关注亿速云行业资讯频道。

推荐阅读:
  1. 模拟实现strcmp
  2. python如何实现模拟数据结构模型?

免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。

r语言

上一篇:JavaScript如何交换变量值

下一篇:JavaScript如何扁平化数组

相关阅读

您好,登录后才能下订单哦!

密码登录
登录注册
其他方式登录
点击 登录注册 即表示同意《亿速云用户服务条款》