概述

margprev是一款专为二元回归模型（如logit、logistic、probit或glm）的后处理设计的Stata社区命令，用于计算人群边际患病率（Marginal Prevalence）。该命令利用margins和nlcom框架，对情景下的边际患病率及其置信区间进行精确估计，是量化公共卫生干预效果的重要统计工具。

依赖

Stata 14.0 及以上

安装

通过Stata Journal安装：

net sj 13-4 st0314
net install st0314.pkg, replace

示例

示例1，分别估计了低出生体重在现有情景（即实际观测数据所反映的真实世界）与一个假设情景（所有母亲均不吸烟）下的边际患病率：

webuse lbwlogit low i.race i.smoke, or robust
margprev
margprev, at(smoke=0)

示例2，演示了margprev命令与可从SSC下载的parmest命令的联合应用：首先，通过margprev（启用post选项）估计低出生体重的边际患病率，并利用parmest将结果（以logit转换后的参数形式）保存至内存中的新数据集（覆盖原始数据集），该数据集包含一条观测记录，对应转换后的参数（命名为Scenario_1），并存储其估计值、置信区间、P值及其他参数属性；随后，使用replace命令将logit尺度下参数的对称置信区间，替换为原始概率尺度（即未转换的边际患病率）下更符合实际分布特征的非对称置信区间，最后对更新后的数据集进行描述并列出其内容。

logit low i.race i.smoke, or robust
margprev, eform post
parmest, norestore
foreach Y of var estimate min* max* {
replace `Y'=invlogit(`Y')
}
describe
list

作者

姓名：Roger Newson
机构：Imperial College London, UK.
邮箱：r.newson@imperial.ac.uk