메뉴 건너뛰기

Sub Promotion
자유게시판

Performance Diagram in R

오석근 2018.10.02 00:43 조회 수 : 3481

아래 샘플은 Performance Diagram을 R로 표출하는 방법입니다.
아래 샘플은, 홍수지역 탐지에 대한 검증때 이용한 것입니다.
안개 검증시 활용하면 좋을 것 같네요

특히, 한가지 사례에대해 검증하더라고 부스트랩기법을 이용하여 1000번 리샘플링을 통한 불확실성 검증결과까지 제시하니
논문 해석시에 유용할 것으로 보입니다.
수정할 부분과 간단한 코멘트는 아래 샘플자료에 넣어두었습니다.
혹, 이용하다고 더 궁금한 사항이 생기면 메일 주세요 (poet1535@gmail.com)

그리고 퍼포먼스 다이어그램에 관련된 논문은 아래 링크에서 확인가능하니
논문쓸때 참고하면 됩니다!!

https://www.nssl.noaa.gov/users/brooks/public_html/feda/papers/roebber2009.pdf

=======
install.packages("verification")
library(verification)
performance.diagram(main = "2009 Red-River spring flood")
performance.diagram(main = " ")
##region1
RB1 <- matrix(c(12444, 7853, 13062, 13441), ncol = 2) ## <----- 2x2 contingency table의 값들입니다. 이곳만 수정하면 됨!!!
##region2
RB2 <- matrix(c(40207, 8590, 79903, 79092), ncol = 2) ## <----- 수정!!! 분석 영역이나, 케이스가 여러개일 경우 이부분을 추가하고 아래 에서 해당하
는 부분을 추가 및 삭제하면 됨

## at point
#region1
pts <- table.stats(RB1)
boot.pts <- table.stats.boot(RB1, R = 1000, alpha = 0.05, fudge = 0.01 )
#region2
pts2 <- table.stats(RB2)
boot.pts2 <- table.stats.boot(RB2, R = 1000, alpha = 0.05, fudge = 0.01 )

## add confidence intervals
#region1
segments(x0=1-pts$FAR, y0=boot.pts["up","pod"],
x1=1-pts$FAR, y1=boot.pts["dw", "pod"], col=2)

segments(x0=1-boot.pts["up","far"], y0=pts$POD,
x1=1-boot.pts["dw","far"], y1=pts$POD, col=2)
points(1 - pts$FAR, pts$POD, col = 2, cex = 3)

#region2
segments(x0=1-pts2$FAR, y0=boot.pts2["up","pod"],
x1=1-pts2$FAR, y1=boot.pts2["dw", "pod"], col=4)

segments(x0=1-boot.pts2["up","far"], y0=pts2$POD,
x1=1-boot.pts2["dw","far"], y1=pts2$POD, col=4)
points(1 - pts2$FAR, pts2$POD, col = 4, cex = 3)
=====
번호 제목 글쓴이 날짜 조회 수
» Performance Diagram in R 오석근 2018.10.02 3481

Meteorological Satellite Research Lab.
Department of Atmospheric Science Konju National University
56, Gongjudaehak-ro, Gongju-si, Chungcheongnam-do, Korea
TEL : +82-41-850-8533 l E-mail : sms4@kongju.ac.kr

© k2s0o1d4e0s2i1g5n. All Rights Reserved