如何用R程式编写有关从二项分布中抽样本？
能否替我看看这程式中有什么问题？

谢谢

prob <- function(n,p)
        {
                u <- runif(n, 0, 1)
                x <- 0
                prob_sum <- 0
       
               
                        for (i in 0:n)
                        {
                                if (x <= n)
                                {
                                        combination <- factorial(n)/(factorial(x)*factorial(n-x))
                                        prob <- combination * (p ^ x) * ((1 - p) ^ (n - x))
                                        prob_sum <- prob_sum + prob
                       
                                        ifelse (prob_sum >= u[i], x, x <- x + 1)
                                        print (x)
                                }
                        }               
        }