Stata中coefplot替换变量估计值为常数的实现方法

当某一模型未包含目标变量（或该变量无有效估计值），需用常数替代展示时，可通过以下几种方法实现：

方法一：用estadd手动补充模型结果 #

先保存各模型的估计结果，再为缺失变量手动添加系数和标准误：

* 拟合并保存第一个包含所有变量的模型
regress y var1 var2
estimates store model1

* 拟合并保存第二个仅包含部分变量的模型
regress y var1
estimates store model2

* 为第二个模型补充缺失变量的系数与标准误
estimates restore model2
estadd scalar b_var2 = 0  // 将var2的系数设为指定常数（此处为0）
estadd scalar se_var2 = 0  // 设定标准误为0，避免显示误差棒；若需误差线可设为自定义值
estimates store model2_updated

* 绘制系数图，指定要展示的变量
coefplot model1 model2_updated, keep(var1 var2) vertical ///
    coeflabels(var1="变量1" var2="变量2") ///
    modelnames("全变量模型" "简化模型")

方法二：构建自定义系数数据集 #

提取所有模型的系数到数据集，手动补充缺失变量的常数后，用coefplot读取该数据集：

* 提取第一个模型的系数到临时数据集
estimates restore model1
estimates stats, clear
rename (b se) (b_model1 se_model1)
keep varname b_model1 se_model1
save model1_coef, replace

* 提取第二个模型的系数并与第一个模型合并
estimates restore model2
estimates stats, clear
rename (b se) (b_model2 se_model2)
merge 1:1 varname using model1_coef, nogen

* 手动添加缺失变量的常数行
if !inlist("var2", varname) {
    set obs `=_N+1'
    replace varname = "var2" in `=_N'
    replace b_model2 = 0 in `=_N'  // 设定常数
    replace se_model2 = 0 in `=_N'
}

* 保存合并后的系数数据集
save all_coef_data, replace

* 基于自定义数据集绘制系数图
coefplot from(all_coef_data), vertical ///
    coeflabels(var1="变量1" var2="变量2") ///
    modelnames("全变量模型" "简化模型") ///
    ytitle("系数估计值")

方法三：用addplot手动添加常数点 #

若仅需快速添加常数点（无需匹配误差棒），可直接在coefplot中用addplot选项插入：

* 先绘制基础系数图
coefplot model1 model2, keep(var1 var2) vertical ///
    coeflabels(var1="变量1" var2="变量2") ///
    modelnames("全变量模型" "简化模型")

* 添加简化模型中var2的常数点（0为常数，2为var2在x轴的位置，需根据变量顺序调整）
addplot(scatteri 0 2, msymbol(circle) mcolor(blue) legend(off))

注意事项 #

• 补充的变量名必须与其他模型中的变量名完全一致（包括大小写、下划线等），否则coefplot无法匹配变量；
• 若需显示误差棒，可将标准误设为合理的自定义值（如变量的标准差），而非0；
• 常数的数值可根据需求调整（如变量均值、特定阈值等），不限于0。

内容的提问来源于stack exchange，提问作者Data_Science_110