search for: tcen

...for simply using optim function and Newton Raphson #### dat2<-c(split(data[1:2],data$Censored==1)) # Split the data(cen/uncen) # tun<-c((dat2$'TRUE')$Ti) # Life times data set for uncensored # Stat.uncen<-c((dat2$'TRUE')$Censored) # uncensored Status data set # tcen<-c((dat2$'FALSE')$Ti) # Life times data set for censored # Stat.cen<-c((dat2$'FALSE')$Censored) # censoring Status data set # ##################### Estimation using optim function ################ ml= function(par){mu=par[1] s=par[2] -sum(dlnorm(tun,mu,s,log=TRUE))...

Hi, I am trying to get the maximum likelihood estimator for lognormal distribution with censored data;when we have left, interval and right censord. I built my code in R, by writing the deriving of log likelihood function and using newton raphson method but my estimators were too high " overestimation", where the values exceed the 1000 in some runing of my code. is there any one can