2015-05-07 28 views
0

我開始使用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) 

或類似的東西?我還沒有找到解決方案。或者有可能做一些事情,比如將輸出保存在一個矩陣中,並用行和列索引來表示它? 我還沒有在文檔中找到解決方案。

+0

請注意網站的政策,保持對帖子的個人評論到最低限度。這就是爲什麼你的帖子被削減。我們很高興相信你很有禮貌,並會很感激。見例如http://meta.stackexchange.com/questions/2950/should-hi-thanks-taglines-and-salutations-be-removed-from-posts –

+0

'return list'顯示可用的內容。它不包括你想要的。不是我的領域,所以我不能評論爲什麼。 –

回答

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 
+0

但是你鬆了rr_mh。我可以做'cs case cases [fw = pop]'dis r(rr)',這是粗糙的,然後與分層一起用'r(rr)'得到rr_mh。但我想同時做到這一點。 – dan

+0

只需創建一個新的標量rr_mh。我已經修改了答案來展示這一點。 –

+0

@dan如果你接受史蒂夫的回答,你將獲得聲望(和史蒂夫一樣)。 –