boot.ppfst about hierfstat

[jgx65]

could you send a (small) example reproducing the behavior? Cheers

[jwwhale]

Hopefully you mean the following...
data <- read.genepop("medtest.gen", ncode = 2L)
Converting data from a Genepop .gen file to a genind object...
File description: Mediana_7773SNPs_12Pops
...done.
medtest <- boot.ppfst(data, nboot = 100, quant = c(0.025, 0.975), diploid = TRUE)
Error in i + 1 : non-numeric argument to binary operator
wc(data)$sigma.loc
SNP lsiga lsigb lsigw
M00001 2.566286e-02 3.802011e-03 0.47457627
M00002 6.928440e-02 4.692108e-02 0.34745763
M00003 3.936859e-02 2.079705e-01 0.13131313
...

Then when converting the genind object

data1 <- genind2hierfstat(data)
class(data1)
[1] "data.frame"
medtest<- boot.ppfst(data1, nboot = 100, quant = c(0.025, 0.975), diploid = TRUE)
Error in i + 1 : non-numeric argument to binary operator

[jgx65]

Sorry, what I meant was a small dataset reproducing the error you are getting. I suspect it is due to some loci being missing for whole populations, but would need to see the data to be sure. Cheers

[jgx65]

alternatively you can share the data via PM if you can't reproduce the behavior

[jgx65]

The issue seems to be with your population naming scheme, I have to dig further (I was able to reproduce your error once, but then the function ran smoothly). A quick fix is to rename populations using the following;

data1<-genind2hierfstat(data)
data1[,1]<-as.integer(data1[,1])
data.ci<-boot.ppfst(data1,...)

Will investigate

[jwwhale]

Awesome!! Thanks Jerome!