search for: double_ep

...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)*x > (a+1)) { flipped =...

...integer without a cast corStruct.c:524: error: expected expression before 'double' corStruct.c:524: warning: initialization makes pointer from integer without a cast corStruct.c:526: warning: incompatible implicit declaration of built-in function 'sqrt' corStruct.c:526: error: 'DOUBLE_EPS' undeclared (first use in this function) corStruct.c:534: error: expected expression before 'double' corStruct.c:534: warning: assignment makes pointer from integer without a cast corStruct.c:558: error: 'dqrdc2' undeclared (first use in this function) corStruct.c:558: error: c...