用序数自变量估计概率的R包？答案

【问题标题】：R package for estimating probit with ordinal independent variables?用序数自变量估计概率的R包？
【发布时间】：2010-07-20 00:00:49
【问题描述】：

我希望估计一个回归模型，其中因变量是一个虚拟变量（编码为 0/1），我有五个或六个序数自变量（我将虚拟出来），再加上一堆其他的东西。任何人都可以推荐一个可以以最少的麻烦或以其他方式处理序数 RHS 变量的软件包吗？谢谢

【问题讨论】：

标签： r

【解决方案1】：

您可以使用内置的glm 函数来完成这一切，此外还可以在公式中应设为虚拟变量的变量周围适当使用factor。

例子：

R> y <- rbinom(100, 1, .5)
R> x1 <- sample(1:5, 100, replace = TRUE)
R> x2 <- sample(1:5, 100, replace = TRUE)
R> m1 <- glm(y ~ factor(x1) + factor(x2), family = binomial(link = "probit"))
R> m1

Call:  glm(formula = y ~ factor(x1) + factor(x2), family = binomial(link = "probit")) 

Coefficients:
(Intercept)  factor(x1)2  factor(x1)3  factor(x1)4  factor(x1)5  factor(x2)2  
      0.335       -0.729       -0.670       -0.639       -0.740        0.327  
factor(x2)3  factor(x2)4  factor(x2)5  
     -0.106        0.624        0.483  

Degrees of Freedom: 99 Total (i.e. Null);  91 Residual
Null Deviance:      138 
Residual Deviance: 129  AIC: 147

您可能还想看看dummies 包。

【讨论】：