[R] Setting max values for rpois

About this list Date view Thread view Subject view Author view Attachment view

From: Bret Collier (bacolli@uark.edu)
Date: Thu 13 May 2004 - 07:11:08 EST


Message-id: <5.2.1.1.0.20040512160017.00b045d0@mail.uark.edu>

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

______________________________________________
R-help@stat.math.ethz.ch 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


About this list Date view Thread view Subject view Author view Attachment view

This archive was generated by hypermail 2.1.3 : Mon 31 May 2004 - 23:05:09 EST