b <- unlist(boeal) # concatenate all samples
ml <- mutestim(b) # maximum likelihood
gf <- mutestim(b,method="GF") # generating function
p0 <- mutestim(b,method="P0") # P0 method
cbind(ml,gf,p0) # compare 3 methods
# test values of mutations and fitness
flan.test(b,alternative=c("greater","less"),mutations0=0.6,fitness0=1)
b1 <- unlist(boeal[1:10]) # first 10 samples
b2 <- unlist(boeal[11:20]) # next 10 samples
flan.test(list(b1,b2)) # test equality
Run the code above in your browser using DataLab