I am a bit stuck on how to automate this programming issue. I am trying to write a macro to perform this but can't seem to figure out a clean way of doing it so it is generalizable.    The input to this macro is going to be a numeric dataset with ~30-50 observations and anywhere from 2 to 10 variables. This dataset will be structure like this: DATA input; do Item = 1 to 30; Score1 = rand('uniform'); Score2 = rand('uniform'); Score3 = rand('uniform'); output; end; run; But the exact numbers will vary for each such input dataset.   What I need to do is evaluate for each Item (i.e. row) whether or not the set of Scores (columns) meets a condition. The condition is: exactly one of the Scores is >= 0.6 and each other score is <=0.4. Now, for the example dataset, this can be coded pretty easily using something like:   DATA output; set input; if Score1 >= 0.6 & Score2 <= 0.4 & Score3 <= 0.4 then condition = 1; else if Score1 <= 0.4 & Score2 >=0.6 & Score3 <= 0.4 then condition = 1; else if Score1 <= 0.4 & Score2 <= 0.4 & Score3 >=0.6 then condition = 1; else condition = 0; run; The problem of course is that this isn't easily generalizable, this only works if there are exactly 3 variables.  How would I automate/generalize this logic in a macro? I'm pretty stuck, I don't really know what to do. I'm guessing there is a way to do this with arrays but I can't quite figure out how to implement that.  ... View more

Below, I have some code that produces some data and generates a plot:   %let Planned_E_Bt = 1.81; %let Current_E_Bt = 0.82; /* Generate lines */ DATA figure1; do t = 0 to 1 by 0.01; planned = &Planned_E_Bt*t; current = &Current_E_Bt*t; if t >= 0.5 then do; CP_current = current; CP_planned = &Current_E_Bt*0.5 + &Planned_E_Bt*(t-0.5); current = .; end; output; end; run; /* Generate markers and labels */ DATA figure1; set figure1; if round(t,0.001)=0.5 then do; marker = CP_current; laby1 = marker - 0.1; label1 = "B(t=0.5)=" || strip(round(marker,0.001)); end; if round(t,0.001)=0.4 then do; laby2 = planned + 0.2; label2= "E{B(t)}=" || strip(round(&Planned_E_Bt,0.01)) || "t"; end; run; PROC SGPLOT data=figure1 noautolegend; scatter x = t y = marker / markerattrs=(color=black); series x = t y = planned / lineattrs=(color=black pattern=1); series x = t y = current / lineattrs=(color=black pattern=2); series x = t y = CP_current / lineattrs=(color=black pattern=1) curvelabel="CP=0.04" curvelabelloc=outside; series x = t y = CP_planned / lineattrs=(color=black pattern=1) curvelabel="CP=0.41" curvelabelloc=outside; yaxis label = "B(t)"; xaxis label = "t" values=(0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0); run;   This all works as intended. However, I am trying to add some labels as text to the plot. As soon as I add them (see modified SGPLOT code chunk below), it changes the x-axis. It adds white space at either end of the axis (<0 and >1).    PROC SGPLOT data=figure1 noautolegend; scatter x = t y = marker / markerattrs=(color=black); series x = t y = planned / lineattrs=(color=black pattern=1); series x = t y = current / lineattrs=(color=black pattern=2); series x = t y = CP_current / lineattrs=(color=black pattern=1) curvelabel="CP=&Current_CP." curvelabelloc=outside; series x = t y = CP_planned / lineattrs=(color=black pattern=1) curvelabel="CP=&Planned_CP." curvelabelloc=outside; text x = t y = laby1 text = label1; text x = t y = laby2 text = label2; yaxis label = "B(t)"; xaxis label = "t" values=(0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0); run;   How do I stop it from adding this blank space at either end? I can't figure out how to force it to not rescale the axis. Curvelabel doesn't work, because that will only let me add the labels to the beginning or end of the line instead of the middle. I've tried Datalabel, but that also adds this blank space at either end. ... View more

I have a dataset full of times, indicating when an individual went to bed at night and when they woke up the next morning. I want to create a new variable that indicates the number of hours that they were asleep. However, SAS doesn't seem to understand how clocks work, and is unable to "wrap" times around the clock to do the calculation in a sensible way.   For example, if someone went to bed at 9pm and woke up at 5am, we would WANT the value of hours slept to be 8. However, SAS will return +/- 16, since it is only able to calculate the difference in a linear fashion. There are a couple of workarounds to this (included in my code snippet below), but they feel very clunky to me. One of them is pulling in the date variables and using the dhms function and then converting the result into hours, the other is converting the bedtime variable into the # of hours before midnight and then adding that to the wake up time variable.    Here's an example:   DATA time; /* Bed time and wakey time */ bed_time = 75600; format bed_time time8.; wake_time = 18000; format wake_time time8.; /* Wrong answer */ wrong_hours_slept = intck('hour',bed_time,wake_time); /* Clunky solution 1 */ day1 = today(); format day1 date8.; day2 = intnx('day',day1,1); format day2 date8.; seconds_slept = dhms(day2,0,0,wake_time) - dhms(day1,0,0,bed_time); right_hours_slept1 = seconds_slept/3600; /* Clunky solution 2 */ hours_before_midnight = (86400 - bed_time)/3600; right_hours_slept2 = hours_before_midnight + wake_time/3600; run; Are there any other solutions? Obviously both of the above work, and ultimately don't require too much extra coding, but it feels rather frustrating to me that SAS doesn't have a built-in way of manipulating clock time. Further, both of these clunky solutions also have points of failure (notably with people who go to sleep after midnight), requiring some if/then logic to correct.         ... View more