@stdlib/stats-base-dists-f-pdf
Version:
F distribution probability density function (PDF).
13 lines (12 loc) • 2.94 kB
JavaScript
;var I=function(e,r){return function(){return r||e((r={exports:{}}).exports,r),r.exports}};var V=I(function(tr,R){
var E=require('@stdlib/math-base-special-gamma-lanczos-sum-expg-scaled/dist'),T=require('@stdlib/math-base-special-expm1/dist'),y=require('@stdlib/math-base-special-log1p/dist'),_=require('@stdlib/math-base-special-sqrt/dist'),z=require('@stdlib/math-base-special-abs/dist'),p=require('@stdlib/math-base-special-exp/dist'),m=require('@stdlib/math-base-special-pow/dist'),W=require('@stdlib/math-base-special-max/dist'),O=require('@stdlib/math-base-special-min/dist'),N=require('@stdlib/math-base-special-ln/dist'),g=require('@stdlib/constants-float64-max-ln/dist'),D=require('@stdlib/constants-float64-min-ln/dist'),G=require('@stdlib/constants-float64-gamma-lanczos-g/dist'),Y=require('@stdlib/constants-float64-e/dist');function Z(e,r,a,t,l){var i,F,q,f,c,o,n,u,v,P,h,w,A,s;if(!l)return m(a,e)*m(t,r);if(A=e+r,f=e+G-.5,c=r+G-.5,o=A+G-.5,i=E(A),i/=E(e)*E(r),i*=_(c/Y),i*=_(f/o),n=(a*r-t*f)/f,u=(t*e-a*c)/c,O(z(n),z(u))<.2)if(n*u>0||O(e,r)<1)z(n)<.1?i*=p(e*y(n)):i*=m(a*o/f,e),z(u)<.1?i*=p(r*y(u)):i*=m(t*o/c,r);else if(W(z(n),z(u))<.5)F=e<r,q=r/e,F&&q*u<.1||!F&&n/q>.1?(v=T(q*y(u)),v=n+v+v*n,v=e*y(v),i*=p(v)):(v=T(y(n)/q),v=u+v+v*u,v=r*y(v),i*=p(v));else if(z(n)<z(u))if(s=e*y(n)+r*N(t*o/c),s<=D||s>=g){if(s+=N(i),s>=g)return NaN;i=p(s)}else i*=p(s);else if(s=r*y(u)+e*N(a*o/f),s<=D||s>=g){if(s+=N(i),s>=g)return NaN;i=p(s)}else i*=p(s);else if(h=a*o/f,w=t*o/c,n=e*N(h),u=r*N(w),n>=g||n<=D||u>=g||u<=D)if(e<r)if(P=m(w,r/e),v=e*(N(h)+N(P)),v<g&&v>D)i*=m(P*h,e);else{if(u+=n+N(i),u>=g)return NaN;i=p(u)}else if(P=m(h,e/r),v=(N(P)+N(w))*r,v<g&&v>D)i*=m(P*w,r);else{if(u+=n+N(i),u>=g)return NaN;i=p(u)}else i*=m(h,e)*m(w,r);return i}R.exports=Z
});var L=I(function(ur,X){
var $=V();function b(e,r,a){var t,l;return t=$(r,a,e,1-e,!0),l=(1-e)*e,t/=l,t}X.exports=b
});var C=I(function(vr,B){
var M=require('@stdlib/math-base-assert-is-nan/dist'),j=require('@stdlib/constants-float64-pinf/dist'),k=L();function x(e,r,a){var t,l,i;return M(e)||M(r)||M(a)||r<=0||a<=0?NaN:e<0||e===j?0:e===0?r<2?j:r===2?1:0:(t=r*e,t>a?(l=a*r/((a+t)*(a+t)),l*k(a/(a+t),a/2,r/2)):(i=a+t,l=(i*r-e*r*r)/(i*i),l*k(t/(a+t),r/2,a/2)))}B.exports=x
});var Q=I(function(nr,K){
var d=require('@stdlib/utils-constant-function/dist'),S=require('@stdlib/math-base-assert-is-nan/dist'),H=require('@stdlib/constants-float64-pinf/dist'),J=L();function rr(e,r){var a,t,l,i;if(S(e)||S(r)||e<=0||r<=0)return d(NaN);return i=e*r,t=e/2,l=r/2,a=0,e<2?a=H:e===2&&(a=1),F;function F(q){var f,c,o;return S(q)?NaN:q<0||q===H?0:q===0?a:(f=e*q,f>r?(c=i/((r+f)*(r+f)),c*J(r/(r+f),l,t)):(o=r+f,c=(o*e-q*e*e)/(o*o),c*J(e*q/(r+f),t,l)))}}K.exports=rr
});var er=require('@stdlib/utils-define-nonenumerable-read-only-property/dist'),U=C(),ir=Q();er(U,"factory",ir);module.exports=U;
/** @license Apache-2.0 */
/** @license Apache-2.0 */
//# sourceMappingURL=index.js.map