Displaying 2 results from an estimated 2 matches for "c_qnbeta".
Did you mean:
c_qbeta
2020 Mar 26
2
unstable corner of parameter space for qbeta?
Despite the need to focus on pbeta, I'm still willing to put in some effort.
But I find it really helps to have 2-3 others involved, since the questions back
and forth keep matters moving forward. Volunteers?
Thanks to Martin for detailed comments.
JN
On 2020-03-26 10:34 a.m., Martin Maechler wrote:
>>>>>> J C Nash
>>>>>> on Thu, 26 Mar 2020
2020 Mar 26
0
unstable corner of parameter space for qbeta?
This is also strange:
qbeta <- function (p, shape1, shape2, ncp = 0, lower.tail = TRUE, log.p = FALSE)
{
if (missing(ncp))
.Call(C_qbeta, p, shape1, shape2, lower.tail, log.p)
else .Call(C_qnbeta, p, shape1, shape2, ncp, lower.tail,
log.p)
}
Since the default value is 0 for non-centrality, it seems like the logic above is wrong. When ncp=0, C_qnbeta would be called rather than C_qbeta.
Am I missing something?
Ravi
________________________________
From: R-devel <r-devel...