Difference between revisions of "SEQLinkage Commands in Exercise"

From Statistical Genetics Courses

Jump to: navigation, search
Line 1: Line 1:
 
__NOTITLE__
 
__NOTITLE__
  
==SEQLinkage Commands in Exercise
+
==SEQLinkage Commands in Exercise==
 +
 
 
<pre>seqlink -h
 
<pre>seqlink -h
 
awk '{if ($0~/^#/) {print} else {{maf=substr($8,12)}; if (maf<=0.01) {print}}}' hearing_with_high_freq.vcf >hearing.vcf
 
awk '{if ($0~/^#/) {print} else {{maf=substr($8,12)}; if (maf<=0.01) {print}}}' hearing_with_high_freq.vcf >hearing.vcf

Revision as of 15:32, 7 June 2018

SEQLinkage Commands in Exercise

seqlink -h
awk '{if ($0~/^#/) {print} else {{maf=substr($8,12)}; if (maf<=0.01) {print}}}' hearing_with_high_freq.vcf >hearing.vcf
bgzip -c hearing.vcf >hearing.vcf.gz 
tabix -p vcf hearing.vcf.gz
seqlink --fam hearing.fam --vcf hearing.vcf.gz --format MERLIN
seqlink --fam hearing.fam --vcf hearing.vcf.gz --format MERLIN --output nocollapsing --bin 1
seqlink --fam hearing.fam --vcf hearing.vcf.gz --format MERLIN --output completecollapsing --bin 0
seqlink --fam hearing.fam --vcf hearing.vcf.gz --freq EXACSASMAF -o linkageanalysis -K 0.002  --moi AR -W 0 -M 1 --run-linkage