similar to: Logistic Regression with 9 classes

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

Hello all. I am doing one part of an evaluation of a mandatory welfare-to-work programme in the UK. As with all evaluations, the problem is to determine what would have happened if the initiative had not taken place. In our case, we have a number of pilot areas and no possibility of random assignment. Therefore we have been given control areas. My problem is to select for survey individuals in

Hi! I'm trying to do some ploytomous logistic regression using multinom() in the nnet package, but am a bit confused about interpretation of the results Is it possible to get the following quantities: I: maximum likelihood estimates to test for fit of model and significance of each predictor (I would like to produce a table of the following type) Analysis of Variance: MLE (values 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

Dear all, I'm interested in extracting the score from multinomial logistic regression models fit using multinom, to assess the stregth of assocation of the parameter with the response (akin to the score from clogit/cox regression). currently I'm using R 1.8.1. Is there a function that will extract the score from a multinom object or how i can get back to it? or from using glm? I

Dear all R-users, I am a new user of R and I am trying to build a discrete choice model (with more than two alternatives A, B, C and D) using logistic regression. I have data that describes the observed choice probabilities and some background information. An example below describes the data: Sex Age pr(A) pr(B) pr(C) pr(D) ... 1 11 0.5 0.5 0 0 1 40 1 0 0 0 0 34 0 0 0 1 0 64 0.1 0.5 0.2 0.2 ...

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

? Hallo, which function can i use to do (baseline) logistic regression + goodness of fit tests? so far i found: # logistic on binary data lrm combined with resid(model,'gof') # logistic on binary data glm with no gof-test # baseline logit on binary data

Dear R helpers! Is there an utility for multiple logistic regression where the response is not ordinal? As the predictors do not satisfy the necessary hypothesis, I cannot use discriminant analysis. Thanks and have a nice R-day Anne [[alternative HTML version deleted]]

Hello, I am using logistic discriminant analysis to check whether a known classification Yobs can be predicted by few continuous variables X. What I do is to predict class probabilities with multinom() in nnet(), obtaining a predicted classification Ypred and then compute the percentage P(obs) of objects classified the same in Yobs and Ypred. My problem now is to figure out whether P(obs) is

Is it possible to use function "glm" in case when my outcome variable has 5 different classes? I have seen examples only when using binomial outcome variable. What about using function "multinom"? How do I to get the signifigance and the confidence levels of the coefficients and the value of goodness of the model with this function? Thank You for Your help! -- View this

Hi, Does anybody know some package or function that computes canonical variates? Luis & Janete -------------------------------------------- sapo.pt/kitadsl -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

Dear R-users, Is there a function in R that classifies data in more than 2 groups using logistic regression/classification? I want to compare the c-indices of earlier research (lrm, binary response variables) with new c-indices obtained from 'multiple' (more response variables) logistic regression. Best regards, Stephanie Delalieux Department Biosystems M?-BIORES Group of Geomatics

Dear Colleagues I have been fitting some multinomial logistic regression models using R (version 1.6.1 on a linux box) and Stata 7. Although the vast majority of the parameter estimates and standard errors I get from R are the same as those from Stata (given rounding errors and so on), there are a few estimates for the same model which are quite different. I would be most grateful if

I was hoping for some advice regarding possible explanations for the fitted probability values I obtained for a multinomial logistic regression. The analysis aims to predict whether Capgras delusions (present/absent) are associated with group (ABH, SV, homicide; values = 1,2,3,), controlling for previous violence. What has me puzzled is that for each combination the fitted probabilities are

Dear All, I am trying to reproduce the example that I found online here http://bit.ly/11VG4ha However, when I run my script (pasted at the end of the email), I notice that there is a factor 2 between the values for the coefficients for the categorical variable female calculated by my script and in the online example. Any idea about where this difference comes from? Besides, how can I

Christoph Lehman had problems with seperated data in two-class logistic regression. One useful little trick is to penalize the logistic regression using a quadratic penalty on the coefficients. I am sure there are functions in the R contributed libraries to do this; otherwise it is easy to achieve via IRLS using ridge regressions. Then even though the data are separated, the penalized

Hello, Does logistic regression really provide better results than lda or qda ? (my purpose is not classification but highlighting of discriminant variables). If this is true, where could I get an R script performing stepwise logistic regression ? Thanks -- Daniel AMORESE Lab. M2C "Morphodynamique Continentale et C?ti?re" UMR CNRS 6143 Caen/Rouen Centre de G?omorphologie UCBN

Hello! I need to test the adjustment of a (Negative) Exponential Distribution to a dataset. The parameter of the distribution is unknown. What is the appropriate test to do? I've tried the ks.test, although I think this isn't the appropriate one, as I don't know the parameter. Can anybody help me? Thanks in advance, Janete -- Janete da Silva Borges janeteborges at gmx.net Ab