R/mahasuhab.r

In ClementCalenge/adehabitat: Analysis of Habitat Selection by Animals

"mahasuhab" <- function(kasc, pts, type=c("distance", "probability"))
  {
      ## Verifications
      x<-pts
      type<-match.arg(type)
      if (!inherits(kasc, "kasc"))
          stop("should be an object of class \"kasc\"")
      if (ncol(x)!=2)
          stop("x should have 2 columns")
      kasc<-managNAkasc(kasc)

      ## Computation of the variance-covariance matrix of the used points:
      hihi<-join.kasc(x, kasc)
      used<-list()

      ## factors are transformed into dummy variables
      for (i in 1:ncol(hihi)) {
          if (is.factor(hihi[,i]))
              used[[i]]<-acm.disjonctif(data.frame(hihi[,i]))[,-1]
          else
              used[[i]]<-hihi[,i]
      }
      used[[i+1]]<-rep(1, nrow(hihi))
      hihi<-as.data.frame(used)
      hihi<-hihi[!is.na(hihi[,1]),]
      mu<-apply(hihi,2, function(x) mean(x, na.rm=TRUE))
      varcov<-t(as.matrix(hihi))%*%as.matrix(hihi)/nrow(hihi)


      ## habitat Availability
      ava<-list()

      ## factors are transformed into dummy variables
      for (i in 1:ncol(kasc)) {
          if (is.factor(kasc[,i]))
              ava[[i]]<-acm.disjonctif(data.frame(kasc[,i]))[,-1]
          else
              ava[[i]]<-kasc[,i]
      }

      ava[[i+1]]<-rep(1, nrow(kasc))
      df<-as.data.frame(ava)

      ## computation of the Mahalanobis distances
      map<-mahalanobis(as.matrix(df), mu, varcov)
      if (type=="probability")
          map<-1-pchisq(map, ncol(hihi)-1)
      mat<-matrix(map, attr(kasc,"ncol"), attr(kasc,"nrow"))
      mat<-getascattr(getkasc(kasc, names(kasc)[1]), mat)

      return(mat)
  }