Is there Harmonic Number in SAS-IML?

Salah — Sat, 31 Jul 2021 13:25:32 GMT

Hello

Does SAS have HarmonicNumber? When I looked this up, I only found "Harmean" function, but this calculates the mean of vector of size n.

I tried to take n=1, unfortunately, this didn't work.

In my problem, I derived the log-likelihood function using Mathematica and then I used the output in my code (see pic below).

Thank you

Re: Is there Harmonic Number in SAS-IML?

Rick_SAS — Sat, 31 Jul 2021 13:41:56 GMT

I hesitate to correct Mathematica, but I thought Harmonic Numbers were only defined for integers. The arguments do not look like integers to me.

I don't think SAS has a built-in function for harmonic numbers. However, the harmonic numbers are just the sum of the reciprocals of integers, so you can define the following function:

proc iml;
start HarmonicNumber(n);
   return sum(1/(1:n));
finish;

/* test it: The 10th harmonic number is 7382/2520 */
H = HarmonicNumber(10);
print H[format=FRACT.];

topic Is there Harmonic Number in SAS-IML? in SAS/IML Software and Matrix Computations

Is there Harmonic Number in SAS-IML?

Re: Is there Harmonic Number in SAS-IML?