Dear R helpers, I am raising one query regarding this "Binomial" thread with the sole intention of learning something more as I understand R forum is an ocean of knowledge.
I was going through all the responses, but wondered that original query was about generating Binomial random numbers while what the R code produced so far generates the Bernoulli Random no.s i.e. 0 and 1. True Binomial distribution is nothing but no of Bernoulli trials. As I said I am a moron and don't understand much about Statistics. Just couldn't stop from asking my stupid question. Regards Sarah --- On Thu, 5/12/11, David Winsemius <dwinsem...@comcast.net> wrote: From: David Winsemius <dwinsem...@comcast.net> Subject: Re: [R] Binomial To: "Alexander Engelhardt" <a...@chaotic-neutral.de> Cc: r-help@r-project.org, "blutack" <x-jess-...@hotmail.co.uk> Date: Thursday, May 12, 2011, 11:08 AM On May 12, 2011, at 5:02 AM, Alexander Engelhardt wrote: > Am 12.05.2011 10:46, schrieb blutack: >> Hi, I need to create a function which generates a Binomial random number >> without using the rbinom function. Do I need to use the choose function or >> am I better just using a sample? >> Thanks. > > I think I remember other software who generates binomial data with e.g. > pi=0.7 by > > pi <- 0.7 I hope Allan knows this and is just being humorous here, but for the less experienced in the audience ... Choosing a different threshold variable name might be less error prone. `pi` is one of few built-in constants in R and there may be code that depends on that fact. > pi [1] 3.141593 > pi <- 0.7 > pi [1] 0.7 > rm(pi) > pi [1] 3.141593 > x <- runif(100)>pi > summary(x) Another method would be: x <- sample(c(0,1) , 100, replace=TRUE, prob=c(0.7, 0.3) ) -- David Winsemius, MD West Hartford, CT ______________________________________________ R-help@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. [[alternative HTML version deleted]]
______________________________________________ R-help@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.
