pan><span class="op">(</span><span class="va">P</span><span class="op">)</span></span>
<span><span class="op">}</span></span>
<span></span>
<span><span class="co">#*** IRF 1PL</span></span>
<span><span class="va">P_1PL</span> <span class="op">&lt;-</span> <span class="kw">function</span><span class="op">(</span> <span class="va">par</span>, <span class="va">Theta</span>, <span class="va">ncat</span><span class="op">)</span><span class="op">{</span></span>
<span>    <span class="va">b</span> <span class="op"