minf=function(x,y){
library(sm)
air3  = cbind(x,y)
result.12 = sm.density(air3,   eval.points = air3,   display = "none")
result.12 = diag(result.12$estimate)

result.1  = sm.density(x,  eval.points = x,  display = "none")$estimate
result.2  = sm.density(y, eval.points = y, display = "none")$estimate

#Mutual Information from the hand out..
#tobs      = mean(log(result.12$estimate /(result.1$estimate * result.2$estimate)))

#Average Mutual Information .. Moon et al. (1995)
tobs      = sum(result.12 * log2(result.12/(result.1 * result.2)))

#Mutual Information
tobs
}