[R] Generating a new matrix using rbinom and a matrix of probabilities.

From: Economics Guy <economics.guy_at_gmail.com>
Date: Tue, 11 Mar 2008 11:15:44 -0400

I am having a little trouble getting R to do something without writing a couple of very awkward loops.

I have a matrix of probabilities and I want to generate a new matrix of ones and zeros where each element in the new matrix is the result of a draw from a binomial distribution where the probability of getting a 1 is the corresponding element in the matrix of probabilities.

Example Code:

## First I generate the matrix of probabilities for example purposes.

probMatrix <- matrix(NA,5,5){

for (i in 1:5)

probVectorI <- runif(5,0,1)

probMatrix[i,] <- probVectorI


# Now I want to take each element in probMatrix and use it as the
probability parameter in rbinom draw and generate a new matrix. Something like this:

binomialMatrix <- rbinom(1,1,probMatrix)

# But that does not work. I know I can run a loop across each vector
of the matrix, but this seems like an bad way to do this.

---------------------------End Code--------------------------------

So any help would be appreciated.



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. Received on Tue 11 Mar 2008 - 15:30:12 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 11 Mar 2008 - 16:30:21 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive