-
Notifications
You must be signed in to change notification settings - Fork 1
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
Error in run_eggla_gwas when use_info=TRUE #80
Comments
I hardcoded the imputation quality field in INFO as "INFO" while it could also be "R2". |
This should be fixed now, can you try to install the development version? remotes::install_github("mcanouil/eggla") PS: do not forget to provide the proper information for the imputation quality metrics by changing the default value for |
Thanks, got the 0.17.4 and running the analysis now. |
Ok, let me know if there is still an issue or if it goes smoothly, and in that case I'll release |
Will do, with use_info=TRUE it took about 2,5 days to finish. |
I'll keep the issue open until you can confirmed it is fixed then. |
This worked! I quickly looked at the results, seems like everything else is OK I just need to add CALL_RATE. |
Great! I'll start the release process for v0.17.4. I don't think the CALL_RATE is really needed though. |
Bug description
I'm trying run_eggla_gwas again, and managed to run it with argument use_info=FALSE. With the following command
I get error and message:
Formatting VCFs ...
[chr1.dose.vcf.gz] Performing PLINK2 regression ...
[chr1.dose.vcf.gz] Extracting INFO ...
[chr1.dose.vcf.gz] Combining results files ...
[chr1.dose.vcf.gz] Writing annotated results file ...
Error in colnamesInt(x, neworder, check_dups = FALSE) :
argument specifying columns specify non existing column(s): cols[21]='INFO'
Calls: run_eggla_gwas ... resolve.list -> signalConditionsASAP -> signalConditions
Execution halted
srun: error: localhost: task 0: Exited with exit code 1
Is there something that is missing in the arguments, or in my vcf files? I don't catch where it's coming from.
eggla
version outputChecklist
eggla
you're running, by pasting the output from runningpackageVersion('eggla')
in the "eggla
version output" text area?The text was updated successfully, but these errors were encountered: