Your problem is nothing but that of drawing a random sample from a multinomial distribution.  Check out the "Multinomial" function.
 
?Multinomial
 
Best,
Ravi.

________________________________

From: r-help-bounces@stat.math.ethz.ch on behalf of vincent@7d4.com
Sent: Tue 10/4/2005 1:59 AM
To: r-help@stat.math.ethz.ch
Subject: Re: [R] sampling vectors



Eric Pante a écrit :

> Hello Listers,
> I am trying to sample a vector to create a new one of sample length,
> witha  sum equal to the sum of the initial vector:
> initial = 10, 30, 10 (sum=50)
> sample example = 5, 35, 10 (sum=50) or 25, 15, 10  (sum=50), etc ...
> My problem is to control the sum, so it stays constant.

f0 = function()
{
s1 = 50;
for (i in 0:s1)
        {
        s2 = s1 - i;
        for (j in 0:s2)
                {
                s3 = s2 - j;
                print (c(i,j,s3));
                }
        }
}

(If I have well understood the question) ?
hih
Vincent

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html



	[[alternative HTML version deleted]]


