I want to calculate the sensitivity and specificity of a diagnostic test for TB. The index test is a numerical variable TBscore with values ranging from 0.86 to 100. However, I also want to adjust for age, sex, cough and BMI, and I have decided to use logistic regression as below
logistic TBDisease TBscore age sex cough bmi
predict p if e(sample), p
roctab TBDisease p, detail
I then get the output below
Detailed report of Sensitivity and Specificity
------------------------------------------------------------------------------
Correctly
Cutpoint Sensitivity Specificity Classified LR+ LR-
------------------------------------------------------------------------------
( >= 1.46e-10 ) 100.00% 0.00% 2.53% 1.0000
( >= 6.42e-08 ) 100.00% 0.04% 2.57% 1.0004 0.0000
( >= .0001824 ) 100.00% 0.07% 2.60% 1.0007 0.0000
( >= .0003155 ) 100.00% 0.11% 2.63% 1.0011 0.0000
( >= .0003179 ) 100.00% 0.14% 2.67% 1.0014 0.0000
( >= .0003397 ) 100.00% 0.18% 2.70% 1.0018 0.0000
( >= .000353 ) 100.00% 0.21% 2.74% 1.0021 0.0000
( >= .0003924 ) 100.00% 0.25% 2.77% 1.0025 0.0000
( >= .0004425 ) 100.00% 0.28% 2.81% 1.0028 0.0000
( >= .0004429 ) 100.00% 0.32% 2.84% 1.0032 0.0000
( >= .0004773 ) 100.00% 0.35% 2.87% 1.0035 0.0000
( >= .0004885 ) 100.00% 0.39% 2.91% 1.0039 0.0000
( >= .0005115 ) 100.00% 0.42% 2.94% 1.0042 0.0000
( >= .0005456 ) 100.00% 0.46% 2.98% 1.0046 0.0000
( >= .0005482 ) 100.00% 0.49% 3.01% 1.0049 0.0000
( >= .0005705 ) 100.00% 0.53% 3.04% 1.0053 0.0000
( >= .0005847 ) 100.00% 0.56% 3.08% 1.0056 0.0000
( >= .0005957 ) 100.00% 0.60% 3.11% 1.0060 0.0000
( >= .0006096 ) 100.00% 0.63% 3.15% 1.0064 0.0000
( >= .0006124 ) 100.00% 0.67% 3.18% 1.0067 0.0000
Question: How do I get the actual scores in the adjusted model (cut-point)? I want to know what value of the actual TB score each of these predicted values correspond to. I can get the values in the unadjusted model by simply using "roctab TBDisease TBScore, detail" but this wont give me the same AUC as well as sensitivity and specificity given by the adjusted model.Is there a way of back transforming the predicted values to the actual score. When I select the cut-off using the predicted values, I want to know what TB score values correspond to it.
Thanks
Humphrey
Related Posts with Predicted probabilities and actual values
Rolling window granger causality or time varying granger causalityHi everyone, I wanted to know if someone would know how to implement a rolling window granger causa…
Model specificationHello STATALISTERS! I have a question regarding model specification. Is it okay if I transformed tw…
Categorical variable per key (no unique keys)Hi guys, currently I am struggling with below table: Array The second column contains a combination…
What is the problem in my code on graphing?gr tw (scatter pgratioregional year, by(region)) (lfit pgratioregional year, by(region)) /// xlabel(…
color of box plotsHi, I created two box plots for the birth weight of infants based on two interventions during pregna…
Subscribe to:
Post Comments (Atom)
0 Response to Predicted probabilities and actual values
Post a Comment