```r true_cdf <- data.frame( delay = 1:26, cdf = pdislnorm(1:26, 2, 0.5) / pdislnorm(26, 2, 0.5), Method = "Data generating" ) ``` It might be worth updating this to be the double censored representation so it is closer to the true data generating dist?