R/boot.mlcm.R

In MLCM: Maximum Likelihood Conjoint Measurement

boot.mlcm <- function(x, nsim, ...){
	d <- as.matrix(x$obj$data[, -1])
	p <- fitted(x)
	rsim <- matrix(rbinom(length(p) * nsim, 1, p), 
		nrow = length(p), ncol = nsim)
    bts.samp <- apply(rsim, 2, function(y, dd) {
        psct <- glm.fit(dd, y, family = binomial(x$link), ...)$coefficients
        names(psct) <- x$stimulus[-1]
        psct
    }, dd = d)
    list(boot.samp = bts.samp, 
    	bt.mean = apply(bts.samp, 1, mean), 
       bt.sd = apply(bts.samp, 1, sd), N = nsim)
	}