Re: [R] rbinom for a matrix

From: Ben Bolker <>
Date: Wed, 09 Jul 2008 15:59:21 -0400

Hash: SHA1

Dylan Beaudette wrote:
| On Wednesday 09 July 2008, Ben Bolker wrote:
|> ACroske <Audy3272 <at>> writes:
|>> I have a large matrix full of probabilities; I would like to convert
|>> probability to a 1 or a 0 using rbinom.
|>> How can I do this on the entire matrix? The matrix was converted from a
|>> raster ArcMap dataset, so the matrix is essentially a map. Because of
|>> this, I have no column headings.
|>> Thanks!
|> How about
|> matrix(rbinom(length(m),prob=m,size=1),nrow=nrow(m))
|> or (perhaps marginally more efficiently?)
|> y <- (runif(m)<m)
|> storage.mode(y) <- "double"
|> Ben Bolker
| Wait a second. Are you trying to convert each probability into the most
| likely '1' or '0' through rounding? The code example above will give
you a
| different answer every time you run it. Is that what you are looking for?
| Just curious,
| Dylan

~ I assumed that since the original poster said "convert each probability to a 1 or 0 using rbinom" that they did indeed want a random assignment, rather than rounding ... -----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.6 (GNU/Linux)
Comment: Using GnuPG with Mozilla -

iD8DBQFIdRiZc5UpGjwzenMRAuLkAKCA6ZhjrVsg5RJGYFGwq+6vz2pWxgCfetYk 9V/ZPGi2jfVGV5jx+LGPs8s=
-----END PGP SIGNATURE----- mailing list PLEASE do read the posting guide and provide commented, minimal, self-contained, reproducible code. Received on Wed 09 Jul 2008 - 20:05:27 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 Thu 10 Jul 2008 - 02:32:02 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.

list of date sections of archive