# NOT RUN {
data(butterflies)
# }
# NOT RUN {
optimal.params(butterflies)
# }
# NOT RUN {
#takes too long without PARI/GP
#Now the one from Etienne 2005, supplementary online info:
zoo <- count(c(pigs=1, dogs=1, cats=2, frogs=3, bats=5, slugs=8))
l <- logkda.R(zoo, use.brob=TRUE) # Use logkda() if pari/gp is available
optimal.params(zoo, log.kda=l) #compare his answer of 7.047958 and 0.22635923.
# }
Run the code above in your browser using DataCamp Workspace