manhattan(d1,chr="Chr",bp="Pos",p="p",snp="Marker",main="Manhattan plot: logistic") tiff("y1-曼哈顿图.tiff") manhattan(d1,chr="Chr",bp="Pos",p="p",snp="Marker",main="Manhattan plot: logistic") dev.off() qq(d1$p,main="Q-Q plot of GWAS p-values : log") tiff("y1-...
MLM的矫正过于严格,会把一些真实相关的SNP标记也过滤掉,因此CMLM模型目的是重新检测到那些假阴性SNP标记。 $ vim cmlm.sh $ perl run_pipeline.pl-Xmx10g-Xms512m-fork1-vcfTroot.vcf-fork2-r trait.txt-fork3-qQ.txt-excludeLastTrait-fork4-kK.txt-combine5-input1-input2-input3-intersect-combine6-i...