R/score.HIT.R

In metno/wgen: Stochastic Weather Generator

#' score.HIT
#'
#' Hit
#' @param data1
#' @param data2
#' @param ths1 threshold#1
#' @param ths2 threshold#2
#' @keywords score
#' @export
#' @examples
#' score.HIT()
score.HIT<-function(data1, #at least one inside the area
                          data2,
						  ths1,ths2) {
   if (length(ths1)==2) temp<-((data1>ths1[1]) & (data1<=ths1[2]))
   if (length(ths1)==1) temp<-(data1>ths1)
   
   if (length(ths2)==2) temp2<-((data2>ths2[1]) & (data2<=ths2[2]))
   if (length(ths2)==1) temp2<-(data2>ths2)
         
   BOOL<-Reduce(`|`, as.data.frame(t(temp)))
   results<-BOOL & temp2
   return(results)
}