3
$\begingroup$

I have some data for numbers that are positive and these numbers are capped above at C (so if a sample from my data should have exceeded C, the data generating process will simply return C).

I know which of the samples were capped and which ones were not.

I am happy to assume that these data points are lognormally distributed and would like to estimate the parameters of the uncapped lognormal distribution. How do I use my uncapped data and apply MLE to achieve this? If lognormal is too hard for some reason, I am open to using normal or other distributions.

$\endgroup$
5
  • $\begingroup$ You can just write a loglikelihood function and optimize it numerically. $\endgroup$ Oct 14, 2015 at 2:48
  • $\begingroup$ 1. Take logs. 2. Fit a truncated normal (truncated at $\log(C)$). The parameters are now ML for your lognormal. $\endgroup$
    – Glen_b
    Oct 14, 2015 at 6:37
  • $\begingroup$ There's an algorithm outlined here that's suitable for the normal case and my previous comment gets you the lognormal from that. $\endgroup$
    – Glen_b
    Oct 14, 2015 at 6:46
  • 1
    $\begingroup$ @Glen_b That algorithm is for right truncated data which is not the same as right censored data. $\endgroup$ Jul 29, 2017 at 22:12
  • $\begingroup$ @Jarle you're correct, I was responding to the title (which previously said "truncated), but the description is indeed of a censored problem. I should have typed censored in my first comment and the one at the link doesn't apply either, as you say. Thanks for pointing this out. $\endgroup$
    – Glen_b
    Jul 30, 2017 at 1:02

1 Answer 1

3
$\begingroup$

Since you know that the event $y>C$ occurred for some observations you have right censored data. If such events instead went unrecorded, the data would have been right truncated. Fitting a lognormal, intercept only, survival regression model as follows gives you MLEs of the location and scale parameter of the lognormal:

> library(survival)
> y <- rlnorm(1000, 0, 1) # Simulated data with location 0 and scale 1
> delta <- y<2 # right censoring at C = 2
> y[!delta] <- 2
> survreg(Surv(y, delta) ~ 1, dist="lognormal")
Call:
survreg(formula = Surv(y, delta) ~ 1, dist = "lognormal")

Coefficients:
 (Intercept) 
-0.008081718 

Scale= 0.9564241 

Loglik(model)= -978.1   Loglik(intercept only)= -978.1
n= 1000 
$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.