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

Popular posts from this blog

java - unable show chart in xls document using jasper reports -

java.lang.RuntimeException: Could not generate dummy secret at sun.security.ssl.RSAClientKeyExchange.<init> -

javascript - How to change image src on success AJAX -