this.j$=this.jStat=function(e){function s(n,t){var r=t<n?n:t;return e.pow(10,17-~~(e.log(0<r?r:-r)*e.LOG10E))}function c(n){return"[object Function]"===f.call(n)}function r(n){return"number"==typeof n&&n==n}function i(){return new i._init(arguments)}function o(){return 0}function a(){return 1}function u(n,t){return n===t?1:0}var t=Array.prototype.concat,n=Array.prototype.slice,f=Object.prototype.toString,l=Array.isArray||function(n){return"[object Array]"===f.call(n)};i.fn=i.prototype,(i._init=function(n){var t;if(l(n[0]))if(l(n[0][0])){for(c(n[1])&&(n[0]=i.map(n[0],n[1])),t=0;t<n[0].length;t++)this[t]=n[0][t];this.length=n[0].length}else this[0]=c(n[1])?i.map(n[0],n[1]):n[0],this.length=1;else if(r(n[0]))this[0]=i.seq.apply(null,n),this.length=1;else{if(n[0]instanceof i)return i(n[0].toArray());this[0]=[],this.length=1}return this}).prototype=i.prototype,(i._init.constructor=i).utils={calcRdx:s,isArray:l,isFunction:c,isNumber:r,toVector:function(n){return t.apply([],n)}},i.extend=function(n){var t,r;if(1===arguments.length){for(r in n)i[r]=n[r];return this}for(t=1;t<arguments.length;t++)for(r in arguments[t])n[r]=arguments[t][r];return n},i.rows=function(n){return n.length||1},i.cols=function(n){return n[0].length||1},i.dimensions=function(n){return{rows:i.rows(n),cols:i.cols(n)}},i.row=function(n,t){return n[t]},i.col=function(n,t){for(var r=new Array(n.length),e=0;e<n.length;e++)r[e]=[n[e][t]];return r},i.diag=function(n){for(var t=i.rows(n),r=new Array(t),e=0;e<t;e++)r[e]=[n[e][e]];return r},i.antidiag=function(n){for(var t=i.rows(n)-1,r=new Array(t),e=0;0<=t;t--,e++)r[e]=[n[e][t]];return r},i.transpose=function(n){var t,r,e,i,o,a=[];for(l(n[0])||(n=[n]),r=n.length,e=n[0].length,o=0;o<e;o++){for(t=new Array(r),i=0;i<r;i++)t[i]=n[i][o];a.push(t)}return 1===a.length?a[0]:a},i.map=function(n,t,r){var e,i,o,a,u;for(l(n[0])||(n=[n]),i=n.length,o=n[0].length,a=r?n:new Array(i),e=0;e<i;e++)for(a[e]||(a[e]=new Array(o)),u=0;u<o;u++)a[e][u]=t(n[e][u],e,u);return 1===a.length?a[0]:a},i.cumreduce=function(n,t,r){var e,i,o,a,u;for(l(n[0])||(n=[n]),i=n.length,o=n[0].length,a=r?n:new Array(i),e=0;e<i;e++)for(a[e]||(a[e]=new Array(o)),0<o&&(a[e][0]=n[e][0]),u=1;u<o;u++)a[e][u]=t(a[e][u-1],n[e][u]);return 1===a.length?a[0]:a},i.alter=function(n,t){return i.map(n,t,!0)},i.create=function(n,t,r){var e,i,o=new Array(n);for(c(t)&&(r=t,t=n),e=0;e<n;e++)for(o[e]=new Array(t),i=0;i<t;i++)o[e][i]=r(e,i);return o},i.zeros=function(n,t){return r(t)||(t=n),i.create(n,t,o)},i.ones=function(n,t){return r(t)||(t=n),i.create(n,t,a)},i.rand=function(n,t){return r(t)||(t=n),i.create(n,t,e.random)},i.identity=function(n,t){return r(t)||(t=n),i.create(n,t,u)},i.symmetric=function(n){var t,r,e=n.length;if(n.length!==n[0].length)return!1;for(t=0;t<e;t++)for(r=0;r<e;r++)if(n[r][t]!==n[t][r])return!1;return!0},i.clear=function(n){return i.alter(n,o)},i.seq=function(n,t,r,e){c(e)||(e=!1);var i,o=[],a=s(n,t),u=(t*a-n*a)/((r-1)*a),f=n;for(i=0;f<=t;f=(n*a+u*a*++i)/a)o.push(e?e(f,i):f);return o};var m=i.prototype;return m.length=0,m.push=Array.prototype.push,m.sort=Array.prototype.sort,m.splice=Array.prototype.splice,m.slice=Array.prototype.slice,m.toArray=function(){return 1<this.length?n.call(this):n.call(this)[0]},m.map=function(n,t){return i(i.map(this,n,t))},m.cumreduce=function(n,t){return i(i.cumreduce(this,n,t))},m.alter=function(n){return i.alter(this,n),this},function(n){for(var t=0;t<n.length;t++)!function(e){m[e]=function(n){var t,r=this;return n?(setTimeout(function(){n.call(r,m[e].call(r))}),this):(t=i[e](this),l(t)?i(t):t)}}(n[t])}("transpose clear symmetric rows cols dimensions diag antidiag".split(" ")),function(n){for(var t=0;t<n.length;t++)!function(e){m[e]=function(n,t){var r=this;return t?(setTimeout(function(){t.call(r,m[e].call(r,n))}),this):i(i[e](this,n))}}(n[t])}("row col".split(" ")),function(n){for(var t=0;t<n.length;t++)r=n[t],m[r]=new Function("return jStat(jStat."+r+".apply(null, arguments));");var r}("create zeros ones rand identity".split(" ")),i}(Math),function(f,m){function d(n,t){return n-t}function p(n,t,r){return m.max(t,m.min(n,r))}var s=f.utils.isFunction;f.sum=function(n){for(var t=0,r=n.length;0<=--r;)t+=n[r];return t},f.sumsqrd=function(n){for(var t=0,r=n.length;0<=--r;)t+=n[r]*n[r];return t},f.sumsqerr=function(n){for(var t,r=f.mean(n),e=0,i=n.length;0<=--i;)e+=(t=n[i]-r)*t;return e},f.sumrow=function(n){for(var t=0,r=n.length;0<=--r;)t+=n[r];return t},f.product=function(n){for(var t=1,r=n.length;0<=--r;)t*=n[r];return t},f.min=function(n){for(var t=n[0],r=0;++r<n.length;)n[r]<t&&(t=n[r]);return t},f.max=function(n){for(var t=n[0],r=0;++r<n.length;)n[r]>t&&(t=n[r]);return t},f.mean=function(n){return f.sum(n)/n.length},f.meansqerr=function(n){return f.sumsqerr(n)/n.length},f.geomean=function(n){return m.pow(f.product(n),1/n.length)},f.median=function(n){var t=n.length,r=n.slice().sort(d);return 1&t?r[t/2|0]:(r[t/2-1]+r[t/2])/2},f.cumsum=function(n){return f.cumreduce(n,function(n,t){return n+t})},f.cumprod=function(n){return f.cumreduce(n,function(n,t){return n*t})},f.diff=function(n){var t,r=[],e=n.length;for(t=1;t<e;t++)r.push(n[t]-n[t-1]);return r},f.mode=function(n){var t,r=n.length,e=n.slice().sort(d),i=1,o=0,a=0,u=[];for(t=0;t<r;t++)e[t]===e[t+1]?i++:(o<i?(u=[e[t]],o=i,a=0):i===o&&(u.push(e[t]),a++),i=1);return 0===a?u[0]:u},f.range=function(n){return f.max(n)-f.min(n)},f.variance=function(n,t){return f.sumsqerr(n)/(n.length-(t?1:0))},f.stdev=function(n,t){return m.sqrt(f.variance(n,t))},f.meandev=function(n){var t,r=0,e=f.mean(n);for(t=n.length-1;0<=t;t--)r+=m.abs(n[t]-e);return r/n.length},f.meddev=function(n){var t,r=0,e=f.median(n);for(t=n.length-1;0<=t;t--)r+=m.abs(n[t]-e);return r/n.length},f.coeffvar=function(n){return f.stdev(n)/f.mean(n)},f.quartiles=function(n){var t=n.length,r=n.slice().sort(d);return[r[m.round(t/4)-1],r[m.round(t/2)-1],r[m.round(3*t/4)-1]]},f.quantiles=function(n,t,r,e){var i,o,a,u,f,s=n.slice().sort(d),c=[t.length],l=n.length;for(void 0===r&&(r=3/8),void 0===e&&(e=3/8),i=0;i<t.length;i++)f=p((a=l*(o=t[i])+(r+o*(1-r-e)))-(u=m.floor(p(a,1,l-1))),0,1),c[i]=(1-f)*s[u-1]+f*s[u];return c},f.percentile=function(n,t){var r=n.slice().sort(d),e=t*(r.length-1),i=parseInt(e),o=e-i;return i+1<r.length?r[i]*(1-o)+r[i+1]*o:r[i]},f.percentileOfScore=function(n,t,r){var e,i,o=0,a=n.length,u=!1;for("strict"===r&&(u=!0),i=0;i<a;i++)e=n[i],(u&&e<t||!u&&e<=t)&&o++;return o/a},f.histogram=function(n,t){var r,e=f.min(n),i=t||4,o=(f.max(n)-e)/i,a=n.length;t=[];for(r=0;r<i;r++)t[r]=0;for(r=0;r<a;r++)t[m.min(m.floor((n[r]-e)/o),i-1)]+=1;return t},f.covariance=function(n,t){var r,e=f.mean(n),i=f.mean(t),o=n.length,a=new Array(o);for(r=0;r<o;r++)a[r]=(n[r]-e)*(t[r]-i);return f.sum(a)/(o-1)},f.corrcoeff=function(n,t){return f.covariance(n,t)/f.stdev(n,1)/f.stdev(t,1)},f.spearmancoeff=function(n,t){sortfunction=function(n,t){return n[1]-t[1]};var r=n.length;if(n.length!=t.length)return 0;if(0==r)return 0;var e,i=new Array(r),o=new Array(r),a=new Array(r),u=new Array(r),f=0;for(e=0;e<r;e++)i[e]=[e,n[e]],o[e]=[e,t[e]];var s=i.sort(sortfunction),c=o.sort(sortfunction);for(e=0;e<r;e++)u[c[a[s[e][0]]=e][0]]=e;for(e=f=0;e<r;e++)f+=(a[e]-u[e])*(a[e]-u[e]);return 1-6*f/(r*(r*r-1))},f.stanMoment=function(n,t){var r=f.mean(n),e=f.stdev(n),o=n.length,a=0;for(i=0;i<o;i++)a+=m.pow((n[i]-r)/e,t);return a/n.length},f.skewness=function(n){return f.stanMoment(n,3)},f.kurtosis=function(n){return f.stanMoment(n,4)-3};var c=f.prototype;(function(n){for(var t=0;t<n.length;t++)!function(o){c[o]=function(n,t){var r=[],e=0,i=this;if(s(n)&&(t=n,n=!1),t)return setTimeout(function(){t.call(i,c[o].call(i,n))}),this;if(1<this.length){for(i=!0===n?this:this.transpose();e<i.length;e++)r[e]=f[o](i[e]);return r}return f[o](this[0],n)}}(n[t])})("cumsum cumprod".split(" ")),function(n){for(var t=0;t<n.length;t++)!function(o){c[o]=function(n,t){var r=[],e=0,i=this;if(s(n)&&(t=n,n=!1),t)return setTimeout(function(){t.call(i,c[o].call(i,n))}),this;if(1<this.length){for("sumrow"!==o&&(i=!0===n?this:this.transpose());e<i.length;e++)r[e]=f[o](i[e]);return!0===n?f[o](f.utils.toVector(r)):r}return f[o](this[0],n)}}(n[t])}("sum sumsqrd sumsqerr sumrow product min max mean meansqerr geomean median diff mode range variance stdev meandev meddev coeffvar quartiles histogram skewness kurtosis".split(" ")),function(n){for(var t=0;t<n.length;t++)!function(u){c[u]=function(){var n=[],t=0,r=this,e=Array.prototype.slice.call(arguments);if(s(e[e.length-1])){var i=e[e.length-1],o=e.slice(0,e.length-1);return setTimeout(function(){i.call(r,c[u].apply(r,o))}),this}function a(n){return f[u].apply(r,[n].concat(e))}i=void 0;if(1<this.length){for(r=r.transpose();t<r.length;t++)n[t]=a(r[t]);return n}return a(this[0])}}(n[t])}("quantiles percentileOfScore".split(" "))}(this.jStat,Math),function(g,v){g.gammaln=function(n){var t,r,e,i=0,o=[76.18009172947146,-86.50532032941678,24.01409824083091,-1.231739572450155,.001208650973866179,-5395239384953e-18],a=1.000000000190015;for(e=(r=t=n)+5.5,e-=(t+.5)*v.log(e);i<6;i++)a+=o[i]/++r;return v.log(2.5066282746310007*a/t)-e},g.gammafn=function(n){var t,r,e,i,o=[-1.716185138865495,24.76565080557592,-379.80425647094563,629.3311553128184,866.9662027904133,-31451.272968848367,-36144.413418691176,66456.14382024054],a=[-30.8402300119739,315.35062697960416,-1015.1563674902192,-3107.771671572311,22538.11842098015,4755.846277527881,-134659.9598649693,-115132.2596755535],u=!1,f=0,s=0,c=0,l=n;if(l<=0){if(!(i=l%1+36e-17))return 1/0;u=(1&l?-1:1)*v.PI/v.sin(v.PI*i),l=1-l}for(r=(e=l)<1?l++:(l-=f=(0|l)-1)-1,t=0;t<8;++t)c=(c+o[t])*r,s=s*r+a[t];if(i=c/s+1,e<l)i/=e;else if(l<e)for(t=0;t<f;++t)i*=l,l++;return u&&(i=u/i),i},g.gammap=function(n,t){return g.lowRegGamma(n,t)*g.gammafn(n)},g.lowRegGamma=function(n,t){var r,e=g.gammaln(n),i=n,o=1/n,a=o,u=t+1-n,f=1/1e-30,s=1/u,c=s,l=1,m=-~(8.5*v.log(1<=n?n:1/n)+.4*n+17);if(t<0||n<=0)return NaN;if(t<n+1){for(;l<=m;l++)o+=a*=t/++i;return o*v.exp(-t+n*v.log(t)-e)}for(;l<=m;l++)c*=(s=1/(s=(r=-l*(l-n))*s+(u+=2)))*(f=u+r/f);return 1-c*v.exp(-t+n*v.log(t)-e)},g.factorialln=function(n){return n<0?NaN:g.gammaln(n+1)},g.factorial=function(n){return n<0?NaN:g.gammafn(n+1)},g.combination=function(n,t){return 170<n||170<t?v.exp(g.combinationln(n,t)):g.factorial(n)/g.factorial(t)/g.factorial(n-t)},g.combinationln=function(n,t){return g.factorialln(n)-g.factorialln(t)-g.factorialln(n-t)},g.permutation=function(n,t){return g.factorial(n)/g.factorial(n-t)},g.betafn=function(n,t){return n<=0||t<=0?void 0:170<n+t?v.exp(g.betaln(n,t)):g.gammafn(n)*g.gammafn(t)/g.gammafn(n+t)},g.betaln=function(n,t){return g.gammaln(n)+g.gammaln(t)-g.gammaln(n+t)},g.betacf=function(n,t,r){var e,i,o,a,u=1e-30,f=1,s=t+r,c=t+1,l=t-1,m=1,d=1-s*n/c;for(v.abs(d)<u&&(d=u),a=d=1/d;f<=100&&(d=1+(i=f*(r-f)*n/((l+(e=2*f))*(t+e)))*d,v.abs(d)<u&&(d=u),m=1+i/m,v.abs(m)<u&&(m=u),a*=(d=1/d)*m,d=1+(i=-(t+f)*(s+f)*n/((t+e)*(c+e)))*d,v.abs(d)<u&&(d=u),m=1+i/m,v.abs(m)<u&&(m=u),a*=o=(d=1/d)*m,!(v.abs(o-1)<3e-7));f++);return a},g.gammapinv=function(n,t){var r,e,i,o,a,u,f=0,s=t-1,c=g.gammaln(t);if(1<=n)return v.max(100,t+100*v.sqrt(t));if(n<=0)return 0;for(r=1<t?(a=v.log(s),u=v.exp(s*(a-1)-c),o=n<.5?n:1-n,r=(2.30753+.27061*(e=v.sqrt(-2*v.log(o))))/(1+e*(.99229+.04481*e))-e,n<.5&&(r=-r),v.max(.001,t*v.pow(1-1/(9*t)-r/(3*v.sqrt(t)),3))):n<(e=1-t*(.253+.12*t))?v.pow(n/e,1/t):1-v.log(1-(n-e)/(1-e));f<12;f++){if(r<=0)return 0;if((r-=e=(i=(g.lowRegGamma(t,r)-n)/(e=1<t?u*v.exp(-(r-s)+s*(v.log(r)-a)):v.exp(-r+s*v.log(r)-c)))/(1-.5*v.min(1,i*((t-1)/r-1))))<=0&&(r=.5*(r+e)),v.abs(e)<1e-8*r)break}return r},g.erf=function(n){var t,r,e,i,o=[-1.3026537197817094,.6419697923564902,.019476473204185836,-.00956151478680863,-.000946595344482036,.000366839497852761,42523324806907e-18,-20278578112534e-18,-1624290004647e-18,130365583558e-17,1.5626441722e-8,-8.5238095915e-8,6.529054439e-9,5.059343495e-9,-9.91364156e-10,-2.27365122e-10,96467911e-18,2394038e-18,-6886027e-18,894487e-18,313092e-18,-112708e-18,381e-18,7106e-18,-1523e-18,-94e-18,121e-18,-28e-18],a=o.length-1,u=!1,f=0,s=0;for(n<0&&(n=-n,u=!0),r=4*(t=2/(2+n))-2;0<a;a--)f=r*(e=f)-s+o[a],s=e;return i=t*v.exp(-n*n+.5*(o[0]+r*f)-s),u?i-1:1-i},g.erfc=function(n){return 1-g.erf(n)},g.erfcinv=function(n){var t,r,e,i,o=0;if(2<=n)return-100;if(n<=0)return 100;for(i=n<1?n:2-n,t=-.70711*((2.30753+.27061*(e=v.sqrt(-2*v.log(i/2))))/(1+e*(.99229+.04481*e))-e);o<2;o++)t+=(r=g.erfc(t)-i)/(1.1283791670955126*v.exp(-t*t)-t*r);return n<1?t:-t},g.ibetainv=function(n,t,r){var e,i,o,a,u,f,s,c,l,m,d=t-1,p=r-1,h=0;if(n<=0)return 0;if(1<=n)return 1;for(f=1<=t&&1<=r?(o=n<.5?n:1-n,f=(2.30753+.27061*(a=v.sqrt(-2*v.log(o))))/(1+a*(.99229+.04481*a))-a,n<.5&&(f=-f),s=(f*f-3)/6,c=2/(1/(2*t-1)+1/(2*r-1)),l=f*v.sqrt(s+c)/c-(1/(2*r-1)-1/(2*t-1))*(s+5/6-2/(3*c)),t/(t+r*v.exp(2*l))):(e=v.log(t/(t+r)),i=v.log(r/(t+r)),n<(a=v.exp(t*e)/t)/(l=a+(u=v.exp(r*i)/r))?v.pow(t*l*n,1/t):1-v.pow(r*l*(1-n),1/r)),m=-g.gammaln(t)-g.gammaln(r)+g.gammaln(t+r);h<10;h++){if(0===f||1===f)return f;if((f-=a=(u=(g.ibeta(f,t,r)-n)/(a=v.exp(d*v.log(f)+p*v.log(1-f)+m)))/(1-.5*v.min(1,u*(d/f-p/(1-f)))))<=0&&(f=.5*(f+a)),1<=f&&(f=.5*(f+a+1)),v.abs(a)<1e-8*f&&0<h)break}return f},g.ibeta=function(n,t,r){var e=0===n||1===n?0:v.exp(g.gammaln(t+r)-g.gammaln(t)-g.gammaln(r)+t*v.log(n)+r*v.log(1-n));return!(n<0||1<n)&&(n<(t+1)/(t+r+2)?e*g.betacf(n,t,r)/t:1-e*g.betacf(1-n,r,t)/r)},g.randn=function(n,t){var r,e,i,o,a;if(t=t||n,n)return g.create(n,t,function(){return g.randn()});for(;r=v.random(),e=1.7156*(v.random()-.5),.27597<(a=(i=r-.449871)*i+(o=v.abs(e)+.386595)*(.196*o-.25472*i))&&(.27846<a||e*e>-4*v.log(r)*r*r););return e/r},g.randg=function(n,t,r){var e,i,o,a,u,f,s=n;if(r=r||t,n=n||1,t)return(f=g.zeros(t,r)).alter(function(){return g.randg(n)}),f;n<1&&(n+=1),e=n-1/3,i=1/v.sqrt(9*e);do{for(;(a=1+i*(u=g.randn()))<=0;);a*=a*a,o=v.random()}while(o>1-.331*v.pow(u,4)&&v.log(o)>.5*u*u+e*(1-a+v.log(a)));if(n==s)return e*a;for(;0===(o=v.random()););return v.pow(o,1/s)*e*a},function(n){for(var t=0;t<n.length;t++)!function(t){g.fn[t]=function(){return g(g.map(this,function(n){return g[t](n)}))}}(n[t])}("gammaln gammafn factorial factorialln".split(" ")),function(n){for(var t=0;t<n.length;t++)!function(n){g.fn[n]=function(){return g(g[n].apply(null,arguments))}}(n[t])}("randn".split(" "))}(this.jStat,Math),function(m,d){(function(n){for(var t=0;t<n.length;t++)!function(o){m[o]=function(n,t,r){return this instanceof arguments.callee?(this._a=n,this._b=t,this._c=r,this):new arguments.callee(n,t,r)},m.fn[o]=function(n,t,r){var e=m[o](n,t,r);return e.data=this,e},m[o].prototype.sample=function(n){var t=this._a,r=this._b,e=this._c;return n?m.alter(n,function(){return m[o].sample(t,r,e)}):m[o].sample(t,r,e)},function(n){for(var t=0;t<n.length;t++)!function(i){m[o].prototype[i]=function(n){var t=this._a,r=this._b,e=this._c;return n||0===n||(n=this.data),"number"!=typeof n?m.fn.map.call(n,function(n){return m[o][i](n,t,r,e)}):m[o][i](n,t,r,e)}}(n[t])}("pdf cdf inv".split(" ")),function(n){for(var t=0;t<n.length;t++)!function(n){m[o].prototype[n]=function(){return m[o][n](this._a,this._b,this._c)}}(n[t])}("mean median mode variance".split(" "))}(n[t])})("beta centralF cauchy chisquare exponential gamma invgamma kumaraswamy lognormal noncentralt normal pareto studentt weibull uniform  binomial negbin hypgeom poisson triangular".split(" ")),m.extend(m.beta,{pdf:function(n,t,r){return 1<n||n<0?0:1==t&&1==r?1:t<512||r<512?d.pow(n,t-1)*d.pow(1-n,r-1)/m.betafn(t,r):d.exp((t-1)*d.log(n)+(r-1)*d.log(1-n)-m.betaln(t,r))},cdf:function(n,t,r){return 1<n||n<0?1*(1<n):m.ibeta(n,t,r)},inv:function(n,t,r){return m.ibetainv(n,t,r)},mean:function(n,t){return n/(n+t)},median:function(n,t){return m.ibetainv(.5,n,t)},mode:function(n,t){return(n-1)/(n+t-2)},sample:function(n,t){var r=m.randg(n);return r/(r+m.randg(t))},variance:function(n,t){return n*t/(d.pow(n+t,2)*(n+t+1))}}),m.extend(m.centralF,{pdf:function(n,t,r){var e;return n<0?0:t<=2?0===n&&t<2?1/0:0===n&&2===t?1:d.sqrt(d.pow(t*n,t)*d.pow(r,r)/d.pow(t*n+r,t+r))/(n*m.betafn(t/2,r/2)):(e=t*n/(r+n*t),t*(r/(r+n*t))/2*m.binomial.pdf((t-2)/2,(t+r-2)/2,e))},cdf:function(n,t,r){return n<0?0:m.ibeta(t*n/(t*n+r),t/2,r/2)},inv:function(n,t,r){return r/(t*(1/m.ibetainv(n,t/2,r/2)-1))},mean:function(n,t){return 2<t?t/(t-2):void 0},mode:function(n,t){return 2<n?t*(n-2)/(n*(t+2)):void 0},sample:function(n,t){return 2*m.randg(n/2)/n/(2*m.randg(t/2)/t)},variance:function(n,t){return t<=4?void 0:2*t*t*(n+t-2)/(n*(t-2)*(t-2)*(t-4))}}),m.extend(m.cauchy,{pdf:function(n,t,r){return r/(d.pow(n-t,2)+d.pow(r,2))/d.PI},cdf:function(n,t,r){return d.atan((n-t)/r)/d.PI+.5},inv:function(n,t,r){return t+r*d.tan(d.PI*(n-.5))},median:function(n,t){return n},mode:function(n,t){return n},sample:function(n,t){return m.randn()*d.sqrt(1/(2*m.randg(.5)))*t+n}}),m.extend(m.chisquare,{pdf:function(n,t){return n<0?0:0===n&&2===t?.5:d.exp((t/2-1)*d.log(n)-n/2-t/2*d.log(2)-m.gammaln(t/2))},cdf:function(n,t){return n<0?0:m.lowRegGamma(t/2,n/2)},inv:function(n,t){return 2*m.gammapinv(n,.5*t)},mean:function(n){return n},median:function(n){return n*d.pow(1-2/(9*n),3)},mode:function(n){return 0<n-2?n-2:0},sample:function(n){return 2*m.randg(n/2)},variance:function(n){return 2*n}}),m.extend(m.exponential,{pdf:function(n,t){return n<0?0:t*d.exp(-t*n)},cdf:function(n,t){return n<0?0:1-d.exp(-t*n)},inv:function(n,t){return-d.log(1-n)/t},mean:function(n){return 1/n},median:function(n){return 1/n*d.log(2)},mode:function(n){return 0},sample:function(n){return-1/n*d.log(d.random())},variance:function(n){return d.pow(n,-2)}}),m.extend(m.gamma,{pdf:function(n,t,r){return n<0?0:0===n&&1===t?1/r:d.exp((t-1)*d.log(n)-n/r-m.gammaln(t)-t*d.log(r))},cdf:function(n,t,r){return n<0?0:m.lowRegGamma(t,n/r)},inv:function(n,t,r){return m.gammapinv(n,t)*r},mean:function(n,t){return n*t},mode:function(n,t){return 1<n?(n-1)*t:void 0},sample:function(n,t){return m.randg(n)*t},variance:function(n,t){return n*t*t}}),m.extend(m.invgamma,{pdf:function(n,t,r){return n<=0?0:d.exp(-(t+1)*d.log(n)-r/n-m.gammaln(t)+t*d.log(r))},cdf:function(n,t,r){return n<=0?0:1-m.lowRegGamma(t,r/n)},inv:function(n,t,r){return r/m.gammapinv(1-n,t)},mean:function(n,t){return 1<n?t/(n-1):void 0},mode:function(n,t){return t/(n+1)},sample:function(n,t){return t/m.randg(n)},variance:function(n,t){return n<=2?void 0:t*t/((n-1)*(n-1)*(n-2))}}),m.extend(m.kumaraswamy,{pdf:function(n,t,r){return 0===n&&1===t?r:1===n&&1===r?t:d.exp(d.log(t)+d.log(r)+(t-1)*d.log(n)+(r-1)*d.log(1-d.pow(n,t)))},cdf:function(n,t,r){return n<0?0:1<n?1:1-d.pow(1-d.pow(n,t),r)},inv:function(n,t,r){return d.pow(1-d.pow(1-n,1/r),1/t)},mean:function(n,t){return t*m.gammafn(1+1/n)*m.gammafn(t)/m.gammafn(1+1/n+t)},median:function(n,t){return d.pow(1-d.pow(2,-1/t),1/n)},mode:function(n,t){return 1<=n&&1<=t&&1!==n&&1!==t?d.pow((n-1)/(n*t-1),1/n):void 0},variance:function(n,t){throw new Error("variance not yet implemented")}}),m.extend(m.lognormal,{pdf:function(n,t,r){return n<=0?0:d.exp(-d.log(n)-.5*d.log(2*d.PI)-d.log(r)-d.pow(d.log(n)-t,2)/(2*r*r))},cdf:function(n,t,r){return n<0?0:.5+.5*m.erf((d.log(n)-t)/d.sqrt(2*r*r))},inv:function(n,t,r){return d.exp(-1.4142135623730951*r*m.erfcinv(2*n)+t)},mean:function(n,t){return d.exp(n+t*t/2)},median:function(n,t){return d.exp(n)},mode:function(n,t){return d.exp(n-t*t)},sample:function(n,t){return d.exp(m.randn()*t+n)},variance:function(n,t){return(d.exp(t*t)-1)*d.exp(2*n+t*t)}}),m.extend(m.noncentralt,{pdf:function(n,t,r){return d.abs(r)<1e-14?m.studentt.pdf(n,t):d.abs(n)<1e-14?d.exp(m.gammaln((t+1)/2)-r*r/2-.5*d.log(d.PI*t)-m.gammaln(t/2)):t/n*(m.noncentralt.cdf(n*d.sqrt(1+2/t),t+2,r)-m.noncentralt.cdf(n,t,r))},cdf:function(n,t,r){var e=1e-14;if(d.abs(r)<e)return m.studentt.cdf(n,t);var i=!1;n<0&&(i=!0,r=-r);for(var o=m.normal.cdf(-r,0,1),a=1+e,u=a,f=n*n/(n*n+t),s=0,c=d.exp(-r*r/2),l=d.exp(-r*r/2-.5*d.log(2)-m.gammaln(1.5))*r;s<200||e<u||e<a;)u=a,0<s&&(c*=r*r/(2*s),l*=r*r/(2*(s+.5))),o+=.5*(a=c*m.beta.cdf(f,s+.5,t/2)+l*m.beta.cdf(f,s+1,t/2)),s++;return i?1-o:o}}),m.extend(m.normal,{pdf:function(n,t,r){return d.exp(-.5*d.log(2*d.PI)-d.log(r)-d.pow(n-t,2)/(2*r*r))},cdf:function(n,t,r){return.5*(1+m.erf((n-t)/d.sqrt(2*r*r)))},inv:function(n,t,r){return-1.4142135623730951*r*m.erfcinv(2*n)+t},mean:function(n,t){return n},median:function(n,t){return n},mode:function(n,t){return n},sample:function(n,t){return m.randn()*t+n},variance:function(n,t){return t*t}}),m.extend(m.pareto,{pdf:function(n,t,r){return n<t?0:r*d.pow(t,r)/d.pow(n,r+1)},cdf:function(n,t,r){return n<t?0:1-d.pow(t/n,r)},inv:function(n,t,r){return t/d.pow(1-n,1/r)},mean:function(n,t){return t<=1?void 0:t*d.pow(n,t)/(t-1)},median:function(n,t){return n*t*d.SQRT2},mode:function(n,t){return n},variance:function(n,t){return t<=2?void 0:n*n*t/(d.pow(t-1,2)*(t-2))}}),m.extend(m.studentt,{pdf:function(n,t){return t=1e100<t?1e100:t,1/(d.sqrt(t)*m.betafn(.5,t/2))*d.pow(1+n*n/t,-(t+1)/2)},cdf:function(n,t){var r=t/2;return m.ibeta((n+d.sqrt(n*n+t))/(2*d.sqrt(n*n+t)),r,r)},inv:function(n,t){var r=m.ibetainv(2*d.min(n,1-n),.5*t,.5);return r=d.sqrt(t*(1-r)/r),.5<n?r:-r},mean:function(n){return 1<n?0:void 0},median:function(n){return 0},mode:function(n){return 0},sample:function(n){return m.randn()*d.sqrt(n/(2*m.randg(n/2)))},variance:function(n){return 2<n?n/(n-2):1<n?1/0:void 0}}),m.extend(m.weibull,{pdf:function(n,t,r){return n<0?0:r/t*d.pow(n/t,r-1)*d.exp(-d.pow(n/t,r))},cdf:function(n,t,r){return n<0?0:1-d.exp(-d.pow(n/t,r))},inv:function(n,t,r){return t*d.pow(-d.log(1-n),1/r)},mean:function(n,t){return n*m.gammafn(1+1/t)},median:function(n,t){return n*d.pow(d.log(2),1/t)},mode:function(n,t){return t<=1?void 0:n*d.pow((t-1)/t,1/t)},sample:function(n,t){return n*d.pow(-d.log(d.random()),1/t)},variance:function(n,t){return n*n*m.gammafn(1+2/t)-d.pow(this.mean(n,t),2)}}),m.extend(m.uniform,{pdf:function(n,t,r){return n<t||r<n?0:1/(r-t)},cdf:function(n,t,r){return n<t?0:n<r?(n-t)/(r-t):1},inv:function(n,t,r){return t+n*(r-t)},mean:function(n,t){return.5*(n+t)},median:function(n,t){return m.mean(n,t)},mode:function(n,t){throw new Error("mode is not yet implemented")},sample:function(n,t){return n/2+t/2+(t/2-n/2)*(2*d.random()-1)},variance:function(n,t){return d.pow(t-n,2)/12}}),m.extend(m.binomial,{pdf:function(n,t,r){return 0===r||1===r?t*r===n?1:0:m.combination(t,n)*d.pow(r,n)*d.pow(1-r,t-n)},cdf:function(n,t,r){var e=[],i=0;if(n<0)return 0;if(n<t){for(;i<=n;i++)e[i]=m.binomial.pdf(i,t,r);return m.sum(e)}return 1}}),m.extend(m.negbin,{pdf:function(n,t,r){return n===n>>>0&&(n<0?0:m.combination(n+t-1,t-1)*d.pow(1-r,n)*d.pow(r,t))},cdf:function(n,t,r){var e=0,i=0;if(n<0)return 0;for(;i<=n;i++)e+=m.negbin.pdf(i,t,r);return e}}),m.extend(m.hypgeom,{pdf:function(n,t,r,e){if(n!=n|0)return!1;if(n<0||n<r-(t-e))return 0;if(e<n||r<n)return 0;if(t<2*r)return t<2*e?m.hypgeom.pdf(t-r-e+n,t,t-r,t-e):m.hypgeom.pdf(e-n,t,t-r,e);if(t<2*e)return m.hypgeom.pdf(r-n,t,r,t-e);if(r<e)return m.hypgeom.pdf(n,t,e,r);for(var i=1,o=0,a=0;a<n;a++){for(;1<i&&o<e;)i*=1-r/(t-o),o++;i*=(e-a)*(r-a)/((a+1)*(t-r-e+a+1))}for(;o<e;o++)i*=1-r/(t-o);return d.min(1,d.max(0,i))},cdf:function(n,t,r,e){if(n<0||n<r-(t-e))return 0;if(e<=n||r<=n)return 1;if(t<2*r)return t<2*e?m.hypgeom.cdf(t-r-e+n,t,t-r,t-e):1-m.hypgeom.cdf(e-n-1,t,t-r,e);if(t<2*e)return 1-m.hypgeom.cdf(r-n-1,t,r,t-e);if(r<e)return m.hypgeom.cdf(n,t,e,r);for(var i=1,o=1,a=0,u=0;u<n;u++){for(;1<i&&a<e;){var f=1-r/(t-a);o*=f,i*=f,a++}i+=o*=(e-u)*(r-u)/((u+1)*(t-r-e+u+1))}for(;a<e;a++)i*=1-r/(t-a);return d.min(1,d.max(0,i))}}),m.extend(m.poisson,{pdf:function(n,t){return d.pow(t,n)*d.exp(-t)/m.factorial(n)},cdf:function(n,t){var r=[],e=0;if(n<0)return 0;for(;e<=n;e++)r.push(m.poisson.pdf(e,t));return m.sum(r)},mean:function(n){return n},variance:function(n){return n},sample:function(n){for(var t=1,r=0,e=d.exp(-n);r++,e<(t*=d.random()););return r-1}}),m.extend(m.triangular,{pdf:function(n,t,r,e){return r<=t||e<t||r<e?NaN:n<t||r<n?0:n<e?2*(n-t)/((r-t)*(e-t)):n===e?2/(r-t):2*(r-n)/((r-t)*(r-e))},cdf:function(n,t,r,e){return r<=t||e<t||r<e?NaN:n<=t?0:r<=n?1:n<=e?d.pow(n-t,2)/((r-t)*(e-t)):1-d.pow(r-n,2)/((r-t)*(r-e))},inv:function(n,t,r,e){return r<=t||e<t||r<e?NaN:n<=(e-t)/(r-t)?t+(r-t)*d.sqrt(n*((e-t)/(r-t))):t+(r-t)*(1-d.sqrt((1-n)*(1-(e-t)/(r-t))))},mean:function(n,t,r){return(n+t+r)/3},median:function(n,t,r){return r<=(n+t)/2?t-d.sqrt((t-n)*(t-r))/d.sqrt(2):(n+t)/2<r?n+d.sqrt((t-n)*(r-n))/d.sqrt(2):void 0},mode:function(n,t,r){return r},sample:function(n,t,r){var e=d.random();return e<(r-n)/(t-n)?n+d.sqrt(e*(t-n)*(r-n)):t-d.sqrt((1-e)*(t-n)*(t-r))},variance:function(n,t,r){return(n*n+t*t+r*r-n*t-n*r-t*r)/18}})}(this.jStat,Math),function(g,h){var i=Array.prototype.push,l=g.utils.isArray;g.extend({add:function(n,e){return l(e)?(l(e[0])||(e=[e]),g.map(n,function(n,t,r){return n+e[t][r]})):g.map(n,function(n){return n+e})},subtract:function(n,e){return l(e)?(l(e[0])||(e=[e]),g.map(n,function(n,t,r){return n-e[t][r]||0})):g.map(n,function(n){return n-e})},divide:function(n,t){return l(t)?(l(t[0])||(t=[t]),g.multiply(n,g.inv(t))):g.map(n,function(n){return n/t})},multiply:function(n,t){var r,e,i,o,a=n.length,u=n[0].length,f=g.zeros(a,i=l(t)?t[0].length:u),s=0;if(l(t)){for(;s<i;s++)for(r=0;r<a;r++){for(e=o=0;e<u;e++)o+=n[r][e]*t[e][s];f[r][s]=o}return 1===a&&1===s?f[0][0]:f}return g.map(n,function(n){return n*t})},dot:function(n,t){l(n[0])||(n=[n]),l(t[0])||(t=[t]);for(var r,e,i=1===n[0].length&&1!==n.length?g.transpose(n):n,o=1===t[0].length&&1!==t.length?g.transpose(t):t,a=[],u=0,f=i.length,s=i[0].length;u<f;u++){for(a[u]=[],e=r=0;e<s;e++)r+=i[u][e]*o[u][e];a[u]=r}return 1===a.length?a[0]:a},pow:function(n,t){return g.map(n,function(n){return h.pow(n,t)})},exp:function(n){return g.map(n,function(n){return h.exp(n)})},log:function(n){return g.map(n,function(n){return h.log(n)})},abs:function(n){return g.map(n,function(n){return h.abs(n)})},norm:function(n,t){var r=0,e=0;for(isNaN(t)&&(t=2),l(n[0])&&(n=n[0]);e<n.length;e++)r+=h.pow(h.abs(n[e]),t);return h.pow(r,1/t)},angle:function(n,t){return h.acos(g.dot(n,t)/(g.norm(n)*g.norm(t)))},aug:function(n,t){for(var r=[],e=0;e<n.length;e++)r.push(n[e].slice());for(e=0;e<r.length;e++)i.apply(r[e],t[e]);return r},inv:function(n){for(var t,r=n.length,e=n[0].length,i=g.identity(r,e),o=g.gauss_jordan(n,i),a=[],u=0;u<r;u++)for(a[u]=[],t=e;t<o[0].length;t++)a[u][t-e]=o[u][t];return a},det:function(n){var t,r=n.length,e=2*r,i=new Array(e),o=r-1,a=e-1,u=o-r+1,f=a,s=0,c=0;if(2===r)return n[0][0]*n[1][1]-n[0][1]*n[1][0];for(;s<e;s++)i[s]=1;for(s=0;s<r;s++){for(t=0;t<r;t++)i[u<0?u+r:u]*=n[s][t],i[f<r?f+r:f]*=n[s][t],u++,f--;u=--o-r+1,f=--a}for(s=0;s<r;s++)c+=i[s];for(;s<e;s++)c-=i[s];return c},gauss_elimination:function(n,t){var r,e,i,o,a=0,u=0,f=n.length,s=n[0].length,c=1,l=0,m=[];for(r=(n=g.aug(n,t))[0].length,a=0;a<f;a++){for(e=n[a][a],o=(u=a)+1;o<s;o++)e<h.abs(n[o][a])&&(e=n[o][a],u=o);if(u!=a)for(o=0;o<r;o++)i=n[a][o],n[a][o]=n[u][o],n[u][o]=i;for(u=a+1;u<f;u++)for(c=n[u][a]/n[a][a],o=a;o<r;o++)n[u][o]=n[u][o]-c*n[a][o]}for(a=f-1;0<=a;a--){for(l=0,u=a+1;u<=f-1;u++)l+=m[u]*n[a][u];m[a]=(n[a][r-1]-l)/n[a][a]}return m},gauss_jordan:function(n,t){for(var r=g.aug(n,t),e=r.length,i=r[0].length,o=0;o<e;o++){for(var a=o,u=o+1;u<e;u++)h.abs(r[u][o])>h.abs(r[a][o])&&(a=u);var f=r[o];r[o]=r[a],r[a]=f;for(u=o+1;u<e;u++){c=r[u][o]/r[o][o];for(var s=o;s<i;s++)r[u][s]-=r[o][s]*c}}for(o=e-1;0<=o;o--){c=r[o][o];for(u=0;u<o;u++)for(s=i-1;o-1<s;s--)r[u][s]-=r[o][s]*r[u][o]/c;r[o][o]/=c;for(s=e;s<i;s++)r[o][s]/=c}return r},lu:function(n,t){throw new Error("lu not yet implemented")},cholesky:function(n,t){throw new Error("cholesky not yet implemented")},gauss_jacobi:function(n,t,r,e){for(var i,o,a,u,f=0,s=0,c=n.length,l=[],m=[],d=[];f<c;f++)for(l[f]=[],m[f]=[],d[f]=[],s=0;s<c;s++)s<f?(l[f][s]=n[f][s],m[f][s]=d[f][s]=0):f<s?(m[f][s]=n[f][s],l[f][s]=d[f][s]=0):(d[f][s]=n[f][s],l[f][s]=m[f][s]=0);for(a=g.multiply(g.multiply(g.inv(d),g.add(l,m)),-1),o=g.multiply(g.inv(d),t),i=r,u=g.add(g.multiply(a,r),o),f=2;h.abs(g.norm(g.subtract(u,i)))>e;)i=u,u=g.add(g.multiply(a,i),o),f++;return u},gauss_seidel:function(n,t,r,e){for(var i,o,a,u,f,s=0,c=n.length,l=[],m=[],d=[];s<c;s++)for(l[s]=[],m[s]=[],d[s]=[],i=0;i<c;i++)i<s?(l[s][i]=n[s][i],m[s][i]=d[s][i]=0):s<i?(m[s][i]=n[s][i],l[s][i]=d[s][i]=0):(d[s][i]=n[s][i],l[s][i]=m[s][i]=0);for(u=g.multiply(g.multiply(g.inv(g.add(d,l)),m),-1),a=g.multiply(g.inv(g.add(d,l)),t),o=r,f=g.add(g.multiply(u,r),a),s=2;h.abs(g.norm(g.subtract(f,o)))>e;)o=f,f=g.add(g.multiply(u,o),a),s+=1;return f},SOR:function(n,t,r,e,i){for(var o,a,u,f,s,c=0,l=n.length,m=[],d=[],p=[];c<l;c++)for(m[c]=[],d[c]=[],p[c]=[],o=0;o<l;o++)o<c?(m[c][o]=n[c][o],d[c][o]=p[c][o]=0):c<o?(d[c][o]=n[c][o],m[c][o]=p[c][o]=0):(p[c][o]=n[c][o],m[c][o]=d[c][o]=0);for(f=g.multiply(g.inv(g.add(p,g.multiply(m,i))),g.subtract(g.multiply(p,1-i),g.multiply(d,i))),u=g.multiply(g.multiply(g.inv(g.add(p,g.multiply(m,i))),t),i),a=r,s=g.add(g.multiply(f,r),u),c=2;h.abs(g.norm(g.subtract(s,a)))>e;)a=s,s=g.add(g.multiply(f,a),u),c++;return s},householder:function(n){for(var t,r,e,i,o=n.length,a=n[0].length,u=0,f=[],s=[];u<o-1;u++){for(t=0,i=u+1;i<a;i++)t+=n[i][u]*n[i][u];for(t=(0<n[u+1][u]?-1:1)*h.sqrt(t),r=h.sqrt((t*t-n[u+1][u]*t)/2),(f=g.zeros(o,1))[u+1][0]=(n[u+1][u]-t)/(2*r),e=u+2;e<o;e++)f[e][0]=n[e][u]/(2*r);s=g.subtract(g.identity(o,a),g.multiply(g.multiply(f,g.transpose(f)),2)),n=g.multiply(s,g.multiply(n,s))}return n},QR:function(n,t){for(var r,e,i,o,a=n.length,u=n[0].length,f=0,s=[],c=[],l=[];f<a-1;f++){for(e=0,r=f+1;r<u;r++)e+=n[r][f]*n[r][f];for(e=(0<n[f+1][f]?-1:1)*h.sqrt(e),i=h.sqrt((e*e-n[f+1][f]*e)/2),(s=g.zeros(a,1))[f+1][0]=(n[f+1][f]-e)/(2*i),o=f+2;o<a;o++)s[o][0]=n[o][f]/(2*i);c=g.subtract(g.identity(a,u),g.multiply(g.multiply(s,g.transpose(s)),2)),n=g.multiply(c,n),t=g.multiply(c,t)}for(f=a-1;0<=f;f--){for(0,r=f+1;r<=u-1;r++)l[r]*n[f][r];l[f]=t[f][0]/n[f][f]}return l},jacobi:function(n){for(var t,r,e,i,o,a,u,f=1,s=n.length,c=g.identity(s,s),l=[];1===f;){for(0,o=n[0][1],i=1,t=e=0;t<s;t++)for(r=0;r<s;r++)t!=r&&o<h.abs(n[t][r])&&(o=h.abs(n[t][r]),e=t,i=r);for(a=n[e][e]===n[i][i]?0<n[e][i]?h.PI/4:-h.PI/4:h.atan(2*n[e][i]/(n[e][e]-n[i][i]))/2,(u=g.identity(s,s))[e][e]=h.cos(a),u[e][i]=-h.sin(a),u[i][e]=h.sin(a),u[i][i]=h.cos(a),c=g.multiply(c,u),n=g.multiply(g.multiply(g.inv(u),n),u),f=0,t=1;t<s;t++)for(r=1;r<s;r++)t!=r&&.001<h.abs(n[t][r])&&(f=1)}for(t=0;t<s;t++)l.push(n[t][t]);return[c,l]},rungekutta:function(n,t,r,e,i,o){var a,u,f;if(2===o)for(;e<=r;)i=i+((a=t*n(e,i))+(u=t*n(e+t,i+a)))/2,e+=t;if(4===o)for(;e<=r;)i=i+((a=t*n(e,i))+2*(u=t*n(e+t/2,i+a/2))+2*(f=t*n(e+t/2,i+u/2))+t*n(e+t,i+f))/6,e+=t;return i},romberg:function(n,t,r,e){for(var i,o,a,u,f,s=0,c=(r-t)/2,l=[],m=[],d=[];s<e/2;){for(f=n(t),a=t,u=0;a<=r;a+=c,u++)l[u]=a;for(i=l.length,a=1;a<i-1;a++)f+=(a%2!=0?4:2)*n(l[a]);f=c/3*(f+n(r)),d[s]=f,c/=2,s++}for(o=d.length,i=1;1!==o;){for(a=0;a<o-1;a++)m[a]=(h.pow(4,i)*d[a+1]-d[a])/(h.pow(4,i)-1);o=m.length,d=m,m=[],i++}return d},richardson:function(n,t,r,e){function i(n,t){for(var r,e=0,i=n.length;e<i;e++)n[e]===t&&(r=e);return r}n.length;for(var o,a,u,f,s,c=h.abs(r-n[i(n,r)+1]),l=0,m=[],d=[];c<=e;)o=i(n,r+e),a=i(n,r),m[l]=(t[o]-2*t[a]+t[2*a-o])/(e*e),e/=2,l++;for(f=m.length,u=1;1!=f;){for(s=0;s<f-1;s++)d[s]=(h.pow(4,u)*m[s+1]-m[s])/(h.pow(4,u)-1);f=d.length,m=d,d=[],u++}return m},simpson:function(n,t,r,e){for(var i,o=(r-t)/e,a=n(t),u=[],f=t,s=0,c=1;f<=r;f+=o,s++)u[s]=f;for(i=u.length;c<i-1;c++)a+=(c%2!=0?4:2)*n(u[c]);return o/3*(a+n(r))},hermite:function(n,t,r,e){for(var i,o=n.length,a=0,u=0,f=[],s=[],c=[],l=[];u<o;u++){for(f[u]=1,i=0;i<o;i++)u!=i&&(f[u]*=(e-n[i])/(n[u]-n[i]));for(i=s[u]=0;i<o;i++)u!=i&&(s[u]+=1/(n[u]-n[i]));c[u]=(1-2*(e-n[u])*s[u])*f[u]*f[u],l[u]=(e-n[u])*f[u]*f[u],a+=c[u]*t[u]+l[u]*r[u]}return a},lagrange:function(n,t,r){for(var e,i,o=0,a=0,u=n.length;a<u;a++){for(i=t[a],e=0;e<u;e++)a!=e&&(i*=(r-n[e])/(n[a]-n[e]));o+=i}return o},cubic_spline:function(n,t,r){for(var e,i,o=n.length,a=0,u=[],f=[],s=[],c=[],l=[],m=[];a<o-1;a++)c[a]=n[a+1]-n[a];for(s[0]=0,a=1;a<o-1;a++)s[a]=3/c[a]*(t[a+1]-t[a])-3/c[a-1]*(t[a]-t[a-1]);for(a=1;a<o-1;a++)u[a]=[],f[a]=[],u[a][a-1]=c[a-1],u[a][a]=2*(c[a-1]+c[a]),u[a][a+1]=c[a],f[a][0]=s[a];for(i=g.multiply(g.inv(u),f),e=0;e<o-1;e++)l[e]=(t[e+1]-t[e])/c[e]-c[e]*(i[e+1][0]+2*i[e][0])/3,m[e]=(i[e+1][0]-i[e][0])/(3*c[e]);for(e=0;e<o&&!(n[e]>r);e++);return t[e-=1]+(r-n[e])*l[e]+g.sq(r-n[e])*i[e]+(r-n[e])*g.sq(r-n[e])*m[e]},gauss_quadrature:function(){throw new Error("gauss_quadrature not yet implemented")},PCA:function(n){var t,r,e,i,o,a=n.length,u=n[0].length,f=0,s=[],c=[],l=[],m=[],d=[],p=[],h=[];for(f=0;f<a;f++)s[f]=g.sum(n[f])/u;for(f=0;f<u;f++)for(d[f]=[],t=0;t<a;t++)d[f][t]=n[t][f]-s[t];for(d=g.transpose(d),f=0;f<a;f++)for(p[f]=[],t=0;t<a;t++)p[f][t]=g.dot([d[f]],[d[t]])/(u-1);for(o=(e=g.jacobi(p))[0],c=e[1],h=g.transpose(o),f=0;f<c.length;f++)for(t=f;t<c.length;t++)c[f]<c[t]&&(r=c[f],c[f]=c[t],c[t]=r,l=h[f],h[f]=h[t],h[t]=l);for(i=g.transpose(d),f=0;f<a;f++)for(m[f]=[],t=0;t<i.length;t++)m[f][t]=g.dot([h[f]],[i[t]]);return[n,c,h,m]}}),function(n){for(var t=0;t<n.length;t++)!function(e){g.fn[e]=function(n,t){var r=this;return t?(setTimeout(function(){t.call(r,g.fn[e].call(r,n))},15),this):"number"==typeof g[e](this,n)?g[e](this,n):g(g[e](this,n))}}(n[t])}("add divide multiply subtract dot pow exp log abs norm angle".split(" "))}(this.jStat,Math),function(s,c){function o(n,t,r,e){if(1<n||1<r||n<=0||r<=0)throw new Error("Proportions should be greater than 0 and less than 1");var i=(n*t+r*e)/(t+e);return(n-r)/c.sqrt(i*(1-i)*(1/t+1/e))}var l=[].slice,a=s.utils.isNumber,r=s.utils.isArray;s.extend({zscore:function(){var n=l.call(arguments);return a(n[1])?(n[0]-n[1])/n[2]:(n[0]-s.mean(n[1]))/s.stdev(n[1],n[2])},ztest:function(){var n,t=l.call(arguments);return r(t[1])?(n=s.zscore(t[0],t[1],t[3]),1===t[2]?s.normal.cdf(-c.abs(n),0,1):2*s.normal.cdf(-c.abs(n),0,1)):2<t.length?(n=s.zscore(t[0],t[1],t[2]),1===t[3]?s.normal.cdf(-c.abs(n),0,1):2*s.normal.cdf(-c.abs(n),0,1)):(n=t[0],1===t[1]?s.normal.cdf(-c.abs(n),0,1):2*s.normal.cdf(-c.abs(n),0,1))}}),s.extend(s.fn,{zscore:function(n,t){return(n-this.mean())/this.stdev(t)},ztest:function(n,t,r){var e=c.abs(this.zscore(n,r));return 1===t?s.normal.cdf(-e,0,1):2*s.normal.cdf(-e,0,1)}}),s.extend({tscore:function(){var n=l.call(arguments);return 4===n.length?(n[0]-n[1])/(n[2]/c.sqrt(n[3])):(n[0]-s.mean(n[1]))/(s.stdev(n[1],!0)/c.sqrt(n[1].length))},ttest:function(){var n,t=l.call(arguments);return 5===t.length?(n=c.abs(s.tscore(t[0],t[1],t[2],t[3])),1===t[4]?s.studentt.cdf(-n,t[3]-1):2*s.studentt.cdf(-n,t[3]-1)):a(t[1])?(n=c.abs(t[0]),1==t[2]?s.studentt.cdf(-n,t[1]-1):2*s.studentt.cdf(-n,t[1]-1)):(n=c.abs(s.tscore(t[0],t[1])),1==t[2]?s.studentt.cdf(-n,t[1].length-1):2*s.studentt.cdf(-n,t[1].length-1))}}),s.extend(s.fn,{tscore:function(n){return(n-this.mean())/(this.stdev(!0)/c.sqrt(this.cols()))},ttest:function(n,t){return 1===t?1-s.studentt.cdf(c.abs(this.tscore(n)),this.cols()-1):2*s.studentt.cdf(-c.abs(this.tscore(n)),this.cols()-1)}}),s.extend({anovafscore:function(){var n,t,r,e,i,o,a,u,f=l.call(arguments);if(1===f.length){for(i=new Array(f[0].length),a=0;a<f[0].length;a++)i[a]=f[0][a];f=i}if(2===f.length)return s.variance(f[0])/s.variance(f[1]);for(t=new Array,a=0;a<f.length;a++)t=t.concat(f[a]);for(r=s.mean(t),a=n=0;a<f.length;a++)n+=f[a].length*c.pow(s.mean(f[a])-r,2);for(n/=f.length-1,a=o=0;a<f.length;a++)for(e=s.mean(f[a]),u=0;u<f[a].length;u++)o+=c.pow(f[a][u]-e,2);return n/(o/=t.length-f.length)},anovaftest:function(){var n,t,r,e,i=l.call(arguments);if(a(i[0]))return 1-s.centralF.cdf(i[0],i[1],i[2]);for(anovafscore=s.anovafscore(i),n=i.length-1,e=r=0;e<i.length;e++)r+=i[e].length;return t=r-n-1,1-s.centralF.cdf(anovafscore,n,t)},ftest:function(n,t,r){return 1-s.centralF.cdf(n,t,r)}}),s.extend(s.fn,{anovafscore:function(){return s.anovafscore(this.toArray())},anovaftes:function(){var n,t=0;for(n=0;n<this.length;n++)t+=this[n].length;return s.ftest(this.anovafscore(),this.length-1,t-this.length)}}),s.extend({normalci:function(){var n,t=l.call(arguments),r=new Array(2);return n=4===t.length?c.abs(s.normal.inv(t[1]/2,0,1)*t[2]/c.sqrt(t[3])):c.abs(s.normal.inv(t[1]/2,0,1)*s.stdev(t[2])/c.sqrt(t[2].length)),r[0]=t[0]-n,r[1]=t[0]+n,r},tci:function(){var n,t=l.call(arguments),r=new Array(2);return n=4===t.length?c.abs(s.studentt.inv(t[1]/2,t[3]-1)*t[2]/c.sqrt(t[3])):c.abs(s.studentt.inv(t[1]/2,t[2].length-1)*s.stdev(t[2],!0)/c.sqrt(t[2].length)),r[0]=t[0]-n,r[1]=t[0]+n,r},significant:function(n,t){return n<t}}),s.extend(s.fn,{normalci:function(n,t){return s.normalci(n,t,this.toArray())},tci:function(n,t){return s.tci(n,t,this.toArray())}}),s.extend(s.fn,{oneSidedDifferenceOfProportions:function(n,t,r,e){var i=o(n,t,r,e);return s.ztest(i,1)},twoSidedDifferenceOfProportions:function(n,t,r,e){var i=o(n,t,r,e);return s.ztest(i,2)}})}(this.jStat,Math),j=jStat,function(){var o=!1,a=/xyz/.test(function(){})?/\b_super\b/:/.*/;this.Class=function(){},Class.extend=function(n){function t(){!o&&this.init&&this.init.apply(this,arguments)}var i=this.prototype;o=!0;var r=new this;for(var e in o=!1,n)r[e]="function"==typeof n[e]&&"function"==typeof i[e]&&a.test(n[e])?function(r,e){return function(){var n=this._super;this._super=i[r];var t=e.apply(this,arguments);return this._super=n,t}}(e,n[e]):n[e];return t.prototype=r,(t.constructor=t).extend=arguments.callee,t}}(),jStat.ONE_SQRT_2PI=.3989422804014327,jStat.LN_SQRT_2PI=.9189385332046728,jStat.LN_SQRT_PId2=.22579135264472744,jStat.DBL_MIN=222507e-313,jStat.DBL_EPSILON=2220446049250313e-31,jStat.SQRT_32=5.656854249492381,jStat.TWO_PI=6.283185307179586,jStat.DBL_MIN_EXP=-999,jStat.SQRT_2dPI=.79788456080287,jStat.LN_SQRT_PI=.5723649429247,jStat.trunc=function(n){return 0<n?Math.floor(n):Math.ceil(n)},jStat.isFinite=function(n){return!isNaN(n)&&n!=Number.POSITIVE_INFINITY&&n!=Number.NEGATIVE_INFINITY};var Range=Class.extend({init:function(n,t,r){this._minimum=parseFloat(n),this._maximum=parseFloat(t),this._numPoints=parseFloat(r)},getMinimum:function(){return this._minimum},getMaximum:function(){return this._maximum},getNumPoints:function(){return this._numPoints},getPoints:function(){for(var n=[],t=this._minimum,r=(this._maximum-this._minimum)/(this._numPoints-1),e=0;e<this._numPoints;e++)n[e]=parseFloat(t.toFixed(6)),t+=r;return n}});Range.validate=function(n){return!(!n instanceof Range)&&!(isNaN(n.getMinimum())||isNaN(n.getMaximum())||isNaN(n.getNumPoints())||n.getMaximum()<n.getMinimum()||n.getNumPoints()<=0)};var ContinuousDistribution=Class.extend({init:function(n){this._name=n},toString:function(){return this._string},getName:function(){return this._name},getClassName:function(){return this._name+"Distribution"},density:function(n){if(isNaN(n)){if(Range.validate(n)){for(var t=n.getPoints(),r=[],e=0;e<t.length;e++)r[e]=parseFloat(this._pdf(t[e]));return r}throw"Invalid parameter supplied to "+this.getClassName()+".density()"}return parseFloat(this._pdf(n).toFixed(15))},cumulativeDensity:function(n){if(isNaN(n)){if(Range.validate(n)){for(var t=n.getPoints(),r=[],e=0;e<t.length;e++)r[e]=parseFloat(this._cdf(t[e]));return r}throw"Invalid parameter supplied to "+this.getClassName()+".cumulativeDensity()"}return parseFloat(this._cdf(n).toFixed(15))},getRange:function(n,t){null==n&&(n=5),null==t&&(t=100);var r=this.getMean()-n*Math.sqrt(this.getVariance()),e=this.getMean()+n*Math.sqrt(this.getVariance());return"GammaDistribution"==this.getClassName()||"LogNormalDistribution"==this.getClassName()?(r=0,e=this.getMean()+n*Math.sqrt(this.getVariance())):"BetaDistribution"==this.getClassName()&&(r=0,e=1),new Range(r,e,t)},getVariance:function(){},getMean:function(){},getQuantile:function(n){var s=this;return function n(t,r){for(var e=t.getPoints(),i=0,o=999,a=0;a<e.length;a++){var u=s.cumulativeDensity(e[a]),f=Math.abs(u-r);f<o&&(i=a,o=f)}return o<=1e-5?e[i]:n(new Range(e[i-1],e[i+1],20),r)}(this.getRange(5,20),n)}}),NormalDistribution=ContinuousDistribution.extend({init:function(n,t){this._super("Normal"),this._mean=parseFloat(n),this._sigma=parseFloat(t),this._string="Normal ("+this._mean.toFixed(2)+", "+this._sigma.toFixed(2)+")"},_pdf:function(n,t){null==t&&(t=!1);var r=this._sigma,e=this._mean;if(!jStat.isFinite(r))return t?Number.NEGATIVE_INFINITY:0;if(!jStat.isFinite(n)&&e==n)return Number.NaN;if(r<=0){if(r<0)throw"invalid sigma in _pdf";return n==e?Number.POSITIVE_INFINITY:t?Number.NEGATIVE_INFINITY:0}return n=(n-e)/r,jStat.isFinite(n)?t?-(jStat.LN_SQRT_2PI+.5*n*n+Math.log(r)):jStat.ONE_SQRT_2PI*Math.exp(-.5*n*n)/r:t?Number.NEGATIVE_INFINITY:0},_cdf:function(n,t,r){null==t&&(t=!0),null==r&&(r=!1);var e,i,o,a=this._mean,u=this._sigma;if(o=t?r?(i=Number.NEGATIVE_INFINITY,0):(i=0,1):r?(i=0,Number.NEGATIVE_INFINITY):(i=1,0),jStat.isFinite(n)||a!=n){if(Number.NaN,u<=0)return u<0?(console.warn("Sigma is less than 0"),Number.NaN):n<a?i:o;if(e=(n-a)/u,!jStat.isFinite(e))return n<a?i:o;var f=function(n,t,r,e,i){var o,a,u,f,s,c,l,m,d,p,h=[2.2352520354606837,161.02823106855587,1067.6894854603709,18154.98125334356,.06568233791820745],g=[47.202581904688245,976.0985517377767,10260.932208618979,45507.78933502673],v=[.39894151208813466,8.883149794388377,93.50665613217785,597.2702763948002,2494.5375852903726,6848.190450536283,11602.65143764735,9842.714838383978,1.0765576773720192e-8],b=[22.266688044328117,235.387901782625,1519.3775994075547,6485.558298266761,18615.571640885097,34900.95272114598,38912.00328609327,19685.429676859992],y=[.215898534057957,.12740116116024736,.022235277870649807,.0014216191932278934,29112874951168793e-21,.023073441764940174],w=[1.284260096144911,.4682382124808651,.06598813786892856,.0037823963320275824,7297515550839662e-20];if(s=.5*jStat.DBL_EPSILON,d=1!=e,p=0!=e,(l=Math.abs(n))<=.67448975){if(s<l)for(a=h[4]*(c=n*n),o=c,m=0;m<3;++m)a=(a+h[m])*c,o=(o+g[m])*c;else a=o=0;u=n*(a+h[3])/(o+g[3]),d&&(t=.5+u),p&&(r=.5-u),i&&(d&&(t=Math.log(t)),p&&(r=Math.log(r)))}else if(l<=jStat.SQRT_32){for(a=v[8]*l,o=l,m=0;m<7;++m)a=(a+v[m])*l,o=(o+b[m])*l;u=(a+v[7])/(o+b[7]),f=(n-(c=jStat.trunc(16*n)/16))*(n+c),i?(t=.5*-c*c+.5*-f+Math.log(u),(d&&0<n||p&&n<=0)&&(r=jStat.log1p(-Math.exp(.5*-c*c)*Math.exp(.5*-f)*u))):r=1-(t=Math.exp(.5*-c*c)*Math.exp(.5*-f)*u),0<n&&(u=t,d&&(t=r),r=u)}else if(i&&l<1e170||d&&-37.5193<n&&n<8.2924||p&&-8.2924<n&&n<37.5193){for(a=y[5]*(c=1/(n*n)),o=c,m=0;m<4;++m)a=(a+y[m])*c,o=(o+w[m])*c;u=c*(a+y[4])/(o+w[4]),u=(jStat.ONE_SQRT_2PI-u)/l,f=(n-(c=jStat.trunc(16*n)/16))*(n+c),i?(t=.5*-c*c+.5*-f+Math.log(u),(d&&0<n||p&&n<=0)&&(r=jStat.log1p(-Math.exp(.5*-c*c)*Math.exp(.5*-f)*u))):r=1-(t=Math.exp(.5*-c*c)*Math.exp(.5*-f)*u),0<n&&(u=t,d&&(t=r),r=u)}else r=0<n?(t=i?0:1,i?Number.NEGATIVE_INFINITY:0):(t=i?Number.NEGATIVE_INFINITY:0,i?0:1);return[t,r]}(n=e,e,void 0,!t,r);return t?f[0]:f[1]}},getMean:function(){return this._mean},getSigma:function(){return this._sigma},getVariance:function(){return this._sigma*this._sigma}});jStat.pnorm=function(n,t,r,e,i){null==t&&(t=0),null==r&&(r=1),null==e&&(e=!0),null==i&&(i=!1);var o=new NormalDistribution(t,r);if(isNaN(n)){if(n.length){for(var a=[],u=0;u<n.length;u++)a.push(o._cdf(n[u],e,i));return a}throw"Illegal argument: x"}return o._cdf(n,e,i)},function(d){function t(n,r,e){(n=n.not(".thead .tr")).each(function(n){var t=d(this);t.closest(".tr").index()==r&&t.val(e)})}function r(n){var t=d(".tr:nth-child("+(n+1)+")").not(".thead .tr"),r=p(t.find(".td-conversions input").val(),t.find(".td-visitors input").val());r+="%",t.find(".td-percent-value").html(r)}function p(n,t){return 0===t.length?0:r=100<(r=100*n/t)?100:isNaN(r)?0:r;var r}function e(){var n,t,r,e,i,o,a=d(".test-list .tbody .tr").slice(1),u=d(".tbody .tr .td-visitors .td-value"),f=0,s=0,c=999999999,l=0,m=p(d(".tbody .tr .td-conversions .td-value").first().val(),d(".tbody .tr .td-visitors .td-value").first().val());a.each(function(n,t){var r=p(d(this).find(".td-conversions .td-value").val(),d(this).find(".td-visitors .td-value").val());r<c&&(c=r,l=n),f<r&&(f=r,s=n)}),l++,s++,f<m?(f=m,s=0):(c=m,l=0),n=u.eq(s).val(),t=u.eq(l).val(),i=h(s+1),o=h(l+1),r=Math.ceil(100*(f-c)/c),e=Math.ceil(function(n,t,r,e){if(0<=n&&0<=t){n/=100,t/=100;var i=Math.sqrt(n*(1-n)/r),o=Math.sqrt(t*(1-t)/e),a=(n-t)/Math.sqrt(Math.pow(i,2)+Math.pow(o,2));return isNaN(a)?0:100*(t<n?jStat.pnorm(a,0,1,null,null):1-jStat.pnorm(a,0,1,null,null))}return 0}(f,c,n,t)),d(".table .tbody").find(".tr-winner").removeClass("tr-winner"),d(".table .tbody .tr:nth-child("+(s+1)+")").addClass("tr-winner"),d(".results-better").html(function(n,t,r,e){isFinite(r)?r+="%":r="";var i='<span class="results-converted"><strong>Test "'+n+'" converted '+r+' better</strong> than Test "'+t+'".</span> <span class="results-certainty">I&nbsp;am <strong>'+e+'%</strong> certain that the changes in Test "'+n+'" will improve your conversion rate.</span>';"br"==lang&&(i='<span class="results-converted"><strong>O Teste '+n+" converteu "+r+" melhor</strong> do que o Teste "+t+'.</span> <span class="results-certainty">Estou <strong>'+e+"%</strong> certo de que as mudanças no Teste "+n+" vão aumentar sua taxa de conversão.</span>");return i}(i,o,r,e)),d(".results-text").html(function(n){if(95<=n)var t=np_localize_ab.ab_calc_significant;else if(90<=n)t=np_localize_ab.ab_calc_significant_questionable;else t=np_localize_ab.ab_calc_significant_not;return t}(e)),d(".tr .td-percent-value").each(function(n){var t=d(this).text();t=Math.round(t.slice(0,-1)),isNaN(t)||d(this).html(t+"%")})}function h(n){return String.fromCharCode(64+n)}d(".td-visitors input").live("change keyup input",function(){var n=d(this).closest(".tr").index();t(d(".td-visitors input"),n,this.value),r(n),e()}),d(".td-conversions input").live("change keyup input",function(){var n=d(this).closest(".tr").index();t(d(".td-conversions input"),n,this.value),r(n),e()}),d(".test-more .btn").click(function(){var n=1+d(".test-list .tbody .tr").length;if(n<=4){var t=h(n);d(".test-list .tbody").append('<tr class="tr"><td class="td td-id">'+t+'</td><td class="td td-visitors"><div class="td-value-wrap"><input class="td-value" type="number" min="0" value="" placeholder="disabled" disabled="disabled" title="form boxes are disabled in cached version of the site" value="Disabled" style="opacity:0.3;" ></div></td><td class="td td-conversions"><div class="td-value-wrap"><input class="td-value" type="number" min="0" value="" placeholder="disabled" disabled="disabled" title="form boxes are disabled in cached version of the site" value="Disabled" style="opacity:0.3;" ></div></td> <td class="td td-gt"><span class="td-gt-value">&gt;</span></td><td class="td td-percent"><span class="td-percent-value"></span></td></tr>'),4==n&&d(".test-more .btn").addClass("btn-disabled")}}),d(document).ready(function(){e()})}(jQuery);