Stata - 队列研究 - 显示原油风险比（如r（rr_crude））

Question

我开始使用Stata 14.我尝试进行一些基本的风险比分析，但我不知道如何提取单个结果。给出以下代码：

clear all
webuse ugdp
cs case exposed [fw=pop], by(age)

我们得到的输出具有四种风险比，适用于年龄类别，粗略类型和M-H类型。与

dis r(rr)

我得到最后一个（？）比率，但可以指定它吗？像

dis r(rr_crude)
dis r(rr_mh)

或类似的东西？我还没有找到解决方案。或者是否可以执行诸如将输出保存在矩阵中并使用行和列索引指示它的操作？ 我还没有在文档中找到解决方案。

Answer 1

修改 只需创建持久存储在内存中的标量

 clear all
 webuse ugdp
 cs case exposed [fw=pop], by(age)
 scalar rr_mh = r(rr)

然后使用glm：

 glm case exposed [fw = pop], family(binomial) link(log)
 scalar rr_crude = exp(_b[exposed])

或

 cs case exposed [fw = pop]
 scalar rr_crude = r(rr)

在任何一种情况下：

 di rr_crude 
 di rr_mh