------------------------------------------------------------------------------- log: C:\WDDtext\log\11.2.Isoproterenol.log log type: text opened on: 7 Dec 2001, 17:51:12 . * 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 Observation 1 id 1 race White fbf0 1 fbf10 1.4 fbf20 6.4 fbf60 19.1 fbf150 25 fbf300 24.6 fbf400 28 Observation 22 id 22 race Black fbf0 2.1 fbf10 1.9 fbf20 3 fbf60 4.8 fbf150 7.4 fbf300 16.7 fbf400 21.2 . generate baseline = fbf0 . * . * Convert data from one record per patient to one record per observation. . * . 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 . * . 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) . graph whitefbf blackfbf wsdbar whitefbf bsdbar blackfbf plotdose, xlog xlabel(10,20,60,100,150,200,3 > 00,400) xtick(6,30,40,50,70,80,90,250,300,350) ylabel(0 5 10 15 20 25) connect(llIIII) symbol(OOiiii) > l1title(Forearm Blood Flow (ml/min/dl)) gap(3) . more . * . * Plot individual responses for white and black patients . * . use C:\WDDtext\11.2.Long.Isoproterenol.dta, clear . sort id plotdose . * . * Plot responses for white patients. . * . graph fbf plotdose if race==1, xlog xlabel(10, 20, 30, 60,100,150,200,300,400) xtick(6,30,40,50,70,8 > 0,90,250,300,350) ylabel(0 5 to 40) connect(L) symbol(O) l1title(Forearm Blood Flow (ml/min/dl)) gap( > 3) . more . * . * Plot responses for black patients. . * . graph fbf plotdose if race==2, xlog xlabel(10, 20, 30, 60,100,150,200,300,400) xtick(6,30,40,50,70,80 > ,90,250,300,350) ylabel(0 5 to 40) connect(L) symbol(O) l1title(Forearm Blood Flow (ml/min/dl)) gap(3 > ) . log close log: C:\WDDtext\log\11.2.Isoproterenol.log log type: text closed on: 7 Dec 2001, 17:51:23 -------------------------------------------------------------------------------------------------------