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

Manhattan Plot #7

Open
dvitale199 opened this issue Sep 22, 2021 · 1 comment
Open

Manhattan Plot #7

dvitale199 opened this issue Sep 22, 2021 · 1 comment
Assignees
Labels
enhancement New feature or request

Comments

@dvitale199
Copy link
Owner

Is your feature request related to a problem? Please describe.
Need plotting method for summary statistics

Describe the solution you'd like
A clean, interactive manhattan plot implemented in plotly. zooming, snp selection would be great!

Describe alternatives you've considered
None

Additional context
None

@dvitale199 dvitale199 added the enhancement New feature or request label Sep 22, 2021
@dvitale199 dvitale199 assigned dvitale199 and nvk23 and unassigned ines570 Dec 16, 2023
@dvitale199
Copy link
Owner Author

create function from the following:

import gwaslab as gl

# read sumstat (you may need to read the gwas data in to a pandas dataframe first!)
EUR_sumstat = gl.Sumstats(
    "gwas_results/logistic.hybrid",
    snpid="ID",
    chrom="#CHROM",
    pos="POS",
    p="P",
    build="38"
)

#get lead variants
EUR_sumstat.get_lead(anno=True, build = "38", windowsizekb=500, sig_level=5e-8)

# annotate the lead variants
annot = ['chr1:5670117', 'chr2:97393816', 'chr16:46387290', 'chr19:44918903']

# manhattan plot
EUR_sumstat.plot_mqq(
    sig_level=5e-8,
    cut=15,
    anno="GENENAME",
    anno_set=annot,
    highlight=annot,
    highlight_windowkb=250,
    title="EUR CC",
    build="38",
    save="EUR_sumstat_gwaslab.manhattan.png"
)

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

No branches or pull requests

3 participants