Displaying 1 result from an estimated 1 matches for "ibeta_larg".
Did you mean:
ibeta_large
1998 Mar 26
1
R-beta: problem with locfit
...s seem to go well (= the previous example works and
give sensible results).
Perhaps, a better possibility is to define igamma and ibeta using
the R pgamma and pbeta function. However, I have not tried it.
guido masarotto
--------------------------------------------------------------------
#define IBETA_LARGE 1.0e30
#define IBETA_SMALL 1.0e-30
#define IGAMMA_LARGE 1.0e30
#define DOUBLE_EPS 2.2204460492503131E-16
double ibeta(x, a, b)
double x, a, b;
{ int flipped = 0, i, k, count;
double I = 0, temp, pn[6], ak, bk, next, prev, factor, val;
if (x <= 0) return(0);
if (x >= 1) r...