# NOT RUN {
data(msats.g)
msats.g <- stratify(msats.g, "fine")
# Just an overall Chi-squared test
ovl <- overallTest(msats.g, stats = "chi2", nrep = 100)
ovl
#' Just a pairwise test for Gst
pws <- pairwiseTest(msats.g, stats = "gst", nrep = 100)
pws
# }
# NOT RUN {
#' Both overall and pairwise tests for Fst and F'st
full <- popStructTest(msats.g, stats = c("fst", "fst.prime"))
print(full$overall)
print(full$pairwise)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab