phpOMS/Business/Recommendation/MemoryCF.php

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

namespace phpOMS\Business\Recommendation;

use phpOMS\Math\Topology\MetricsND;

/**
 * Memory based collaborative filtering
 *
 * Items or potential customers are found based on how much they like certain items.
 *
 * This requires a item/product rating of some sort in the backend.
 * Such a rating could be either manual user ratings or a rating based on how often it is purchased or how long it is used.
 * Most likely a combination is required.
 *
 * @package phpOMS\Business\Recommendation
 * @license OMS License 2.0
 * @link    https://jingga.app
 * @see     https://realpython.com/build-recommendation-engine-collaborative-filtering/
 * @since   1.0.0
 */
final class MemoryCF
{
    /**
     * All rankings
     *
     * @var array<array>
     * @since 1.0.0
     */
    private array $rankings = [];

    /**
     * Constructor.
     *
     * @param array<array> $rankings Array of item ratings by users (or reverse to find users)
     *
     * @since 1.0.0
     */
    public function __construct(array $rankings)
    {
        $this->rankings = $this->normalizeRanking($rankings);
    }

    /**
     * Normalize all ratings.
     *
     * This is necessary because some users my give lower or higher ratings on average (bias).
     *
     * @param array<array> $rankings Item ratings/rankings
     *
     * @return array<array>
     *
     * @since 1.0.0
     */
    private function normalizeRanking(array $rankings) : array
    {
        foreach ($rankings as $idx => $items) {
            $avg = \array_sum($items) / \count($items);

            foreach ($items as $idx2 => $_) {
                $rankings[$idx][$idx2] -= $avg;
            }
        }

        return $rankings;
    }

    /**
     * Euclidean distance between users
     *
     * @param array        $ranking  Rating to find the distance for
     * @param array<array> $rankings All ratings to find the distance to
     *
     * @return float[]
     *
     * @since 1.0.0
     */
    private function euclideanDistance(array $ranking, array $rankings) : array
    {
        $distances = [];
        foreach ($rankings as $idx => $r) {
            $distances[$idx] = \abs(MetricsND::euclidean($ranking, $r));
        }

        return $distances;
    }

    /**
     * Cosine distance between users
     *
     * @param array        $ranking  Rating to find the distance for
     * @param array<array> $rankings All ratings to find the distance to
     *
     * @return float[]
     *
     * @since 1.0.0
     */
    private function cosineDistance(array $ranking, array $rankings) : array
    {
        $distances = [];
        foreach ($rankings as $idx => $r) {
            $distances[$idx] = \abs(MetricsND::cosine($ranking, $r));
        }

        return $distances;
    }

    /**
     * Assign a item rank/rating based on the distance to other items
     *
     * @param string       $itemId    Id of the item to rank
     * @param array        $distances Distance to other users
     * @param array<array> $users     All user ratings
     * @param int          $size      Only consider the top n distances (best matches with other users)
     *
     * @return float Estimated item rank/rating based on similarity to other users
     *
     * @since 1.0.0
     */
    private function weightedItemRank(string $itemId, array $distances, array $users, int $size) : float
    {
        $rank  = 0.0;
        $count = 0;
        foreach ($distances as $uId => $_) {
            if ($count >= $size) {
                break;
            }

            if (!isset($users[$itemId])) {
                continue;
            }

            ++$count;
            $rank += $users[$uId][$itemId];
        }

        return $count === 0 ? 0.0 : $rank / $count;
    }

    /**
     * Find similar users
     *
     *
     * @param array $ranking Array of item ratings (e.g. products, movies, ...)
     *
     * @return array
     *
     * @since 1.0.0
     */
    public function bestMatchUser(array $ranking, int $size = 10) : array
    {
        $ranking = $this->normalizeRanking([$ranking]);
        $ranking = $ranking[0];

        $euclidean = $this->euclideanDistance($ranking, $this->rankings);
        // $cosine    = $this->cosineDistance($ranking, $this->rankings);

        \asort($euclidean);
        // \asort($cosine);

        $size    = \min($size, \count($this->rankings));
        $matches = [];

        $distancePointer = \array_keys($euclidean);
        // $anglePointer    = \array_keys($cosine);

        // Inspect items of the top n comparable users
        for ($i = 0; $i < $size; ++$i) {
            // $uId = $i % 2 === 0 ? $distancePointer[$i] : $anglePointer[$i];
            $uId = $distancePointer[$i];

            if (!\in_array($uId, $matches)) {
                $matches[] = $uId;
            }
        }

        return $matches;
    }

    /**
     * Find potential users which are a good match for a user.
     *
     * @param array $ranking Array of item ratings (e.g. products, movies, ...)
     *
     * @return array
     *
     * @since 1.0.0
     */
    public function bestMatchItem(array $ranking, int $size = 10) : array
    {
        $ranking = $this->normalizeRanking([$ranking]);
        $ranking = $ranking[0];

        $euclidean = $this->euclideanDistance($ranking, $this->rankings);
        $cosine    = $this->cosineDistance($ranking, $this->rankings);

        \asort($euclidean);
        \asort($cosine);

        $size    = \min($size, \count($this->rankings));
        $matches = [];

        $distancePointer = \array_keys($euclidean);
        $anglePointer    = \array_keys($cosine);

        // Inspect items of the top n comparable users
        for ($i = 0; $i < $size; ++$i) {
            $uId       = $i % 2 === 0 ? $distancePointer[$i] : $anglePointer[$i];
            $distances = $i % 2 === 0 ? $euclidean : $cosine;

            foreach ($this->rankings[$uId] as $iId => $_) {
                // Item is already in dataset or in historic dataset (we are only interested in new)
                if (isset($matches[$iId]) || isset($ranking[$iId])) {
                    continue;
                }

                // Calculate the expected rating the user would give based on what the best comparable users did
                $matches[$iId] = $this->weightedItemRank((string) $iId, $distances, $this->rankings, $size);
            }
        }

        \asort($matches);

        return \array_reverse($matches, true);
    }
}