在Stata中保存置信限
运行在Stata中保存置信限,stata,confidence-interval,Stata,Confidence Interval,运行glm后,我可以键入矩阵列表r(表)并查看我所有结果的表。如果我愿意,我可以将斜率和SEs写入变量,例如,gen B=\u B[x1]或gen se=\u se[x1]。然而,这不适用于置信限ll和ul。我如何才能以类似的方式访问它们?我不确定\u b[]和\u se[]结果是否与r(table)相关——我认为它们是e(b)和e(V)的产品 无论如何,由于您有r(table),您可以将结果保存到另一个矩阵中,然后使用常规矩阵运算将下限和上限放入新矩阵中。如果出于某种原因需要转换为变量(例如,
glm
后,我可以键入矩阵列表r(表)
并查看我所有结果的表。如果我愿意,我可以将斜率和SEs写入变量,例如,gen B=\u B[x1]
或gen se=\u se[x1]
。然而,这不适用于置信限ll和ul。我如何才能以类似的方式访问它们?我不确定\u b[]
和\u se[]
结果是否与r(table)
相关——我认为它们是e(b)
和e(V)
的产品
无论如何,由于您有r(table)
,您可以将结果保存到另一个矩阵中,然后使用常规矩阵运算将下限和上限放入新矩阵中。如果出于某种原因需要转换为变量(例如,绘图),则始终存在-svmat-
sysuse auto,clear
glm price mpg foreign, f(gaussian)
mat r=r(table)
matrix ll=r["ll",....]' // see -help matrix extraction-; transposed for svmat
svmat ll,names(ll) // lower bounds are in variable ll1
B=\u B[x1]
(是否有gen
缺失?)不完全是“将斜率和SEs写入变量”,因为\u B[x1]
中只有一个值。也许你的意思是将结果写入标量或宏?@AspenChen是的,正如我所说明的变量赋值,而不是编程命令。我同意这一点。我将尝试下面您建议的矩阵方法。谢谢。您还可以查看SSC的parmest
。