data(copepod)
data(copepod) optimize(f=theta.likelihood,interval=c(10,100), maximum=TRUE, S=nrow(copepod), J=sum(copepod$count), give.log=TRUE) data(ostracod) preston(ostracod)
Run the code above in your browser using DataLab