全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
5678 3
2010-12-09
版内不知为何都没有关于运用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)
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2010-12-9 05:32:14
LZ的图片是叉烧包 看不到
glst没听说过……不过rr和or 包括95CI是有很多方法得到的 最简单的就是cc和cs命令了
只要定义好exposure(0=unexposed 1=exposure)和case(0=control 1=case)就可以直接通过这两个命令得到rr和or
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2010-12-9 18:26:17
谢谢指点,我输入以下的命令还要加什么命令可以出现剂量效应图呢?
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)
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-8-21 15:03:55
你用的stata是什么版本?可以共享给我不。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群