simulate.blim {pks} | R Documentation |
Simulates responses from the distribution corresponding to a fitted
blim
model object.
## S3 method for class 'blim' simulate(object, nsim = 1, seed = NULL, ...)
object |
an object of class |
nsim |
currently not used. |
seed |
currently not used. |
... |
further arguments passed to or from other methods. None are used in this method. |
Responses are simulated in two steps: First, a knowledge state is drawn
with probability P.K
. Second, responses are generated by applying
rbinom
with probabilities computed from the model
object's beta
and eta
components.
A named vector of frequencies of response patterns.
data(DoignonFalmagne7) blim1 <- blim(DoignonFalmagne7$K, DoignonFalmagne7$N.R) simulate(blim1) ## Parametric bootstrap disc <- replicate(200, blim(blim1$K, simulate(blim1))$discrepancy) hist(disc, col = "lightgray", border = "white", freq = FALSE, breaks = 20, main = "Parametric bootstrap", xlim = c(.05, .3)) abline(v = blim1$discrepancy, lty = 2)