Top Award and First Place Best Presentation of Data Lan Tran-La. Scios Nova, Inc. BLOOD PRESSURE AND HEART RATE vs TIME

Transcription

1 Top Award and First Place Best Presentation of Data Lan Tran-La Scios Nova, Inc. BLOOD PRESSURE AND HEART RATE vs TIME Vital signs were collected before, during, and after the infusion of Drug A. At the end of the study, individual patlenf plots for the vital signs were produced simultaneously with Drug A infusion. The Idea was to show how the vital signs change over time relstlve to the drug infusion. This program used two procedures in SASIGRAPH. PROC GPLOT was used to generate both the upper and lower plots. Then, PROC GREPLAY was implemented to combine the 2 plots together. Postscript 825 printer is the device being used. 1441

2 BLOOD PRESSURE AND HEART RATE vs TIME :E 11 0 :E D Q 0 0,.., ~*flft "" BLOOD PRESSURE.t. "" HEART RATE " " 80 'd 3: ;::: 0.20 H ~ DO.08 Q r TIME J 9 12 ( HOURS) SINCE START OF DRUG A 1442

3 DATA IN VITALS DATA SET DATA IN VITALS DATA SET OBS SBP DBP PULSE HOUR OBS SBP DBP PULSE HOUR DATA IN DOSE DATA SET OBS DOSE HOUR ", h , ,;-.: r l' <, " f ~ ~ '} '. j-, ~' I c ~; , ;; ~ :f f , i~ I f.'

4 /*****************************************************************************/ f* *f f* Create 2 plots: the upper plot SBP, DBP, and heart rate vs. time *f f* the lower plot Drug A vs. time *f f* *f f* Coded by Lan Tran-La *f f* *f f* Input data: vitals *f f* dose *f f* *f /*****************************************************************************1 ** include the macro to get the min and max for the axis ranges ** ; %include "axislim.sas" ; X "purgex vitals.gsf" ; x "renamex vitals.gsf; vitals.gsfbak" filename gsasfile "vitals. gsf" ; ** initialize the graphic output **; goptions reset (axis goptions legend note pattern symbol footnote title) hsize 6 in vsize 9 in vpos = 90 hpos = 72 horigin = 1.5 in vorigin = 1 in autofeed gaccess = gsasfile gsfmode = append nodi splay f* Postscript 825 printer *f device = qmsps820 ** copy vitals and dose data sets into work directory ** libnarne ds I' ; proc copy in = ds out = work memtype = data, select vitals dose ; proc sort data by hour ; run; vitals proc sort data = dose by hour dose ** create data set for finding max ranges ** data minmax ; set vitals dose ** determine the range of the axes for the 2 plots ** proc means noprint data = minmax var hour sbp dbp pulse dose ; 1444

5 output out = range min = minhour minsbp mindbp minpu1se mindose max = maxhour maxsbp maxdbp maxpu1se maxdose data null set range; minbp = min(minsbp, mindbp) ; maxbp = max(maxsbp, maxdbp) ; %axis1im(minhour, maxhour, -6, 42, 3) ; %axis1im(minbp, maxbp, 20, 160, 10) ; %axis1im(minpu1se, maxpu1se, 40, 140, 10) %axislim(mindose, maxdose, 0,.28,.04) ; call symput("minhour", put (minhour, 8.)) call symput ("maxhour", put (maxhour, 8.)) ; call symput("minbp", put(minbp, 8.)) ; call symput ("maxbp", put (maxbp, 8.)) ; call symput ("minpu1se", put (minpulse, 8.)) call symput("maxpu1se", put(maxpulse, 8.)) call symput ("mindose", put (mindose, 8.2)) ; call symput("maxdose", put (maxdose, 8.2)) ; ** define symbols ** ; symbol1 color = yellow line = 1 i = join ; symbo12 color = yellow line = 1 i = join ; symbo13 color = black line 1 i join font special value L height 1.5 symbo14 color = black line = 1 i = stepj ; ** create annotate dataset to shade the blood pressure area %annomac ; data ascend (keep = hour bp) set vitals bp = sbp ; **., data descend (keep set vitals bp=dbp; proc sort data descend by descending hour ; data shade set ascend descend hour bp) length xsys ysys hsys $1 function $8 x y 8 size line angle rotate 8 style color $8 position $1 %system(2, 2, 2) x = hour ; y = bp ; if n eq 1 then do ; %poly(x, y, yellow, solid, 1) 1445

6 else do %polycont(x, y, yellow) end ******************************************** ** sets the origin and axes for the plots ** ******************************************** ** x-axis of the upper plot ** ; axis1 label none order (&minhour to &maxhour by 3) major (width = 4) minor none offset = (0.1 in) value = (height = 1) origin (0.75 in, 4 in) length = 4.5 in ** left y-axis of the upper plot ** axisll label = ( angle = 90 justify = center "BLOOD PRESSURE (MM HG)" ) order (&minbp to &maxbp by 10) ~jor (width = 4) ffilnor (number = 1 width = 3) offset = (0.1 in) value = (height = 1) origin (0.75 in, 4 in) length = 4 in ** right y-axis of the upper plot ** ; axis12 label = ( angle = -90 justify = center "HEART RATE (BPM) " ) order = (&minpulse to &maxpulse by 10) ~jor (width = 4) ffilnor = (number = 1 width = 3) r offset (0.1 in),; value = (height = 1) origin (5.25 in, 4 in) length = 4 in ** x-axis of the lower plot ** ; {- axis2 label (justify = center "TIME (HOURS) SINCE START -OF DRUG A") order (&minhour to &maxhour by 3) major (width = 4) minor none offset = (0.1 in) value = (height = 1) origin (0.75 in,lin) length = 4.5 in ** left y-axis of the lower plot ** axis21 label ( angle = 90 justify = center "DRUG A (UNIT) " ) 1446

7 order (&mindose to &maxdose by 0.04) ~jor (width = 4) minor (number = 1 width = 3) offset = (0.1 in) value = (height = 1) origin = (0.75 in,lin) length = 2 in title1 height title2 height generate upper plot proc gplot data = vitals gout = work run plot sbp * dbp * plot2 pulse note move = height move = font = height font 2 " " ; 3 "B100D PRESSURE AND HEART RATE vs TIME" hour 1 hour 2 / annc shade haxis axis1 vaxis axisll overlay ; * hour = 3 / vaxis axis12, (1.25 in, in) = 2.5 color = yellow font = marker "U" (1.5 in, 8 in) none height = 1 color = black " = BLOOD PRESSURE = 1. 5 font = special "1" none height = 1 n = HEART RATE" " ** generate lower plot ** proc gplot data = dose gout = work plot dose * hour = 4 / haxis vaxis axis2 axis21 ** combine plots ** goption display proc greplay igout = work tc = work nofs ; /* set up template tdef newplot des 1 / def ; template newp10t treplay 1 gplot 1 : gplotl quit ; slzlng */ 'New Plot' 1447

