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

我开始使用Stata 14.我试图做一些基本的风险比分析，但我不知道如何提取单个结果。考虑下面的代码：Stata - 队列研究 - 显示原油风险比（如r（rr_crude））

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

我们得到四个风险比率输出，两个年龄组，粗一个和M-H之一。用

dis r(rr)

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

dis r(rr_crude) 
dis r(rr_mh)

或类似的东西？我还没有找到解决方案。或者有可能做一些事情，比如将输出保存在一个矩阵中，并用行和列索引来表示它？我还没有在文档中找到解决方案。

来源

2015-05-07 dan

请注意网站的政策，保持对帖子的个人评论到最低限度。这就是为什么你的帖子被削减。我们很高兴相信你很有礼貌，并会很感激。见例如http://meta.stackexchange.com/questions/2950/should-hi-thanks-taglines-and-salutations-be-removed-from-posts –

'return list'显示可用的内容。它不包括你想要的。不是我的领域，所以我不能评论为什么。 –

编辑： 只要创建标量，其坚持的内存

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

来源

2015-05-11 21:46:47

但是你松了rr_mh。我可以做'cs case cases [fw = pop]'dis r（rr）'，这是粗糙的，然后与分层一起用'r（rr）'得到rr_mh。但我想同时做到这一点。 – dan

只需创建一个新的标量rr_mh。我已经修改了答案来展示这一点。 –

@dan如果你接受史蒂夫的回答，你将获得声望（和史蒂夫一样）。 –

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

回答

相关问题