美文网首页R数量遗传或生统
GWAS分析-说人话(9)初尝关联分析禁果

GWAS分析-说人话(9)初尝关联分析禁果

作者: 医学小蛋散 | 来源:发表于2019-11-22 11:25 被阅读0次

前言

OK,这里只是初尝关联分析,有很多参数上的设置,及说明

请参考其他文章(有空会在最下面补上)。


书接上一回(GWAS分析-说人话(4)合并数据),当我们整理好了一个合并后的对照组和实验组数据后,就可以“发动咒语,启动大招了”!

/Users/seedson/Downloads/plink_mac_20190617/plink --noweb --bfile /Users/seedson/mergefilesclcvscontrol --allow-no-sex --assoc --out sccontrol

说人话:

告诉电脑Plink位置

/Users/seedson/Downloads/plink_mac_20190617/plink 

不要检查基于Web的版本

--noweb 

“施法对象”,明显是一个二进制文件(bfile),b格高的文件?b是二进制(binary)啦~

--bfile /Users/seedson/mergefilesclcvscontrol 

我的数据没有性别数据,不用看性别这一栏

--allow-no-sex 

#法术名称:关联分析!

--assoc 

#输出结果(就是你分析完后,得有个记录啊,分析后的数据要找个地方写下来啊!)

--out sccontrol

结束画面,看到Done 是不是呗爽?我是不会告诉你我一开始没有写--allow-no-sex报错了几次的~

系统自动生成文件

系统相应位置自动生成关联分析文件

好的,结果出来了。大白话就是:小白可能还不知道看哪个结果...

sccontrol.assoc 

有“.assoc”这个后缀的文件就是关联分析结果啦!

好了,怎样看这个文件?

不要告诉我双击就好!(GWAS分析-说人话(3)打开文件

其实重点是,我希望输出这个文件达到双击的目的,怎样做?

其实就是要后续分析啦,不是单纯为了看这个那么low的目的了。

sort -o

sort -o sort1.txt sccontrol.assoc 

获得:sort1.txt 文件(TXT!!!!!是不是很熟悉的面孔啊!)

右键Excel打开就好了~

看到这个P值栏,有没有很感动?!

Excel结果

这里就获得放进R或者python等作图工具的数据啦~

以上!~

相关文章

网友评论

    本文标题:GWAS分析-说人话(9)初尝关联分析禁果

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