Message-ID: <5073CFD4.7040207@sapo.pt>
Date: 2012-10-09T07:18:44Z
From: Rui Barradas
Subject: simulating a toss of a coin
In-Reply-To: <FD95C2D9-D66D-47A1-94AA-75A571B4C585@suddenlink.net>
Hello,
Try the following.
fun <- function(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
}
Runs <- 1e1
t(replicate(Runs, fun()))
Hope this helps,
Rui Barradas
Em 09-10-2012 06:04, David Arnold escreveu:
> 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
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.