search for: l_converg

Displaying 3 results from an estimated 3 matches for "l_converg".

Did you mean: l_convert
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 res...
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...
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) =