Skip to contents

The function gg_Volcano() creates volcano plots from GAPIT GWAS results.

Specifying a folder and trait is all that is needed to create manhattan plots.

# Plot
mp <- gg_Volcano(
  # Specify a folder with GWAS results
  folder = "GWAS_Results/", 
  # Select a trait to plot
  trait = "DTF_Nepal_2017" )
# Save
ggsave("figures/gg_Volcano_01.png", mp, width = 12, height = 4)


Customized Plot

# Plot
mp <- gg_Volcano(
  # Specify a folder with GWAS results
  folder = "GWAS_Results/", 
  # Select a trait to plot
  trait = "DTF_Nepal_2017",
  # Plot only certain GWAS models
  models = c("MLM","MLMM","FarmCPU","BLINK"),
  # Highlight specific markers
  markers = c("Lcu.1GRN.Chr2p44545877", 
              "Lcu.1GRN.Chr5p1658484"), 
  # Create alt labels for the markers
  labels = c("44545877","1658484") 
  )
# Save
ggsave("figures/gg_Volcano_02.png", mp, width = 12, height = 4)