Displaying 3 results from an estimated 3 matches for "l_converged".
Did you mean:
converged
2003 Sep 22
2
PR#2894
...rrors, for example
>
> r = 1 - pp;
> t = 1 - qq;
These can also be removed by changing qbeta.c:156-157 (R-1.7.1)
to
y = (y-a) * xinbta * (1.0-xinbta) *
exp (logbeta - pp * log(xinbta) - qq * log1p(-xinbta));
I don't understand why you have qbeta.c:167
if (fabs(y)<=acu) goto L_converged
which will fail for certain values of p & q, when x is close to zero as
the beta density tends to infinity. Why not use an exit condition based on
how far away from 'a' you are?
qbeta.c:174 (prevent excess precision)
Might this not just get optimized out? I doubt modern compilers respe...
2003 Nov 24
0
PR#2894
...rrors, for example
>
> r = 1 - pp;
> t = 1 - qq;
These can also be removed by changing qbeta.c:156-157 (R-1.7.1)
to
y = (y-a) * xinbta * (1.0-xinbta) *
exp (logbeta - pp * log(xinbta) - qq * log1p(-xinbta));
I don't understand why you have qbeta.c:167
if (fabs(y)<=acu) goto L_converged
which will fail for certain values of p & q, when x is close to zero as
the beta density tends to infinity. Why not use an exit condition based on
how far away from 'a' you are?
qbeta.c:174 (prevent excess precision)
Might this not just get optimized out? I doubt modern compilers r...
2001 Dec 09
1
error in qbeta (PR#1201)
Full_Name: Ziheng Yang
Version: 1.3.1
OS: Windows 98
Submission from: (NULL) (172.136.54.89)
I noticed that qbeta is sometimes wrong and the error is not even due to the
beta parameters being too extreme. I am calculating the quantiles corresponding
to cdf = 0.05, 0.15, ..., 0.95. The value corresponding to cdf=0.25 is wrong
while all other values are correct.
qbeta(0.05, 0.143891, 0.05) =