首先根据Krinsky and Robb提出的限参Bootstrap方法在随机参数Logit模型估计结果的基础上模拟了2000个样本的偏好参数,并据此计算出2000个边际价值。同时，采用组间T检验法和Poe et al.提出的“非参数完全组合方法”对边际价值进行组间均值检验。
上述情况应用stata软件，应该如何执行呢？相关stata代码又是如何？
* 读取数据
use data.dta
* 估计随机参数Logit模型
// 这里假设模型的自变量为X，因变量为Y，随机效应为Z
// 请根据具体情况修改
mixed Y X || Z:
* 进行限参Bootstrap并得到2000个边际价值
// 这里假设边际价值为M，Bootstrap重抽样次数为B=2000
// 请根据具体情况修改
bootstrap M = r(M), reps(2000) saving(bootstrap.dta, replace): margins X
* 进行组间均值检验
// 这里采用Poe et al.提出的“非参数完全组合方法”
// 请根据具体情况修改
teffects ipw (Y X, Z), atmeans
teffects combperm (M), reps(2000) saving(combinatorial.dta, replace)