Stata coefplot中的编辑置信区间

Stata coefplot中的编辑置信区间,stata,coefplot,Stata,Coefplot,我正在Stata中使用coefplot命令绘制多元回归模型的系数和置信区间。我正在绘制4种不同型号规格的相同系数(X) 有一个模型规范(替代标准误差),我无法确定如何在Stata中进行估计,但我能够使用R进行估计。这意味着对于一个模型规范,我在R中有标准误差,但在Stata中没有 有没有一种简单的方法可以手动更改coefplot中的标准错误 我的代码是: coefplot A B C D, drop(_cons) xline(0) keep(X) 如何将模型D中系数X的标准误差添加到该代码中?

我正在Stata中使用
coefplot
命令绘制多元回归模型的系数和置信区间。我正在绘制4种不同型号规格的相同系数(X)

有一个模型规范(替代标准误差),我无法确定如何在Stata中进行估计,但我能够使用R进行估计。这意味着对于一个模型规范,我在R中有标准误差,但在Stata中没有

有没有一种简单的方法可以手动更改coefplot中的标准错误

我的代码是:

coefplot A B C D, drop(_cons) xline(0) keep(X)

如何将模型D中系数X的标准误差添加到该代码中?您可以手动编辑
e(V)
(方差-协方差矩阵)和
e(b)
向量。为此,请定义一个程序:

est restore estimates1

 capture program drop changeeV
 program changeeV, eclass
   tempname b V 
   matrix `b' = e(b)
   matrix `V' = e(V)
   matrix `V'[1,1] = 1.1 // Add values of new variance-covariance matrix
   matrix `b'[1,1] = 10 // Add new coefficient vector
   ereturn post `b' `V' // Repost new vectors 
   ereturn local cmd "reg outcome treatment covariates"
          // Repost initial command (required)
 end
changeeV // Execute program

est store eaX  // Store new generated estimtes

请注意,要获得协方差矩阵,需要从
R
中的输出中获取标准误差的平方。祝你好运

我可以问一下您使用的是哪种“替代标准错误”吗?另一种解决方法是使用
two-way
组合
rcap
绘图和
scatter
绘图,以生成类似于
coefplot
生成的绘图。缺点是您可能需要使用
input
手动指定所需的系数和上下限。很高兴提供一个对您有用的示例。Coefplot有
se
ci
选项,您可以使用ci/se创建一个矩阵,并将其用作输入,而不是默认的
e(V)
。谢谢,Jonathan!我相信这会奏效,但我还是有点不知所措。我找到了一个合理的解决方案,使用ci(Y),我绘制了Stata图,并绘制了超大置信区间(97.5%而不是95%),因此绘制的Stata置信区间与用R估计的95%置信区间相匹配。不精确,但足够接近我正在绘制的图。