Bret, the following should do:
x <- sum(pmin(rpois(recruit.f, l),3)) Giovanni > Date: Wed, 12 May 2004 16:11:08 -0500 > From: Bret Collier <[EMAIL PROTECTED]> > Sender: [EMAIL PROTECTED] > Precedence: list > > R-users, > I am simulating a birth process for 4 classes of individuals with > l[i] being the average No. fetuses per individual. However, I need to > bound the resulting values for each generated rpois to be <=3 (no > individual can have > 3 offspring). I have not been able to figure out how > to incorporate this into the below example. Any suggestions on integrating > would be appreciated. > > > recruit.f <- c(12, 12, 25, 51) #No. females in each age class > l <- c(.05, 1.22, 1.6, 1.8) #mean No. fetuses for each age class > x <- sapply(lapply(1:4, function(i) rpois(recruit.f[i], l[i])), sum) > > TIA, > > Bret A. Collier > Arkansas Cooperative Fish and Wildlife Research Unit > University of Arkansas > > ______________________________________________ > [EMAIL PROTECTED] mailing list > https://www.stat.math.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html > > -- __________________________________________________ [ ] [ Giovanni Petris [EMAIL PROTECTED] ] [ Department of Mathematical Sciences ] [ University of Arkansas - Fayetteville, AR 72701 ] [ Ph: (479) 575-6324, 575-8630 (fax) ] [ http://definetti.uark.edu/~gpetris/ ] [__________________________________________________] ______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html