Difference between revisions of "GSCAN"
Jump to navigation
Jump to search
| Line 7: | Line 7: | ||
==UK Biobank== | ==UK Biobank== | ||
| − | + | <sytaxhighlight lang="bash"> | |
| − | |||
| − | < | ||
#!/bin/sh | #!/bin/sh | ||
#SBATCH -A S00000524 | #SBATCH -A S00000524 | ||
| Line 29: | Line 27: | ||
assertNotEmpty f | assertNotEmpty f | ||
| − | vcftools --gzvcf chr20impv1.vcf.gz --positions $f --recode --stdout | | + | vcftools --gzvcf chr20impv1.vcf.gz --positions $f --recode --stdout | bgzip -c > $f.chr20impv1.vcf.gz |
| − | |||
</syntaxhighlight> | </syntaxhighlight> | ||
| Line 36: | Line 33: | ||
<syntaxhighlight lang="bash"> | <syntaxhighlight lang="bash"> | ||
| − | cat list_of_files_with_position_chunks.txt | | + | cat list_of_files_with_position_chunks.txt | parallel sbatch --qos janus --reservation=janus-serial chunk.sh |
| − | |||
</syntaxhighlight> | </syntaxhighlight> | ||
Revision as of 04:59, 18 March 2016
GSCAN GWAS
We track studies in two ways. First, we track analyses with this Google doc. Second, we track study contact info with this Google sheet. These documents are only available to study personnel.
We also track our list of dbGaP studies on this Airtable.
UK Biobank
<sytaxhighlight lang="bash">
- !/bin/sh
- SBATCH -A S00000524
- SBATCH --mem=20gb
- SBATCH --ntasks-per-node=12
- SBATCH --nodes=1
- SBATCH --qos janus
- SBATCH --time=0:30:00
assertNotEmpty() {
: "${!1:? "$1 is empty, aborting."}"
}
- variable
- ${f:=$1}
export f assertNotEmpty f
vcftools --gzvcf chr20impv1.vcf.gz --positions $f --recode --stdout | bgzip -c > $f.chr20impv1.vcf.gz </syntaxhighlight>
Then run the command with
cat list_of_files_with_position_chunks.txt | parallel sbatch --qos janus --reservation=janus-serial chunk.sh