jb-leger / blockmodels Goto Github PK

    B <- matrix(0, nrow(adj), ncol(adj))
    for (k in 1:length(covariates)) {
        B <- B + model_parameters[[Q]]$beta[k] * covariates[[k]]
    }
    if (membership_name == "LBM") {
        return(sigmo(memberships[[Q]]$Z1 %*% model_parameters[[Q]]$m %*% 
            t(memberships[[Q]]$Z2) + B))
    }
    else {
        return(sigmo(memberships[[Q]]$Z %*% model_parameters[[Q]]$m %*% 
            t(memberships[[Q]]$Z) + B))