我是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_sold
和sold
?