美文网首页GWASGWAS学习
2018-11-02 GWAS实战(五)plink 进阶之数据过

2018-11-02 GWAS实战(五)plink 进阶之数据过

作者: 小郑的学习笔记 | 来源:发表于2018-11-02 19:12 被阅读54次

    之前讲了太多的很细的数据分析和过滤,在实战中,一般可以把命令写在一起。

    total

    这张图片对比了数据summary和数据过滤的命令之间的差异。

    对于我们自然群体,没有家系结构,所以主要是
    --mind N 删个体(基于缺失率)
    --geno N 删snp(基于缺失率)
    --maf N 删snp (基于最小等位基因频率)
    --hwe N 删snp (基于A和a不符合哈德温伯格平衡)
    --LD(这个一般不做双基因互作的话不考虑,详见我上一篇)

    plink --bfile core_v0.7 --mind 0.1 --maf 0.05 --geno 0.1 --hwe 0.01 --make-bed --out clean
    

    这样可以生成过滤好的文件,进行下一步分析

    生成的新文件

    clean

    相关文章

      网友评论

        本文标题:2018-11-02 GWAS实战(五)plink 进阶之数据过

        本文链接:https://www.haomeiwen.com/subject/dzqntqtx.html