Displaying 1 result from an estimated 1 matches for "holling2".
2011 Oct 17
1
simultaneously maximizing two independent log likelihood functions using mle2
...d AIC.
Can anyone give advice on how to do this?
My likelihood functions are long so I'll use the tadpole predation
example from Ben Bolker's book, Ecological Data and Models in R (p.
268-270).
library(emdbook)
data(ReedfrogFuncresp)
attach(ReedfrogFuncresp)
# Holling Type II Equation
holling2.pred = function(N0, a, h, P, T) {
a * N0 * P * T/(1 + a * h * N0)
}
# Negative log likelihood function
NLL.holling2 = function(a, h, P = 1, T = 1) {
-sum(dbinom(Killed, prob = a * T * P/(1 + a * h * Initial),
size = Initial, log = TRUE))
}
# MLE statement
FFR.holling2 = mle2(NLL.holling2,...