pdfskill: Calculates the overlap between two PDFs.

In sethmcg/climod: Bias correction and other tools for climate model output

Description

Given two datasets (observations and model), this function calculates a skilll score as the area of overlap between their estimated PDFs. This metric comes from Perkins, et al. 2007 in the Journal of Climate (doi: 10.1175/JCLI4253.1).

Usage

pdfskill(obs, mod, ...)

Arguments

obs	A vector of observed data.
mod	A vector of model data.
...	Arguments to the density estimation function.

Details

This implementation uses binned kernel density estimation from the KernSmooth package to estimate the PDFs.

Value

A value between 0 and 1.

Examples

obs <- rnorm(1000, mean=0, sd=1)
mod <- rnorm(1000, mean=0.3, sd=0.7)
pdfskill(obs, mod)

sethmcg/climod documentation built on Nov. 19, 2021, 11:12 p.m.