maxLik 函數中的對數似然

Question

我有一個 PDF 文件，其中我從逆方法中得到了它。 和來自 x<-norm(20,0,1) 的 x。 這是來自指數並行的模型，用於使用協變量進行生存分析，其中我讓 lambda= b0+b1*x

PDF = 2 * (b0+b1*x) * exp(-(b0+b1*x)*t) * (1-exp(-(b0+b1*x)*t))

問題是當我執行 maxLik 函數時，maxLik 中的標准誤差值返回 NA 值。 我將 z1 指定為對數似然函數。

#Likelihood
library(maxLik)
LLF <- function(para){
set.seed(1)

b0 = para[1]
b1 = para[2]

#n = 1

z1 = (n*log(2)) + (n*log(b0+b1*xsum)) - ((b0+b1*xsum)*tsum) + (n*log(1-exp((-(b0 + b1*xsum)*tsum))))

return(z1)
}

mle <- maxLik(LLF, start = c(2,4))
summary(mle)

Maximum Likelihood estimation
Newton-Raphson maximisation, 3 iterations
Return code 1: gradient close to zero
Log-Likelihood: -22.7055 
2  free parameters
Estimates:
     Estimate Std. error t value Pr(> t)
[1,]    1.986         NA      NA      NA
[2,]    3.986         NA      NA      NA

先感謝您。

Answer 1

參見返回碼 1：梯度接近於零。 標准。 可能使用梯度/Hessian 計算錯誤。 如果這些不正確，那么這可能是罪魁禍首。

嘗試使用參數 method='Nelder-Mead' 運行您的函數，看看這是否可以緩解問題。 如果是這樣，您可以使用此方法，或提供您自己的 Jacobian 和 Hessian 函數。