dat <- data.frame(freq = c(243, 64, 58, 55)) # from Riefer et al. (2002)
mpt1 <- mpt(freq ~ list(
c*r,
(1 - c)*u^2,
2*(1 - c)*u*(1- u),
c*(1 - r) + (1 - c)*(1 - u)^2),
dat)
logLik(mpt1)
AIC(mpt1)
Run the code above in your browser using DataLab