GWAS分析-说人话(21)-上位效应(epistasis)分析

前言

今天说一下如何“上位”~


1.大神橙子牛奶糖已经写得很清楚了,不多写:

https://www.cnblogs.com/chenwenyan/p/13813946.html

2.唯一真的要说的就是这个“二分类表型文件phenocc.txt”

只要看过GWAS分析-说人话(2)认识文件名,应该也问题不大

(https://www.jianshu.com/p/343ad060cc99)

就找到对应的列,复制黏贴到1个3列的txt文件处就行(Excel也能做,后面保存成txt就行),

想着都要用linux处理的话,就想多了。

3.如果不想做全基因组范围的,

 利用--extract snps.txt 准备基因型文件即可。

(本小结文件格式可以是ped/map或者bed/bim/fam,大神的格式为ped/map,用的是--file),

如使用bed/bim/fam,注意文件是用-bfile 即可

plink --noweb --bfile SNPs --pheno phenocc.txt --epistasis --epi1 1 --out test

后记:

Plink还是挺“银杏”的~

点赞超过100的话才写epistasis的理论东西吧~

以上~

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

推荐阅读更多精彩内容