Displaying 1 result from an estimated 1 matches for "igamma_larg".
Did you mean:
igamma_large
1998 Mar 26
1
R-beta: problem with locfit
...le 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) return(1);
/* use ibeta(x,a,b) = 1-ibeta(1-x,b,z) */
if ((a+b+1...