Simulate Responses from EBA Models

Simulates responses from the distribution corresponding to a fitted eba model object.

# S3 method for eba
simulate(object, nsim, seed, pool = TRUE, ...)

Arguments

object	an object of class eba, typically the result of a call to eba
nsim	currently not used
seed	currently not used
pool	logical, if TRUE (default), pooled responses (summed across respondents) are returned
...	further arguments passed to or from other methods; none are used in this method

Details

Responses are simulated by rbinom with sizes taken from the original sample and probabilities computed from the model object.

Value

A paired-comparison matrix of (pooled) responses.

See also

eba, boot.

Examples

data(celebrities)                     # absolute choice frequencies
A <- list(c(1,10), c(2,10), c(3,10),
          c(4,11), c(5,11), c(6,11),
          c(7,12), c(8,12), c(9,12))  # the structure of aspects
eba1 <- eba(celebrities, A)           # fit elimination-by-aspects model

## Parametric bootstrap of goodness-of-fit test
LR.stat <- replicate(200, deviance(eba(simulate(eba1), A)))

hist(LR.stat, col="lightgray", border="white", freq=FALSE, breaks=20,
     xlim=c(0, 60), main="Parametric bootstrap")
curve(dchisq(x, df=eba1$goodness.of.fit["df"]), add=TRUE)
abline(v=deviance(eba1), lty=2)