minf=function(x,y){
library(sm)
air3  = cbind(x,y)

xx = sm.density(air3, eval.points = air3, display = "none",verbose=0)
result.12=diag(xx$estimate)
#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 /(result.1 * result.2)))

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

#Mutual Information
tobs
}