I'm trying to write a function in R which produces a set of n i.i.d Bernoulli random variables X1,X2,...,Xn ∼ B(1,π).
To generate M = 1000 random samples. When n = 10 and pi = 0.4.
I'm having trouble starting the code. If anyone could help out that would be greatly appreciated.
Here is what I have run to this point:
Taking what you've done so far:
set.seed(101) ## for reproducibility rr <- replicate(1000,rbinom(n=10,prob=.04,size=1))
creates 1000 samples.
dim(rr) ##  10 1000 mean(rr) ##  0.0411
i.e. this is 10 rows x 1000 columns; you can use
t() to transpose it if you like. If you specify
simplify=FALSE in the
replicate() call you'd get a list of vectors of length 10 instead.
(Your question is inconsistent as to whether you want probability of 0.04 or 0.4, but you should be able to sort that out for yourself.)
For this particular question you could also use