R/GetDC.R

In CompoundEvents: A package for compound events modeling and assessment

#'@title  Occurrence of compound dry-cold events
#'
#'@description  Take the precipitation and temperature and extract dry-cold occurrence based on thresholds.The binary variable of the dry and cold (DC) event can be otained.
#'
#'@param mp Precipitation
#'@param mt Temperature
#'@param threp Threshold of precipitation (e.g., 20th percentile)
#'@param thret Threshold of temperature (e.g., 20th percentile)
#'@usage GetDC(mp,mt,threp,thret)
#'@references Hao, Z. et al (2013). Changes in concurrent monthly precipitation and temperature extremes.Environ. Res. Lett., 8(3): 034014.
#'@return The occurrence of compound wet-hot event (0-1 binary variable)
#'@export
#'@examples
#'mp=matrix(rnorm(120,0,1),ncol=1)
#'mt=matrix(rnorm(120,0,1),ncol=1)
#'threp=20
#'thret=20
#'DC<-GetDC(mp,mt,threp,thret)

GetDC<- function(mp,mt,threp,thret)
{
  n <- length(mp)
  y <- matrix(data=0, nrow = n, ncol = 1)
  # Define the matrix
  y1 <- matrix(data=0, nrow = n, ncol = 1)
  y2 <- matrix(data=0, nrow = n, ncol = 1)

  # Spefify the Quantile (from percentile threp, threp)

  p0 <- quantile(mp,threp/100)
  t0 <- quantile(mt,thret/100)

  for (i in 1:n)
  {
    if (mp[i]<=p0)
    {
      y1[i]= 1
    }

    if (mt[i]<=t0)
    {
      y2[i]= 1
    }

  }

  y <- y1*y2

  return (y)

}