210 |
{ |
{ |
211 |
q = M_PI_2 * (atan(y) + ck * y / b); |
q = M_PI_2 * (atan(y) + ck * y / b); |
212 |
} |
} |
213 |
if ( tail == GSL_CDF_UPPER ) |
if ( ((tail == GSL_CDF_UPPER ) && (t > 0.0))||((tail==GSL_CDF_LOWER)&&(t<=0.0))) |
214 |
{ |
{ |
215 |
result->val = (1.0 - q)/2; |
result->val = (1.0 - q)/2; |
216 |
} |
} |
228 |
*/ |
*/ |
229 |
u = cornish_fisher ( t, nu ); |
u = cornish_fisher ( t, nu ); |
230 |
q = gsl_cdf_gauss ( u, GSL_CDF_UPPER ); |
q = gsl_cdf_gauss ( u, GSL_CDF_UPPER ); |
231 |
if ( tail == GSL_CDF_UPPER ) |
if ( ((tail == GSL_CDF_UPPER )&&(t>0.0))||((tail==GSL_CDF_LOWER)&&(t<=0.0))) |
232 |
{ |
{ |
233 |
result->val = q; |
result->val = q; |
234 |
} |
} |
282 |
diff = (lg2.val) - (lg1.val); |
diff = (lg2.val) - (lg1.val); |
283 |
p *= pow(y,nu) * exp (diff) / sqrt (M_PI); |
p *= pow(y,nu) * exp (diff) / sqrt (M_PI); |
284 |
|
|
285 |
if ( tail == GSL_CDF_UPPER ) |
if ( ((tail == GSL_CDF_UPPER )&&(t>0.0))||((tail==GSL_CDF_LOWER)&&(t<=0.0))) |
286 |
{ |
{ |
287 |
result->val = p; |
result->val = p; |
288 |
} |
} |