Below is the command I am using to calculate the AUC for my prediction models in imputed datasets.
I am searching for a STATA command similar to "roctab: depvar p, detail" for imputed datasets.
Any suggestions for STATA command to calculate the sensitivity and specificity (at different cut-offs) for imputed datasets?
cap program drop eroctab
program eroctab, eclass properties(mi)
version 12.0
args refvar
capture drop p
logistic `refvar' x1 x2 x3
predict p
roctab `refvar' p
tempname b V
mat `b' = r(area)
mat `V' = r(se)^2
local N = r(N)
mat colnames `b' = AUC
mat colnames `V' = AUC
mat rownames `V' = AUC
ereturn post `b' `V', obs(`N')
ereturn local cmd "eroctab"
ereturn local title "ROC area"
end
mi estimate, cmdok: eroctab refvar
Thank you!!!
Related Posts with Multiple imputation STATA command
IV Quantile Regression and over-identificationDear Statalist members, I want to estimate an IV quantile regression but my model is over-identifie…
How to write a joint maximum likelihood estimation?Dear statalist: I want to write a joint maximum likelihood estimation(JMLE) for wage and unemploymen…
How to divide the samples into deciles?Hi! I am new to Stata I need to divide my sample into deciles in each year and industry. So I have …
scatter plot with line of best fitDear Stata users, I am interested in investigating the association between specific measurements (b…
equality of coefficients test with reghdfe: subsample analysis or interaction terms?I am trying to test the impact that margin trading status has on stock liquidity. I regress liquidit…
Subscribe to:
Post Comments (Atom)
0 Response to Multiple imputation STATA command
Post a Comment