R/computePredictor.HmscRandomLevel.R

In hmsc-r/HMSC: Hierarchical Model of Species Communities

computePredictor.HmscRandomLevel = function(rL,Eta,Lambda,pi=c(1:nrow(Eta)),dfpi=NULL){
   if(rL$xDim == 0){
      LRan = Eta[pi,]%*%Lambda
   } else{
      LRan = matrix(0,length(pi),ncol(Lambda))
      for(k in 1:rL$xDim)
         LRan = LRan + (Eta[pi,]*rL$x[as.character(dfpi),k]) %*% Lambda[,,k]

      # np = length(pi)
      # nf = ncol(Lambda)
      # EtaArray = array(Eta[pi,], c(np,nf,rL$xDim))
      # XArray = array(rL$x[rep(as.character(dfpi),nf),], c(np,nf,rL$xDim))
   }
   return(LRan)
}