R/suff.stat.R

In gIPFrm: Generalized Iterative Proportional Fitting for Relational Models

suff.stat <-
function(ModelMatrix, Table)
{
    nr <- nrow(ModelMatrix)
    S <- rep(0,nr)  
    for( h in 1:nr)
    {
       S[h] <- (Table) %*% ModelMatrix[h,]  
    }
    return(S)
}