Displaying 1 result from an estimated 1 matches for "lt_lognormal".
2009 Aug 07
0
Fitting Truncated Distribution
...ay data below/above a
certain threshold are removed from the observations and I do not know
any longer that there are any observations ruled out by a threshold).
(1) Consider the following snippet (along the lines of the suggestion in
the link above)
rm(list=ls())
library(MASS)
set.seed(1234)
lt_lognormal <- function(x, meanlog, sdlog ){
dlnorm(x, meanlog , sdlog )/plnorm(0.5, meanlog , sdlog )
}
my_seq <- rlnorm(10000)
my_fit <- fitdistr(my_seq,"lognormal")
cut_low <- my_seq[which(my_seq>0.5)]
my_fit_low <- fitdistr(cut_low,lt_lognormal,list(meanlog=0.2,sdlog=0.7)...