----------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\MyDocs\MPH\LectureNotes\ClassDoLogData\Analysis of Variance\11.2.Isoproterenol.log
  log type:  text
 opened on:   7 Feb 2010, 11:36:39

. *  11.2.Isoproterenol.log
. *
. *  Plot mean forearm blood flow by race and log dose of isoproternol 
. *  using the data of Lang et al. (1995).  Show standard deviation for 
. *  each race at each drug level.
. *
. use C:\wddtext\11.2.Isoproterenol.dta, clear

. table race, row

----------------------
     Race |      Freq.
----------+-----------
    White |         13
    Black |          9
          | 
    Total |         22
----------------------

. list if id==1 | id ==22

     +----------------------------------------------------------------------+
     | id    race   fbf0   fbf10   fbf20   fbf60   fbf150   fbf300   fbf400 |
     |----------------------------------------------------------------------|
  1. |  1   White      1     1.4     6.4    19.1       25     24.6       28 |
 22. | 22   Black    2.1     1.9       3     4.8      7.4     16.7     21.2 |
     +----------------------------------------------------------------------+

. generate baseline = fbf0

. *
. * Convert data from one record per patient to one record per observation.
. *
. more

. reshape long fbf, i(id) j(dose)
(note: j = 0 10 20 60 150 300 400)

Data                               wide   ->   long
-----------------------------------------------------------------------------
Number of obs.                       22   ->     154
Number of variables                  10   ->       5
j variable (7 values)                     ->   dose
xij variables:
                  fbf0 fbf10 ... fbf400   ->   fbf
-----------------------------------------------------------------------------

. list if id==1 | id ==22

     +-------------------------------------+
     | id   dose    race    fbf   baseline |
     |-------------------------------------|
  1. |  1      0   White      1          1 |
  2. |  1     10   White    1.4          1 |
  3. |  1     20   White    6.4          1 |
  4. |  1     60   White   19.1          1 |
  5. |  1    150   White     25          1 |
     |-------------------------------------|
  6. |  1    300   White   24.6          1 |
  7. |  1    400   White     28          1 |
148. | 22      0   Black    2.1        2.1 |
149. | 22     10   Black    1.9        2.1 |
150. | 22     20   Black      3        2.1 |
     |-------------------------------------|
151. | 22     60   Black    4.8        2.1 |
152. | 22    150   Black    7.4        2.1 |
153. | 22    300   Black   16.7        2.1 |
154. | 22    400   Black   21.2        2.1 |
     +-------------------------------------+

. generate delta_fbf = fbf-baseline
(4 missing values generated)

. label variable delta_fbf "Change in Forearm Blood Flow"

. label variable dose "Isoproterenol Dose (ng/min)"

. generate plotdose = dose

. replace plotdose=6    if dose == 0
(22 real changes made)

. label variable plotdose "Isoproterenol Dose (ng/min)"

. generate logdose = log(dose)
(22 missing values generated)

. label variable logdose "Log Isoproterenol Dose"

. *
. *  Save long format of data for subsequent analyses
. *
. save C:\wddtext\11.2.Long.Isoproterenol.dta, replace
file C:\wddtext\11.2.Long.Isoproterenol.dta saved

. *
. *  Generate Figure 11.1
. *
. more

. collapse (mean) fbfbar=fbf (sd) sd=fbf, by(race plotdose)

. generate blackfbf = .
(14 missing values generated)

. generate whitefbf = .
(14 missing values generated)

. generate whitesd = .
(14 missing values generated)

. generate blacksd = .
(14 missing values generated)

. replace whitefbf = fbfbar if race==1
(7 real changes made)

. replace blackfbf = fbfbar if race==2
(7 real changes made)

. replace blacksd = sd if race==2
(7 real changes made)

. replace whitesd = sd if race==1
(7 real changes made)

. label variable whitefbf "Forearm Blood Flow (ml/min/dl)"

. label variable blackfbf "Forearm Blood Flow (ml/min/dl)"

. generate wsdbar = whitefbf-whitesd
(7 missing values generated)

. generate bsdbar = blackfbf-blacksd
(7 missing values generated)

. replace wsdbar = whitefbf+ whitesd if plotdose < 20
(2 real changes made)

. more

. twoway connected whitefbf plotdose, color(red)          ///
>     || rcap whitefbf wsdbar plotdose, color(red)        ///
>     || connected blackfbf plotdose, color(blue)         ///
>     || rcap blackfbf bsdbar plotdose, color(blue)       ///
>     ||, ytitle(Forearm Blood Flow (ml/min/dl))                           ///
>         legend(ring(0) position(11) col(1) order(1 "Whites" 3 "Blacks")) ///
>         xtitle(Isoproterenol Dose (ng/min)) xscale(log)                  ///
>         xlabel(6 "0" 10 20 30 60 100 150 200 300 400, angle(45))         ///
>         xmtick(40(10)90 250 350) 

. more

. *
. *       Plot individual responses for white and black patients
. *
. use C:\wddtext\11.2.Long.Isoproterenol.dta, clear

. sort id plotdose

. twoway connected fbf plotdose, connect(L)  xscale(log)                          ///
>     xlabel(6 "0" 10 20 30 60 100 150 200 300 400, angle(45))            ///
>     xmtick(40(10)90 250 350) ylabel(0(5)40, angle(0))                   ///
>     ytitle(Forearm Blood Flow (ml/min/dl)) by(race)

. log close
      name:  <unnamed>
       log:  C:\MyDocs\MPH\LectureNotes\ClassDoLogData\Analysis of Variance\11.2.Isoproterenol.log
  log type:  text
 closed on:   7 Feb 2010, 11:37:28
----------------------------------------------------------------------------------------------------------