最近有需求,對WGS測序獲得SNP信息進行篩減,可問題是測序個體少,call rate,maf,hwe,等條件過濾后,snp數量還是千萬級別,所以后面利用plink工具根據LD信息來濾除大量SNP標記。 工具版本:PLINK v1.90b4.6 64-bit (15 Aug 2017) 一、格式 ...
遺傳學研究中,標記間的變異連鎖不平衡程度直接影響Mapping的精確度,剖析相關群體中標記間LD情況可直接估計Mapping精度,具體的計算如下: plink bfile chr ld window ld window kb r ld window r . out chr bfile: 讀取的基因型Bim文件 ld window : 表達結果文件中輸出變異最大數量 ld window kb :計算 ...
2020-05-25 08:53 2 554 推薦指數:
最近有需求,對WGS測序獲得SNP信息進行篩減,可問題是測序個體少,call rate,maf,hwe,等條件過濾后,snp數量還是千萬級別,所以后面利用plink工具根據LD信息來濾除大量SNP標記。 工具版本:PLINK v1.90b4.6 64-bit (15 Aug 2017) 一、格式 ...
plink --bfile file --extract all.snp --r2 --ld-window-kb 1000 --ld-window-r2 0.8 --ld-snp-list all.snp --out all.snp.08.inter all.snp文件如下所示: ...
PLINK提供了“--ld”的參數計算兩個SNP位點的連鎖不平衡值。 命令如下: plink --file file --ld rs123 rs134 --out rs123_rs134 生成如下數據: --ld rs123 rs134: R-sq ...
獲取到的所有文件名放到數組中時必須加上引號,不然 for 循環時會被解析成通配符,或者使用 shell 字典,同樣也需要引號。 shell 字典示例 ...
Plink是一個基於Flink的流處理平台,旨在基於 [Apache Flink]封裝構建上層平台。 提供常見的作業管理功能。如作業的創建,刪除,編輯,更新,保存,啟動,停止,重啟,管理,多作業模板配置等。 Flink SQL 編輯提交功能。如 SQL 的在線開發,智能提示,格式化,語法校驗 ...
用--recode命令,--out表示轉化的文件的名字,本例已經命名為“filter” ...
VCFtools can convert VCF files into formats convenient for use in other programs. One such example is the ability to convert into PLINK format. ...
搜索“資源文件key空格”,搜出來了好幾篇博文,我甚至不知道原創是誰的,真的有必要這樣復制粘貼嗎? ...