search for: yr89

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,...