GWAS分析-说人话(3)网络数据下载后的处理

前言

上天会告诉你那个文件才是需要用的吗?

不会!

但我会!(点个赞,大个赏呗~)


数据下载后,相信我,小白连看哪个数据都会不知道的!

当然,数据分析的第0章,数据下载,也不是普通人可以handle的,

我可以说一天!~

首先我们需要申请dbGAP的权限,通过特殊的工具下载数据到硬盘中(或者服务器中)。

数据长这个样子:

我们要的数据在Genotype文件夹里面

一直切换到“.CADM”,然后你就会看到“map”和“ped”文件。(GWAS分析-说人话(2)认识文件名

首先,这里有一个经验性的注意地方:

硬盘可能不能读写数据,意味着输入指令后,不能够输出任何结果到硬盘上(小白就会一直报错,一直不知道为什么“我又错了”!)

所以要在自己的电脑,想要的地方,创建一个文件夹

(当然,你可以右建新建,不过这个做法是会被鄙视的!)

没错,这就是老屁股对新人的态度~

让我来告诉小白如何辗压“老屎忽”:

1.首先,cd到想要的路径

2.然后,midir XXX

(大神请忽略下面复习内容:

#复习:

#(删除文件夹实例:

#rm -rf/User/Dhyana/desktop (当然也可以cd到该文件夹,rm -rf 文件夹名称,删除)

#将会删除 /User/Dhyana/desktop目录以及其下所有文件、文件夹

#删除文件使用实例:

#rm -f /User/Dhyana/desktop/test.py

#将会强制删除/User/Dhyana/desktop/test.py这个文件)

#使用这个rm -rf的时候一定要格外小心,linux没有回收站的,删除之后再想找回就很难了)

说人话:

可能不能直接在硬盘上写数据,换个自己的工作目录!~

接着,开始下面的plink操作:

把map和ped文件转换成2进制文件(当然是为了后续分析用啊,不然才难得搞呢!)

(为什么要有这个转换?天下武功唯快不破,二进制文件电脑看得高兴啊!~)

(大神请忽略下面复习内容:

#复习一下plink的基本格式:

#plink --filechr2--make-bed --outchr2

调动Plink --处理得文件 --执行的操作 --给个名字

#把map和ped文件转换成2进制文件的代码:

#plink --filechr2--make-bed --outchr2

#这个基本的操作,就是要转化chr2(想要的系列文件,plink都是一组组数据一起处理的)成一个bed文件

#即:PED/MAP转为二进制格式)

好吧,反正有事网上“不华不实”的参考代码!!!!!

说人话:

现实是这样子的:

/Users/seedson/Downloads/plink_mac_20190617/plink --file /Volumes/Seagate\ Backup\ Plus\ Drive/70389_LungSomke/CGEMS/GENEVA_LungCancer/phs000093v2/p2/genotype/phg000206v1/phg000206.v1.GENEVA_LungCancer.genotype-imputed-data.c1.CADM/chr2 --make-bed --out chr2

我会给大家说人话的:

#调动Plink(查找你的plink在哪里了:

/Users/seedson/Downloads/plink_mac_20190617/)

#输入文件( 查找你的文件在哪里,拖进terminal:

/Volumes/Seagate\ Backup\ Plus\ Drive/70389_LungSomke/CGEMS/GENEVA_LungCancer/phs000093v2/p2/genotype/phg000206v1/phg000206.v1.GENEVA_LungCancer.genotype-imputed-data.c1.CADM/chr2 

#产生2进制文件

--make-bed 

#命名一个名字

--out chr2 

这是在个人的mac上运算的,考虑计算能力,只能一个一个做了

(大神可以在服务器上直接全跑了,求不要鄙视~)

运行时的样子:

会显示运行进度

运行后,在指定位置(刚刚cd到的地方)产生的数据:

产生的数据

#本例处理的是染色体3,后面提取的染色体,走的一样的程序,因为我要提取不同的染色体,所以其实就是不断改名字的反复操作了(如2,改3,3改5,根据染色体而定)如:

染色体3:

/Users/seedson/Downloads/plink_mac_20190617/plink --file /Volumes/Seagate\ Backup\ Plus\ Drive/70389_LungSomke/CGEMS/GENEVA_LungCancer/phs000093v2/p2/genotype/phg000206v1/phg000206.v1.GENEVA_LungCancer.genotype-imputed-data.c1.CADM/chr3 --make-bed --out chr3 

染色体4:

/Users/seedson/Downloads/plink_mac_20190617/plink --file /Volumes/Seagate\ Backup\ Plus\ Drive/70389_LungSomke/CGEMS/GENEVA_LungCancer/phs000093v2/p2/genotype/phg000206v1/phg000206.v1.GENEVA_LungCancer.genotype-imputed-data.c1.CADM/chr4 --make-bed --out chr4 

等等.......(我已经粗体需要修改的部分了,不用我把23个染色体都写出来吧?)

以上

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 203,324评论 5 476
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,303评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 150,192评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,555评论 1 273
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,569评论 5 365
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,566评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,927评论 3 395
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,583评论 0 257
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,827评论 1 297
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,590评论 2 320
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,669评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,365评论 4 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,941评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,928评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,159评论 1 259
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,880评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,399评论 2 342

推荐阅读更多精彩内容