\Header{fanny}{Fuzzy Analysis Clustering}
\keyword{cluster}{fanny}
\begin{Description}\relax
Computes a fuzzy clustering of the data into \code{k} clusters.
\end{Description}
\begin{Usage}
\begin{verbatim}
fanny(x, k, diss = inherits(x, "dist"), metric = "euclidean", stand = FALSE)
\end{verbatim}
\end{Usage}
\begin{Arguments}
\begin{ldescription}
\item[\code{x}] data matrix or data frame, or dissimilarity matrix, depending on the
value of the \code{diss} argument.

In case of a matrix or data frame, each row corresponds to an observation,
and each column corresponds to a variable. All variables must be numeric.
Missing values (NAs) are allowed.

In case of a dissimilarity matrix, \code{x} is typically the output
of \code{\Link{daisy}} or \code{\Link{dist}}.  Also a vector of
length n*(n-1)/2 is allowed (where n is the number of observations),
and will be interpreted in the same way as the output of the
above-mentioned functions.  Missing values (NAs) are not allowed.

\item[\code{k}] integer giving the desired number of clusters.  It is
required that \eqn{0 < k < n/2}{} where \eqn{n}{} is the number of
observations.
\item[\code{diss}] logical flag: if TRUE (default for \code{dist} or
\code{dissimilarity} objects), then \code{x} is assumed to be a
dissimilarity matrix.  If FALSE, then \code{x} is treated as
a matrix of observations by variables.

\item[\code{metric}] character string specifying the metric to be used for calculating
dissimilarities between observations.
The currently available options are "euclidean" and "manhattan".
Euclidean distances are root sum-of-squares of differences, and
manhattan distances are the sum of absolute differences.
If \code{x} is already a dissimilarity matrix, then this argument will
be ignored.

\item[\code{stand}] logical; if true, the measurements in \code{x} are
standardized before calculating the dissimilarities.  Measurements
are standardized for each variable (column), by subtracting the
variable's mean value and dividing by the variable's mean absolute
deviation.  If \code{x} is already a dissimilarity matrix, then this
argument will be ignored.
\end{ldescription}
\end{Arguments}
\begin{Details}\relax
In a fuzzy clustering, each observation is ``spread out'' over the various
clusters. Denote by u(i,v) the membership of observation i to cluster v.
The memberships are nonnegative, and for a fixed observation i they sum to 1.
The particular method \code{fanny} stems from chapter 4 of
Kaufman and Rousseeuw (1990).\\
Compared to other fuzzy clustering methods, \code{fanny} has the following
features: (a) it also accepts a dissimilarity matrix; (b) it is
more robust to the \code{spherical cluster} assumption; (c) it provides
a novel graphical display, the silhouette plot (see
\code{\Link{plot.partition}}).

Fanny aims to minimize the objective function
\deqn{\sum_{v=1}^k
\frac{\sum_{i=1}^n\sum_{j=1}^n u_{iv}^2 u_{jv}^2 d(i,j)}{
2 \sum_{j=1}^n u_{jv}^2}}{SUM_[v=1..k] (SUM_(i,j) u(i,v)^2 u(j,v)^2 d(i,j)) / (2 SUM_j u(j,v)^2)}
where \eqn{n}{} is the number of observations, \eqn{k}{} is the number of
clusters and \eqn{d(i,j)}{} is the dissimilarity between observations
\eqn{i}{} and \eqn{j}{}.
\end{Details}
\begin{Value}
an object of class \code{"fanny"} representing the clustering.
See \code{\Link{fanny.object}} for details.
\end{Value}
\begin{SeeAlso}\relax
\code{\Link{agnes}} for background and references;
\code{\Link{fanny.object}}, \code{\Link{partition.object}},
\code{\Link{plot.partition}}, \code{\Link{daisy}}, \code{\Link{dist}}.
\end{SeeAlso}
\begin{Examples}
\begin{ExampleCode}
## generate 25 objects, divided into two clusters, and 3 objects lying
## between those clusters.
x <- rbind(cbind(rnorm(10, 0, 0.5), rnorm(10, 0, 0.5)),
           cbind(rnorm(15, 5, 0.5), rnorm(15, 5, 0.5)),
           cbind(rnorm( 3,3.5,0.5), rnorm( 3,3.5,0.5)))
fannyx <- fanny(x, 2)
fannyx
summary(fannyx)
plot(fannyx)

data(ruspini)
## Plot similar to Figure 6 in Stryuf et al (1996)
plot(fanny(ruspini, 5))
\end{ExampleCode}
\end{Examples}

