Difference between revisions of "GSCAN"
Jump to navigation
Jump to search
| Line 7: | Line 7: | ||
==UK Biobank== | ==UK Biobank== | ||
| + | |||
| + | Script to submit many jobs in parallel. | ||
<syntaxhighlight lang="bash"> | <syntaxhighlight lang="bash"> | ||
#!/bin/sh | #!/bin/sh | ||
| Line 27: | Line 29: | ||
assertNotEmpty f | assertNotEmpty f | ||
| − | vcftools --gzvcf chr20impv1.vcf.gz --positions $f --recode --stdout | bgzip -c > $f.chr20impv1.vcf.gz | + | vcftools --gzvcf chr20impv1.vcf.gz --positions $f --recode --stdout | \ |
| + | bgzip -c > $f.chr20impv1.vcf.gz | ||
</syntaxhighlight> | </syntaxhighlight> | ||
| Line 33: | Line 36: | ||
<syntaxhighlight lang="bash"> | <syntaxhighlight lang="bash"> | ||
| − | cat list_of_files_with_position_chunks.txt | parallel sbatch --qos janus --reservation=janus-serial chunk.sh | + | cat list_of_files_with_position_chunks.txt | \ |
| + | parallel sbatch --qos janus --reservation=janus-serial chunk.sh | ||
</syntaxhighlight> | </syntaxhighlight> | ||
Revision as of 04:50, 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
Script to submit many jobs in parallel.
#!/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
Then run the command with
cat list_of_files_with_position_chunks.txt | \
parallel sbatch --qos janus --reservation=janus-serial chunk.sh