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

비밀번호

댓글 쓴 사람
 




1  2  3  4  [5]  6  7  8  9  10  11  12  13  14  15  ...
NoWriterDateCnt.TitleFile(s)
12108정성태1/10/2020151오류 유형: 587. Kernel Driver 시작 시 127(The specified procedure could not be found.) 오류 메시지 발생
12107정성태1/10/2020266.NET Framework: 877. C# - 프로세스의 모든 핸들을 열람 - 두 번째 이야기
12106정성태1/8/2020300VC++: 136. C++ - OSR Driver Loader와 같은 Legacy 커널 드라이버 설치 프로그램 제작 [1]
12105정성태1/8/2020220디버깅 기술: 153. C# - PEB를 조작해 로드된 DLL을 숨기는 방법
12104정성태1/9/2020435DDK: 9. 커널 메모리를 읽고 쓰는 NT Legacy driver와 C# 클라이언트 프로그램 [2]
12103정성태1/12/2020607DDK: 8. Visual Studio 2019 + WDK Legacy Driver 제작- Hello World 예제 [1]파일 다운로드2
12102정성태1/6/2020314디버깅 기술: 152. User 권한(Ring 3)의 프로그램에서 _ETHREAD 주소(및 커널 메모리를 읽을 수 있다면 _EPROCESS 주소) 구하는 방법
12101정성태1/8/2020268.NET Framework: 876. C# - PEB(Process Environment Block)를 통해 로드된 모듈 목록 열람
12100정성태1/3/2020214.NET Framework: 875. .NET 3.5 이하에서 IntPtr.Add 사용
12099정성태1/3/2020297디버깅 기술: 151. Windows 10 - Process Explorer로 확인한 Handle 정보를 windbg에서 조회
12098정성태1/2/2020292.NET Framework: 874. C# - 커널 구조체의 Offset 값을 하드 코딩하지 않고 사용하는 방법
12097정성태1/2/2020189디버깅 기술: 150. windbg - Wow64, x86, x64에서의 커널 구조체(예: TEB) 구조체 확인
12096정성태1/2/2020345디버깅 기술: 149. C# - DbgEng.dll을 이용한 간단한 디버거 제작
12095정성태12/27/2019339VC++: 135. C++ - string_view의 동작 방식
12094정성태12/26/2019353.NET Framework: 873. C# - 코드를 통해 PDB 심벌 파일 다운로드 방법
12093정성태12/26/2019531.NET Framework: 872. C# - 로딩된 Native DLL의 export 함수 목록 출력파일 다운로드1
12092정성태12/25/2019343디버깅 기술: 148. cdb.exe를 이용해 (ntdll.dll 등에 정의된) 커널 구조체 출력하는 방법
12091정성태12/25/2019487디버깅 기술: 147. pdb 파일을 다운로드하기 위한 symchk.exe 실행에 필요한 최소 파일
12090정성태12/24/2019353.NET Framework: 871. .NET AnyCPU로 빌드된 PE 헤더의 로딩 전/후 차이점
12089정성태12/23/2019260디버깅 기술: 146. gflags와 _CrtIsMemoryBlock을 이용한 Heap 메모리 손상 여부 체크
12088정성태12/23/2019206Linux: 28. Linux - 윈도우의 "Run as different user" 기능을 shell에서 실행하는 방법
12087정성태12/21/2019303디버깅 기술: 145. windbg/sos - Dictionary의 entries 배열 내용을 모두 덤프하는 방법 (do_hashtable.py)
12086정성태12/20/2019367디버깅 기술: 144. windbg - Marshal.FreeHGlobal에서 발생한 덤프 분석 사례
12085정성태12/20/2019304오류 유형: 586. iisreset - The data is invalid. (2147942413, 8007000d) 오류 발생 - 두 번째 이야기 [1]
12084정성태12/21/2019346디버깅 기술: 143. windbg/sos - Hashtable의 buckets 배열 내용을 모두 덤프하는 방법 (do_hashtable.py)
12083정성태12/17/2019528Linux: 27. linux - lldb를 이용한 .NET Core 응용 프로그램의 메모리 덤프 분석 방법 [1]
1  2  3  4  [5]  6  7  8  9  10  11  12  13  14  15  ...