Nothing
#' Gibbs update for multiplicative effects covariance
#'
#' @usage rSuv_fc(U,V, Suv0=NULL,kappa0=NULL)
#' @param U matrix of row random effects
#' @param V matrix of row random effects
#' @param Suv0 prior (inverse) scale matrix for the prior distribution
#' @param kappa0 prior degrees of freedom for the prior distribution
#' @author Peter Hoff
#' @export rSuv_fc
#'
rSuv_fc<-function(U,V,Suv0=NULL,kappa0=NULL)
{
if(is.null(Suv0)){ Suv0<-diag(2*ncol(U)) }
if(is.null(kappa0)){ kappa0<-2+nrow(Suv0) }
UV<-cbind(U,V)
solve(rwish(solve(kappa0*Suv0+t(UV)%*%UV),nrow(U)+kappa0))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.