similar to: ordered probit and cauchit

Dear all, A am trying to fit a generalized linear mixed effects model with a binomial link function, my response data is binary, using the lme4 R package, for the glmer model but with the cauchit link function (CDF of Cauchy distribution), under the package this has not yet been coded and was wondering if anyone knew a way in which I could incorporate this link function into the code. Thankyou

> 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

Full_Name: Henric Nilsson Version: 2.3.0 Patched (2006-05-09 r38014) OS: Windows 2000 SP4 Submission from: (NULL) (83.253.9.137) When supplying an unavailable link to `quasibinomial', the error message looks strange. E.g. > quasibinomial("x") Error in quasibinomial("x") : 'x' link not available for quasibinomial family, available links are "logit",

Just a quick point which may be easy to correct. Whilst typing the wrong thing into R 2.2.1, I noticed the following error messages, which seem to have some stray quotation marks and commas in the list of available families. Perhaps they have been corrected in the latest version (sorry, I don't want to upgrade yet, but it should be easy to check)? > glm(1 ~ 2,

Hello, The following code: > binomial(identity) Generates an error message: Error in binomial(identity) : link "identity" not available for binomial family; available links are ?logit?, ?probit?, ?cloglog?, ?cauchit?, ?log? While : > binomial("identity") Yields an identity-binomial object that works as expected with stats::glm The error in the first example mislead

this may be a better question for r-devel, but ... Is there a particular reason (and if so, what is it) that the inverse link is not in the list of allowable link functions for the binomial family? I initially thought this might have something to do with the properties of canonical vs non-canonical link functions, but since other link functions (probit, cloglog, cauchit, log) are allowed, I

Help needed with lm function: Dear R's, Could anyone tell me how to replace the link function (probit logit, loglog etc.) in lm with an abitrary user-defined function? The task is to perform ML Estimation of betas for a dichotome target variable. Maybe there is already a package for this (I did not find one). Any hints or a code excerpt would be welcome! Thank you -Jeff jeff.pr2 (at)

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

Hi, I'm getting an error message using polr(): Error in optim(start, fmin, gmin, method = "BFGS", hessian = Hess, ...) : initial value in 'vmmin' is not finite The outcome variable is ordinal and factored, and the independant variable is continuous. I've checked the source code for both polr() and optim() and can't find any variable called

Hello, I'm working on the multiple categorical data (5-points scale) using linear mixed effect model and wondering if anyone knows about or works on the linear mixed effect model with ordered logit or probit link. I found that the "lmer" function in R is very flexible and supports various models, but not ordered logit/probit models. I may conduct my analysis by turning my DVs

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

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

Hello, I searched everywhere but I didn't find what I want, that is why I as the question here. Threads discussing this issue on this mailing list are already quite old. Does anybody know of a function in R which allows to estimate ordered probit/logit model with random coefficients. The only mixed effect model I found was clmm of the ordinal package but it only provides random intercepts. I

Hi, I am new in R. I would like to do a ordered probit regression using clm2 (in the ordinal package). My dependent variable y is the way of payment in M&A: y=0 if the deal is financed by stock only, y=1 if the deal is financed by a mix of cash and stock and y=2 if it is by cash only. My independent variables are CollateralB, Cashavailable and Leverage. This is the code I wrote: >

Hello, First time poster here so let me know if you need any more information. I am trying to run an ordered probit with maximum likelihood model in R with a very simple model (model <- econ3 ~ partyid). Everything looks ok until i try to run the optim() command and that's when I get " Error in rep(1, nrow(x)) : invalid 'times' argument". I had to adapt the code from a 4

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

Dear all, I have a very simple set of data and I would like to analyze them with probit analysis. dose event trial 0.0 3 15 1.1 4 15 1.3 4 15 2.0 3 15 2.2 5 15 2.8 4 15 3.7 5 15 3.9 9 15 4.4 8 15 4.8 11 15 5.9 12 15 6.8 13 15 The dose should be transformed with log10(). I use glm(y ~ log10(dose), family=binomial(link=probit)) to do probit analysis, however, I have to exclude the

R users, I got a problem to analyze with probit model. What package contains the algorithm to do probit model. Lawrence N.M Kazembe Mathematical Sciences Department Chancellor College University of Malawi P.O. Box 280 Zomba Malawi Tel: (265) 524 222 ext 284 Fax: (265) 524 046 e-mail: lkazembe at chirunga.sdnp.org.mw url: kazembe.cjb.net kazembe.tsx.org

Not real data. It was gererated randomly. The original codes are the following: par(mfrow=c(2,1)) n <- 500 ######################### #DATA GENERATING PROCESS# ######################### x1 <- rnorm(n,0,1) x2 <- rchisq(n,df=3,ncp=0)-3 sigma <- 1 u1 <- rnorm(n,0,sigma) ylatent1 <-x1+x2+u1 y1 <- (ylatent1 >=0) # create the binary indicator ####################### #THE

Hi Folks, I am interested in extending the repertoire of link functions in glm(Y~X, family=binomial(link=...)) to include a "tan" link: eta = (4/pi)*tan(mu) i.e. this link bears the same relation to the Cauchy distribution as the probit link bears to the Gaussian. I'm interested in sage advice about this from people who know their way aroung glm. >From the surface, it looks