meta.t1.results<-burdenMeta(c1,c2,wts=function(maf){maf<0.01},SNPInfo=SNPInfo) 代码语言:javascript 复制 meta.t1.results<-burdenMeta(c1,c2,wts=1,mafRange=c(0,0.01),SNPInfo=SNPInfo) 通过上面两种方式会得到同样的结果,所有我们只需要设置一个参数即可。 4. skatOMeta 最佳SKAT分析方法。结合SKAT...
GWAS(Genome-wideassociation study),即全基因组关联分析,是指在人类全基因组范围内找出存在的序列变异,即单核苷酸多态性(SNP),从中筛选出与疾病相关的SNPs(摘自百度百科)。如何让GWAS的结果可视化,我们就用到了曼哈顿图来展示其结果。那么在R语言中当然也有研究者开发了相关的R包“qqman”。 首先我们看下函数构成,...