<div dir="ltr"><div><br></div>hi all;<div><br></div><div>I want to write the density function for the doubly non-central F distribution using Rcpp. The algorithm (from Paolella's Intermediate Probability, listing 10.8) computes the weighted sum of the PDF function of the singly non-central F distribution. This is available in R via df; the internal df function is apparently </div><div><div><br></div><div>function (x, df1, df2, ncp, log = FALSE) <br></div><div><div>{</div><div>    if (missing(ncp)) </div><div>        .External(C_df, x, df1, df2, log)</div><div>    else .External(C_dnf, x, df1, df2, ncp, log)</div><div>}</div></div><div><br></div><div>Is it possible to call this C_dnf from Rcpp? One thought I had was to use Rinside, but I want to avoid that overhead if possible.</div><div><br></div><div>cheers,</div><div><br></div><div><br></div>-- <br><div class="gmail_signature"><br><div>--sep</div></div>
</div></div>