similar to: p value for coefficients in multinomial model

Dear all, I search the mail list about this topic and learn that no simple way is available to get "lsmeans" in R as in SAS. Dr.John Fox and Dr.Frank E Harrell have given very useful information about "lsmeans" topic. Dr. Frank E Harrell suggests not to think about lsmeans, just to think about what predicted values wanted and to use the predict

Dear all, When I use aggregate function as: attach(warpbreaks) aggregate(warpbreaks[, 1], list(wool = wool, tension = tension), sum) The results are right but I get a warning message: "number of items to replace is not a multiple of replacement length." BTW: I use R version 2.4.1 in Ubuntu 7.04. Your kind solutions will be great appreciated. Best wishes Yours, sincerely, Xingwang

The partial R-square (or coefficient of partial determination, or squared partial correlation coefficients) measures the marginal contribution of one explanatory variable when all others are already included in multiple linear regression model. The following link has very clear explanations on partial and semi-partial correlation: http://www.psy.jhu.edu/~ashelton/courses/stats315/week2.pdf In

{again a message that was sent to owner-r-help (which is me, currently) why on earth ???!??!? reply to R-help or the original sender Brett Presnell; } I'm teaching a categorical data analysis course this term, and a minor "problem" has resurfaced that I have often thought about before. This applies equally to Splus I suppose, but my undergrads aren't using Splus. It

Dear All, Could we specify the multinomial distribution in R when doing Bayesian data analysis through R2WinBUGS? See below for my issues. t2[i,1:2] ~ dmulti(p2[i,1:2],n2[i]) #t2[] is a matrix with two columns, the problem is how to specify it in R p2[i,1] <- (p[i] * s1[t[i]] * s2[t[i]] + (1 - p[i]) * (1 - c1[t[i]]) * (1 - c2[t[i]]))/p1[i] p2[i,2] <- (p[i] * s1[t[i]] * (1 - s2[t[i]]) +(1 -

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 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

Good morning, I'd like to estimate a simple multinomial logit model in R (not a McFadden conditional logit). For instance, I'd like to estimate the probability of someone having one of eight titles in a company with the independent variables being the company characteristics. A binary logit is well documented. What about the multinomial? Thanks, Walt Paczkowski

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

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]]

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

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 __________________________________

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.

I have had plenty of succes generating one dimensional variables and plotting them, but what do i do for more (specifically 2) dimensional multinomial variables? I figure i have to create a vector consisting of two 1 dim normallly distributed variables, that way i can also control the correlation. But how do i do this :) ? Many thanks -- View this message in context:

Hi everyone? Does anybody know how can I estimate a Latent class multinomial (or conditional) logit using R? I have tried flexmix, poLCA, and they do not seem to support this model. thanks in advance adan -- View this message in context: http://r.789695.n4.nabble.com/Latent-class-multinomial-or-conditional-logit-using-R-tp4230083p4230083.html Sent from the R help mailing list archive at

Hi, I am trying to obtain p-values for coefficient estimates in a multinomial logit model. Although I am able to test for significance using other methods (e.g., Wald statistics), I can't seem to get R to give me simple p-values. I am sure there is a very simple solution to this, but the R archives seem to have nothing on this issue. I would appreciate any help. Thanks in advance! Best,

I have to the treatment effect base on the observational data.And the treatment variable is multinomial rather than binary.Because the treatment assignment is not random,so the selection-bias exists.Under this condition,what's the best way to estimate the treatment effect? I know that if the treatment is binary,I can use propensity score matching using MatchIt package.But what about

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

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

I'm interested in doing multinomial logistic regression with equality constraints on some of the parameter values. For example, with categorical outcomes Y_1 (baseline), Y_2, and Y_3, and covariates X_1 and X_2, I might want to impose the equality constraint that \beta_{2,1} = \beta_{3,2} that is, that the effect of X_1 on the logit of Y_2 is the same as the effect of X_2 on the