Hi there,

I'm just using your diff_stats function to calculate phi_st values. I've tried using a SNP data set of ~18000 loci (a genind object filtered for monomorphs). It keeps producing the error "Error in dimnames(x) <- dn : length of 'dimnames' [1] not equal to array extent". If I drop the number of loci by heavily filtering with a MAF of 0.5 (leaves 18 loci) the function works fine. I was just wondering if this problem has been encountered before and if you might know a solution? Or perhaps it's not optimised to deal with large amounts of data?

Cheers