使用 mlogit 构建错误组件 logit

Question

我想使用 R 的 mlogit 库构建一个错误组件 logit。

我已将我的数据集视为面板数据集（即每一行表示一个替代项），然后构建一个错误组件 logit 模型。

虽然我知道为了构建混合 logit 模型，我需要在 rpar 命令中添加协变量列表。但是，我不想估计协变量的随机参数，而是估计截距项。

Answer 1

在多项式 logit 模型中，您可以估计 J-1 替代特定常数（截距）。使它们随机化的最简单方法是创建替代的特定指标。

例如，假设您有三个备选方案：1、2 和 3，并且这些备选方案存储在变量 alt 中。现在您可以创建两个名为alt_1 和alt_2 的新变量，对于备选方案1 和2，它们分别等于1，否则为0。

data$alt_1 <- ifelse(data$alt == 1, 1, 0)
data$alt_2 <- ifelse(data$alt == 2, 1, 0)

现在使用mlogit.data() 函数。

然后，在您的模型中，您将指定 alt_1 和 alt_2 为遵循某种分布的随机参数。现在你有随机的替代特定常数，你估计平均值和标准偏差。如果您希望它们是具有零均值和单位标准差的简单误差分量，您可以将截距的均值和 sd 参数分别固定为 0 和 1。