topic Re: How to calculate geodesic distance in SAS in SAS Procedures

How to calculate new variable n using the new variable n-1??

OceanDream — Mon, 02 May 2016 00:05:50 GMT

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.

Re: How to calculate geodesic distance in SAS

ChrisHemedinger — Mon, 02 May 2016 00:17:39 GMT

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.

Re: How to calculate geodesic distance in SAS

OceanDream — Mon, 02 May 2016 00:24:56 GMT

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

Re: How to calculate new variable n using the new variable n-1??

ChrisHemedinger — Mon, 02 May 2016 00:35:33 GMT

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.

Re: How to calculate geodesic distance in SAS

stat_sas — Mon, 02 May 2016 00:36:11 GMT

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;

Re: How to calculate geodesic distance in SAS

OceanDream — Mon, 02 May 2016 02:43:12 GMT

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 😞

Re: How to calculate geodesic distance in SAS

Reeza — Mon, 02 May 2016 02:47:48 GMT

Post sample data.

However, I think all you need is a retain.

Add

retain y;


***rest of code;

y=y*0.75+x/2;

Re: How to calculate geodesic distance in SAS

stat_sas — Mon, 02 May 2016 05:23:45 GMT

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;