similar to: Robust clustered errors for probit ordinal regression analysis

If I wanted to fit a logit model and account for clustering of observations, I would do something like: library(Design) f <- lrm(Y1 ~ X1 + X2, x=TRUE, y=TRUE, data=d) g <- robcov(f, d$st.year) What would I do if I wanted to do the same thing with a probit model? ?robcov says the input model must come from the Design package, but the Design package appears not to do probit? Thanks very

A student brought this question to me and I can't find any articles or examples that are directly on point. Suppose there are 2 ordinal logistic regression models, and one wants to set them into a simultaneous equation framework. Y1 might be a 4 category scale about how much the respondent likes the American Flag and Y2 might be how much the respondent likes the Republican Party in America.

Hello, all. I'm investigating the rate at which skeletal joint surfaces pass through a series of ordered stages (changes in morphology). Current statistical methods in this type of research use various logit or probit regression techniques (e.g., proportional odds logit/probit, forward/backward continuation ratio, or restricted/unrestricted cumulative probit). Data typically include the

Hi everyone, I am running an ordinal probit using the Bayesian MCMCpack and I am getting an error saying "attempt for find suitable starting values failed" Here is my code: > posterior <- MCMCoprobit(y ~ x1 + x2 + x3 + x4 + x5 + x6 + x7 + x8 +x9 + x10 + x11 + x12 +x13 , beta.start=c(-10, 0.05, 0.02, 0.04, 0.98, 0.61, -0.29, 0.91, -0.82, 1.34, 0.68, 0.57, 0.09, 0.5), mcmc=10000)

Dear All, I have been struggling to understand why for the housing data in MASS library R and stata give coef. estimates that are really different. I also tried to come up with many many examples myself (see below, of course I did not have the set.seed command included) and all of my `random' examples seem to give verry similar output. For the housing data, I have changed the data into numeric

Hello everyone I am trying to fit an ordered probit/logit model for bank rating prediction. Besides polr() in MASS package which is not written especially for this as far as I know, do you know how else I can do this? I already found the modified polr () version on the Valentin STANESCU Enrst and Young Tel. 402 4000 ---------------------------------------------------------- The information

Hello All, I'm trying to run a multinomial probit on a dataset with 28 data points and five levels (0,1,2,3,4) in the latent choice involving response variable. I downloaded the latest mnp package to run the regression. It starts the calculation and then crashes the rpogram. I wish I could give the error message but it literally shuts down R without a warning. I'm using the R

Hi, I have a problem with the ordered probit model -polr function (library MASS). My independent variables are countinuos. I am not able to understand two main points: a) how to calculate marginal effects b) how to calculate the relative importance of each independent variables If required i will attach my model output. Thanks Franco

Hello, currently I am estimating an ordered probit model with the function polr (MASS package). Is there a simple way to obtain values for the prediction of the index function ($X*\hat{\beta}$)? (E..g. in the GLM function there is the linear.prediction value for this purpose). If not, is there another function / package where this feature is implemented? Thank you very much for

> I make a study in health econometrics and have a categorical > dependent variable (take value 1-5). I would like to fit an ordered > probit or ordered logit but i didn't find a command or package who > make that. Does anyone know if it's exists ? R is very fancy. You won't get mundane things like ordered probit off the shelf. (I will be very happy if someone will show

Dear all, I have two questions that are almost completely related to how to do things in R. I am running an ordinal probit regression analysis in R. The dependent variable has three levels (0=no action; 1=warning; 2=sanction). I use the lrm command in the rms package: print( res1<- lrm(Y ~ x1+x2+x3+x4+x5+x6, y=TRUE, x=TRUE, data=mydata)) I simply couldn't make any sense of the

Hi all, I am trying to do a ordered probit regression using polr(), replicating a result from SAS. >polr(y ~ x, dat, method='probit') suppose the model is y ~ x, where y is a factor with 3 levels and x is a factor with 5 levels, To get coefficients, SAS by default use the last level as reference, R by default use the first level (correct me if I was wrong), The result I got is a

I have been asked to analyse some questionnaire data- which is not data I'm that used to dealing with. I'm hoping that I can make use of the nabble expertise (again). The questionnaire has a section which contains a particular issue and then questions which are related to this issue (and potentially to each other): 1) importance of the issue (7 ordinal categories from -3 to +3) 2) impact

What is the current state of the R-art for ordered probit models, and more esoterically is there any available R strategy for ordered cauchit models, i.e. ordered multinomial alternatives with a cauchy link function. MCMC is an option, obviously, but for a univariate latent variable model this seems to be overkill... standard mle methods should be preferable. (??) Googling reveals that spss

hello everybody solved the problem with summary, now I have another one eg I estimate > try.op <- polr( > as.ordered(sod.sit.ec.fam) ~ > log(y) + > log(1 + nfiglimin) + > log(1 + nfiglimagg) + > log(ncomp - nfiglitot) + > eta + > I(eta^2) + >

Dear Jim, wow! It worked! Thanks a lot. I did as you suggested and it worked well with the real data. Although it gave me this error: Error in if (!is.na(x$Y[i])) { : argument is of length zero. For some reason the X1 produced less observations than it is in the data. But it's not a big deal - I identified those cases and simply deleted from the data (it was countries that only appeared

Hi - I am using polr. I can get a result from polr fit by calling result.plr <- polr(formula, data=mydata, method="probit"); However, from the 'result.plr', how can I access standard error of the estimated coefficients as well as the t statistics for each one of them? What I would like to do ultimately is to see which coefficients are not significant and try to refit the

Hi all. I'm just starting to explore ordinal multinomial regression. My dataset is 300,000 rows, with an outcome (ordinal factor from 1 to 9) and five independent variables (all continuous). My first stab at it was this: pomod <- polr(Npf ~ o_stddev + o_skewness + o_kurtosis + o_acl_1e + dispersal, rlc, Hess=TRUE) And that worked; I got a good model fit. However, a variety of other

hello everybody I'm a newbie with ordered probit and with polr too. The problem is that I have a dependent variable I need to explain with an ordered probit that is > head(dfscale$sod.sit.ec.fam,100) > [1] 5 7 5 6 5 5 6 8 6 8 8 8 6 6 6 5 0 5 NA 6 > [21] 7 NA NA 0 0 2 5 3 6 6 7 6 NA 8 6 6 7 NA NA NA > [41] 4 5 NA 10 3 9 10 10 7 5 5 5 NA

Hi all, I am using polr(). The resulting coefficients of first levels are always 0. What to do if I wnat to get the coefficients of the last level 0. For example, suppose x has 3 levels, 1, 2, 3 probit <- plor(y ~ x, data1, method='probit') will get coefficients of level 2, 3 of x, but I want coefficients of level 1, 2 Thank you, Tian [[alternative HTML version deleted]]