## global variables ##
popname <- c("cool","mild","warm") # population names
n_loci <- 100 # number of loci 
shape1_pops <- c(2,3,5) # shape1 for beta distribution
shape2_pops <- c(50,50,50) # shape2 for beta distribution

# simulate data
test1 <- get_allele.freq(popname, n_loci, shape1_pops, shape2_pops)

# statistical test
freq_aov(x=test1)
##      F_value      P_value 
## 8.342161e+01 1.770953e-29 
## 
##  difference in average allele frequency is significant
# plots
freq_plot(x=test1)