C# を使用してクラスター分析アルゴリズムを作成する方法

using System;
using System.Collections.Generic;
using System.Linq;

public class KMeans
{
    public List<List<double>> Cluster(List<List<double>> data, int k, int maxIterations)
    {
        // 初始化聚类中心
        List<List<double>> centroids = InitializeCentroids(data, k);
        
        for (int i = 0; i < maxIterations; i++)
        {
            // 创建临时的聚类结果
            List<List<List<double>>> clusters = new List<List<List<double>>>();
            for (int j = 0; j < k; j++)
            {
                clusters.Add(new List<List<double>>());
            }
            
            // 将数据样本分配到最近的聚类中心
            foreach (var point in data)
            {
                int nearestCentroidIndex = FindNearestCentroidIndex(point, centroids);
                clusters[nearestCentroidIndex].Add(point);
            }
            
            // 更新聚类中心
            List<List<double>> newCentroids = new List<List<double>>();
            for (int j = 0; j < k; j++)
            {
                newCentroids.Add(UpdateCentroid(clusters[j]));
            }
            
            // 判断聚类结果是否变化，若不再变化则停止迭代
            if (CentroidsNotChanged(centroids, newCentroids))
            {
                break;
            }
            
            centroids = newCentroids;
        }
        
        return centroids;
    }

    private List<List<double>> InitializeCentroids(List<List<double>> data, int k)
    {
        List<List<double>> centroids = new List<List<double>>();
        Random random = new Random();

        for (int i = 0; i < k; i++)
        {
            int randomIndex = random.Next(data.Count);
            centroids.Add(data[randomIndex]);
            data.RemoveAt(randomIndex);
        }

        return centroids;
    }

    private int FindNearestCentroidIndex(List<double> point, List<List<double>> centroids)
    {
        int index = 0;
        double minDistance = double.MaxValue;

        for (int i = 0; i < centroids.Count; i++)
        {
            double distance = CalculateDistance(point, centroids[i]);
            if (distance < minDistance)
            {
                minDistance = distance;
                index = i;
            }
        }

        return index;
    }

    private double CalculateDistance(List<double> PointA, List<double> PointB)
    {
        double sumSquaredDifferences = 0;
        for (int i = 0; i < PointA.Count; i++)
        {
            sumSquaredDifferences += Math.Pow(PointA[i] - PointB[i], 2);
        }

        return Math.Sqrt(sumSquaredDifferences);
    }

    private List<double> UpdateCentroid(List<List<double>> cluster)
    {
        int dimension = cluster[0].Count;
        List<double> centroid = new List<double>();

        for (int i = 0; i < dimension; i++)
        {
            double sum = 0;
            foreach (var point in cluster)
            {
                sum += point[i];
            }
            centroid.Add(sum / cluster.Count);
        }

        return centroid;
    }

    private bool CentroidsNotChanged(List<List<double>> oldCentroids, List<List<double>> newCentroids)
    {
        for (int i = 0; i < oldCentroids.Count; i++)
        {
            for (int j = 0; j < oldCentroids[i].Count; j++)
            {
                if (Math.Abs(oldCentroids[i][j] - newCentroids[i][j]) > 1e-6)
                {
                    return false;
                }
            }
        }

        return true;
    }
}

class Program
{
    static void Main(string[] args)
    {
        // 假设我们有以下数据样本
        List<List<double>> data = new List<List<double>>()
        {
            new List<double>() {1, 1},
            new List<double>() {1, 2},
            new List<double>() {2, 1},
            new List<double>() {2, 2},
            new List<double>() {5, 6},
            new List<double>() {6, 5},
            new List<double>() {6, 6},
            new List<double>() {7, 5},
        };

        KMeans kmeans = new KMeans();
        List<List<double>> centroids = kmeans.Cluster(data, 2, 100);

        Console.WriteLine("聚类中心：");
        foreach (var centroid in centroids)
        {
            Console.WriteLine(string.Join(", ", centroid));
        }
    }
}