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]
[no subject]
Thanks for your help, Carole
On Tue, Mar 25, 2014 at 11:35 AM, Joe Canner <[email protected]> wrote:
> 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/
*
* 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/