版内不知为何都没有关于运用stata软件进行剂量效应分析的讨论,看了一堆文献云里雾里的,glst命令怎么用,还请高手们给予指导!需要文章提供哪些数据呢?剂量,cases,control,OR,RR,95%CI?表格中该如何输入?怎么才能出现如下图(不知道为什么上传不了图片,这是图片的链接http://i.0dxy.cn/upload/2010/12/08/51812658.jpg)的剂量效应图呢?
我用以下命令结果出来了,但是没有图,怎么才能显示出图呢?
. use http://nicolaorsini.altervista.org/stata/data/dose.dta, clear
(Dose-response data)
. gen double logor = log(adjor)
. gen double logorlb = log(lb)
. gen double logorub = log(ub)
. gen double se = ((logorub - logorlb)/(2*invnorm(.975)))
. vwls logor dose in 2/4, sd(se) nocons
. mat list e(V)
. glst logor dose, se(se) cov(N case) cc
. mat list e(Sigma)