Microsoft MVP성태의 닷넷 이야기
글쓴 사람
홈페이지
첨부 파일

MathNET + OxyPlot을 이용한 간단한 통계 정보 처리 - Histogram

MathNET Numerics 라이브러리의,

C# - MathNet.Numerics의 Matrix(행렬) 연산
; https://www.sysnet.pe.kr/2/0/11910

통계 관련한 간단한 사용법을 볼까요? ^^ 우선, 다양한 평균에 대한 메서드입니다.

int [] heightInfo = LoadData("data.txt");

double [] dblHeights = (from elem in heightInfo
                    select (double)elem).ToArray();

Console.WriteLine($"Mean: {Statistics.Mean(dblHeights)}"); // 산술평균
Console.WriteLine($"GeometricMean: {Statistics.GeometricMean(dblHeights)}"); // 기하평균
Console.WriteLine($"HarmonicMean: {Statistics.HarmonicMean(dblHeights)}"); // 조화평균
Console.WriteLine($"RootMeanSquare: {Statistics.RootMeanSquare(dblHeights)}"); // 제곱평균

자동 히스토그램은 구간 수만 정해 주면 됩니다.

Histogram hist = new Histogram(dblHeights, 6);

또는, 상/하위 구간을 지정해 줄 수도 있습니다.

Histogram hist = new Histogram(dblHeights, 6, 141, 170);

아니면, 구간 자체를 사용자 정의하는 것도 가능합니다.

Histogram hist = new Histogram();

hist.AddBucket(new Bucket(141, 145));
hist.AddBucket(new Bucket(145, 150));
hist.AddBucket(new Bucket(150, 155));
hist.AddBucket(new Bucket(155, 160));
hist.AddBucket(new Bucket(160, 165));
hist.AddBucket(new Bucket(165, 170));

hist.AddData(dblHeights);

보는 바와 같이 계급 구간을 Bucket 타입으로 지정하는데,

public Bucket(double lowerBound, double upperBound);

상/하한 기준 값이 double이라는 점과 데이터 비교 기준이 다음과 같다는 점을 염두에 두어야 합니다.

(lowerBound, upperBound]

lowerBound < x <= upperBound

따라서 145 값을 주면 0번 Bucket에 계산되지만, 145.1을 주면 1번 Bucket에 계산됩니다.

상/하한 기준 값이 double이라는 점 때문에 계급값을 범용적으로 계산하는 것도 좀 이상해지긴 하지만, 그런 탓인지 Bucket 타입에는 계급값에 대한 속성이 없습니다. 또한, Histogram의 경우 계급값과 상대도수를 이용한 평균을 내주는 기능도 없으므로 이런 부분들은 임의로 여러분들이 다음과 같은 확장 메서드를 만들어 주어야 합니다.

public static class BucketExtension
{
    public static double GetClassMark(this Bucket bucket)
    {
        return Math.Round((bucket.LowerBound + bucket.UpperBound) / 2, MidpointRounding.AwayFromZero);
    }

    public static double GetRelativeFrequency(this Bucket bucket, Histogram histogram)
    {
        return bucket.Count / histogram.DataCount;
    }
}

public static class HistogramExtension
{
    public static void Output(this Histogram histogram)
    {
        for (int i = 0; i < histogram.BucketCount; i++)
        {
            Console.WriteLine($"{histogram[i].GetClassMark()} == {histogram[i].Count}");
        }
    }

    public static double Mean(this Histogram histogram)
    {
        double mean = 0.0;

        for (int i = 0; i < histogram.BucketCount; i++)
        {
            Bucket bucket = histogram[i];

            double classMark = bucket.GetClassMark();
            double relativeFrequency = bucket.GetRelativeFrequency(histogram);

            mean += (classMark * relativeFrequency);
        }

        return mean;
    }
}

그래서 다음과 같은 정도로 사용해 줄 수 있습니다.

Histogram hist = new Histogram();

hist.AddBucket(new Bucket(141, 145));
hist.AddBucket(new Bucket(145, 150));
hist.AddBucket(new Bucket(150, 155));
hist.AddBucket(new Bucket(155, 160));
hist.AddBucket(new Bucket(160, 165));
hist.AddBucket(new Bucket(165, 170));

hist.AddData(dblHeights);
hist.Output();
Console.WriteLine($"Mean of histogram: {hist.GetMean()}");

/*
출력 결과:
143 == 1
148 == 6
153 == 19
158 == 30
163 == 18
168 == 6
Mean of histogram: 157.75
*/




그래도 히스토그램은 역시 그래프로 봐야 제맛이죠? ^^

OxyPlot을 이용해,

C# Plotting 라이브러리 OxyPlot
; https://www.sysnet.pe.kr/2/0/10973

C# - OxyPlot 라이브러리로 복소수 표현
; https://www.sysnet.pe.kr/2/0/10974

Docs / Model / Series
; http://docs.oxyplot.org/en/latest/models/series/index.html

다음과 같은 정도로 그려낼 수 있습니다.

Histogram hist = GetHistogram(dblHeights);
plotView1.Model = GetHistogramModel(hist);

