2018-10-30 GWAS實戰(zhàn)(二)plink 入門

有很多時候,看網上的教程,不如仔細閱讀工具的使用說明,我準備系統學習一下plink的用法。

按照官網的教程,下載下來的那個包里面有example

plink_mac

其實,看再多,不如一條命令跑一下來的實在

plink --file toy --freq --out toy_analysis
simple

好了跑完了,我們來分析一下

  1. 命令都是plink開頭的
  2. --file 代表的輸入文件 只要所有輸入文件的前綴就可以
  • Input data: '--file toy' tells PLINK to use the genomic data in the text files toy.ped and toy.map.
  1. --freq 這一層是表示要進行運算,后續(xù)有很多命令,換個名字而已,當然還有各種參數,這里freq指的是統計等位基因的頻率
  • Calculation(s)2 to perform: --freq tells PLINK to generate an allele frequency report. The full range of supported calculations is summarized under 'Main functions' in the sidebar.
  1. --out 指定輸出文件,只要指定前綴就可以,建議不要重復,因為文件多了會抹掉前面的文件。如果你不設置的話,應該是會報錯的

By default, the output files generated by PLINK all have names of the form 'plink.xyz', where '.xyz' is one of these extensions. This is fine for a single run, but as soon as you make more use of PLINK, you'll start causing results from previous runs to be overwritten.Therefore, you usually want to choose a different output file prefix for each run. --out causes 'plink' to be replaced with the prefix you provide. E.g. in the example above, '--out toy_analysis' caused PLINK to create a file named toy_analysis.frq instead of plink.frq.

我們來看生成的文件


result

以上就是對簡單的plink程序使用介紹,其實要學會使用plink做GWAS,還得有統計遺傳的基礎知識,要熟悉各種輸入的格式才行

?著作權歸作者所有,轉載或內容合作請聯系作者
【社區(qū)內容提示】社區(qū)部分內容疑似由AI輔助生成,瀏覽時請結合常識與多方信息審慎甄別。
平臺聲明:文章內容(如有圖片或視頻亦包括在內)由作者上傳并發(fā)布,文章內容僅代表作者本人觀點,簡書系信息發(fā)布平臺,僅提供信息存儲服務。

相關閱讀更多精彩內容

  • mean to add the formatted="false" attribute?.[ 46% 47325/...
    ProZoom閱讀 3,224評論 0 3
  • rljs by sennchi Timeline of History Part One The Cognitiv...
    sennchi閱讀 7,872評論 0 10
  • 清真寺買生牛肉,做干煸牛肉 打掃家里衛(wèi)生 運動1小時 看學習視頻1個
    喵喵喵嗚啦啦閱讀 136評論 0 0
  • 他說要和她一起看畫展,他坐14路公交車去找她。公車顛簸動蕩,走過一段泥濘的道路。這個孤獨的小城,似乎永遠都不會發(fā)展...
    Infinite清歡閱讀 269評論 2 1
  • 《賦能:打造應對不確定性的敏捷團隊》 斯坦利·麥克里斯特爾 美國陸軍四星上將,美軍駐阿富汗以及國際安全援助部隊的指...
    王俊婷閱讀 385評論 0 0

友情鏈接更多精彩內容