Hi all,
I'm trying to compare sensitivity and specificity for two diagnostic tests with SROC curves. I am using the metandi & metandiplot commands. My question is how can I combine the SROC curves for test 1 & test2 into a single graph instead of generating separate graphs for both? These are my commands for test 1 & 2, with failp & curep corresponding to true & false positives respectively and failn & curen corresponding to false & true negatives respectively:
metandi failp3 curep3 failn3 curen3 if test3==3 & test1==1,nolog
metandiplot failp3 curep3 failn3 curen3 if test3==3 & test1==1 [aw=1], predopts(off) graphregion(fcolor(white)) graphregion(margin(none)) graphregion(margin(none)lcolor(white)) legend(off) xtit("Specificity", size(medium)) ytit("Sensitivity", size(medium)) summopts(msymb(S)msize(1.5)mcolor(blue)) confopts(lcolor(blue)) curveopt(lcolor(blue)lwidth(0.4)) studyopts(mlcolor(blue) msym(oh)msize(1.5))
metandi failp1 curep1 failn1 curen1 if test3==3 & test1==1,nolog
metandiplot failp1 curep1 failn1 curen1 if test3==3 & test1==1 [aw=1], predopts(off) graphregion(fcolor(white)) graphregion(margin(none)) graphregion(margin(none)lcolor(white)) legend(off) xtit("Specificity", size(medium)) ytit("Sensitivity", size(medium)) summopts(msymb(S)msize(1.5)mcolor(red)) confopts(lcolor(red)) curveopt(lcolor(red)lwidth(0.4)) studyopts(mlcolor(red) msym(oh)msize(1.5))
Thanks in advance!
I'm trying to compare sensitivity and specificity for two diagnostic tests with SROC curves. I am using the metandi & metandiplot commands. My question is how can I combine the SROC curves for test 1 & test2 into a single graph instead of generating separate graphs for both? These are my commands for test 1 & 2, with failp & curep corresponding to true & false positives respectively and failn & curen corresponding to false & true negatives respectively:
metandi failp3 curep3 failn3 curen3 if test3==3 & test1==1,nolog
metandiplot failp3 curep3 failn3 curen3 if test3==3 & test1==1 [aw=1], predopts(off) graphregion(fcolor(white)) graphregion(margin(none)) graphregion(margin(none)lcolor(white)) legend(off) xtit("Specificity", size(medium)) ytit("Sensitivity", size(medium)) summopts(msymb(S)msize(1.5)mcolor(blue)) confopts(lcolor(blue)) curveopt(lcolor(blue)lwidth(0.4)) studyopts(mlcolor(blue) msym(oh)msize(1.5))
metandi failp1 curep1 failn1 curen1 if test3==3 & test1==1,nolog
metandiplot failp1 curep1 failn1 curen1 if test3==3 & test1==1 [aw=1], predopts(off) graphregion(fcolor(white)) graphregion(margin(none)) graphregion(margin(none)lcolor(white)) legend(off) xtit("Specificity", size(medium)) ytit("Sensitivity", size(medium)) summopts(msymb(S)msize(1.5)mcolor(red)) confopts(lcolor(red)) curveopt(lcolor(red)lwidth(0.4)) studyopts(mlcolor(red) msym(oh)msize(1.5))
Thanks in advance!
