新手小白学习STATA网状meta分析的指令分享
跟着网状meta分析公开课的老师认真学习了一周的时间,基本上跑图的流程都摸索了个大概,给大家分享一下:
1.首先提取数据的格式:这一步我自己不能像老师一样直接把excel表给import进去,所以我是直接打开STATA手动录入的。其中:
二分类变量录入格式:
连续性变量录入格式:
2.再确认一下你的各种包是不是都装全了,我做的过程中经常出现各种红字然后丁香园和百度里各种搜索对侧完善安装各种包才好的,出了mvmeta和network包以外,以下的包建议都装上免得又出问题的红字:
ssc install metareg
ssc install metan
ssc install sencode
3.录入完数据后就是对录入数据的各种处理,具体原理老师讲的我也没吃透,就是死记指令,(其中set matsize 5000 是因为我纳入研究太多矩阵不够所以加的扩展矩阵的。)
二分类数据处理:
network setup d n, studyvar(study) trtvar(trt) ref(1)
set matsize 5000
network meta c
network rank max
network meta i
network sidesplit all
连续性数据处理:
network setup mean sd n,studyvar(study) trtvar(trt) ref(1)
set matsize 5000 network
meta c network rank max
network meta i
network sidesplit all
4.然后可以把图上对应的干预措施ABCDEFG替换本来的名字(这是我自己的例子,根据你具体的干预措施替换中文文字)指令如下:
replace _t1="安慰剂" if _t1=="A" .
replace _t1="干扰素" if _t1=="B" .
replace _t1="芬戈莫德" if _t1=="C" .
replace _t2="安慰剂" if _t2=="A" .
replace _t2="干扰素" if _t2=="B" .
replace _t2="芬戈莫德" if _t2=="C"
5.以下是各种图的指令:其中“network convert pairs ”和“network convert augment ”就是两种数据的格式根据图的要求相互转换就可以了。
网状关系图
network convert pairs
networkplot _t1 _t2
贡献图
network convert pairs
netweight _y _stderr _t1 _t2,color(navy) symbol(circle) title(无复发人数)
不一致性检验图
network convert pairs
ifplot _y _stderr _t1 _t2 study, tau2(loop) eform summary
漏斗图
network convert pairs
netfunnel _y _stderr _t1 _t2,bycomparison ytitle(standard error of logor)
倒三角图
network convert augment
set matsize 5000
network meta c
netleague,eform lab(安慰剂 干扰素 芬戈莫德)sort(干扰素 芬戈莫德 安慰剂) nokeep
预测区间图
network convert augment
set matsize 5000 network meta c
intervalplot,eform pred null(1)lab(安慰剂 干扰素 芬戈莫德 )
SUCRA疗效比较图
network convert augment
set matsize 5000
network meta c
network rank max,zero all reps(10000) gen(prob)
sucra prob*,labels(安慰剂 干扰素 芬戈莫德) lcol(red black)
以上是我的一些学习心得,总体来说也算基本入门,通过以上指令能做一些基础的网状meta分析的图了,当然网课老师讲的不止这些,具体想要了解的童鞋们还是值得花钱学习咱们丁香园的公开课滴,收货非常大。
附上以上指令的TXT文本供大家学习
最后编辑于 2022-10-09 · 浏览 2.0 万