I think the most common reason to see different parameter estimates with ordinal regression is that programs set up the model differently.
For example, check out library(MASS) ?polr We see polr uses: logit P(Y <= k | x) = zeta_k - eta and notes that other software packages may use the opposite sign for eta. Also check out that the ordered variables have the same order (reference category). I sometimes have messed that up. On Mon, Sep 27, 2010 at 2:53 PM, Ben Hunter <bjameshun...@gmail.com> wrote: > I am learning R via a textbook that performs analysis with SPSS and SAS. In > trying to reproduce the results for an ordinal logit model, I get very > similar point estimates for my cut-off points, but the parameters for the > covariate q60 do not match. The estimate for q51 also matches. Is this > because I need to change a base case for the ordered covariate q60? Can this > be done in or is it always the first case? > > mod<-polr(as.ordered(q43j)~as.ordered(q60)+q51, method="logistic") > > Perhaps a book using R would be a better idea, but it's the content and > price (free) of this book that I like. > > Thanks a lot, > > Ben > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.