这是一个更部分的实际答案,但在深入理论之前做一些练习对我有用。
这个ats.ucla.edu链接是一个参考,可能有助于以更实用的方式开始了解多项逻辑回归(如 Bill 所指出的)。
它提供了可重现的代码来理解包中的功能multinom,并简要介绍了输出解释。nmetR
考虑这段代码:
va = c('cat','dog','dog','goat','cat','goat','dog','dog')
# cat will be the outcome baseline
vb = c(1,2,1,2,1,2,1,2)
vc = c('blue','red','blue','red','red','blue','yellow','yellow')
# blue will be the vc predictor baseline
set.seed(12)
vd = round(rnorm(8),2)
data = data.frame(cbind(va,vb,vc,vd))
library(nnet)
fit <- multinom(va ~ as.numeric(vb) + vc + as.numeric(vd), data=data)
# weights: 18 (10 variable)
initial value 8.788898
iter 10 value 0.213098
iter 20 value 0.000278
final value 0.000070
converged
fit
Call:
multinom(formula = va ~ as.numeric(vb) + vc + as.numeric(vd),
data = data)
Coefficients:
(Intercept) as.numeric(vb) vcred vcyellow as.numeric(vd)
dog -1.044866 120.3495 -6.705314 77.41661 -21.97069
goat 47.493155 126.4840 49.856414 -41.46955 -47.72585
Residual Deviance: 0.0001656705
AIC: 20.00017
这是您可以解释对数线性拟合多项逻辑模型的方式:
ln(P(va=cat)P(va=dog))ln(P(va=cat)P(va=goat))=b10+b11vb+b12(vc=red)+b13(vc=yellow)+b14vd =b20+b21vb+b22(vc=red)+b23(vc=yellow)+b24vd
以下是有关如何解释模型参数的摘录:
- 变量 vd 增加一个单位与“狗”对“猫”的对数几率减少 21.97069 (b14)。
第二行的逻辑相同,但考虑到“山羊”与“猫”与 (b24=-47.72585)。
- 如果从 vc="blue" 移动到 vc="red"(b12)。
......
文章还有很多,但我认为这部分是核心。
参考:
R 数据分析示例:多项 Logistic 回归。加州大学洛杉矶分校:统计咨询集团。
来自http://www.ats.ucla.edu/stat/r/dae/mlogit.htm(2013年 11 月 5 日访问)。