public PlotModel GetHistogramModel(Histogram hist)
{
    PlotModel model = new PlotModel();

    CategoryAxis xaxis = new CategoryAxis();
    xaxis.Position = AxisPosition.Bottom;
    xaxis.MajorGridlineStyle = LineStyle.Solid;
    xaxis.MinorGridlineStyle = LineStyle.Dot;

    for (int i = 0; i < hist.BucketCount; i++)
    {
        xaxis.Labels.Add(hist[i].GetClassMark().ToString());
    }

    LinearAxis yaxis = new LinearAxis();
    yaxis.Position = AxisPosition.Left;
    yaxis.MajorGridlineStyle = LineStyle.Dot;
    xaxis.MinorGridlineStyle = LineStyle.Dot;

    ColumnSeries series = new ColumnSeries();
    series.IsStacked = true;
    series.LabelPlacement = LabelPlacement.Base;
    series.LabelFormatString = "{0}";

    for (int i = 0; i < hist.BucketCount; i++)
    {
        ColumnItem item = new ColumnItem(hist[i].Count);
        series.Items.Add(item);
    }

    model = new PlotModel();
    model.Title = "Histogram";

    model.Axes.Add(xaxis);
    model.Axes.Add(yaxis);
    model.Series.Add(series);

    return model;
}

plotview_1.png

(첨부 파일은 이 글의 소스 코드를 포함합니다.)




[이 글에 대해서 여러분들과 의견을 공유하고 싶습니다. 틀리거나 미흡한 부분 또는 의문 사항이 있으시면 언제든 댓글 남겨주십시오.]

[연관 글]





[최초 등록일: ]
[최종 수정일: 5/24/2019 ]

Creative Commons License
이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-변경금지 2.0 대한민국 라이센스에 따라 이용하실 수 있습니다.
by SeongTae Jeong, mailto:techsharer@outlook.com

비밀번호

댓글 쓴 사람
 




... [16]  17  18  19  20  21  22  23  24  25  26  27  28  29  30  ...
NoWriterDateCnt.TitleFile(s)
11889정성태5/8/20191249개발 환경 구성: 438. mstest, QTAgent의 로그 파일 설정 방법
11888정성태5/8/20192364.NET Framework: 830. C# - 비동기 호출을 취소하는 CancellationToken의 간단한 예제 코드파일 다운로드1
11887정성태5/8/20191556.NET Framework: 829. C# - yield 문을 사용할 수 있는 메서드의 조건
11886정성태5/8/20191583오류 유형: 534. mstest.exe 실행 시 "Visual Studio Enterprise is required to execute the test." 오류 [2]
11885정성태5/7/20191028오류 유형: 533. mstest.exe 실행 시 "File extension specified '.loadtest' is not a valid test extension." 오류 발생
11884정성태5/5/20191845.NET Framework: 828. C# DLL에서 Win32 C/C++처럼 dllexport 함수를 제공하는 방법 - 두 번째 이야기
11883정성태5/3/20192259.NET Framework: 827. C# - 인터넷 시간 서버로부터 받은 시간을 윈도우에 적용하는 방법파일 다운로드1
11882정성태5/9/20191215.NET Framework: 826. (번역글) .NET Internals Cookbook Part 11 - Various C# riddles파일 다운로드1
11881정성태4/28/20191805오류 유형: 532. .NET Core 프로젝트로 마이그레이션 시 "CS0579 Duplicate 'System.Reflection.AssemblyCompanyAttribute' attribute" 오류 발생
11880정성태4/25/20191107오류 유형: 531. 이벤트 로그 오류 - Task Scheduling Error: m->NextScheduledSPRetry 1547, m->NextScheduledEvent 1547
11879정성태1/27/20201528.NET Framework: 825. (번역글) .NET Internals Cookbook Part 10 - Threads, Tasks, asynchronous code and others파일 다운로드2
11878정성태5/9/20191399.NET Framework: 824. (번역글) .NET Internals Cookbook Part 9 - Finalizers, queues, card tables and other GC stuff파일 다운로드1
11877정성태5/9/20191463.NET Framework: 823. (번역글) .NET Internals Cookbook Part 8 - C# gotchas파일 다운로드1
11876정성태5/9/20191223.NET Framework: 822. (번역글) .NET Internals Cookbook Part 7 - Word tearing, locking and others파일 다운로드1
11875정성태4/21/20191109오류 유형: 530. Visual Studo에서 .NET Core 프로젝트를 열 때 "One or more errors occurred." 오류 발생
11874정성태5/9/20191264.NET Framework: 821. (번역글) .NET Internals Cookbook Part 6 - Object internals파일 다운로드1
11873정성태7/9/20201081.NET Framework: 820. (번역글) .NET Internals Cookbook Part 5 - Methods, parameters, modifiers파일 다운로드1
11872정성태5/9/20191331.NET Framework: 819. (번역글) .NET Internals Cookbook Part 4 - Type members파일 다운로드1
11871정성태5/9/20191590.NET Framework: 818. (번역글) .NET Internals Cookbook Part 3 - Initialization tricks [3]파일 다운로드1
11870정성태4/16/20191276.NET Framework: 817. Process.Start로 실행한 콘솔 프로그램의 출력 결과를 얻는 방법파일 다운로드1
11869정성태5/9/20191338.NET Framework: 816. (번역글) .NET Internals Cookbook Part 2 - GC-related things파일 다운로드1
11868정성태4/15/20191155.NET Framework: 815. CER(Constrained Execution Region)이란?파일 다운로드1
11867정성태4/15/20191236.NET Framework: 814. Critical Finalizer와 SafeHandle의 사용 의미파일 다운로드1
11866정성태4/9/20192519Windows: 159. 네트워크 공유 폴더(net use)에 대한 인증 정보는 언제까지 유효할까요?
11865정성태4/9/20191225오류 유형: 529. 제어판 - C:\ProgramData\Microsoft\Windows\Start Menu\Programs\Administrative Tools is not accessible.
11864정성태4/9/20191375오류 유형: 528. '...' could be '0': this does not adhere to the specification for the function '...'
... [16]  17  18  19  20  21  22  23  24  25  26  27  28  29  30  ...