# Function that computes integrated classification likelihood based on stochastic one-mode and linked block modeling. If `clu` is a list, the method for linked/multilevel networks is applied. The support for multirelational networks is not tested. ```r ICLStochBlock( M, clu, weights = NULL, uWeights = NULL, diagonal = c("ignore", "seperate", "same"), limitType = c("none", "inside", "outside"), limits = NULL, weightClusterSize = 1, addOne = TRUE, eps = 0.001 ) ``` ## Arguments - `M`: A matrix representing the (usually valued) network. For multi-relational networks, this should be an array with the third dimension representing the relation. - `clu`: A partition. Each unique value represents one cluster. If the nework is one-mode, than this should be a vector, else a list of vectors, one for each mode. Similarly, if units are comprised of several sets, clu should be the list containing one vector for each set. - `weights`: The weights for each cell in the matrix/array. A matrix or an array with the same dimmensions as `M`. - `uWeights`: The weights for each unin. A vector with the length equal to the number of units (in all sets). - `diagonal`: How should the diagonal values be treated. Possible values are: * ignore - diagonal values are ignored * seperate - diagonal values are treated seperately * same - diagonal values are treated the same as all other values - `limitType`: Type of limit to use. Forced to 'none' if `limits` is `NULL`. Otherwise, one of either `outer` or `inner`. - `limits`: If `diagonal` is `"ignore"` or `"same"`, an array with dimensions equal to: * number of clusters (of all types) * number of clusters (of all types) * number of relations * 2 - the first is lower limit and the second is upper limit If `diagonal` is `"seperate"`, a list of two array. The first should be as described above, representing limits for off diagonal values. The second should be similar with only 3 dimensions, as one of the first two must be omitted. - `weightClusterSize`: The weight given to cluster sizes (logprobabilites) compared to ties in loglikelihood. Defaults to 1, which is "classical" stochastic blockmodeling. - `addOne`: Should one tie with the value of the tie equal to the density of the superBlock be added to each block to prevent block means equal to 0 or 1 and also "shrink" the block means toward the superBlock mean. Defaults to TRUE. - `eps`: If addOne = FALSE, the minimal deviation from 0 or 1 that the block mean/density can take. ## Returns The value of ICL ## Examples ```r # Create a synthetic network matrix set.seed(2022) library(blockmodeling) k<-2 # number of blocks to generate blockSizes<-rep(20,k) IM<-matrix(c(0.8,.4,0.2,0.8), nrow=2) clu<-rep(1:k, times=blockSizes) n<-length(clu) M<-matrix(rbinom(n*n,1,IM[clu,clu]),ncol=n, nrow=n) clu<-sample(1:2,nrow(M),replace=TRUE) plotMat(M,clu) # Have a look at this random partition ICL_pre<-ICLStochBlock(M,clu) # Calculate its ICL ICL_pre res<-stochBlock(M,clu=clu) # Optimizing the partition plot(res) # Have a look at the optimized partition ICL_post<-res$ICL # Calculate its ICL ICL_post # We expect the ICL pre-optimisation to be smaller: ICL_pre

ICLStochBlock function