如何用R生成服从任意有限离散分布的随机数

myfun2=function(n,x,prob)
#x=(x_1,x_2,...)
#n:sample size
#prob=(p_1,p_2,...)
  { result=rep(0,n)
    p=cumsum(prob)
   
    for(i in 1:n)
    {
      u=runif(1)
      kk=min(which(p-u>0))
      result[i]=x[kk]
    }
  
    return(result)
}

qoiqpwqr 发表于 2011-5-22 19:12
drnd

为什么我复制粘贴过去  代码有问题呢，求大神讲解