Bookmark and Share

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/


© Copyright 1996–2018 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   Site index