# data
data(data.JSM)
cat(data.JSM$genepop, file="JSM_genepop.txt", sep="\n")
cat(data.JSM$popname, file="JSM_popname.txt", sep=" ")
popdata <- read.genepop(genepop="JSM_genepop.txt", popname="JSM_popname.txt")
# Gst estimation
result.gst <- GstN(popdata)
write.csv(result.gst, "result_Gst.csv", na="")
Run the code above in your browser using DataLab