1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
|
#!/bin/bash
if [ ! -f beagle.24Mar20.5f5.jar ]; then
echo
echo "Downloading beagle.24Mar20.5f5.jar"
wget http://faculty.washington.edu/browning/beagle/beagle.24Mar20.5f5.jar
fi
if [ ! -f bref3.24Mar20.5f5.jar ]; then
echo
echo "Downloading bref3.24Mar20.5f5.jar"
wget http://faculty.washington.edu/browning/beagle/bref3.24Mar20.5f5.jar
fi
echo
if [ ! -f test.24Mar20.5f5.vcf.gz ]; then
echo
echo "*** Downloading some 1000 Genomes Project data to file: test.24Mar20.5f5.vcf.gz ***"
wget http://faculty.washington.edu/browning/beagle/test.24Mar20.5f5.vcf.gz
fi
echo
echo "*** Creating test files: ref.24Mar20.5f5.vcf.gz target.24Mar20.5f5.vcf.gz ***"
echo
zcat test.24Mar20.5f5.vcf.gz | cut -f1-190 | tr '/' '|' | gzip > ref.24Mar20.5f5.vcf.gz
zcat test.24Mar20.5f5.vcf.gz | cut -f1-9,191-200 | gzip > target.24Mar20.5f5.vcf.gz
echo
echo "*** Running test analysis with \"gt=\" argument ***"
echo
java -jar beagle.24Mar20.5f5.jar gt=test.24Mar20.5f5.vcf.gz out=out.gt
echo
echo "*** Running test analysis with \"ref=\" and \"gt=\" arguments ***"
echo
java -jar beagle.24Mar20.5f5.jar ref=ref.24Mar20.5f5.vcf.gz gt=target.24Mar20.5f5.vcf.gz out=out.ref
echo
echo "*** Making \"bref3\" file ***"
echo
java -jar bref3.24Mar20.5f5.jar ref.24Mar20.5f5.vcf.gz > ref.24Mar20.5f5.bref3
echo
echo "*** Running test analysis with \"bref3\" file ***"
echo
java -jar beagle.24Mar20.5f5.jar ref=ref.24Mar20.5f5.bref3 gt=target.24Mar20.5f5.vcf.gz out=out.bref3
|
