概要
marglmean是一款专用于回归模型后对特定情景下结局变量的边际均值(即调整协变量后的平均预测值)进行精确统计推断的Stata社区命令。该命令适用于预测值被解释为非负结局变量条件算术均值的模型(如logit、logistic、probit、poisson或glm等广义线性模型),有效弥补了常规margins命令仅提供对称置信区间的局限。
依赖
Stata 12.0及以上
安装
通过Stata Journal安装:
net sj 13-4 st0314 net install st0314.pkg, replace
示例
示例1,分别估计了在真实世界与一个假设情景(该情景中所有汽车车型均为美国制造,但车重与真实世界保持一致)下,汽车里程(即每加仑行驶英里数)的边际均值或其对数值。
sysuse auto glm mpg weight foreign, fam(gamma) link(log) robust eform marglmean, eform marglmean, at(foreign=0) eform marglmean marglmean, at(foreign=0)
示例2,演示了marglmean命令与可从SSC下载的parmest命令的联合使用。通过marglmean(配合post选项)估计里程的边际均值,并利用parmest将结果保存至内存中的数据集,覆盖原有数据集;该数据集包含一条记录,对应未转换的参数(命名为Scenario_1),并包含估计值、置信区间上下限、P值及其他参数属性信息。随后,对新生成的数据集进行描述和列出其内容:
glm mpg weight foreign, fam(gamma) link(log) robust eform marglmean, eform post parmest, eform norestore describe list
更多解读可以参考作者的会议发言幻灯[1]和Stata Journal上的论文[2]。
作者
姓名:Roger Newson
机构:Imperial College London, UK.
邮箱:r.newson@imperial.ac.uk