# NOT RUN {
## Storage-retrieval model for pair clustering (Riefer & Batchelder, 1988)
data(retroact)
spec <- mptspec(
c*r,
(1 - c)*u^2,
2*(1 - c)*u*(1 - u),
c*(1 - r) + (1 - c)*(1 - u)^2,
u,
1 - u
)
m <- mpt(spec, retroact[retroact$lists == 0, ])
summary(m) # parameter estimates, goodness of fit
plot(m) # residuals versus predicted values
confint(m) # approximate confidence intervals
plot(coef(m), axes = FALSE, ylim = 0:1, pch = 16, xlab = "",
ylab="Parameter estimate (MPT model, 95% CI)")
axis(1, 1:3, names(coef(m))); axis(2)
arrows(1:3, plogis(confint(m))[, 1], 1:3, plogis(confint(m))[, 2],
.05, 90, 3)
## See data(package = "mpt") for application examples.
# }
Run the code above in your browser using DataCamp Workspace