probdist/distributions/f.js

Probability Distributions

var distribution = require('../helpers/distribution'),
	beta = require('../helpers/beta');
	
module.exports = function(v1, v2) {
	var beta_v1_v2 = beta(v1 / 2, v2 / 2);
	
	return distribution({
		pdf: function(x) {
			return x >= 0 ? 
			Math.sqrt(Math.pow(v1 * x, v1) * Math.pow(v2, v2) / 
				Math.pow(v1 * x + v2, v1 + v2)) /
			(x * beta_v1_v2) : 0;
		},
		
		mean: v2 > 2 ? v2 / (v2 - 2) : undefined,
		
		variance: v2 > 4 ? 
			2 * Math.pow(v2, 2) * (v1 + v2 - 2) /
			(v1 * Math.pow(v2 - 2, 2) * (v2 - 4)) :
			undefined
	});
};