phpOMS/tests/Math/Stochastic/Distribution/ChiSquaredDistributionTest.php

<?php
/**
 * Orange Management
 *
 * PHP Version 7.2
 *
 * @package    tests
 * @copyright  Dennis Eichhorn
 * @license    OMS License 1.0
 * @version    1.0.0
 * @link       http://website.orange-management.de
 */

namespace phpOMS\tests\Math\Stochastic\Distribution;

use phpOMS\Math\Stochastic\Distribution\ChiSquaredDistribution;

class ChiSquaredDistributionTest extends \PHPUnit\Framework\TestCase
{
    public function testHypothesisFalse() : void
    {
        $p = [0.6, 0.25, 0.15];
        $a = 0.05;

        $total    = 470;
        $observed = [255, 125, 90];
        $expected = [$total * $p[0], $total * $p[1], $total * $p[2]];

        $test = ChiSquaredDistribution::testHypothesis($observed, $expected, $a);

        self::assertEquals(8.46, $test['Chi2'], '', 0.1);
        self::assertNotEquals(0, $test['P']);
        self::assertFalse($test['H0']);
        self::assertEquals(2, $test['df']);
    }

    public function testDegreesOfFreedom() : void
    {
        self::assertEquals(2, ChiSquaredDistribution::getDegreesOfFreedom([1, 2, 3]));
        self::assertEquals(6, ChiSquaredDistribution::getDegreesOfFreedom([
            [1, 2, 3, 4],
            [1, 2, 3, 4],
            [1, 2, 3, 4],
        ]));
    }

    public function testMode() : void
    {
        self::assertEquals(max(5 - 2, 0), ChiSquaredDistribution::getMode(5));
    }

    public function testMean() : void
    {
        $df = 5;

        self::assertEquals($df, ChiSquaredDistribution::getMean($df));
    }

    public function testVariance() : void
    {
        $df = 5;

        self::assertEquals(2 * $df, ChiSquaredDistribution::getVariance($df));
    }

    public function testMedian() : void
    {
        $df = 5;

        self::assertEquals($df * (1 - 2 / (9 * $df)) ** 3, ChiSquaredDistribution::getMedian($df));
    }

    public function testSkewness() : void
    {
        $df = 5;

        self::assertEquals(sqrt(8 / $df), ChiSquaredDistribution::getSkewness($df));
    }

    public function testExKurtosis() : void
    {
        $df = 5;

        self::assertEquals(12 / $df, ChiSquaredDistribution::getExKurtosis($df));
    }

    public function testMgdf() : void
    {
        $df = 5;
        $t  = 0.3;

        self::assertEquals((1 - 2 * $t) ** (-$df / 2), ChiSquaredDistribution::getMgf($df, $t));
    }

    public function testHypothesisSizeException() : void
    {
        self::expectedException(\Exception::class);

        ChiSquaredDistribution::testHypothesis([1, 2], [2]);
    }

    public function testHypothesisDegreesOfFreedomException() : void
    {
        self::expectedException(\Exception::class);

        ChiSquaredDistribution::testHypothesis([], []);
    }

    public function testPdfOutOfBoundsException() : void
    {
        self::expectedException(\OutOfBoundsException::class);

        ChiSquaredDistribution::getPdf(-1, 0);
    }

    public function testMgfOutOfBoundsException() : void
    {
        self::expectedException(\OutOfBoundsException::class);

        ChiSquaredDistribution::getMgf(1, 0.6);
    }
}