/* fig6_1.do

   Figure 6.1 (Example 6.1)

   test distributions for hearing vs. not, by frequency/intensity

   data required: dp2.dta
      Stover/Gorga DPOAE test data

   last update:  23 Apr 2003
*/
version 7
set more off
cap semt_profile
cap log close

use ${semt_data}dp2,clear

egen int fl = group(f l), l
compress
sort f l d y
keep f l fl d y
order f l fl d y
log using ${semt_log}fig6_1, replace

* categories corresponding to 9 plots:

ta f l

la l fl

log close

qui sum y ,meanonly
local xmin = r(min)
local xmax = r(max)
local xlab "-40(20)20"
local xtick "-30(10)10"
local ymax = .068

tempfile f1 f2 f3 f4 f5 f6 f7 f8 f9

preserve
set textsize 200

forvalues i = 1/9 {
    #delimit ;
    kdensity y if d==0 & fl == `i', width(1.8)
      gen(x0 fx0) nogr ;
    kdensity y if d==1 & fl == `i', width(2.5) /* 2.5 */
      gen(x1 fx1) nogr  ;
    keep if x0~=. & x1~=. ;
    stack x0 fx0 x1 fx1, into(x y) wide clear;
    gr fx0 fx1 x, s(..) c(ss)
      xscale(`xmin', `xmax') xlab(`xlab')
      /* t1(`titl2') */
      b2(" ") /* ysca(0,`ymax') */
      title(" ")
      saving(`f`i'') ;
    #delimit cr
    restore, preserve
}

gprefs set window xsize 9.0
gprefs set window ysize 6.5
gprefs set window usegphsize off
#delimit ;

gr using `f1' `f2' `f3' `f4' `f5' `f6' `f7' `f8' `f9', margin(10) t1("  ") t2(" ")
saving(${semt_log}fig6_1,replace)
;
#delimit cr
/*  regenerate and save density vars for plotting externally */

#delimit ;

forvalues i = 1/9 { ;
    kdensity y if d==0 & fl == `i', width(1.8) n(150)
      gen(x0_`i' fx0_`i') nogr ;
    kdensity y if d==1 & fl == `i', width(2.5) n(150)
      gen(x1_`i' fx1_`i') nogr  ;
} ;
/* 11/11/02 change: use -SNR rather than SNR, s.t. normals have lower values: */

foreach var of varlist x* {
    replace `var' = `var'
}
d
pause
#delimit cr
keep if x0_1 ~= .
keep x* fx*
save ${semt_log}fig6_1,replace