GWAS大家都很熟悉了，随着目前越来越多的重测序文章的发表，可见如今SNP数据是非常的丰富，以及GWAS也发展了15年了，如今到了一个后GWAS时代。
GWAS我也做了很多了，经常有朋友问我，GWAS阈值如何确定的问题，其实这个问题很简单，基本上都是以Genetic Type I error方法计算独立的SNP即SNP（Me），然后以显著水平除以这个数就可以了。那么问题来了，独立的SNP怎么算呢？下面这个软件可以轻松快速上手：
Genetic type 1 E rror C alculator (GEC)是Miao-Xin Li老师写的软件，基于java，无需编译直接使用，操作很简单
#软件支持很多格式，我以我最常用的plink进行展示
java -jar /public/home/hguo/packages/gec/gec.jar -Xmx1g --effect-number --plink-binary gwas_all_chr.filter --genome --out ./leaf
#没错就这么简单，这个软件会把SNP 阈值以及SNP（Me）都算出来，也支持算某些区域的
基本上就这样，五分钟写一个帖子
网友评论