Nakagami分布の対数尤度は無限大ですR

Question

[60,80] $の中の私のデータセット$ x \の正規化されたヒストグラムをNakagami分布に当てはめています。まず私は、次のMLEコードをVGAMパッケージのdnakaを使用して、スケールと形状パラメータを推定した：

：次に

ll <- function(par) { 
    if(par[1]>0 & par[2]>0) {return(-sum(log(dnaka(x, scale = par[1], shape = par[2]))))} # m=shape, ohm or spread = scale 
    else return(Inf) 
} 
mle = optim(c(1000,1), ll) 

、私は以下のコードを推定したパラメータに基づいて対数尤度値を推定しています

lik = sum(log(dnaka(x, shape = mle$par[1], scale = mle$par[2]))) 

ただし、対数尤度値likは-Infです。この無限の価値は、PDFの中カ分布の方程式のexp（。）項によるものであることを理解します。 [60,80] $でデータセット$ x \の中上分布の有限log-likelihood値を推定する方法はありますか？ありがとうございました。

Answer 1

元の質問に私のコメントを見てください。ここで

はscale = 1.5とshape = 1

set.seed(2017); 
x <- rnaka(10^4, scale = 1.5, shape = 1); 

ll <- function(par) { 
    if (par[1] >= 0.5 && par[2] > 0) { 
     return(-sum(log(dnaka(x, scale = par[1], shape = par[2])))); 
    } 
    else return(Inf); 
} 

mle <- optim(c(0.5, 1), ll); 

mle$par; 
#[1] 1.4833965 0.9938022 

ll(mle$par); 
#[1] 7946.478