[R] simulating a toss of a coin

[David Arnold]

All,

I've devised a way to toss a coin three times and record the number of heads, then use the table command to summarize the number of time 0, 1, 2, or 3 heads were obtained.

n <- 100
x <- sample(c("H","T"),3*n,replace=TRUE)
dim(x) <- c(3,n)
num_heads <- apply(x,2,function(x) sum(x=="H"))
table(num_heads)/n

What I'd like to do next is run this code 10 times, and gather the results in a matrix with 10 rows and 4 columns, where the first row records the output of the first run of the above code, the second row records the second run of the above code, etc.

Looking for some simple ideas on how to accomplish this.

Thanks.

David