-
Notifications
You must be signed in to change notification settings - Fork 9
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
Association output only on SNPs #16
Comments
Hi, Thank you for your interest in HATK. It's highly likely that the bMarkerGenerator step failed. You should check If the problem isn't solved, then could you consider sending your input HLA type and SNP data? (Anonymize samples' label if necessary) I'll try to work on them. |
Hello! Thanks a lot for your reply. From our log file, it seems IMGT2Seq and NomenCleaner have run successfully (?) as no error was reported.
though we did not run IMGT2Seq and NomenCleaner separately - we ran the pipeline following the example in the documentation (see below), which says 'This command will implement (1) IMGT2Seq, (2) NomenCleaner, (3) bMarkerGenerator, (4) HLA_Analyzer(Association Test - logistic regression), (5) Manhattan Plot and (6) Heatmap Plot, which are the minimal components for HLA fine-mapping analysis.'
So is there any difference by running the above script vs running IMGT2Seq and NomenCleaner separately? The complete log file is also attached. Would appreciate your suggestions! And we are happy to send anonymized data for further examination. Thanks! |
Hi, schen2643. Thank you for your interest in HATK. No. there is no difference. The 'bMarkerGenerator' is practically the main process for HLA fine-mapping while it requires some preprocessings for data from the IPD-IMGT/HLA database and raw HLA type data. The IMGT2Seq and NomenCleaner are for this. In other words, Once those required preprocessings are done, users don't have to repeat them and can save time. The IMGT2Seq can be skipped with the '--dict-AA' and '--dict-SNPS' arguments and the NomenCleaner can be skipped with the '--chped' argument, where those arguments essentailly take output from those preprocessing steps. If Logistic regression result and Manhattan plot were generated well, then there would be no significant problem. In your log file, plotting heatmap for 8 HLA genes failed anyway. If you need help related to this, Please consider sending the anonymized genotype and HLA type data. |
@schen2643 & @ciaran-campbell Were you able to identify the issue. I am running into the same issue. And in the association file I see results for SNPs but NA for HLA markers and AAs. Any help is much appreciated. |
Hi there,
When I run HATK on one dataset, the association output has only worked for SNPs with rsids.
The following error is displayed and no heatmaps are produced?
Do you know what could be causing this?
Thanks for the help
The text was updated successfully, but these errors were encountered: