R/getMI2vs1.R

In th86/gislkit: gislkit

getMI2vs1 <- function(x, y, z, bin=6, so=3, normalize=TRUE){
  n <- length(x)
  if(length(y) != n){stop("legnth of the vectors are different!")}
  if(length(z) != n){stop("legnth of the vectors are different!")}
  if(so >= bin){stop("spline order must be less than bin")}
  out <- .C("mi2vs1R", x = as.double(x), y=as.double(y), z=as.double(z), n=as.integer(n), bin=as.integer(bin), so = as.integer(so), miOut = 0, norm = as.integer(normalize))  
  mi <- out$miOut
  return (mi)
  
}