Geometric mean ≈ MLE of median [Regulatives / Guidelines]

posted by Helmut Homepage – Vienna, Austria, 2024-11-14 10:54 (190 d 18:43 ago) – Posting: # 24275
Views: 5,523

Hi BEQool,

❝ ❝ Well, the geometric mean is the maximum likelihood estimator of the median.

❝ Is geometric mean always equal to median when we have "perfect" log-normal distribution of data set (very large sample)?

Yes. Try a simple simulation:

set.seed(20241114)
mue  <- 1
CV   <- 0.2
n    <- 1e8
x    <- rlnorm(n = n, meanlog = log(mue), sdlog = sqrt(log(CV^2 + 1)))
# instead of prod(x)^(1/n) we use exp(mean(log(x))) because with
# very large n we will underrun the numeric resolution of the machine

y    <- data.frame(estimator = c("arithmetic mean", "geometric mean", "median"),
                   estimate  = c(mean(x), exp(mean(log(x))), median(x)))
y$RE <- sprintf("%+.5f%%", 100 * (y$estimate - mue) / mue)
print(y, row.names = FALSE)

       estimator estimate        RE
 arithmetic mean 1.019806 +1.98057%
  geometric mean 1.000005 +0.00046%
          median 1.000012 +0.00115%


Dif-tor heh smusma 🖖🏼 Довге життя Україна! [image]
Helmut Schütz
[image]

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,681 registered users;
40 visitors (0 registered, 40 guests [including 29 identified bots]).
Forum time: 06:38 CEST (Europe/Vienna)

There are many questions which fools can ask
that wise men cannot answer.    George Pólya

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5