/* fig5_9.do

   fig 5.9 from Example 5.7

   ROC-GLM (dfroc) curves based on model for t<.20
   (generate plot coordinate data for graphics pkg import)

   programs required
       dfroc.ado  (v 1.3.0 or later)

   empirical ROC curves for this figure were obtained from fig 5.4 (fig5_4.do)

   last update:  30 Apr 2003

*/
version 7
set more off
cap semt_profile
cap log close
pause on

clear
use ${semt_data}wiedat2b
log using ${semt_log}fig5_9, replace

dfroc y1 d, par(20) nobs
sca a1 = r(a)
sca b1 = r(b)
pause q <return> to continue . . .

dfroc y2 d, par(20) nobs
sca a2 = r(a)
sca b2 = r(b)

drop _all
set obs 200
range t 0 .2

gen roc1 = normprob(a1 + b1 * invnorm(t))
replace roc1 = 0 if t==0  /* otherwise missing/undefined */

gen roc2 = normprob(a2 + b2 * invnorm(t))
replace roc2 = 0 if t==0  /* otherwise missing/undefined */

sort t
save ${semt_log}fig5_9,replace

qui log off