draw many elements julia

 # a vector of frequencies

p=rand(Beta(2,2),10)

# a vector of "distributions"

a=Binomial.(2,p)

# a vector of vectors

 b=rand.(a,1)

# collapsed (row vector)

reduce(hcat,b)

# draw 5 animals at once

b=rand.(a,5)

# collapse

bb=reduce(hcat,b)

#compute frequencies

pest = mean(bb,dims=1)/2