\HeaderA{interaction}{Compute Factor Interactions}{interaction}
\keyword{category}{interaction}
\begin{Description}\relax
\code{interaction} computes a factor which represents the interaction
of the given factors.  The result of \code{interaction} is always unordered.
\end{Description}
\begin{Usage}
\begin{verbatim}
interaction(..., drop = FALSE, sep = ".", lex.order = FALSE)
\end{verbatim}
\end{Usage}
\begin{Arguments}
\begin{ldescription}
\item[\code{...}] the factors for which interaction is to be computed, or a
single list giving those factors.
\item[\code{drop}] if \code{drop} is \code{TRUE}, unused factor levels
are dropped from the result.  The default is to retain all
factor levels.
\item[\code{sep}] string to construct the new level labels by joining the
constituent ones.
\item[\code{lex.order}] logical indicating if the order of factor concatenation
should be lexically ordered.
\end{ldescription}
\end{Arguments}
\begin{Value}
A factor which represents the interaction of the given factors.
The levels are labelled as the levels of the individual factors joined
by \code{sep} which is \code{.} by default.

By default, when \code{lex.order = FALSE}, the levels are ordered so
the level of the first factor varies fastest, then the second and so
on.  This is the reverse of lexicographic ordering (which you can get
by  \code{lex.order = TRUE}), and differs from
\code{\LinkA{:}{:}}.  (It is done this way for compatibility with S.)
\end{Value}
\begin{References}\relax
Chambers, J. M. and Hastie, T. J. (1992)
\emph{Statistical Models in S}.
Wadsworth \& Brooks/Cole.
\end{References}
\begin{SeeAlso}\relax
\code{\LinkA{factor}{factor}};
\code{\LinkA{:}{:}} where \code{f:g} is similar to
\code{interaction(f, g, sep=":")} when \code{f} and  \code{g} are factors.
\end{SeeAlso}
\begin{Examples}
\begin{ExampleCode}
a <- gl(2, 4, 8)
b <- gl(2, 2, 8, labels = c("ctrl", "treat"))
s <- gl(2, 1, 8, labels = c("M", "F"))
interaction(a, b)
interaction(a, b, s, sep = ":")
stopifnot(identical(a:s,
                    interaction(a, s, sep = ":", lex.order = TRUE)),
          identical(a:s:b,
                    interaction(a, s, b, sep = ":", lex.order = TRUE)))
\end{ExampleCode}
\end{Examples}

