*! rocmicd.ado 1.0.0 17th March 2009, R.Froud.
* This program reports minimally important change (MIC) thresholds for deterioration using 45 degree tangent line intersection
* and the smallest residual sum of sensitivity ans specificity.
capture program drop rocmicd
version 10.1
program define rocmicd, rclass
syntax varlist(min=2 max=2 numeric)[, scale(real 0)]
return local varname `varlist'
local ref: word 1 of `varlist'
local dv: word 2 of `varlist'
quietly {
capture drop se1 sp1 gap1 gap thresh emg emg1
logistic `ref' `dv'
lsens, gensens(se1) genspec(sp1)
gen gap1=se1-sp1
gen gap=sqrt(gap1^2)
gen thresh= `dv' - `scale'
sort gap
return scalar mic = thresh in 1/1
gen emg1 = ((1-se1)+(1-sp1))
gen emg = sqrt(emg1^2)
sort emg
return scalar emgo = thresh in 1/1
}
lroc
display in smcl as text "the MIC calculated using a 45 degree tangent line intersection is " as result %05.3f return(mic)
display in smcl as text "the MIC calculated using the smallest sum of residual sensitivity and specificity is " as result %05.3f return(emgo)
display in smcl as text "use bootstrap for Conf. Intervals"
roctab `ref' `dv'
capture drop se1 sp1 gap1 gap thresh emg emg1
end