R/posterior.inverse.R

In MCMCglmm: MCMC Generalised Linear Mixed Models

"posterior.inverse"<-function(x){

   if(is.null(dim(x))[1]){
     x<-as.matrix(x)
   }
   if(sqrt(dim(x)[2])%%1!=0){
    stop("cannot coerce rows of x into a matrix - the number of columns should be a triangular number")
   }
   coerce.inverse<-function(x){solve(matrix(x, sqrt(length(x)), sqrt(length(x))))}

   if(is.mcmc(x)==FALSE){
     warning("posterior.inverse expecting mcmc object")
     t(apply(x, 1, coerce.inverse))
   }else{   
     as.mcmc(t(apply(x, 1, coerce.inverse)))
   }  
}