In queezzz/qzmle: Maximum likelihood estimation

devtools::install_github("queezzz/qzmle")

rfp <- transform(emdbook::ReedfrogPred, 
                 killed=density-surv, 
                 nsize=as.numeric(size)-1)

form <- killed ~ dbinom(size = density, 
                        prob = exp(log_a)/(1 + exp(log_a) * h * density))

## devtools::install_github("queezzz/qzmle")
library(qzmle)

qzmle::mle(form, 
           start=list(h=4,log_a=c(1,2)), 
           parameters=list(log_a~1+nsize), 
           links=list(a="log"), data=rfp)

qzmle::mle(form, start=list(h=4,log_a=2), 
           parameters=list(log_a~1+nsize),
           links=list(a="log"), 
           data=rfp, 
           method="TMB")

(emdbook::ReedfrogPred)

queezzz/qzmle documentation built on Nov. 24, 2021, 6:34 p.m.