ppsbm (version 1.0.0)

sortIntensities: Sort intensities

Description

Sort intensities associated with the estimated clustering \(\hat z\) "in the same way" as the original intensities associated with true clustering \(z\) by permutation of rows.

Usage

sortIntensities(intensities, z, hat.z, directed)

Arguments

intensities

Matrix whose rows contain piecewise constant intensities \(\alpha^{(q,l)}\), given as a vector of values on a regular partition of the time interval.

z

Matrix of size \(Q \times n\).

hat.z

Matrix of size \(Q \times n\).

directed

Boolean for directed (TRUE) or undirected (FALSE) case.

References

HUBERT, L. & ARABIE, P. (1985). Comparing partitions. J. Classif. 2, 193–218.

MATIAS, C., REBAFKA, T. & VILLERS, F. (2018). A semiparametric extension of the stochastic block model for longitudinal networks. Biometrika. 105(3): 665-680.

Examples

# True and estimated clusters for n=6 nodes clustered into Q=3 groups
z <- matrix(c(1,1,0,0,0,0, 0,0,1,1,0,0, 0,0,0,0,1,1), nrow = 3, byrow = TRUE)
hat.z <- matrix(c(0,0,1,1,0,0, 1,1,0,0,0,0, 0,0,0,0,1,1), nrow = 3, byrow = TRUE)

# Set constant intensities for each directed pair (q,l)
intens <- matrix(c(1,1,1,2,2,2,3,3,3),9)

# Permute the rows according to the permutation that "matches" hat.z with z
sortIntensities(intens,z,hat.z, TRUE)