Displaying 1 result from an estimated 1 matches for "origstepaic".
2011 Nov 29
0
Any function\method to use automatically Final Model after bootstrapping using boot.stepAIC()
...a, B = 50)
bglmfit
In the summary of Bootstrapping the 'stepAIC()' procedure, Following
information is listed:
Initial Model:
y1 ~ x1 + x2 + x3 + x4 + x5
Final Model:
y1 ~ x1 + x5
Is there any function or method for using Final Model by Bootstrapping the
'stepAIC()' procedure, like OrigstepAIC model as shown below.
n <- 200
x1 <- runif(n, -3, 3)
x5 <- factor(sample(letters[1:2], n, rep = TRUE))
eta <- 0.1 + 1.6 * x1 - 2.5 * as.numeric(as.character(x5) == levels(x5)[1])
data1 <- data.frame(x1, x5)
data1$probscore <- predict(bglmfit$OrigStepAIC , data1)
Is there any way...