BookmarkSubscribeRSS Feed
OceanDream
Fluorite | Level 6

Hi guys,

 

I'm trying to calculate a new variable WANT(n), based on the previous observation's variable WANT(n-1).

 

For example, if n=1, WANT=HAVE/2;

else if n>1 then WANT(n)=WANT(n-1)*0.75+HAVE/2;

 

I have no idea which code to start. Could anyone give me some idea? Thanks in advance.

7 REPLIES 7
ChrisHemedinger
Community Manager

I'm not sure if your subject (geodesic distance) is what you meant in your question.  If so: 

 

Check out this communities article: How to calculate geodesic distance in SAS by @MikeZdeb!

 

But it sounds like you need an approach that uses the LAG function (to store the value from a current obs to calculate a new value in the next record).  There are lots of papers and articles about this -- here is one.

It's time to register for SAS Innovate! Join your SAS user peers in Las Vegas on April 16-19 2024.
OceanDream
Fluorite | Level 6

Chris,

 

Thank you for your reply. I changed my title. It was automatically filled in by sas when I typed my title.

 

My question is not only a LAG function question. I think it's more of a LOOP function question, but I don't know how to use it, since the lagged variable I need to use to construct my new variable also need to be constructed the same way, because it's the same new variable, if this makes sense.

 

Regards,

Haimeng

ChrisHemedinger
Community Manager

I think that the implicit loop in the DATA step takes care of this for you.  Here's a simple example.

 

data test;
  set sashelp.class(rename=(weight=have) keep=weight);  
  prev_have=lag(have);
  if _n_ = 1 then want=have/2;
  else want=prev_have*0.75 + have/2;
run;

Run that and check the math to see if it's what you want.  Note that you must always load up prev_have using the lag function, unconditionally .  That makes the prev obs value available for calculation in the next iteration.

It's time to register for SAS Innovate! Join your SAS user peers in Las Vegas on April 16-19 2024.
stat_sas
Ammonite | Level 13

Please provide a sample data set. Based on your requirement this may be you are looking for

 

data have;
input x;
datalines;
2
4
6
8
;

 

data want;
set have;
y=lag(x)*0.75+x/2;
if _n_=1 then y=x/2;
run;

OceanDream
Fluorite | Level 6

stat_sas, what I need is Y=lag(Y)*0.75+x/2, not y=lag(x)*0.75+x/2. That's my problem 😞

Reeza
Super User

Post sample data.

 

However, I think all you need is a retain.

 

Add

 

retain y;


***rest of code;

y=y*0.75+x/2;

 

stat_sas
Ammonite | Level 13

As mentioned earlier by Reeza, retain statement will be required to get previous value of y for calcuation. Please try the following code:

 

data want;
set have;
retain y;
if _n_=1 then y=x/2;
else y=y*0.75+x/2;
run;

sas-innovate-2024.png

Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.

Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.

 

Register now!

What is Bayesian Analysis?

Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.

Find more tutorials on the SAS Users YouTube channel.

Click image to register for webinarClick image to register for webinar

Classroom Training Available!

Select SAS Training centers are offering in-person courses. View upcoming courses for:

View all other training opportunities.

Discussion stats
  • 7 replies
  • 1123 views
  • 4 likes
  • 4 in conversation