##############################################################
######################## GWAS analyses #######################
##############################################################

### BGENIE v1.2, see https://jmarchini.org/bgenie/
# Bycroft et al. (2018) Genome-wide genetic data on ~500,000 UK Biobank participants

for i in {1..22}; do

~/bgenie_v1.2_static1 \
--bgen ~/chr$i\_v3_MAF1_INFO4.bgen \
--pheno ~/phenotye_adj \ # the phenotype was adjusted for six ancestry principal components, batch and centre effects
--miss NA \
--scale_phenotypes \
--scale_genotypes \
--pvals \
--out chr$i \
--thread 6

done

##############################################################
######################## LDSC analyses #######################
##############################################################

# Bulik-Sullivan, B., et al. An Atlas of Genetic Correlations across Human Diseases and Traits. Nature Genetics, 2015. Preprint available on bioRxiv doi: http://dx.doi.org/10.1101/014498

# https://github.com/bulik/ldsc

## Munge data
~/munge_sumstats.py \
--sumstats trait1_results_no_MHC \ # summary statistics obtained from the GWAS excluding MHC region for trait1
--N 25285 \ # n of individuals
--out trait1_noMHC \
--merge-alleles w_hm3.snplist

## calculate genetic correlation between trait1 and trait2:
~/ldsc.py \
--rg trait1_noMHC.sumstats.gz,trait2_noMHC.sumstats.gz \
--ref-ld-chr eur_w_ld_chr/ \
--w-ld-chr eur_w_ld_chr/ \
--out trait_1-2_rg