similar to: Likelihood Function for Multinomial Logistic Regression and its partial derivatives

Hi there I am trying to learn how to compute mle in R for a multinomial negative log likelihood function. I am using for this the book by B. Bolker "Ecological models and data in R", chapter 6: "Likelihood an all that". But he has no example for multinomial functions. What I did is the following: I first defined a function for the negative log likelihood:

R Experts: I am conducting a meta-analysis where the effect measures to be pooled are simple proportions. For example, consider this data from Fleiss/Levin/Paik's Statistical methods for rates and proportions (2003, p189) on smokers: Study N Event P(Event) 1 86 83 0.965 2 93 90 0.968 3 136 129 0.949 4 82 70 0.854 Total

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 R-group, Could somebody recommend a package that can deal with a multinomial response variable (choice of breeding tactic in mice, which has four unordered levels), multiply-imputed data (generated using the Amelia package) and two non-nested random effects: individual identity (133 individuals made up to four choices each) and year (for which there are six levels and sample size varies

Hello, I''ve spent a lot of time during the past month trying to get negative multinomial regression models for clustered event counts as described in (Guang Guo. 1996. "Negative Multinomial Regression Models For Clustered Event Counts." Sociological Methodology 26: 113-132., abstract at http://depts.washington.edu/socmeth2/4abst96.htm) implemented in R. A FORTRAN version of the

I have two questions: 1. If I want to do a binomial logit, how to handle the categorical response variable? Data for the response variables are not numerical, but text. 2. What if I want to do a multinomial logit, still with categorical response variable? The variable has 5 non-numerical response levels, I have to do it with a multinomial logit. Any input is highly appreciated! Thanks! Ed

I want to get the logLik to calculate McFadden.R2 ,ML.R2 and Cragg.Uhler.R2, but the value from multinom does not have logLik.So my quetion is : is logLik meaningful to multinomial logistic model from multinom?If it does, how can I get it? Thank you! ps: I konw VGAM has function to get the multinomial logistic model with logLik, but I prefer use the function from "official" R

Dear all, I have found some difference of the results between multinom() function in R and multinomial logistic regression in SPSS software. The input data, model and parameters are below: choles <- c(94, 158, 133, 164, 162, 182, 140, 157, 146, 182); sbp <- c(105, 121, 128, 149, 132, 103, 97, 128, 114, 129); case <- c(1, 3, 3, 2, 1, 2, 3, 1, 2, 2); result <- multinom(case ~ choles

I tried to use the effect() to get predicted marginals for multinomial logistic as I did for general logistic regression, but failed. Is there anyway to do that? Thx! -- View this message in context: http://www.nabble.com/How-to-get-predicted-marginal-%28aka-predicted-mean%29-after-multinomial-logistic--tp17200114p17200114.html Sent from the R help mailing list archive at Nabble.com.

Hi, how can I do multinomial logistic regression in R? I think glm() can only handle binary response variable, and polr() can only handle ordinal response variable. how to do logistic regression with multinomial response variable? Thanks __________________________________

Dear helpers, I am using R version 2.5.1 to estimate a multinomial logit model using my own maximum likelihood function (I work with share data and the default function of R cannot deal with that). However, the computer (I have an Athlon XP 3200+ with 512 GB ram) takes quite a while to estimate the model. With 3 categories, 5 explanatory variables and roughly 5000 observations it takes 2-3 min.

Sir, I want to fit a multinomial logistic regression in R.I think mlogit() is the function for doing this. mlogit () is in packege globaltest.But, I can not install this package. I use the following: install.packages("globaltest") Can you help me? Regards, Suman Dhara [[alternative HTML version deleted]]

Hello, after calculating a multinomial logit regression on my data, I compared the output to an output retrieved with SPSS 18 (Mac). The coefficients appear to be the same, but the logLik (and therefore fit) values differ widely. Why? The regression in R: set.seed(1234) df <- data.frame( "y"=factor(sample(LETTERS[1:3], 143, repl=T, prob=c(4, 1, 10))), "a"=sample(1:5,

Dear Sirs, Could you please be so kind as to send us some information on residuals in multinomial logistic models? Is it possible to use R software? We thank you in advance. Sincerely yours Luciana Alves,MSc Beatriz Leimann, MD -- Luciana Correia Alves Doutoranda em Sa??de P??blica ENSP - Fiocruz

I am attempted to examine the temporal independence of my data set and think I need an unordered multinomial logistic regression (or logit model) with repeated measures to do so. The data in question is location of chickens. Chickens could be in any one of 5 locations when a snapshot sample was taken. The locations of chickens (bird) in 8 pens (pen) were scored twice a day (AMPM) for 20 days

Dear All I am trying to fit a multinomial logistic regression to a data set with a size of 94279 by 14 entries. The data frame has one "sample" column which is the categorical variable, and the number of different categories is 9. The size of the data set (as a csv file) is less than 10 MB. I tried to fit a multinomial logistic regression, either using vglm() from the VGAM package or

Thanks for your answers Stephen and Ben, I hope I am posting on the correct list now. I managed so far to run the multinomial model with random effect with the following command: MCMCglmm(fixed=cbind(Apsy,Mygl,Crle,Crru,Miag,empty) ~ habitat:trait,random=~idh(trait):mesh,family="multinomial12", data=dataA,rcov=~trait:units) (where multiple responses are different species, Habitat

I am looking into using the multinomial logistic regression option in the nnet library and have two questions about formatting the data. 1. Can data be analysed in the following format or does it need to be transformed into count data, such as the housing data in MASS? Id Crime paranoia hallucinate toc disorg crimhist age 1 2 1 0 1 0 1 25 2 2 0 1 1 1 1 37 3 1 1 0 1 1 0 42 4 3 0

On May 24, 2011; 11:06pm Belle wrote: > Does anyone know how to run Multinomial logistical Model in R in order to > get predicted probability? Yes. I could stop there but you shouldn't. The author of the package provides plenty of examples (and two good vignettes) showing you how to do this. Suggest you do some work in that area. Look especially at how model formulas are

Hi all, I'm working on a multinomial (or "polytomous") logistic regression using R and have made great progress using multinom() from the nnet library. My response variable has three categories, and there are two different possible predictors. I'd like to use the likelihoods of certain models (ie, saturated, fitteds, and null) to calculate Nagelkerke R-squared values for