Nothing
#' Simulate a and Sab from full conditional distributions under bin likelihood
#'
#' Simulate a and Sab from full conditional distributions under bin likelihood
#'
#'
#' @usage raSab_bin_fc(Z, Y, a, b, Sab, Sab0=NULL, eta0=NULL, SS = round(sqrt(nrow(Z))))
#' @param Z a square matrix, the current value of Z
#' @param Y square binary relational matrix
#' @param a current value of row effects
#' @param b current value of column effects
#' @param Sab current value of Cov(a,b)
#' @param Sab0 prior (inverse) scale matrix for the prior distribution
#' @param eta0 prior degrees of freedom for the prior distribution
#' @param SS number of iterations
#' @return \item{Z}{new value of Z} \item{Sab}{new value of Sab} \item{a}{new
#' value of a}
#' @author Peter Hoff
#' @export raSab_bin_fc
raSab_bin_fc <-
function(Z,Y,a,b,Sab,Sab0=NULL,eta0=NULL,SS=round(sqrt(nrow(Z))))
{
if(is.null(Sab0)){ Sab0<-diag(2) }
if(is.null(eta0)){ eta0<-4 }
E<-Z-a%*%t(rep(1,nrow(Z)))
MEL<-MEU<- -E
MEL[!is.na(Y) & Y==0]<- -Inf
MEU[!is.na(Y) & Y==1]<-Inf
MEL[is.na(Y)]<- -Inf ; MEU[is.na(Y)]<- Inf
# diag(MEU)<-Inf;diag(MEL)<- -Inf
lba<-apply(MEL,1,max)
lba[is.na(lba)]<- -Inf
uba<-apply(MEU,1,min)
uba[is.na(uba)]<- Inf
for(ss in 1:SS)
{
ea<-b*Sab[1,2]/Sab[2,2]
sa<-sqrt(Sab[1,1]-Sab[1,2]^2/Sab[2,2])
a<-ea+sa*qnorm(runif(nrow(Z),pnorm((lba-ea)/sa),pnorm((uba-ea)/sa)))
Sab<-solve(rwish(solve(eta0*Sab0+crossprod(cbind(a,b))),eta0+nrow(Z)))
}
list(Z=E+a%*%t(rep(1,nrow(Z))),a=a,Sab=Sab)
}
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.