stochastic-process - 错误：找不到对象

Question

我是 R 新手，我正在模拟一个实验来应用一些理论结果。实验是这样的：在拉里仓库，顾客的数量是泊松分布的，参数为lambda。顾客购买西装的概率为p。k出售套房的概率是多少？

这是我建议的代码，但有问题。请帮忙！

 Larry<-function(reps=10000, p=1/5, lambda=5, k=2){
 sold<-c(rep(0, reps))
 k_sold<-0
 r<-rpois(reps,lambda)
 for(i in 1:reps){
   for(j in 1:r[i]){
   sold[i]<-sold[i]+rbinom(1,1,p)   
                   }
   k_sold<- if(sold[i]==k){ k_sold+1}   
                }  
    }
  cat("The probability that there are", k ,"clothes sold in Larry's store is", k_sold/reps, "\n")
  return(list(sold))
   }

为什么 R 说对象k_sold并sold没有找到？