stata做完混合logit模型后,如何计算选择概率,使用什么命令,比如mixlpred,举例说明
时间: 2023-10-02 21:06:51 浏览: 111
在 Stata 中,可以使用 `mixlpred` 命令来计算混合 logit 模型的选择概率。具体使用方法如下:
1. 首先需要使用 `mixlogit` 命令估计混合 logit 模型,比如:
```
mixlogit y x1 x2, group(id) nologit
```
这里假设因变量为 `y`,自变量为 `x1` 和 `x2`,`group(id)` 表示使用 `id` 变量来指定混合 logit 模型中的组别。
2. 然后使用 `mixlpred` 命令来计算选择概率,比如:
```
mixlpred, draw(1000) saving(pred)
```
这里使用 `draw(1000)` 参数来指定模拟抽样的次数,`saving(pred)` 参数表示将结果保存到名为 `pred` 的文件中。
`mixlpred` 命令会自动读取 `mixlogit` 命令估计的混合 logit 模型,并根据 `draw` 参数进行模拟抽样,然后计算每个个体在每个选择情况下的选择概率,并将结果保存到指定的文件中。
例如,假设 `pred` 文件中保存了以下内容:
```
id y1 y2 y3
---------------------------------------
1 0.7307 0.1801 0.0892
2 0.1973 0.6895 0.1132
3 0.8097 0.0775 0.1128
4 0.2255 0.2749 0.4996
5 0.6541 0.2463 0.0996
6 0.1394 0.4736 0.3870
7 0.4600 0.1181 0.4219
8 0.2849 0.5842 0.1309
9 0.7493 0.1050 0.1457
10 0.5177 0.1682 0.3141
```
其中 `id` 表示个体编号,`y1`、`y2`、`y3` 分别表示选择情况 1、2、3 的选择概率,则可以看出第一个个体在选择情况 1、2、3 分别的选择概率分别为 0.7307、0.1801、0.0892。