Displaying 1 result from an estimated 1 matches for "prdata".
Did you mean:
pdata
2011 Oct 19
1
hypothetical prediction after polr
...###########
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="probs")
################################################################
I do not have any problem estimating the model, but when it comes to
the last time (predict), I got the fol...