2018-11-05 GWAS实战(六)plink 进阶之主成分分析和群体分层

要学习这一节,首先是要学习什么是主成分,什么是聚类
这个我一直有点搞混,可以百度谷歌查一查
http://www.cnblogs.com/pinard/p/6239403.html
https://www.jianshu.com/p/794e91f60170

聚类这一块我做了,似乎结果不理想,以后参考这个http://www.cog-genomics.org/plink/1.9/strat#cluster

这里我着重介绍如何用plink做主成分http://www.cog-genomics.org/plink/1.9/strat#pca

PLINK 1.9 provides two dimension reduction routines: --pca, for principal components analysis (PCA) based on the variance-standardized relationship matrix, and --mds-plot, for multidimensional scaling (MDS) based on raw Hamming distances. Top principal components are generally used as covariates in association analysis regressions to help correct for population stratification, while MDS coordinates help with visualizing genetic distances.

感觉这一块学的不是很懂,我先继续走下去,后面再来补充

似乎如果已经定义了种群和聚类,可以在各个类中继续做PCA

这里我先不管,继续

plink --bfile clean --pca 5 --out clean

过程如下

pac.png

产生两个文件
eigenval eigenvec

clean.eigenval.png

这个文件是看每个主成分可以解释多少基因的变异

clean.eigenvec.png

这个文件是每个个体在每个主成分上的投影

这一章的内容我是为了做EigenGWAS 做准备

©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容

  • 凝眸如水 晓露痕, 犹带薄暮泪融姿。 皓肌如莹, 若花娇。 弱心笑颜光荡漾。
    叶靖枫阅读 224评论 0 2
  • 那笑容不是笑,仿佛一片羽毛轻飘飘 文化程度不高的人们总会道听途说,以己之见去对待别人 学习真的会改变一个人 人无信仰
    小刺猬Mandy阅读 162评论 0 0
  • 共读营排长碎碎念 1/10 我九排的小伙伴们。早上好今天我的城市在下雨,你的城市呢?希望你的城市能有个好天气。 1...
    王人三阅读 171评论 0 4
  • 无力感、恐慌感、自卑感 讲道理,四月比三月轻松许多了,可是我并没有觉得有多容易。反倒是,矫情劲儿上来,逃都没处逃。...
    便利店鲸鱼阅读 188评论 0 0