Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Trying to run PRSice on files received by a collaborator but BGEN file is not being opened #348

Open
varsh19 opened this issue Nov 14, 2023 · 1 comment

Comments

@varsh19
Copy link

varsh19 commented Nov 14, 2023

I used the following command to calculate PRS of a sequenced file coming from a collaborator. I imputed the vcf file which gave me separate vcf files for each chromosome. I then converted them to bgen and generated bgi and sample files too.

Rscript ../PRSice/PRSice.R
--base ../ukb_prs/PGS000785.txt
--target ../tempus_files/202030905/207532750086_R08C02/#,../tempus_files/202030905/207532750086_R08C02/1.sample
--type bgen
--snp SNP
--chr CHR
--bp BP
--A1 A1
--A2 A2
--stat BETA
--pvalue P
--binary-target T
--no-clump
--fastscore
--ignore-fid
--pheno ../tempus_files/202030905/207532750086_R08C02/207532750086_R08C02_pheno.txt
--out 207532750086_R08C02_results
--prsice ../PRSice/PRSice_linux

But I get the following error: Error: Cannot open bgen file ../tempus_files/202030905/207532750086_R08C02/1.bgen

The paths are all correct, but I cannot identify what the problem is. Please help me figure this out.

@choishingwan
Copy link
Owner

Did you check if you can see the file "../tempus_files/202030905/207532750086_R08C02/1.bgen" in the folder that your run the script?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants