1

確率密度関数 f(x)、累積分布関数 F(x)、および分位関数 Q(u) を使用したワイブル分布のプログレッシブ タイプ II 打ち切りスキームに基づくシミュレーション研究を行っていました。

ここに画像の説明を入力 すべての m 個のプログレッシブ タイプ II 右打ち切り順序統計の同時確率密度関数は、次の形式です。 ここに画像の説明を入力

N Balakrishnan と Rita Aggarwala による本「Progressive Censoring Theory, Methods, and Applications」から使用したアルゴリズムは次のとおりです。 ここに画像の説明を入力

データ生成用の R コードを作成し、次のように "nlm" 関数を使用してその最尤推定値を見つけようとしました。

alpha = 3
beta = 2
m = 20
R = c(replicate(20,1))
c = 40^20
w = runif(m,0,1)
v = c(rep(0,m))
for(i in 1:m){
 sm = 0
  for (j in (m-i+1):m){
    sm = sm+R[j]
  }
 v[i] =w[i]^(1/(i+sm))
 i = i+1
}
u = rep(0,m)
for (i in 1:m ) {
prod = 1
for (j in (m-i+1):m){
  prod = prod*v[j]
  z = prod
}
u[i] = 1-z
i = i+1
}
x = ((-log(1-u))^(1/alpha))/beta
x

L = function(par){
alpha = par[1]; beta = par[2];
NLL = log(c)+m*log(alpha)+alpha*m*log(beta)+(alpha-1)*sum(log(x))- beta^(alpha)*sum(x^(alpha))-(beta^(alpha))*sum(R*x^(alpha))
return(-NLL)
}
initial < -c(alpha,beta)
est = nlm(L,initial,hessian=TRUE)
est

est$estimate[1]
est$estimate[2]

このコードは、言及されたいくつかの特別な条件のアルゴリズムを満たしているため、混乱しています。Morover私は次のようにnlm関数のエラーに直面しています:

"Warning messages: 1: In log(beta) : NaNs produced 2: In nlm(L,
initial, hessian = TRUE) :   NA/Inf replaced by maximum positive
value" ```

誰かがこの問題で私を助けてくれれば、それは素晴らしいことです.

ありがとう

4

0 に答える 0