I keep getting this error when trying to run OutFLANK:
Error in pi0est(p, ...) :
ERROR: The estimated pi0 <= 0. Check that you have valid p-values or use a different range of lambda.
I have biallelic SNP data and seem to be able to make the diploid Fst input file with no errors. I have included my code and data below - thanks in advance for your help.
Hi Katie,
I keep getting this error when trying to run OutFLANK: Error in pi0est(p, ...) : ERROR: The estimated pi0 <= 0. Check that you have valid p-values or use a different range of lambda.
I have biallelic SNP data and seem to be able to make the diploid Fst input file with no errors. I have included my code and data below - thanks in advance for your help.
load("1223ind_86snps_38pop.RData")
locus_names = locNames(data_filt)
poplist = as.vector(data_filt$pop)
data = read.table("1223ind_86snps_38pop.str.lfmm")
library(OutFLANK)
inputfile = MakeDiploidFSTMat(data, locusNames = locus_names, popNames = poplist)
outflank = OutFLANK(inputfile, LeftTrimFraction = 0.05, RightTrimFraction = 0.05, Hmin = 0.10, NumberOfSamples=38, qthreshold = 0.05)
Data_Files.zip