1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/fun.by.blocks.default.R,
% R/fun.by.blocks.opt.more.par.R, R/funByBlocks.R
\name{funByBlocks.default}
\alias{funByBlocks.default}
\alias{funByBlocks.optMorePar}
\alias{funByBlocks.opt.more.par}
\alias{funByBlocks}
\alias{fun.by.blocks}
\title{Computation of function values by blocks}
\usage{
\method{funByBlocks}{default}(
x = M,
clu,
M = x,
ignore.diag = "default",
sortNames = TRUE,
FUN = "mean",
...
)
\method{funByBlocks}{optMorePar}(x, which = 1, orderClu = FALSE, sortNames = NULL, ...)
\method{funByBlocks}{opt.more.par}(x, which = 1, orderClu = FALSE, sortNames = NULL, ...)
funByBlocks(x, ...)
fun.by.blocks(x, ...)
}
\arguments{
\item{x}{An object of suitable class or a matrix/array representing the (usually valued) network. For multi-relational networks, this should be an array with the third dimension representing the relation.
The network can have one or more modes (different kinds of units with no ties among themselves.
If the network is not two-mode, the matrix must be square.}
\item{clu}{A partition. Each unique value represents one cluster.
If the network is one-mode, then this should be a vector, else a list of vectors, one for each mode.}
\item{M}{A matrix representing the (usually valued) network. For multi-relational networks, this should be an array with the third dimension representing the relation.
The network can have one or more modes (different kinds of units with no ties among themselves.
If the network is not two-mode, the matrix must be square.}
\item{ignore.diag}{Should the diagonal be ignored.}
\item{sortNames}{Should the rows and columns of the matrix be sorted based on their names.}
\item{FUN}{The function to be computed over the blocks.}
\item{\dots}{Further arguments to \code{funByBlocks.default}.}
\item{which}{Which (if several) of the "best" solutions should be used.}
\item{orderClu}{Should the partition be ordered before computing. \code{FALSE} by default. If \code{TRUE}, \code{\link{orderClu}} is used (using default arguments) to order the clusters in a partition in "decearsing" (see \code{\link{orderClu}} for interpretation) order. If \code{TRUE}, \code{sortNames} is set to \code{FALSE}.}
}
\value{
A numerical matrix of \code{FUN} values by blocks, induced by a partition \code{clu}.
}
\description{
Computes a value of a function over blocks of a matrix, defined by a partition.
}
\examples{
n <- 8 # If larger, the number of partitions increases dramatically,
# as does if we increase the number of clusters
net <- matrix(NA, ncol = n, nrow = n)
clu <- rep(1:2, times = c(3, 5))
tclu <- table(clu)
net[clu == 1, clu == 1] <- rnorm(n = tclu[1] * tclu[1], mean = 0, sd = 1)
net[clu == 1, clu == 2] <- rnorm(n = tclu[1] * tclu[2], mean = 4, sd = 1)
net[clu == 2, clu == 1] <- rnorm(n = tclu[2] * tclu[1], mean = 0, sd = 1)
net[clu == 2, clu == 2] <- rnorm(n = tclu[2] * tclu[2], mean = 0, sd = 1)
# Optimizing 10 random partitions with optRandomParC
res <- optRandomParC(M = net, k = 2, rep = 10, approaches = "hom", homFun = "ss", blocks = "com")
plot(res) # Hopefully we get the original partition
funByBlocks(res)
# Computing mean by blocks, ignoring the diagonal (default)
}
\references{
\enc{Žiberna, A.}{Ziberna, A.} (2007). Generalized Blockmodeling of Valued Networks. Social Networks, 29(1), 105-126. doi: 10.1016/j.socnet.2006.04.002
\enc{Žiberna, A.}{Ziberna, A.} (2008). Direct and indirect approaches to blockmodeling of valued networks in terms of regular equivalence. Journal of Mathematical Sociology, 32(1), 57-84. doi: 10.1080/00222500701790207
}
\seealso{
\code{\link{optRandomParC}}, \code{\link{optParC}}
}
\author{
\enc{Aleš Žiberna}{Ales Ziberna}
}
\keyword{cluster}
\keyword{math}
|
