Hello, I am trying to estimate global FST and 95% CI for a data set. The original data was imported as genind object and converted to hierfstat format using genind2hierfstat() function. To estimate CI, I did the following:
library(hierfstat) trees.hf <- genind2hierfstat(trees.genind) df.bootfst <- boot.vc(trees.hf[,1], trees.hf[,-1], nboot=1000)
Error in 1:dim(data[dum, ])[1] : argument of length 0 The data formatting is not an issue because I was able to run basic.stats(trees.hf) successfully. I have also tried specifying full dimensions for the loci as follows without success:
df.bootfst <- boot.vc(trees.hf[,1], trees.hf[,2:500], nboot=1000)
I get the same error when using the varcomp.glob() function. Thank you for any help. Vikram