similar to: factors in multinom function (nnet)

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 the following table that I want to analyze using multinom function freq segments sample 4271 Seg1 tumour 4311 Seg2 tumour 3515 Seg1 normal 3561 Seg2 normal I want to compare model with both factors to the one where only sample is present. model1=multinom(freq~segments+sample,data=table) model2=multinom(freq~ sample,data=table)

Dear R users, The nnet package includes the multinom method for the confint function. The R Help file (?confint) for the generic function in the stats package and the help files for the glm and nls methods in the MASS package indicate that one can use the "parm" argument as "a specification of which parameters are to be given confidence intervals, either a vector of numbers or

Hello, I apologize if this seems like an obvious question, but I have been looking everywhere and have yet to find an answer. I am doing a multinomial regression with multinom() in the nnet package. I have a 3 level ordered response (ordered()) variable and 4 predictors, 3 of which are numerical and one which is an ordered factor (also ordered()) with 5 levels (a, b, c, d, e). My question is in

Dear R users, Does anyone know if it is possible to use multinom to do a polychotomous fit using one categorical and one numeric variable as response. The doc. for multinom states that for formula , response can be K>2 classes. Is this 2 and more, or as I have understood it only greater than 2. I have tried fitting my data, but have only encountered error messages. On another note, Is it

Is anyone aware of a way to check whether multinom has converged by checking a component of the output ? I'm not familar with nnet but, since multinom calls nnet , maybe there is an extra argument once can send to multinom to capture this information. Thanks.

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

Hi, I am trying to specify a multinomial logit model using the multinom function from the nnet package. Now I add another independent variable and it halves the AIC as given by summary(multinom()). But when I call Anova(multinom()) from the car package, it tells me that this added variable is insignificant (Pr(>Chisq)=0.39). Thus, the improved AIC suggests to keep the variable but the Anova

Hi, I found that using different contrasts (e.g. contr.helmert vs. contr.treatment) will generate different fitted probabilities from multinomial logistic regression using multinom(); while the fitted probabilities from binary logistic regression seem to be the same. Why is that? and for multinomial logisitc regression, what contrast should be used? I guess it's helmert? here is an example

Hi, I am having difficulty entering a 'programmable' argument into the multinom function from the nnet package. Interactively, I can get the function to work fine by calling it this way: z1=multinom(formula = class.ind(grp[-outgroup])~ (PC1 + PC2 + PC3), data=data.frame(scores)) However I need to be able to change the number of variables I am looking for in 'scores' and so am

Full_Name: Rob Foxall Version: 2.2.0 OS: Windows XP Submission from: (NULL) (149.155.96.5) I was using multinom from nnet package, when I did something stupid -- I entered in an incorrect factor variable as response. This factor had only one level. Instead of R telling me not to be so dumb, it crashed, clicking on debug coming up with the message "An exception 'Unhandled Win32

Best R Help, I like to estimate a Multinomial Logit Model with 10 Classes. The problem is that the number of observations differs a lot over the 10 classes: Class | num. Observations A | 373 B | 631 C | 171 D | 700 E | 87 F | 249 G | 138 H | 133 I | 162 J | 407 Total: 3051 Where my data looks like: x1 x2 x3 x4 Class 1 1,02 2 1 A 2 7,2 1 5 B 3 4,2 1 4 H 1 4,1 1 8 F 2 2,4 3 7 D 1 1,2 0 4 J 2 0,9

Hi: I'm working with multinomial models with library nnet, and I'm trying to get the explained deviance (pseudo R^2) of my models. I am assuming that: pseudo R^2= 1 - dev(model) / dev (null) where dev(model) is the deviance for the fitted model and dev(null) is the deviance for the null model (with the intercept only). library(nnet) full.model<- multinom(cbind(factor1,

Running a binary logit model on the data df <- data.frame(y=sample(letters[1:3], 100, repl=T), x=rnorm(100)) reveals some residual deviance: summary(glm(y ~ ., data=df, family=binomial("logit"))) However, running a multinomial model on that data (multinom, nnet) reveals a residual deviance: summary(multinom(y ~ ., data=df)) On page 203, the MASS book says that "here the

DeaR useRs: I am looking for a function which fits a multinomial model and in Baron?s page I find the function "multinom" in package "nnet" but this package is deprecated. I suppose that this function is now in other package but I can't find it. Can you help me? Thanks.

I would like to perform a multinomial logistic regression on a large data set, but do not know how. I've only thought of a few possibilities and write to seek advice and guidance on them or deepening or expanding my search. On smaller data sets, I have successfully loaded the data and issued commands such as: length(levels(factor(data$response))) [1] 6 # implies polychotomy library(nnet)

Hello folks, Let me first apologize: I'm not a professional nor a mathematician, just an ordinary guy, fooling around with the excellent R-package. I know the basic principles behind statistics, but haven't read anything more advanced than the ordinary first probability and statistics courses. Enough disclaimers? Good! I was examining the multinom-function (in the nnet-package) the other

I see that "case weights" can be optioned in multinom(). I wanted to make sure I understand what weights= is expecting. My weights (not really mine but I'm stuck with them) are noninteger, are not scaled to sum to the sample size, and larger weights are intended to increase influence. The description of various types of weights is a perennial confusion for me; sorry. STS

Good Day: When I used: multinom(formula = Y ~ X1 + X2 + X3 + X1:X2 + X1:X3 + X3:X2 + X1^2 + X2^2 + X3^2, data = DATASET), I get estimates and AIC for the model containing main effects and interactions only (no squared terms)...and FYI, all predictors are continuous. Is this "normal" behavior? If I run this in S-Plus I get estimates and AIC for the model containing all terms(including

I have RedHat linux 9 with R 1.8. I'm estimating models with multinom with a dependent variable that has 3 different values. Sometimes the models run fine and I can understand the results. Sometimes when I put in another variable, I see an indication that the estimation did work, but then I can't get the summary method to work. It's like this: > votemn1 <-