R 我想计算密度的最大似然估计量，但我的输出产生NaN'；当我的方程式在数学上是正确的时候

我认为我的代码“几乎正确”，但为什么我会得到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

N=100
est1=功能（m）{
tau=rep（0，m）
b0=代表（0，m）
b1=代表（0，m）
tau[1]=1
b0[1]=1
b1[1]=1
x我认为它们在数学上是不正确的，因为你从来没有定义过
y
！
你应该a）更全面地解释你想做什么，b）特别是解释
y
应该是什么（如果不是打字错误的话）。