r语言 - 我想计算密度的最大似然估计器，但是当我的方程在数学上正确时，我的输出会产生 NaN

Question

我认为我的代码"几乎正确"，但是为什么我会得到NaN？

N=100

est1 = function(m){

tau = rep(0,m)
b0 = rep(0,m)
b1 = rep(0,m)
tau[1] = 1
b0[1] = 1
b1[1] = 1
x<- runif(100,0,1)

for(i in 2:m){
tau[i] = N/tau[i-1] + sum(log(y)) -(1 + exp(b0[i-1]+b1[i-1]*x) )*sum(y^tau[i-1]*log(y)/(1+y^tau[i-1]))  
b0[i]= N - log(  1 +y^tau[i-1] )*sum( exp(b0[i-1]+b1[i-1]*x))
b1[i] = (sum(x) )-log(1+y^tau[i-1])*sum( x*exp(b0[i-1]+ b1[i-1]*x ))
}
cbind(tau,b0,b1)
}
output = est1(100)
output

Answer 1

我认为它们在数学上是不正确的，因为你从来没有定义过y！