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

<?php
/**
 * Jingga
 *
 * PHP Version 8.1
 *
 * @package   tests
 * @copyright Dennis Eichhorn
 * @license   OMS License 2.0
 * @version   1.0.0
 * @link      https://jingga.app
 */
declare(strict_types=1);

namespace phpOMS\tests\Math\Stochastic\Distribution;

use phpOMS\Math\Stochastic\Distribution\ParetoDistribution;

/**
 * @internal
 */
final class ParetoDistributionTest extends \PHPUnit\Framework\TestCase
{
    /**
     * @covers phpOMS\Math\Stochastic\Distribution\ParetoDistribution
     * @group framework
     */
    public function testPdf() : void
    {
        self::assertEqualsWithDelta(0.263374485596, ParetoDistribution::getPdf(3, 2, 4), 0.001);
    }

    /**
     * @covers phpOMS\Math\Stochastic\Distribution\ParetoDistribution
     * @group framework
     */
    public function testCdf() : void
    {
        self::assertEqualsWithDelta(0.8024691358, ParetoDistribution::getCdf(3, 2, 4), 0.001);
    }

    /**
     * @covers phpOMS\Math\Stochastic\Distribution\ParetoDistribution
     * @group framework
     */
    public function testMean() : void
    {
        self::assertEqualsWithDelta(8 / 3, ParetoDistribution::getMean(2, 4), 0.001);
        self::assertEquals(\PHP_FLOAT_MAX, ParetoDistribution::getMean(2, 1));
    }

    /**
     * @covers phpOMS\Math\Stochastic\Distribution\ParetoDistribution
     * @group framework
     */
    public function testVariance() : void
    {
        self::assertEqualsWithDelta(2, ParetoDistribution::getVariance(3, 4), 0.001);
        self::assertEqualsWithDelta(\PHP_FLOAT_MAX, ParetoDistribution::getVariance(3, 2), 0.001);
    }

    /**
     * @covers phpOMS\Math\Stochastic\Distribution\ParetoDistribution
     * @group framework
     */
    public function testStandardDeviation() : void
    {
        self::assertEqualsWithDelta(\sqrt(2), ParetoDistribution::getStandardDeviation(3, 4), 0.001);
    }

    /**
     * @covers phpOMS\Math\Stochastic\Distribution\ParetoDistribution
     * @group framework
     */
    public function testExKurtosis() : void
    {
        self::assertEqualsWithDelta(35.666666666666664, ParetoDistribution::getExKurtosis(6), 0.001);
        self::assertEquals(0.0, ParetoDistribution::getExKurtosis(4));
    }

    /**
     * @covers phpOMS\Math\Stochastic\Distribution\ParetoDistribution
     * @group framework
     */
    public function testSkewness() : void
    {
        self::assertEqualsWithDelta(3.810317377662722, ParetoDistribution::getSkewness(6), 0.001);
        self::assertEquals(0.0, ParetoDistribution::getSkewness(3));
    }

    /**
     * @covers phpOMS\Math\Stochastic\Distribution\ParetoDistribution
     * @group framework
     */
    public function testMedian() : void
    {
        self::assertEquals(3 * \pow(2, 1 / 4), ParetoDistribution::getMedian(3, 4));
    }

    /**
     * @covers phpOMS\Math\Stochastic\Distribution\ParetoDistribution
     * @group framework
     */
    public function testMode() : void
    {
        self::assertEquals(3, ParetoDistribution::getMode(3));
    }

    /**
     * @covers phpOMS\Math\Stochastic\Distribution\ParetoDistribution
     * @group framework
     */
    public function testEntropy() : void
    {
        self::assertEquals(
            \log(3 / 4 * \exp(1 + 1 / 4)),
            ParetoDistribution::getEntropy(3, 4)
        );
    }

    /**
     * @covers phpOMS\Math\Stochastic\Distribution\ParetoDistribution
     * @group framework
     */
    public function testFisherInformation() : void
    {
        self::assertEquals(
            [
                [4 / (3 ** 2), -1 / 3],
                [-1 / 3, 1 / (4 ** 2)],
            ],
            ParetoDistribution::getFisherInformation(3, 4)
        );
    }
}