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
Placebo test using an RD designHi all, I have a question regarding the performance of a placebo test when using a regression disco…
too many signficant variablesHello I am doing a multilevel regression using PASEC2014 data at the first level I included several …
Placebo test using an RD designHi all, I have a question regarding the performance of a placebo test when using a regression disco…
Apply values across all in the same groupHi all, I have a dataset of about 10k participants across 600 groups (ParentSiteID in the dataset). …
Plotting weights in meta-analysisHello, I need to plot the study weights in forest plot. I am using the following command to create …
Subscribe to:
Post Comments (Atom)
0 Response to Multiple imputation STATA command
Post a Comment