2018-11-06 GWAS實戰(zhàn)(七)plink 進階之關聯(lián)分析

這一節(jié)說的是association analysis 也是我們分析單標記回歸得到的結果(P值)的一步,這一步的結果可以用來與我第一講 聯(lián)合起來肮柜,形成一個閉環(huán)。跑跑標準流程。

Association analysis 可以有很多用處朴上,比如:

The basic association test is for a disease trait and is based on comparing allele frequencies between cases and controls (asymptotic and empirical p-values are available). Also implemented are the Cochran-Armitage trend test, Fisher’s exact test, di?erent genetic models (dominant, recessive and general), tests for stratified samples (e.g. Cochran-Mantel-Haenszel, Breslow-Day tests), a test for a quantitative trait; a test for dif- ferences in missing genotype rate between cases and controls; multilocus tests, using either Hotelling’s T(2) statistic or a sum-statistic approach (evaluated by permutation) as well as haplotype tests. The basic tests can be performed with permutation, described in the following section to provide empirical p-values, and allow for dierent designs (e.g. by use of structured, within-cluster permutation).

我這了主要介紹一個Linear and logistic models

These two features allow for multiple covariates when testing for both quantitative trait and disease trait SNP association, and for interactions with those covariates. The covariates can either be continuous or binary (i.e. for categorical covariates, you must first make a set of binary dummy variables).

這個主要是可以加入?yún)f(xié)變量作為控制,很靈活卒煞,但是可能速度會慢一點

說明上說最基礎的用法是這樣:


basic

但是這里我遇到一個問題
我bed bim fam 是不包含表型數(shù)據(jù)的痪宰,所以我要自己重新定義一個表型文件

我這里用quantitative traits作為例子

一般來說,就是自己設置一個文件
然后使用 --pheno 指定這個文件

--pheno causes phenotype values to be read from the 3rd column of the specified space- or tab-delimited file, instead of the .fam or .ped file. The first and second columns of that file must contain family and within-family IDs, respectively.

總共三列
前兩列是family and within-family IDs 第三列是表型

我這里用第一主成分作為表型

pheno

咱們來試試看

失敗了畔裕,程序運行錯誤

Warning: Skipping --linear since # variables >= # samples.
error.png

記住要加一個 --allow-no-sex

--allow-no-sex is now required if you want to retain phenotype values for missing-sex samples. This is a change from PLINK 1.07; we believe it would be more confusing to continue treating regular and --pheno phenotypes differently, and apologize for any temporary inconvenience we've caused.

程序就是這么設定的

plink --bfile clean --linear --pheno clean_one.eigenvec --allow-no-sex
done.png

成功啦

生成一個文件 assoc.linear

assoc.png

這個就可以用來畫圖啦

畫圖回到一

這里我解釋一下每一列的意義吧

  1. 染色體
  2. snp 名字
  3. base-pair 物理位置
  4. Tested allele (minor allele by default)
  5. Code for the test 估計就是模型
  6. Number of non-missing individuals included in analysis 個體數(shù)目
  7. Regression coefficient (--linear) or odds ratio (--logistic) 也就是beta值(回歸系數(shù))
  8. Coefficient t-statistic (beta除以standard error, 越大越顯著)
    9 Asymptotic p-value for t-statistic P值 看顯著

這個我就簡單介紹到這里

這里我還遇到一些實戰(zhàn)的問題
比如有些數(shù)據(jù)衣撬,我算出來極顯著,P值等于0 扮饶,這樣后面畫圖 ylim不能為無窮大具练,會產(chǎn)生問題,還有就是我的圖莫名其妙變瘦了甜无,很奇怪扛点。

