
Create Manhattan Plot for Pleiotropic SNPs
Source:R/plot_pleiotropy_manhattan.R
plot_pleiotropy_manhattan.RdGenerates a Manhattan plot highlighting pleiotropic SNPs across chromosomes. The plot shows the distribution of significant associations and can highlight specific SNPs of interest.
Usage
plot_pleiotropy_manhattan(
pleio_data,
highlight_snp = NULL,
significance_line = 5e-08,
title = "Manhattan Plot of Pleiotropic SNPs"
)Arguments
- pleio_data
A data.frame containing pleiotropy analysis results.
- highlight_snp
Character. SNP identifier to highlight (default: NULL).
- significance_line
Numeric. Genome-wide significance threshold as a p-value, drawn as a red line (default: 5e-8).
- title
Character. Plot title (default: "Manhattan Plot of Pleiotropic SNPs").
Examples
data(gwas_subset)
pleio_results <- detect_pleiotropy(gwas_subset)
if (nrow(pleio_results) > 0) {
p <- plot_pleiotropy_manhattan(pleio_results, highlight_snp = "rs814573")
print(p)
}