similar to: logit link + alternatives

I am trying to figure out how to apply a loglog link to a binomial model (dichotomous response variable with far more zeros than ones). I am aware that there are several relevant posts on this list, but I am afraid I need a little more help. The two suggested approaches seem to be: 1) modify the make.link function in GLM, or 2) use the loglog or cloglog functions in the VGAM package.

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

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

> 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

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

if I have the following in my routes.rb: map.resources :product_categories it provides urls such as: /product_categories /product_categories/1 /product_categories/1;edit etc. and uses the ProductCategories controller. how can I keep the same controller (and model) but set up the resources so that I can have any abitrary url point to the existing controller, much like specifying the controller

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

Dear all, can somebody please help me how to calculate "McFadden R-square" for a LOGIT model? Corresponding definition can be found here: http://publib.boulder.ibm.com/infocenter/spssstat/v20r0m0/index.jsp?topic=%2Fcom.ibm.spss.statistics.help%2Falg_plum_statistics_rsq_mcfadden.htm Here is my data: Data <- structure(c(1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1,

Dear R-Help, I'm hoping to find a Multinomial Nested Logit package in R. It would be great to find something analogous to "PROC MDC" in SAS: > The MDC (Multinomial Discrete Choice) procedure analyzes models > where the > choice set consists of multiple alternatives. This procedure > supports conditional logit, > mixed logit, heteroscedastic extreme value,

Dear everybody! Are there algorithms for multinominal logit/probit available for R? Is it my fault that I cannot find these in CRAN? Has somebody programmed these? with best wishes Ott Toomet Ott.Toomet at mail.ee -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

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

Does anybody know how to do multiple logit/probit analysis with R? Thanks in advance!

Hi Folks, I wonderif someone who is familiar with the details of vglm in the VGAM package can assist me. I'm new to using it, and there doesn;t seem much in the documentation that's relevant to the question below. Say I have a vector x of 0/1 responses and another vector y of 0/1 responses, these in fact being a bivariate set of 0/1 responses equivalent to cbind(x,y). E.g.

I would like to do an R glm() with family = binomial(link="loglog") Right now, the cloglog link exists, which is nice when the data have a heavy tail to the left. I have the opposite case and the loglog link is what I need. Can someone suggest how to add the loglog link onto glm()? It would be lovely to have it there by default, and it certainly makes sense to have the two opposite

Hi, I've tried two SIP clients so far and both have unusable outgoing audio quality. Skype app sounds fine, and recording the same mic sounds fine, so I can only assume there is an issue with the clients themselves. Both clients allow you to register and make calls via SIP with any abitrary provider and credentials, so they'll work with Asterisk. I've tried them with two good

Dear List I am trying to fit a simutaneous equation logit model. i.e., the response variables of the structured equations are binomial, I am not sure if systemfit can do this job. A google search doesn't yield too much helpful information. Your knowledge on any other packages or codes are appreciated. Thanks will

Hello all, I am wondering if there is a way to do a spatial error probit/logit model in R? I can't seem to find it in any of the packages. I can do it in MATLAB with Gibbs sampling, but would like to confirm the results. Ideally I would like to use this model to predict probability of parcel conversion in a future time period. This seems especially difficult in a binary outcome model

Hello! One note/question hier about specification of control-parameters in the lme(...,control=list(...)) function call: i tried to specify tne number of iteration needed via lme(....,control=list(maxIter=..., niterEM=...,msVerbose=TRUE)) but every time i change the defualt values maxIter (e.g. maxIter=1, niterEM=0) on ones specified by me, the call returns all the iterations needed until

Hi, I have recently been attempting to find the LD50 from two predicted fits (For male and females) in a Generalised linear model which models the effect of both sex + logdose (and sex*logdose interaction) on proportion survival (formula = y ~ ldose * sex, family = "binomial", data = dat (y is the survival data)). I can obtain the LD50 for females using the dose.p() command in the MASS

Hi! I was trying to implement a probit model on a dichotomous outcome variable and found that the predictions were outside the (0,1) interval that one should get. I later tried it with some simulated data with a similar result. Here is a toy program I wrote and I cant figure why I should be getting such odd predictions. x1<-rnorm(1000) x2<-rnorm(1000) x3<-rnorm(1000)