This function generates a Manhattan plot using the statistical summary data
obtained from the manhattan()
function. Top 3 genes (based on p-value) are annotated per chromosome.
Arguments
- data
Data frame containing the necessary columns from
manhattan()
output for plotting: -variant_position
: Variant position -variant_chromosome
: Variant chromosome -pval
: P-value -variant_id
: Variant ID -best_locus2genes_score
: Best locus2genes score -best_locus2genes_gene_symbol
: Best locus2genes gene symbol
Examples
if (FALSE) {
p <- manhattan(study_id = "GCST003044") %>% plot_manhattan()
p
}