Displaying 2 results from an estimated 2 matches for "yr89".
Did you mean:
r89
2011 Oct 20
1
effect function in the effects package
...fect plot after the proportional
odds logistic regression model. 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))))
##############################################################################...
2011 Oct 19
1
hypothetical prediction after polr
...somehow this does not
work. I guess I must have missed something 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,...