/* ms2_fig6.do

  do-file to generate figure 6 for SJ manuscript:
    "Accomodating covariates in ROC analysis".
       Janes H, Longton G, and Pepe M.

  last update: 4/22/2008

  uses Norton neonatal audoilogy dataset:  nnhs2.dta
        -comproc- version 1.1.7 or later
    and -roccurve- version 1.3.9 or later

  Figure 6
    incremental value
        Y = DPOAE
        Z = gender, age

  ROC curves for logistic predicted values
  and comparison with -comproc-

*/
version 10
set more off
cap log close

set scheme sj

use http://www.fhcrc.org/science/labs/pepe/book/data/nnhs2, clear

set seed 73289755

log using sjms2_fig6.txt, replace

logit d currage gender
predict p1
logit d currage gender y1
predict p2

la var p2 "(age, gender, DPOAE)"
la var p1 "(age, gender)"

comproc d p2 p1, cl(id) auc pauc(.20) roc(.20) noccsamp

qui log close

roccurve d p2 p1, roc(.2) noccsamp cl(id)  ///
    name(g1, replace) title("")