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

What arguments should I pass to pass to run in exhaustive mode? #19

Open
sahwa opened this issue Aug 24, 2022 · 1 comment
Open

What arguments should I pass to pass to run in exhaustive mode? #19

sahwa opened this issue Aug 24, 2022 · 1 comment

Comments

@sahwa
Copy link

sahwa commented Aug 24, 2022

Hi @WansonChoi,

Thanks for your piece of software.

I'd like to run the omnibus test in exhaustive mode. However, it's not clear to me what needs passing to the --exhaustive argument on the command line. Trying --exhaustive True or any other string after results in the error HLAassoc: error: unrecognized arguments: True and leaving the argument empty --exhaustive results in HLAassoc/src/run_omnibus_test_LOGISTIC.R: error: argument --exhaustive: expected one argument Execution halted in the log file.

What should I pass to this argument to run in exhaustive mode?

My overall commmand is like this

python \
        -m HLAassoc OMNIBUS_LOGISTIC \
        --phased chr6.dose.G_group.filtered.${RC}.bgl.gz \
        --bim chr6.dose.G_group.filtered.${RC}.bim \
        --fam chr6.dose.G_group.filtered.${RC}.fam \
        --out output/hep_b_python_HLA_snps_logistic.${RC}.exhaustive \
        --pheno hep_b_phenos_omnibus.txt  \
        --covars hep_b_covars_omnibus.${RC}.txt \
        --exhaustive True \
        --aa-only \
        --maf-threshold 0

Thanks,

Sam

@yluo86
Copy link
Contributor

yluo86 commented Nov 16, 2022

sorry for the lack of documentation on this the --exhaustive argument should follow by the gene that one wants to perform exhaustive search on. You can also provide further argument to further restrict the search. For example, the following arguments will allow you three amino acid exhaustive search in HLA-B while fixing position 67 and 97

--exhaustive B --exhaustive-aa-pos 67 97 --exhaustive-min-aa 3 --exhaustive-max-aa 3 --exhaustive-no-filter

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