最近有需求,对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空格”,搜出来了好几篇博文,我甚至不知道原创是谁的,真的有必要这样复制粘贴吗? ...