org.jgap.gp.function.statistics.Kurtosis.java Source code

Introduction

Here is the source code for org.jgap.gp.function.statistics.Kurtosis.java

Source

/*
 * This file is part of JGAP.
 *
 * JGAP offers a dual license model containing the LGPL as well as the MPL.
 *
 * For licensing information please see the file license.txt included with JGAP
 * or have a look at the top of class org.jgap.Chromosome which representatively
 * includes the JGAP license policy applicable for any file delivered with JGAP.
 */
package org.jgap.gp.function.statistics;

import org.apache.commons.math.stat.descriptive.*;
import org.jgap.*;
import org.jgap.gp.*;
import org.jgap.gp.impl.*;
import org.jgap.util.*;

/**
 * Computes the Kurtosis of the available values.
 * <p>
 * We use the following (unbiased) formula to define kurtosis:</p>
 *  <p>
 *  kurtosis = { [n(n+1) / (n -1)(n - 2)(n-3)] sum[(x_i - mean)^4] / std^4 } - [3(n-1)^2 / (n-2)(n-3)]
 *  </p><p>
 *  where n is the number of values, and std is the
 * Standard Deviation</p>
 *
 * @author Luis Garcia
 * @since 3.5
 */
public class Kurtosis extends CommandDynamicArity implements ICloneable {
    public Kurtosis(GPConfiguration a_conf, int a_arityInitial, int a_arityMin, int a_arityMax, Class a_returnType)
            throws InvalidConfigurationException {
        super(a_conf, a_arityInitial, a_arityMin, a_arityMax, a_returnType);
    }

    @Override
    public String toString() {
        String s = "Kurtosis(";
        int size = size();
        for (int i = 0; i < size; i++) {
            if (i > 0) {
                s += ";";
            }
            s += "&" + (i + 1);
        }
        return s + ")";
    }

    @Override
    public double execute_double(ProgramChromosome c, int n, Object[] args) {
        int size = size();
        DescriptiveStatistics stats = new DescriptiveStatistics();
        for (int i = 0; i < size; i++) {
            stats.addValue(c.execute_double(n, i, args));
        }
        return stats.getKurtosis();
    }

    @Override
    public float execute_float(ProgramChromosome c, int n, Object[] args) {
        int size = size();
        DescriptiveStatistics stats = new DescriptiveStatistics();
        for (int i = 0; i < size; i++) {
            stats.addValue(c.execute_float(n, i, args));
        }
        return (float) stats.getKurtosis();
    }

    @Override
    public String getName() {
        return "Kurtosis";
    }

    @Override
    public Object clone() {
        try {
            Kurtosis result = new Kurtosis(getGPConfiguration(), getArity(null), getArityMin(), getArityMax(),
                    getReturnType());
            return result;
        } catch (Exception ex) {
            throw new CloneException(ex);
        }
    }
}