R/summary.hmsc.R

In HMSC: Hierarchical modelling of species community

summary.hmsc<-function(object,...){
#### F. Guillaume Blanchet - June 2013
##########################################################################################
	nsp<-dim(object$model$paramX)[1]
	nsites<-nrow(object$data$X)
	
	modelSummary<-array(dim=c(nsites,nsp,5))
	dimnames(modelSummary)[[1]]<-rownames(object$data$Y)
	dimnames(modelSummary)[[2]]<-colnames(object$data$Y)
	dimnames(modelSummary)[[3]]<-c("means","sd","0.025","0.5","0.975")
	
	modelSummary[,,1:2]<-basicStatHMSC(object)
	modelSummary[,,3:5]<-quantile(object)
	
	modelSummary
}