Scripts for the LARGE-PD PRS manuscript
Scripts for calculating and evaluating PRS models
- calc_PRS.R: script for calculating the GWAS-significant PRS
- PRS_functions.R: contains functions for calculating a PRS
- test_PRS.R: script for testing a PRS using logistic regression and 10 folds
- survival_analysis.R: script for generating K-M curves and Cox regression
- PRS_plots.R: script for generating AUC plots and OR plots
- plot_PRS_dist.R: script for plotting distribution of PD PRS
- PRS_wilcox.R: script for characterizing the PD PRS by LARGE-PD recruitment site
- PD_PRS_1KG.R: script for charaterizing the PD PRS in 1000 Genomes
- AUC_PVAL.R: script for getting pvalues for the AUC of different PRS models
- PRS_PCA_plots.R: script for plotting PCA and PRS distributions
- PRS.sh: sample bash script for calculating and testing PRS
- run_PRSice.sh: bash script for running PRSice-2
Scripts for performing haplotype analysis in SNCA
- hap_blocks.sh: bash script for using PLINK to estimate haplotype blocks
- haplotpes.R: script for extracting haplotypes from a phased VCF file
- hap_functions.R: functions for extracting haplotypes
- plot_haplotypes.R: script for plotting haplotypes
- test_haplotypes.R: script for testing haplotyeps for association with PD
- haplotype_descriptive_stats.R: script for getting descriptive stats from haplotype data
- convert_nexus.R: script for converting haplotpye data to nexus format for visualizing in PopArt
Data management scripts
- align_merge.sh: script for aligning and merging all datasets used in study
- subset_SNCA.sh: script for subsetting region around rs356182 in SNCA for datasets used in study
- merge_imputed.sh: script for merging imputed data from LARGE-PD and Luo et al.
- concat.sh: script for concatennating LARGE-PD and Luo et al. data
- align: contains scipts for aligning genotype data, with align.R being the primary script.
miscellaneous scripts used in this project
- filter_rels.R: script for resolving relative pairs
- get_PCSs.R: script for calculating PCs and a kinship matrix using PC-AiR and PC-Relate
- merged_GWAS.R: script for running a GWAS usng GENESIS on the merged LARGE-PD/Luo et al. data
- phase_beagle.sh : script for phasing using BEAGLE