From: Hanke, Alex <HankeA_at_mar.dfo-mpo.gc.ca>

Date: Tue, 08 Jul 2008 09:19:49 -0300

Peng Jiang 江鹏 ,Ph.D. Candidate

Antai College of Economics & Management 安泰经济管理学院

Department of Mathematics

数学系

Shanghai Jiaotong University (Minhang Campus) 800 Dongchuan Road

200240 Shanghai

P. R. China

R-help_at_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.

R-help_at_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. Received on Tue 08 Jul 2008 - 12:23:04 GMT

Date: Tue, 08 Jul 2008 09:19:49 -0300

Hi Peng,

Does this help?

sum(diff(c(0,sort(sample(seq(1,99,1),50,replace=T)),100)))
Regards

Alex

-----Original Message-----

From: r-help-bounces_at_r-project.org [mailto:r-help-bounces_at_r-project.org] On Behalf Of Peng Jiang
Sent: July 8, 2008 3:57 AM

To: m_olshansky_at_yahoo.com

Cc: r-help_at_stat.math.ethz.ch

Subject: Re: [R] Sum(Random Numbers)=100

> If they are really random you can not expect their sum to be 100.

*> However, it is not difficult to get that given that the sum of n
**> independent Poisson random variables equals N, any individual one has
**> the conditional binomial distribution with size = N and p = 1/n, i.e.
**> P(Xi=k/Sn=N) = (N over k)*(1/n)^k*((n-1)/n)^(N-k).
**> So you can generate X1 binomial with size = 100 and p = 1/50; if X1 =
**> k1 then the sum of the rest 49 must equal 100 - k1, so now you
**> generate X2 binomial with size = 100-k1 and p = 1/49; if X2 = k2 then
**> generate X3 binomial with size = 100 -(k1+k2) and p = 1/48, etc.
**>
**> Why do you need this?
**>
**>
**> --- On Tue, 8/7/08, Shubha Vishwanath Karanth
**> <shubhak_at_ambaresearch.com
**> > wrote:
**>
**>> From: Shubha Vishwanath Karanth <shubhak_at_ambaresearch.com>
**>> Subject: [R] Sum(Random Numbers)=100
**>> To: r-help_at_stat.math.ethz.ch
**>> Received: Tuesday, 8 July, 2008, 3:58 PM Hi R,
**>>
**>>
**>>
**>> I need to generate 50 random numbers (preferably poisson),
**>> such that
**>> their sum is equal to 100. How do I do this?
**>>
**>>
**>>
**>>
**>>
**>> Thank you,
**>>
**>> Shubha
**>>
**>>
**>>
**>> This e-mail may contain confidential and/or privileged
**>> i...{{dropped:13}}
**>>
**>> ______________________________________________
**>> R-help_at_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.
**>
**> ______________________________________________
**> R-help_at_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.
*

Peng Jiang 江鹏 ,Ph.D. Candidate

Antai College of Economics & Management 安泰经济管理学院

Department of Mathematics

数学系

Shanghai Jiaotong University (Minhang Campus) 800 Dongchuan Road

200240 Shanghai

P. R. China

R-help_at_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.

R-help_at_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. Received on Tue 08 Jul 2008 - 12:23:04 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 Tue 08 Jul 2008 - 12:31:40 GMT.

*
Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help.
Please read the posting
guide before posting to the list.
*