est_between {bigergm} | R Documentation |
Estimate between-block parameters
Description
Function to estimate the between-block model by relying on the maximum likelihood estimator.
Usage
est_between(
formula,
network,
add_intercepts = TRUE,
clustering_with_features = FALSE
)
Arguments
formula |
An R |
network |
a network object with one vertex attribute called 'block' representing which node belongs to which block |
add_intercepts |
Boolean value to indicate whether adequate intercepts should be added to the provided formula so that the model in the first stage of the estimation is a nested model of the estimated model in the second stage of the estimation |
clustering_with_features |
Boolean value to indicate if the clustering
was carried out making use of the covariates or not (only important if |
Value
'ergm' object of the estimated model.
References
Morris M, Handcock MS, Hunter DR (2008). Specification of Exponential-Family Random Graph Models: Terms and Computational Aspects. Journal of Statistical Software, 24.
Examples
adj <- c(
c(0, 1, 0, 0, 1, 0),
c(1, 0, 1, 0, 0, 1),
c(0, 1, 0, 1, 1, 0),
c(0, 0, 1, 0, 1, 1),
c(1, 0, 1, 1, 0, 1),
c(0, 1, 0, 1, 1, 0)
)
adj <- matrix(data = adj, nrow = 6, ncol = 6)
rownames(adj) <- as.character(1001:1006)
colnames(adj) <- as.character(1001:1006)
# Use non-consecutive block names
block <- c(50, 70, 95, 50, 95, 70)
g <- network::network(adj, matrix.type = "adjacency")
g %v% "block" <- block
est <- est_between(
formula = g ~ edges,network = g,
add_intercepts = FALSE, clustering_with_features = FALSE
)