File tree Expand file tree Collapse file tree 1 file changed +6
-6
lines changed
03.selection_signature/XP-EHH Expand file tree Collapse file tree 1 file changed +6
-6
lines changed Original file line number Diff line number Diff line change @@ -65,28 +65,28 @@ function main() {
6565
6666mkdir XP-EHH.progress
6767
68- # extract sample
68+ # extract sample
6969$vcftools --gzvcf $vcf --keep $ref --recode --recode-INFO-all --out ./XP-EHH.progress/01.ref
7070$vcftools --gzvcf $vcf --keep $tag --recode --recode-INFO-all --out ./XP-EHH.progress/01.tag
7171
7272cd XP-EHH.progress
7373for (( k= 1 ; k<= $chr ; k++ )) ;
7474do
75- # splite chr for ref and tag
75+ # splite chr for ref and tag
7676$vcftools --vcf 01.ref.recode.vcf --recode --recode-INFO-all --chr ${k} --out ref.chr${k}
7777$vcftools --vcf 01.tag.recode.vcf --recode --recode-INFO-all --chr ${k} --out tag.chr${k}
78- # calculate map distance
78+ # calculate map distance
7979$vcftools --vcf ref.chr${k} .recode.vcf --plink --out chr${k} .MT
8080awk ' BEGIN{OFS=" "} {print 1,".",$4/1000000,$4}' chr${k} .MT.map > chr${k} .MT.map.distance
8181done
8282
8383for (( k= 1 ; k<= $chr ; k++ )) ;
8484do
85- # XP-EHH
85+ # XP-EHH
8686$selscan --xpehh --vcf tag.chr${k} .recode.vcf --vcf-ref ref.chr${k} .recode.vcf --map chr${k} .MT.map.distance --threads $thread --out chr${k} .ref_tag
87- # norm
87+ # norm
8888$norm --xpehh --files chr${k} .ref_tag.xpehh.out --bp-win --winsize $win
89- # 加窗口步长
89+ # add win and step
9090python ../XPEHH_Win_step.py --file chr${k} .ref_tag.xpehh.out.norm --chr $k --window $win --step $step
9191done
9292cat {1.." $chr " }.XPEHH > ../${output} .XPEHH
You can’t perform that action at this time.
0 commit comments