topic Geometric mean in log base 10 in SAS Procedures

Geometric mean in log base 10

amoreau1 — Thu, 07 Jun 2012 07:33:29 GMT

I need to calculate the geometric mean of data.

These data are first transformed in log10.

Then I do the arithmetic mean of these log10 values and then i do the exponential of the arithmetic mean.

Do you confirm I am right?

(I find very different values from arithmetic mean so I am not sure I do the calculation correclty).

Moreover, someome told me that I should do GM = exponential of the arithmetic mean * log(10) (log function is the natural logarithm in SAS) and not just GM = exponential of the arithmetic mean

Thanks.for your help

Re: Geometric mean in log base 10

SteveDenham — Thu, 07 Jun 2012 11:30:20 GMT

The backtransformation is not exp(), but rather 10**(), when the transformation is log10. The expression log(10)*exp() will give the same value.

Steve Denham

Re: Geometric mean in log base 10

Ksharp — Fri, 08 Jun 2012 01:52:30 GMT

SAS has a function GEOMEAN() to do this way . Do you need it ?

Re: Geometric mean in log base 10

amoreau1 — Fri, 08 Jun 2012 08:56:59 GMT

Dear Steve,

Thanks for your reply.

Just to clarify, if I do 10**(mean_of_log_values), it is not equal to log(10)*exp(mean_of_log_values).

I found the same result if I do exp(mean_of_log_values * log(10).

So, what is the good formula?

Thanks.

Re: Geometric mean in log base 10

SteveDenham — Fri, 08 Jun 2012 09:48:03 GMT

My mistake.

Let X = mean of log10 values.

Then geometric mean = 10**X

OR geometric mean = exp(log(10)*X), where log is the natural log of 10 = 2.303 (approximately).

If the variables in question are x1 to xN, then Ksharps recommendation of geomean(x1, x2,..., xN) is the way to go. Note that this is invalid for any x less than zero, and returns zero if any of the x's are equal to zero. This may require transposing the data if the x values are on separate records.

Steve Denham