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
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/nkpar.R, R/nkpartitions.R
\name{nkpar}
\alias{nkpar}
\alias{nkpartitions}
\title{Functions for listing all possible partitions or just counting the number of them}
\usage{
nkpar(n, k)
nkpartitions(n, k, exact = TRUE, print = FALSE)
}
\arguments{
\item{n}{Number of units/objects.}
\item{k}{Number of clusters/groups.}
\item{exact}{Search for partitions with exactly \code{k} or at most \code{k} clusters.}
\item{print}{Print results as they are found.}
}
\value{
The matrix or number of possible partitions.
}
\description{
The function \code{nkpartitions} lists all possible partitions of n objects in to k clusters.
}
\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)
# Computation of criterion function with the correct partition
nkpar(n = n, k = length(tclu)) # Computing the number of partitions
all.par <- nkpartitions(n = n, k = length(tclu)) # Forming the partitions
all.par <- lapply(apply(all.par, 1, list), function(x) x[[1]])
# to make a list out of the matrix
res <- critFunC(M = net, clu = clu, approaches = "val",
blocks = c("nul", "com"), preSpecM = 4)
plot(res) # We get the original partition
}
\author{
Chris Andrews
}
\keyword{cluster}
|
