SAS suggestions for sens-spec CI: 24170 - Estimating sensitivity, specificity, positive and negative predictive values, and other statistics
... View more
Yes, thank you, I completely agree with you. I cannot understand why SAS propose to calculate confidence limits with proc freq IF it's clear that the results would be different. And I cannot understand how can I now obtain my CLs for specificity and sensitivity. Should I/Could I use an online-calculator? I don't know.
... View more
Ok, I have found how the SAS calculates a classificaiton table: SAS/STAT(R) 9.2 User's Guide, Second Edition Now I need to know how it calculates confidence intervals for sensitivity ans specificty, LR- and LR- using the proc logistic.
... View more
They are pretty close. From 147 subjects (47 in desease) 26 of diases and 92 with no disease were correctly classified by proc logistic (true positive/true negative), while 27 and 95 respectively were correctly classified with proc freq. It takes me a difference of 3% for sensitivity and for specificity,
... View more
We could create a classification table in two ways: 1. Using proc logistic with ctable pprob=xxx Example: proc logistic desc data=mmse ; model fn= lhippoc lmidtemp eicv c_age_a c_age_b ss/ ctable pprob=0.32; run; 2. Using output and manipulating with data: proc logistic desc data=mmse ; model fn= lhippoc lmidtemp eicv c_age_a c_age_b ss/ ctable pprob=0.32; output out=ci_with_outl p=rsk; run; data ci_with_outl; set ci_with_outl; if rsk >=0.32 then pos=1; else pos=0; run; proc sort; by descending pos descending fn; proc freq order=data; table pos*fn; run; The problem is following. I have received two different classification table, with different numbers of true/faux positives and negatives. The question is: what is an algorithm of calculation of true/faux positives and negatives in proc logistic ctable?
... View more