<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<base href="x-msg://434/">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">
You could try this:
<div><br>
</div>
<div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; ">dvector rmultinom(<span style="color: rgb(187, 44, 162); ">const</span> <span style="color: rgb(187, 44, 162); ">int</span>& seed, <span style="color: rgb(187, 44, 162); ">const</span> <span style="color: rgb(187, 44, 162); ">int</span>&
size,<span style="color: rgb(187, 44, 162); ">const</span> dvector& prob)</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; color: rgb(0, 132, 0); ">
<span style="color: rgb(0, 0, 0); ">{<span style="white-space: pre-wrap; "> </span>
</span>//Returns a multinomial sample, of size n, based on sampling probabilities p.</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; color: rgb(0, 132, 0); ">
<span style="color: rgb(0, 0, 0); "><span style="white-space: pre-wrap; "></span></span>//p is normalized internally, based on the same method employed in R</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><span style="white-space: pre-wrap; "></span>random_number_generator rng(seed);</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><span style="white-space: pre-wrap; "></span><span style="color: rgb(187, 44, 162); ">int</span> i,n,lb,ub;</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; color: rgb(187, 44, 162); ">
<span style="color: rgb(0, 0, 0); "><span style="white-space: pre-wrap; "></span></span>float<span style="color: rgb(0, 0, 0); "> p;</span></div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><span style="white-space: pre-wrap; "></span>lb=prob.indexmin(); ub=prob.indexmax();</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><span style="white-space: pre-wrap; "></span>dvector freq(lb,ub); freq.initialize();</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><span style="white-space: pre-wrap; "></span>dvector P=prob;</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><span style="white-space: pre-wrap; "></span>P/=sum(P);</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><span style="white-space: pre-wrap; "></span>dvector bisamp(<span style="color: rgb(39, 42, 216); ">1</span>,size); bisamp.fill_randbi(P[lb],rng);</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><span style="white-space: pre-wrap; "></span>freq[lb]=sum(bisamp);</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><span style="white-space: pre-wrap; "></span><span style="color: rgb(187, 44, 162); ">for</span>(i=lb+<span style="color: rgb(39, 42, 216); ">1</span>;i<=ub;i++)</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><span style="white-space: pre-wrap; "></span>{</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><span style="white-space: pre-wrap; "></span>n=size-sum(freq);</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><span style="white-space: pre-wrap; "></span>p=P[i]/(<span style="color: rgb(39, 42, 216); ">1.</span>-sum(P(lb,i-<span style="color: rgb(39, 42, 216); ">1</span>)));<span style="white-space: pre-wrap; ">
</span><span style="color: rgb(0, 132, 0); ">//Corrected version</span></div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; min-height: 13px; ">
<br>
</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; color: rgb(0, 132, 0); ">
<span style="color: rgb(0, 0, 0); "><span style="white-space: pre-wrap; "></span></span>//cout<<ub-i<<endl;</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><span style="white-space: pre-wrap; "></span>dvector bisamp(<span style="color: rgb(39, 42, 216); ">1</span>,n); bisamp.fill_randbi(p,rng);</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><span style="white-space: pre-wrap; "></span>freq[i]=sum(bisamp);</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><span style="white-space: pre-wrap; "></span><span style="color: rgb(187, 44, 162); ">if</span>(sum(freq)==size) <span style="color: rgb(187, 44, 162); ">break</span>;</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><span style="white-space: pre-wrap; "></span>}</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><span style="white-space: pre-wrap; "></span><span style="color: rgb(187, 44, 162); ">return</span> (freq);</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; ">}</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><br>
</div>
<div style="margin: 0px; font-size: 11px; font-family: Menlo; "><br>
</div>
<div>
<div>On 2013-08-01, at 5:35 PM, Peter Stephenson <<a href="mailto:Peter.Stephenson@fish.wa.gov.au">Peter.Stephenson@fish.wa.gov.au</a>></div>
<div> wrote:</div>
<br class="Apple-interchange-newline">
<blockquote type="cite">
<div lang="EN-AU" link="blue" vlink="purple" style="font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; ">
<div class="WordSection1" style="page: WordSection1; ">
<div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; ">
I am wanting to draw a multinomial sample of size say 300 with a vector of known proportions at age.<o:p></o:p></div>
<div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; ">
randbi does not seem to work and the fill_multinomial does not give me what I am looking for.<o:p></o:p></div>
<div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; ">
The .tpl and .dat files are attached. Any suggestions. Cheers Peter<o:p></o:p></div>
</div>
<span><MultinomialSIM.tpl></span><span><MultinomialSIM.dat></span>_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@admb-project.org" style="color: purple; text-decoration: underline; ">Users@admb-project.org</a><br>
<a href="http://lists.admb-project.org/mailman/listinfo/users" style="color: purple; text-decoration: underline; ">http://lists.admb-project.org/mailman/listinfo/users</a><br>
</div>
</blockquote>
</div>
<br>
</div>
</body>
</html>