Stata coefplot中的编辑置信区间
我正在Stata中使用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的标准误差添加到该代码中?
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%置信区间相匹配。不精确,但足够接近我正在绘制的图。