跳至内容
售后
用户工具
登录
站点工具
搜索
工具
显示页面
修订记录
反向链接
最近更改
媒体管理器
网站地图
登录
>
最近更改
媒体管理器
网站地图
您的足迹:
xp_clr
编辑本页后请点击“保存”。请参阅
syntax
了解维基语法。只有在您能
改进
该页面的前提下才编辑它。如果您想尝试一些东西,请先到
playground
热身。
媒体文件
======简介====== XP-CLR (cross-population composite likelihood ratio):跨群体混合似然比法是利用群体间的等位基因频率的差异,对覆盖扩展连锁区域末端等位基因变异频率进行扫描,高XP-CLR的临近的窗口被聚集到同一种特征下,每一种特征代表一个单一选择扫描效应从而有效定位受正选择的位点。 ======功能====== 检测基因组选择信号分析 ======数据准备====== 单个样本的vcf文件 genotype矩阵,基因组索引等文件 ======数据分析====== 1 示例路径====== /TJNAS01/AFS_RESEQ/Share/GWAS/XPCLR/01.test ====== 2 分析脚本 <code> geno=/TJNAS01/AFS_RESEQ/Share/GWAS/XPCLR/01.test/test.chrz.geno.gz list=/TJNAS01/AFS_RESEQ/Share/GWAS/XPCLR/01.test/FT.vs.NFT chr=/TJNAS01/AFS_RESEQ/Share/GWAS/XPCLR/01.test/chr.list genofile=/TJNAS01/AFS_RESEQ/Share/GWAS/XPCLR/01.test/geno_file xpclr=/TJNAS01/AFS_RESEQ/Share/GWAS/XPCLR/00.bin/get_xpclr.pl fai=/TJNAS01/AFS_RESEQ/Share/GWAS/XPCLR/01.test/duck_chr_new.fa.fai chr=/TJNAS01/AFS_RESEQ/Share/GWAS/XPCLR/01.test/chr.list genegenenames=/TJNAS01/AFS_RESEQ/Share/GWAS/XPCLR/01.test/FTduck.gene.name geneID=/TJNAS01/AFS_RESEQ/Share/GWAS/XPCLR/01.test/FTduck.gene_ID outpre=FT.vs.NFT queue="afsreseq.q" perl $xpclr -outDir ./ -list $list -geno $geno -chr $chr -geno_file $genofile -gd 2.8 -slidingwindowsize 0.1 -mostwindowSNP 300 -gridsize 10000 -cutoff 0.99 -fai $fai -chr $chr -genegenenames $genegenenames -geneID $geneID -outpre $outpre -mis 2 -queue $queue </code> 3 流程参数 <code> -outDir 输出结果目录 -list 用来计算的群体信息,每个群体是一列 -geno 群体的基因型信息 -mis 过滤的mis值,不过滤填2(若基因型文件是已过滤的vcf生成的,则不必再过滤) -chr 用来计算的染色体,XPCLR分染色体计算(染色体id必须为数字) -geno_file 输出的单个染色体的基因型文件夹 (若存在则直接给目录) -GD 基因组大小/物种遗传距离(单位Mb/cM) (根据物种不同自行配置) -top 取top -mostwindowSNP 窗口内最大SNP数(过高会影响计算速度,过低影响准确性) -slidingwindowsize 窗口大小(单位为cM,用该值乘以GD即可得到窗口大小(物理)) -gridsize 步长 -cutoff 置信区间(一般0.95即可) </code> ======交付结果====== 选择信号分析的相关表格和图片
保存
预览
取消
编辑摘要
当您选择开始编辑本页,即寓示你同意将你贡献的内容按下列许可协议发布:
CC Attribution-Share Alike 4.0 International
xp_clr.txt
· 最后更改: 2023/03/27 08:07 由
hongxiang
页面工具
显示页面
修订记录
反向链接
回到顶部