Obtaining polygenic risk scores
Commands used to generate the best output (in Linux terminal):
Rscript PRSice/PRSice.R
--base PGS000074.txt
--target /home/sharedFolder/referenceData/ukb/imputed_genotypes/ukb_imp_chr#_v3,/home/sharedFolder/referenceData/ukb/imputed_genotypes/ukb_imp_chr1_v3.sample
--type bgen
--snp SNP --chr CHR --bp BP --A1 A1 --A2 A2 --stat BETA --pvalue P
# --base-maf allelefrequency_effect:0.2
# --maf 0.01
--binary-target T
--no-clump
--fastscore
--ignore-fid
--pheno /home/vsrinivasan75/cases_controls_new/pheno_white.txt
--cov /home/vsrinivasan75/cases_controls_new/covariates_white.txt
--cov-factor SEX
--quantile 100 --quant-break 10,20,30,40,50,60,70,80,90,100 --quant-ref 50
--out /ukb_prs/colorectal_cancer/using_prsice/PGS000074/PGS000074_results
--prsice PRSice/PRSice_linux
Rscript PRSice/PRSice.R
--base PGS000785.txt
--target /home/sharedFolder/referenceData/ukb/imputed_genotypes/ukb_imp_chr#_v3,/home/sharedFolder/referenceData/ukb/imputed_genotypes/ukb_imp_chr1_v3.sample
--type bgen
--snp SNP --chr CHR --bp BP --A1 A1 --A2 A2 --stat BETA --pvalue P
# --base-maf allelefrequency_effect:0.2
# --maf 0.01
--binary-target T
--no-clump
--fastscore
--ignore-fid
--pheno /home/vsrinivasan75/cases_controls_new/pheno_white.txt
--cov /home/vsrinivasan75/cases_controls_new/covariates_white.txt
--cov-factor SEX
--quantile 100 --quant-break 10,20,30,40,50,60,70,80,90,100 --quant-ref 50
--out /ukb_prs/colorectal_cancer/using_prsice/PGS000785/PGS000785_results
--prsice PRSice/PRSice_linux
Last updated