R/barplot_N_alleles.R

In nikostourvas/PopGenUtils: A collection of useful functions to deal with genetic data in R

barplot.N.alleles <- function(obj) {
div <- summary(obj)

all <- as.data.frame(div$pop.n.all)
all <- all %>%
        mutate(population = rownames(all)) %>%
        rename(count = "div$pop.n.all")
all$maf <- maf

all$count <- as.integer(all$count)

# str(all)

p <- ggplot(all, aes(x=population, y=count) ) +
        geom_bar(stat = "identity",
                 color = "black", fill = "white",
                 width = .8) +
        geom_text(aes(label=count, vjust = 1.6)) +
        ggtitle("Total Number of alleles per population")


all_list <- list("data" = all, "plot" = p)
return(all_list)
}