Displaying 1 result from an estimated 1 matches for "bootfn".
Did you mean:
bootfs
2003 Oct 11
1
boot statictic fn for dual estimation of 2 stats?
...4,
1, 1, 1, 1, 3, 1, 1, 1, 4, 1, 1, 1, 1, 2, 3, 1, 1, 1, 2, 2, 2, 1, 2, 1, 1,
1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 3, 3, 1, 3, 2, 3, 2, 3, 1, 1, 3, 2,
2, 3, 2, 1, 4, 1, 3)
sanction <- data.frame(cost=cost, mil=mil, coop=coop)
est <- polr(as.factor(cost) ~ mil + coop, data = sanction)
bootfn <- function(data, i, object) {
d <- data[i,]
fit <- update(object, data = d)
c(fit$coef, fit$zeta)
}
res <- boot(sanction, bootfn, R = 10, object = est)