R/rb_fc.R

In eigenmodel: Semiparametric Factor and Regression Models for Symmetric Relational Data

#' Sample from the full conditional distribution of the regression coefficients
#' 
#' Sample from the full conditional distribution of the regression coefficients
#' in an eigenmodel
#' 
#' 
#' @param E a symmetric matrix
#' @return a p x 1 vector
#' @author Peter Hoff
#' @keywords multivariate models
"rb_fc" <-
function( E=Z-ULU(UL) ) {

## sample b from fc
## needs xtx, tx, pm_b and pp_b 

  if(length(pm_b)>0) {
    vr<- solve(  xtx  +  pp_b   )
    mn<- vr%*%( pp_b%*%pm_b +  tx%*%c(E[upper.tri(E)]) )
    b<-rmvnorm(mn,vr)
                      } 
  b 
                   }