Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
st: RE: P-values for the difference in sensitivity in metandi
From
Joe Canner <[email protected]>
To
"[email protected]" <[email protected]>
Subject
st: RE: P-values for the difference in sensitivity in metandi
Date
Tue, 25 Mar 2014 18:35:41 +0000
Carole,
You say you want the "difference in sensitivity/specificity". What are you calculating the difference *between*? In your -metandi- command all I see is the components for computing a single sensitivity and specificity, not for comparing the sensitivity/specificity between two different modalities or two different thresholds.
Regards,
Joe Canner
Johns Hopkins University School of Medicine
-----Original Message-----
From: [email protected] [mailto:[email protected]] On Behalf Of Carole Lunny
Sent: Tuesday, March 25, 2014 2:02 PM
To: [email protected]
Subject: st: P-values for the difference in sensitivity in metandi
Dear Stata users,
I am using metandi to get results of a diagnostic accuracy review I am
doing. I am using the HSROC model, and would like to know how to get
the P-values for the difference in sensitivity/specificity in metandi
(or another stata command).
These are the results from metandi:
. metandi tp fp fn tn
Refining starting values:
Iteration 0: log likelihood = -233.08848
Iteration 1: log likelihood = -230.8551
Iteration 2: log likelihood = -230.79615
Iteration 3: log likelihood = -230.74257
Performing gradient-based optimization:
Iteration 0: log likelihood = -230.74257
Iteration 1: log likelihood = -230.74239
Iteration 2: log likelihood = -230.74239
Meta-analysis of diagnostic accuracy
Log likelihood = -230.74239 Number of studies = 31
------------------------------------------------------------------------------
| Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
Bivariate |
E(logitSe) | 1.552802 .1604973 1.238233 1.867371
E(logitSp) | 2.490983 .1836023 2.131129 2.850837
Var(logitSe) | .6461271 .2113009 .3403707 1.226546
Var(logitSp) | .8975503 .2730492 .4944354 1.629326
Corr(logits) | .1192086 .209656 -.2886231 .4904251
-------------+----------------------------------------------------------------
HSROC |
Lambda | 3.980269 .2660604 3.4588 4.501737
Theta | -.3043524 .2406498 -.7760173 .1673124
beta | .1643365 .2221673 0.74 0.459 -.2711034 .5997763
s2alpha | 1.704627 .5034386 .9555173 3.041027
s2theta | .3353756 .1085868 .1778002 .632602
-------------+----------------------------------------------------------------
Summary pt. |
Se | .8253181 .0231386 .7752564 .8661538
Sp | .9235073 .01297 .8938922 .9453619
DOR | 57.04186 14.58478 34.55841 94.15288
LR+ | 10.7895 1.883955 7.662542 15.19251
LR- | .1891505 .0254589 .1452912 .2462498
1/LR- | 5.286795 .7115817 4.060917 6.882731
------------------------------------------------------------------------------
Covariance between estimates of E(logitSe) & E(logitSp) .0029529
Any help would be appreciated. Thank you, Carole
*
* For searches and help try:
* http://www.stata.com/help.cgi?search
* http://www.stata.com/support/faqs/resources/statalist-faq/
* http://www.ats.ucla.edu/stat/stata/
*
* For searches and help try:
* http://www.stata.com/help.cgi?search
* http://www.stata.com/support/faqs/resources/statalist-faq/
* http://www.ats.ucla.edu/stat/stata/