omega: This function finds the bivariate joint probability or the binary correlation from the corresponding Gaussian correlation `x`

Description

This function finds the bivariate joint probability or the binary correlation from the corresponding Gaussian correlation x

Usage

omega(x = 0.5, p0_v1 = 0.5, p0_v2 = NA, correlation = FALSE)

Arguments

value of expected correlation between the corresponding Gaussian-distributed variables

p0_v1, p0_v2

probability of no precipitation occurrences for the v1 and v2 time series respectively. See Notes.

correlation

logical numeric value. Default is FALSE. If TRUE the function returns the binary correlation like eq. 6 of Mhanna, et al.,2011.

Value

probability of no precipitation occurrence in both v1 and v2 simultaneously. It is a matrix if x is a matrix.

References

D.S. Wilks (1998), Multisite Generalization of a Daily Stochastic Precipitation Generation Model, Journal of Hydrology, Volume 210, Issues 1-4, September 1998, Pages 178-191, https://www.sciencedirect.com/science/article/pii/S0022169498001863

Muamaraldin Mhanna and Willy Bauwens (2011) A Stochastic Space-Time Model for the Generation of Daily Rainfall in the Gaza Strip, International Journal of Climatology, Volume 32, Issue 7, pages 1098-1112, 10.1002/joc.2305, https://rmets.onlinelibrary.wiley.com/doi/abs/10.1002/joc.2305

Examples

Run this code

# NOT RUN {
rho <- 0.4
p00 <- omega(x=rho,p0_v1=0.5,p0_v2=0.5)
cor00 <- omega(x=rho,p0_v1=0.5,p0_v2=0.5,correlation=TRUE)
# }