/* * Apache License * Version 2.0, January 2004 * http://www.apache.org/licenses/ * * Copyright 2013 Aurelian Tutuianu * Copyright 2014 Aurelian Tutuianu * Copyright 2015 Aurelian Tutuianu * Copyright 2016 Aurelian Tutuianu * * Licensed under the Apache License, Version 2.0 (the "License"); * you may not use this file except in compliance with the License. * You may obtain a copy of the License at * * http://www.apache.org/licenses/LICENSE-2.0 * * Unless required by applicable law or agreed to in writing, software * distributed under the License is distributed on an "AS IS" BASIS, * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. * See the License for the specific language governing permissions and * limitations under the License. * */ package rapaio.core.distributions; import static rapaio.sys.WS.formatFlex; /** * Continuous uniform distribution * * @author <a href="mailto:padreati@yahoo.com">Aurelian Tutuianu</a> */ public class Uniform implements Distribution { private static final long serialVersionUID = -6077483164719205038L; private final double a; private final double b; public Uniform(double a, double b) { this.a = a; this.b = b; } public double a() { return a; } public double b() { return b; } @Override public String name() { return "Uniform(a=" + formatFlex(a) + ",b=" + formatFlex(b) + ")"; } @Override public boolean discrete() { return false; } @Override public double pdf(double x) { if (x < a || x > b) { return 0; } if (a == b) { return 1; } return 1 / (b - a); } @Override public double cdf(double x) { if (x < a) { return 0; } if (x > b) { return 1; } return (x - a) / (b - a); } @Override public double quantile(double p) { if (p < 0 || p > 1) { throw new IllegalArgumentException("probability value should lie in [0,1] interval"); } return a + p * (b - a); } @Override public double min() { return a; } @Override public double max() { return b; } @Override public double mean() { return a + (b - a) / 2.0; } @Override public double mode() { return mean(); } @Override public double var() { return Math.pow(b - a, 2) / 12.0; } @Override public double skewness() { return 0; } @Override public double kurtosis() { return -6.0 / 5.0; } @Override public double entropy() { return Math.log(b - a); } }