topic Formula calculation in SAS Procedures

Formula calculation

user1 — Sat, 24 Mar 2018 22:08:27 GMT

Hi,

Can anyone help the calculation of Y below using SAS code?

Y=A*(B)^-1+ A*(B)^-2 + A*(B)^-3 + A*(B)^-4 +……+ A*(B)^-N

Thanks a lot!

Re: Formula calculation

Astounding — Sat, 24 Mar 2018 23:43:22 GMT

Within a DATA step, and assuming you have variables A, B, and N:

y=0;

do k=1 to n;

y + a*b**(-k);

end;

If B is actually an expression rather than a variable, put it in parentheses:

y + a*(b)**(-k);

Re: Formula calculation

PGStats — Sun, 25 Mar 2018 03:31:40 GMT

Or this, which should be faster:

    u = 0;
    v = 1/b;
    do k = 1 to n;
        u = (1 + u) * v;
        end;
    y = a*u;

Re: Formula calculation

user1 — Tue, 27 Mar 2018 00:23:01 GMT

Thanks for you both for providing the answer, much appreciated! PG Stats, since the power is negative, should the sas code be revised?

as u = 0;
    v = 1/b;
    do k = -1 to -n; u = (1 + u) * v; end;
    y = a*u;

Thanks and regards.

Re: Formula calculation

PGStats — Tue, 27 Mar 2018 20:26:38 GMT

No. k is only a counter, not the exponent. The exponent comes from the number of multiplications by v (= 1/b, = b**-1), i.e. the negative exponent results from multiplying by a negative power of b.

Re: Formula calculation

user1 — Wed, 28 Mar 2018 02:33:22 GMT

Thanks again, PG Stats!!

Re: Formula calculation

user1 — Wed, 28 Mar 2018 03:14:01 GMT

Thanks! How about if N has one decimal place, e.g. 14.1?

Re: Formula calculation

Astounding — Wed, 28 Mar 2018 09:32:36 GMT

What is the first term in the equation if N has one decimal place?

Re: Formula calculation

user1 — Wed, 28 Mar 2018 19:56:46 GMT

In the equation, everything is the same, only if N is not an integer, the example is below and I need to calculate for each record (Y₁........Y_N)

Y₁=A*(B)^-1+ A*(B)^-2 + A*(B)^-3 + A*(B)^-4 +……+ A*(B)^-10.2

Y₂=A*(B)^-1+ A*(B)^-2 + A*(B)^-3 + A*(B)^-4 +……+ A*(B)^-20.3

^…….

Y_N=A*(B)-1+ A*(B)-2 + A*(B)-3 + A*(B)-4 +……+ A*(B)^-N.5

Re: Formula calculation

PGStats — Wed, 28 Mar 2018 20:06:08 GMT

Not clear. What is the term before the last?

Re: Formula calculation

user1 — Thu, 29 Mar 2018 02:31:34 GMT

Thanks very much, PG!

the term before the last is: A*(B)^-N+1

For example,

Y₁=A*(B)^-1+ A*(B)^-2 + A*(B)^-3 + A*(B)^-4 +……+ A*(B)^-9 +A*(B)^-10.2

When I run the above sas code, the decimal 0.2 from the last term is missing.

Re: Formula calculation

user1 — Thu, 29 Mar 2018 02:35:27 GMT

sorry, it is not -N+1, it should be integer -N+1, for example if N is 10.2, then the term before the last one should be A*(B)^-9 , not 9.2.