checking goodness of fit for Gompertz distribution in R -
i know how check weibull , exponential using fitdistr , ks.test.
but don't work gompertz.
i think need use maxlik , rgompertz , ks.gompertz don't understand how use maxlik find estimates of alpha , theta.
i found in maxlik pdf:
estimate parameter of exponential distribution
t <- rexp(100, 2) loglik <- function(theta) log(theta) - theta*t gradlik <- function(theta) 1/theta - t hesslik <- function(theta) -100/theta^2 estimate analytic gradient , hessian <- maxlik(loglik, gradlik, hesslik, start=1) summary( )
how change gompertz?
Comments
Post a Comment