8 /***********************************************i***** ************************/ /* Macro AXISLIM.SAS */ /* Macro to set the min, max, and tick mark range for the axis. */ / * parameters: /* varmin The variable for which you are finding the minimum value. */ /* varmax The variable for which you are finding the maximum value. */ /* min The minimum value as defined in the specs. If the actual */ /* data value is greater than this given value, this value */ /* will be the minimum value ; otherwise, the minimum value */ /* will be calculated depending on the increment of the tick */ /* mark. */ /* max The maximum value as defined in the specs. If the actual */ /* data value is less then this given value, this value will */ /* be the maximum value ; otherwise, the maximum value will */ /* be calculated depending on the increment of the tick mark. */ /* tickrang = Constant tick mark or variable for which you are finding */ /* the tick mark. */ /*****************************************************************************/ %macro axislim(varmin, varmax, min, max, tickrang) ; if (&varmin ne.) and (&varmax ne.) then do if (&varmax lt &max) then do ; if &varmin It &min then range = &max - &varmin else range = &max - &min else do ; if &varmin gt &min then range = &varmax - &min else range = &varmax - &varmin ; minmax = &max - &min ; %if %substr(&tickrang, 1, 1 ) ge A %then %do ; &tickrang = ceil(range / minmax) * &tickrang % if (&varmin ge &min) then &varmin = &min ; else do ; if (mod (&varmin, &tickfang) ne 0) then do ; &varmin = &varmin - mod(&varmin, &tickrang) - &tickrang if (&varmax 1t &max) then &varmax = &max ; else do ; if (mod (&varmax, &tickrang) ne 0) then do ; &varmax = &varmax - mod(&varmax, &tickrang) + &tickrang if (mod (&varmin, &tickrang) ne 0) then do ; &varmin = floor(&varmin / &tickrang) * &tickrang if (mod (&varmax, &tickrang) ne 0) then do ; &varmax = ceil(&varmax / &tickrang) * &tickrang else do ; put "WAR" "NING missing &varmax and/or &varmin &varmin &min &varmax = &max ; %mend axislim ; II * / 1448