错误:找不到对象

时间:2017-02-13 07:12:12

标签: stochastic-process

我是R的新手,我正在模拟一个实验来应用一些理论结果。实验是这样的:在Larry仓库,客户数量是泊松分布参数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_soldsold

0 个答案:

没有答案