如何用R软件做多类别Logistic回归

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

要求用R 软件做多类别logistic回归（y有值：1，2，3，4）
我只会做二项的binomal的，求多类别的，谢谢大家

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：logistic回归 logistic ogistic logisti logist 如何 软件

1. # Logistic Regression
   
2. # where F is a binary factor and
   
3. # x1-x3 are continuous predictors
   
4. fit <- glm(F~x1+x2+x3,data=mydata,family=binomial())
   
5. summary(fit) # display results
   
6. confint(fit) # 95% CI for the coefficients
   
7. exp(coef(fit)) # exponentiated coefficients
   
8. exp(confint(fit)) # 95% CI for exponentiated coefficients
   
9. predict(fit, type="response") # predicted values
   
10. residuals(fit, type="deviance") # residuals

复制代码

ryusukekenji 发表于 2013-6-17 00:02

谢谢你的代码，不过你做的也是y取0,1两个值，可是我们的数据里y有1,2,3,4,四个值，这个该怎么做啊

library(nnet)
multinom函数
详见MASS 7.3节
也可以试试library(mlogit)

hugebear 发表于 2013-6-19 09:39
library(nnet)
multinom函数
详见MASS 7.3节

谢谢你，这个函数我也听说了，不过惭愧的是，不太会使，分析出来的结果看太懂啊

找不到mlogit这个函数呀

谢谢！！！！找了好久！！！

multinom可以做多类别的啊

  result <- multinom(Y ~ X+age, data = tdata)
  z <- summary(result) $ coefficients/summary(result) $ standard.errors
  p <- (1-pnorm(abs(z),0,1))*2
  out <- data.frame(cbind(names(tdata)[q], t(summary(result) $ coefficients[,2]), t(p[,2])))
  timo <- rbind(timo,out)