我知道如何使用fitdistr和ks.test来检查Weibull和指数。
但这些不适用于Gompertz。
我想我需要使用maxlik和rgompertz以及ks.gompertz,但我不明白如何使用maxlik来查找alpha和theta的估计值。
我在maxlik pdf中找到了这个:
估计指数分布的参数
t <- rexp(100, 2)
loglik <- function(theta) log(theta) - theta*t
gradlik <- function(theta) 1/theta - t
hesslik <- function(theta) -100/theta^2
Estimate with analytic gradient and hessian
a <- maxLik(loglik, gradlik, hesslik, start=1)
summary( a )
如何更改Gompertz?