stata反分析代码剂命令.pdfVIP

Fruit:

gendoublelogrr=log(adjrr)

gendoubleloglb=log(lb)

gendoublelogub=log(ub)

gendoublese=((logub-loglb)/(2*invnorm(.975)))

glstlogrrdose,se(se)cov(ncase)pfirst(idstudy)

（发现存在异质性，此时的P值可作为最后是否有意义的P值）

glstlogrrdose,se(se)cov(ncase)pfirst(idstudy)random

（不可以一开始直接就用随机）

capturedropdoses*

_pctiledose,percentile(5356595)

retlist

mksplinedoses=dose,knots(`=r(r1)`=r(r2)`=r(r3)`=r(r4))cubicdisplayknots

glstlogrrdoses*,se(se)cov(ncase)pfirst(idstudy)random

testparmdoses2doses3

glstlogrrdose,se(se)cov(ncase)pfirst(idstudy)ts(f)

predictnllrr_lin=_b[dose]*dose

genrr_lin=exp(lrr_lin)

glstlogrrdoses*,se(se)cov(ncase)pfirst(idstudy)random

predictnllogrrwithref=_b[doses1]*doses1+_b[doses2]*doses2+_b[doses3]*doses3,ci(lohi)

genrrwithref=exp(logrrwithref)

genlbwithref=exp(lo)

genubwithref=exp(hi)

levelsofdose,local(level)

xblcdoses*,c(dose)at(`r(levels))ref(0.5)eform

twoway(linelbwithrefubwithrefrrwithrefdose,sortlp(longdashlongdashl)lc(blackblackblack))

(linerr_lindose,sortlp(shortdash)lc(black)),scheme(s1mono)ylabel(.8511.151.3,angle(horiz)

format(%3.2fc))xlabel(0.5(1)4.5)legend(off)ytitle(RelativeRisk,margin(right))xtitle(Fruit

intake,servings/day,margin(top_bottom))name(figure1,replace)yscale(log)

plotregion(style(none))

Fruitandvegetables:

gendoublelogrr=log(adjrr)

gendoubleloglb=log(lb)

gendoublelogub=log(ub)

gendoublese=((logub-loglb)/(2*invnorm(.975)))

glstlogrrdose,se(se)cov(ncase)pfirst(idstudy)

（未发现存在异质性，此时的P值可作为最后是否有意义的P值）

capturedropdoses*

_pctiledose,perce