turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

- Home
- /
- SAS Programming
- /
- Base SAS Programming
- /
- adapting moving average code for additional dimens...

Topic Options

- RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

03-28-2017 03:57 PM

Looking at the code taken from here to calculate the moving average for one dimension/category patient:

```
data ds1;
do patient='A','B','C';
do month=1 to 7;
num=int(ranuni(0)*10);
output;
end;
end;
run;
proc sort;
by patient;
run;
%let n = 4;
data ds2;
set ds1;
by patient;
retain num_sum 0;
if first.patient then do;
count=0;
num_sum=0;
end;
count+1;
last&n=lag&n(num);
if count gt &n then num_sum=sum(num_sum,num,-last&n);
else num_sum=sum(num_sum,num);
if count ge &n then mov_aver=num_sum/&n;
else mov_aver=.;
run;
title 'Moving average within BY-Group';
proc print;
run;
```

I am wondering if this could be adapted for an additional dimenension/category Cat1. A starting point would be:

data ds1; do patient='A','B','C'; do Cat1='C1','C2'; do month=1 to 7; num=int(ranuni(0)*10); output; end; end; end; run; proc sort; by patient Cat1; run;

Any help would be very much appreciated. Thanks a lot.

- Mark as New
- Bookmark
- Subscribe
- RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

Posted in reply to csetzkorn

03-28-2017 04:08 PM

1. Change BY statement in SORT (done & shown) and Data Step (DS2)

2. Any time you see a reference to first.patient change it to first.cat1

I think that's all...but untested.

csetzkorn wrote:

Looking at the code taken from here to calculate the moving average for one dimension/category patient:

`data ds1; do patient='A','B','C'; do month=1 to 7; num=int(ranuni(0)*10); output; end; end; run; proc sort; by patient; run; %let n = 4; data ds2; set ds1; by patient; retain num_sum 0; if first.patient then do; count=0; num_sum=0; end; count+1; last&n=lag&n(num); if count gt &n then num_sum=sum(num_sum,num,-last&n); else num_sum=sum(num_sum,num); if count ge &n then mov_aver=num_sum/&n; else mov_aver=.; run; title 'Moving average within BY-Group'; proc print; run;`

I am wondering if this could be adapted for an additional dimenension/category Cat1. A starting point would be:

data ds1; do patient='A','B','C'; do Cat1='C1','C2'; do month=1 to 7; num=int(ranuni(0)*10); output; end; end; end; run; proc sort; by patient Cat1; run;Any help would be very much appreciated. Thanks a lot.

- Mark as New
- Bookmark
- Subscribe
- RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

Posted in reply to Reeza

03-28-2017 04:20 PM

Thanks. Are you sure the sort has only to be done for one dimension? Please note that the data generation code is just an example and the data can be in any order. Ultimately the MA is to be applied at the time dimension month.

- Mark as New
- Bookmark
- Subscribe
- RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

Posted in reply to csetzkorn

03-28-2017 04:34 PM

What do you mean the sort is by one dimension? Your proc sort showed two variables in the BY statement. That's what you need. In both BY statements, the data step and sort.

And the fastest way is to answer this is to test it out