Displaying 2 results from an estimated 2 matches for "ordwarm2".
2011 Oct 19
1
hypothetical prediction after polr
...ething here. I first used the polr
function in the MASS package, and I create a data frame and supply it
to the predict function (see below):
###############################################################
myologit <- polr(factor(warm) ~ yr89 + male + white + age + ed + prst,
data=ordwarm2, method=c("logistic"))
yr89 <- c(1)
male <- c(1)
white <- c(1)
age <- c(mean(ordwarm2$age))
ed <- c(mean(ordwarm2$ed))
prst <- c(mean(ordwarm2$prst))
prdata <- data.frame(yr89, male, white, age, ed, prst)
prdata$rankR <-predict(myologit,newdata=prdata,type...
2011 Oct 20
1
effect function in the effects package
...del. There is no problem for me to estimate
the model, but when it comes to the graphing, I was stuck. see the
codes below:
##############################################################################
myologit <- polr(factor(warm) ~ yr89 + male + white + age + ed + prst,
+ data=ordwarm2, method=c("logistic"))
require(effects)
plot(effect("age", myologit, xlevles=list(age=seq(20, 80, 5),
given.values(male=1, yr89=1))))
##############################################################################
I got the following error message:
Warning message:
package ...