library(adegenet); library(hierfstat) setwd("~/Desktop/ddRAD_epiRAD_powan/demultiplexed_reads/EpiRAD/Intermediate_files3/06.Population_genetics/combined/") genind1 <- read.structure("combined.recode.str", n.ind = 187, n.loc = 5116, onerowperind = FALSE, col.lab = 1, NA.char = "-9", ask = FALSE, row.marknames = 1, quiet = FALSE,col.pop = 2) whitefish.hfstat <- genind2hierfstat(genind1, pop = genind1$pop) AR <- allelic.richness(whitefish.hfstat, min.n = 16) lom_epi.rich <- mean(AR$Ar[,10]) lom_ddr.rich <- mean(AR$Ar[,11]) eck_epi.rich <- mean(AR$Ar[,3]) eck_ddr.rich <- mean(AR$Ar[,4]) slo_epi.rich <- mean(AR$Ar[,8]) slo_ddr.rich <- mean(AR$Ar[,9])