I run a multinomial regression on a data set with an outcome that has three values. First, I build an initial model, b.mod. Then I run a loop to bootstrap the coefficients. For the initial model, using "predict()", I can print the wrong/false predictions table. But how do I get this table (that is, the predictions) for the bootstrapped model? Thanks for hints, *S*

df <- data.frame(
  "y"=factor(sample(LETTERS[1:3], 200, repl=T)),
  "a"=rnorm(200),
  "b"=rnorm(200)
)
library(nnet)
b.mod <- multinom(y ~ ., data=df, trace=F)
i <- 1; k <- 50
c.mat.1 <- matrix(nrow=k, ncol=length(names(df)), dimnames=list(NULL, c("cons", names(df)[-1]))) c.mat.2 <- matrix(nrow=k, ncol=length(names(df)), dimnames=list(NULL, c("cons", names(df)[-1])))
set.seed(123)
while(i <= k){
  l <- sample(1:length(df[, 1]), replace=T)
  m <- update(b.mod, . ~ ., data=df[l, ], trace=F)
  c.mat.1[i, ] <- coef(m)[1, ]
  c.mat.2[i, ] <- coef(m)[2, ]
  i <- i + 1
}
(coefs1 <- apply(c.mat.1, 2, mean))
(coefs2 <- apply(c.mat.2, 2, mean))
summary(predict(b.mod)==model.frame(b.mod)$y)

--
Sascha Vieweg, saschav...@gmail.com

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to