[程式] STATA meta analysis一個問題
Generalized least squares for trend estimation
of summarized dose–response data
這是我看的paper 網路上可以免費下載
http://www.stata-journal.com/article.html?article=st0096
在做 dose–response
我看不懂這個pdf中的一個操作方法
4.4 Meta-analysis: Lactose intake and ovarian cancer risk
這裡他有6個case control study
3個cohort study
他一開始先用把這9個研究一起分析
指令是
glst logrr dose, se(se) cov(n case) pfirst(id study) eform
這裡可以理解
但是作者說 因為這兩種研究方法不同
所以要分開分析 這也可以理解
但是它分開的方法就難以理解
首先將case control研究和cohort研究定義為0和1兩種type
接著製造了一個交互作用項
叫做doseXtypes
再把這個doseXtypes 一起放入model中跑
所以指令是
glst logrr dose doseXtypes, se(se) cov(n case) pfirst(id study)
接著 如果想看case control 研究的結果
就只要鍵入
lincom dose + doseXtypes*0, eform
如果要看cohort的結果
就輸入
lincom dose + doseXtypes*1, eform
這就是難以理解的地方
為什麼創造一個交互作用項
一起跑個model
再分別輸入這兩個指令 就能夠看到兩個不同研究方式的結果呢?
有沒有人了解這個分析方法可以說明
拜託了謝謝您
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 125.227.132.216
※ 文章網址: http://www.ptt.cc/bbs/Statistics/M.1418119496.A.884.html