Microsoft MVP성태의 닷넷 이야기
Math: 51. MathNET + OxyPlot을 이용한 간단한 통계 정보 처리 - Histogram [링크 복사], [링크+제목 복사],
조회: 21204
글쓴 사람
정성태 (techsharer at outlook.com)
홈페이지
첨부 파일
(연관된 글이 2개 있습니다.)

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 at outlook.com

비밀번호

댓글 작성자
 




1  2  3  [4]  5  6  7  8  9  10  11  12  13  14  15  ...
NoWriterDateCnt.TitleFile(s)
13879정성태2/5/20255295오류 유형: 946. Ubuntu - N: Updating from such a repository can't be done securely, and is therefore disabled by default.
13878정성태2/3/20255248오류 유형: 945. Windows - 최대 절전 모드 시 DRIVER_POWER_STATE_FAILURE 발생 (pacer.sys)
13877정성태1/25/20254482닷넷: 2315. C# - PCI 장치 열거 (레지스트리, SetupAPI)파일 다운로드1
13876정성태1/25/20255522닷넷: 2314. C# - ProcessStartInfo 타입의 Arguments와 ArgumentList파일 다운로드1
13875정성태1/24/20254108스크립트: 69. 파이썬 - multiprocessing 패키지의 spawn 모드로 동작하는 uvicorn의 workers
13874정성태1/24/20255321스크립트: 68. 파이썬 - multiprocessing Pool의 기본 프로세스 시작 모드(spawn, fork)
13873정성태1/23/20253944디버깅 기술: 217. WinDbg - PCI 장치 열거파일 다운로드1
13872정성태1/23/20253557오류 유형: 944. WinDbg - 원격 커널 디버깅이 연결은 되지만 Break (Ctrl + Break) 키를 눌러도 멈추지 않는 현상
13871정성태1/22/20254051Windows: 278. Windows - 윈도우를 다른 모니터 화면으로 이동시키는 단축키 (Window + Shift + 화살표)
13870정성태1/18/20254924개발 환경 구성: 741. WinDbg - 네트워크 커널 디버깅이 가능한 NIC 카드 지원 확대
13869정성태1/18/20254257개발 환경 구성: 740. WinDbg - _NT_SYMBOL_PATH 환경 변수에 설정한 경로로 심벌 파일을 다운로드하지 않는 경우
13868정성태1/17/20253893Windows: 277. Hyper-V - Windows 11 VM의 Enhanced Session 모드로 로그인을 할 수 없는 문제
13867정성태1/17/20255395오류 유형: 943. Hyper-V에 Windows 11 설치 시 "This PC doesn't currently meet Windows 11 system requirements" 오류
13866정성태1/16/20255613개발 환경 구성: 739. Windows 10부터 바뀐 device driver 서명 방법
13865정성태1/15/20255169오류 유형: 942. C# - .NET Framework 4.5.2 이하의 버전에서 HttpWebRequest로 https 호출 시 "System.Net.WebException" 예외 발생
13864정성태1/15/20255163Linux: 114. eBPF를 위해 필요한 SELinux 보안 정책
13863정성태1/14/20254078Linux: 113. Linux - 프로세스를 위한 전용 SELinux 보안 문맥 지정
13862정성태1/13/20254585Linux: 112. Linux - 데몬을 위한 SELinux 보안 정책 설정
13861정성태1/11/20254694Windows: 276. 명령행에서 원격 서비스를 동기/비동기로 시작/중지
13860정성태1/10/20254357디버깅 기술: 216. WinDbg - 2가지 유형의 식 평가 방법(MASM, C++)
13859정성태1/9/20255135디버깅 기술: 215. Windbg - syscall 이후 실행되는 KiSystemCall64 함수 및 SSDT 디버깅
13858정성태1/8/20255124개발 환경 구성: 738. PowerShell - 원격 호출 시 "powershell.exe"가 아닌 "pwsh.exe" 환경으로 명령어를 실행하는 방법
13857정성태1/7/20255410C/C++: 187. Golang - 콘솔 응용 프로그램을 Linux 데몬 서비스를 지원하도록 변경파일 다운로드1
13856정성태1/6/20254317디버깅 기술: 214. Windbg - syscall 단계까지의 Win32 API 호출 (예: Sleep)
13855정성태12/28/20245963오류 유형: 941. Golang - os.StartProcess() 사용 시 오류 정리
13854정성태12/27/20245854C/C++: 186. Golang - 콘솔 응용 프로그램을 NT 서비스를 지원하도록 변경파일 다운로드1
1  2  3  [4]  5  6  7  8  9  10  11  12  13  14  15  ...