I am trying to simulate a series of ones and zeros (1 or 0) and I am using "rbinom" but realizing that the number of successes expected is not accurate. Any advice out there.

This is the example:

N<-500
status<-rbinom(N, 1, prob = 0.15)
count<-sum(status)

15 percent of 500 should be 75 but what I obtain from the "count" variable is 77 that gives the probability of success to be 0.154. Not very good.

Is there another way beyond using "sample" and "rep" together?


A Smile costs Nothing  

     But Rewards Everything

Happiness is not perfected until it is shared
                                                              -Jane Porter        


       
	[[alternative HTML version deleted]]

