R/entropy.hand.R

In metno/wgen: Stochastic Weather Generator

#' entropy.hand
#'
#' H(X,Y)
#' @param X
#' @param Y
#' @keywords entropy
#' @export
#' @examples
#' entropy.hand()
entropy.hand<-function(X,Y) {

   tmp<-entropy.Xsumup(X)
   nb_classX<-tmp$nb_classX
   val_classX<-tmp$val_classX
   
   tmp<-entropy.Xsumup(Y)
   nb_classY<-tmp$nb_classX
   val_classY<-tmp$val_classX
   
   hand_tmp<-matrix(0,nr=nb_classX,nc=nb_classY) 
   Pij<-entropy.Pij(X,Y)
   
   for (i in 1:nb_classX) {
      for (j in 1:nb_classY) {
	     if (Pij[i,j]!=0) {
            hand_tmp[i,j]<--(Pij[i,j])*log(Pij[i,j])/log(2)	
		} else hand_tmp[i,j]<-0 								  
      }
   }	
   hand<-sum(hand_tmp,na.rm=TRUE)
   return(hand)
}