?著作權歸作者所有,轉載或內容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市岂丘,隨后出現(xiàn)的幾起案子陵究,更是在濱河造成了極大的恐慌,老刑警劉巖奥帘,帶你破解...
    沈念sama閱讀 217,084評論 6 503
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件铜邮,死亡現(xiàn)場離奇詭異,居然都是意外死亡,警方通過查閱死者的電腦和手機松蒜,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,623評論 3 392
  • 文/潘曉璐 我一進店門扔茅,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人秸苗,你說我怎么就攤上這事召娜。” “怎么了惊楼?”我有些...
    開封第一講書人閱讀 163,450評論 0 353
  • 文/不壞的土叔 我叫張陵萤晴,是天一觀的道長。 經(jīng)常有香客問我胁后,道長店读,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 58,322評論 1 293
  • 正文 為了忘掉前任攀芯,我火速辦了婚禮屯断,結果婚禮上,老公的妹妹穿的比我還像新娘侣诺。我一直安慰自己殖演,他們只是感情好,可當我...
    茶點故事閱讀 67,370評論 6 390
  • 文/花漫 我一把揭開白布年鸳。 她就那樣靜靜地躺著趴久,像睡著了一般。 火紅的嫁衣襯著肌膚如雪搔确。 梳的紋絲不亂的頭發(fā)上彼棍,一...
    開封第一講書人閱讀 51,274評論 1 300
  • 那天,我揣著相機與錄音膳算,去河邊找鬼座硕。 笑死,一個胖子當著我的面吹牛涕蜂,可吹牛的內容都是我干的华匾。 我是一名探鬼主播,決...
    沈念sama閱讀 40,126評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼机隙,長吁一口氣:“原來是場噩夢啊……” “哼蜘拉!你這毒婦竟也來了?” 一聲冷哼從身側響起有鹿,我...
    開封第一講書人閱讀 38,980評論 0 275
  • 序言:老撾萬榮一對情侶失蹤旭旭,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后印颤,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體您机,經(jīng)...
    沈念sama閱讀 45,414評論 1 313
  • 正文 獨居荒郊野嶺守林人離奇死亡穿肄,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內容為張勛視角 年9月15日...
    茶點故事閱讀 37,599評論 3 334
  • 正文 我和宋清朗相戀三年年局,在試婚紗的時候發(fā)現(xiàn)自己被綠了际看。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 39,773評論 1 348
  • 序言:一個原本活蹦亂跳的男人離奇死亡矢否,死狀恐怖仲闽,靈堂內的尸體忽然破棺而出,到底是詐尸還是另有隱情僵朗,我是刑警寧澤赖欣,帶...
    沈念sama閱讀 35,470評論 5 344
  • 正文 年R本政府宣布,位于F島的核電站验庙,受9級特大地震影響顶吮,放射性物質發(fā)生泄漏。R本人自食惡果不足惜粪薛,卻給世界環(huán)境...
    茶點故事閱讀 41,080評論 3 327
  • 文/蒙蒙 一悴了、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧违寿,春花似錦湃交、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,713評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至掂咒,卻和暖如春才沧,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背绍刮。 一陣腳步聲響...
    開封第一講書人閱讀 32,852評論 1 269
  • 我被黑心中介騙來泰國打工糜工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人录淡。 一個月前我還...
    沈念sama閱讀 47,865評論 2 370
  • 正文 我出身青樓捌木,卻偏偏與公主長得像,于是被迫代替她去往敵國和親嫉戚。 傳聞我的和親對象是個殘疾皇子刨裆,可洞房花燭夜當晚...
    茶點故事閱讀 44,689評論 2 354

推薦閱讀更多精彩內容

  • 騎馬思維 病危的國王為挑選繼承人,給兩個兒子出了道題:“給你們兩匹馬彬檀,白馬給老大帆啃,黃馬給老二,你們騎馬到清泉邊去飲...
    祁國瑜閱讀 2,059評論 0 1
  • 伴隨冬季來的不僅只有逼人的寒氣窍帝,還有烤地瓜誘人的香氣努潘。 冬天和好友一起逛街淘吃,相對于奶茶、漢堡疯坤、烤串报慕、夾餅等之類...
    一粒麥子_6fd0閱讀 290評論 1 1
  • 汝為門中客,知我相思苦压怠。 青梅與君共眠冈,紅豆為誰煮? 蠟臺翻秋影菌瘫,素帳攏鴛鴦蜗顽。 垂淚相銅鏡,碧痕夜夜殤雨让。
    長安舊人閱讀 489評論 2 4
  • 無聊雇盖。未來,天天想栖忠,不去做刊懈,白搭。 房子認知太淺娃闲。
    九點三閱讀 152評論 0 0