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
|
\name{Crossings}
\alias{Crossings}
%- Also NEED an '\alias' for EACH other topic documented here.
\title{
Crossings Interaction of Factors
}
\description{
Given two ordered factors in a square, n x n frequency table, \code{Crossings}
creates an n-1 column matrix corresponding to different degrees of difficulty
in crossing from one level to the next, as described by Goodman (1972).
}
\usage{
Crossings(...)
}
%- maybe also 'usage' for other objects documented here.
\arguments{
\item{\dots}{
Two factors
}
}
%\details{
%% ~~ If necessary, more details than the description above ~~
%}
\value{
For two factors of \code{n} levels, returns a binary indicator matrix
of \code{n*n} rows and \code{n-1} columns.
}
\references{
Goodman, L. (1972). Some multiplicative models for the analysis of cross-classified data.
In: \emph{Proceedings of the Sixth Berkeley Symposium on Mathematical Statistics and Probability},
Berkeley, CA: University of California Press, pp. 649-696.
}
\author{
Michael Friendly and Heather Turner
}
%\note{
%% ~~further notes~~
%}
%% ~Make other sections like Warning with \section{Warning }{....} ~
\seealso{
\code{\link[stats]{glm}},
\code{\link[gnm]{gnm}} for model fitting functions for frequency tables
\code{\link[gnm]{Diag}},
\code{\link[gnm]{Mult}},
\code{\link[gnm]{Symm}},
\code{\link[gnm]{Topo}}
for similar extensions to terms in model formulas.
}
\examples{
data(Hauser79)
# display table
structable(~Father+Son, data=Hauser79)
hauser.indep <- gnm(Freq ~ Father + Son, data=Hauser79, family=poisson)
hauser.CR <- update(hauser.indep, ~ . + Crossings(Father,Son))
LRstats(hauser.CR)
hauser.CRdiag <- update(hauser.indep, ~ . + Crossings(Father,Son) + Diag(Father,Son))
LRstats(hauser.CRdiag)
}
% Add one or more standard keywords, see file 'KEYWORDS' in the
% R documentation directory.
\keyword{models}
\keyword{manip}% __ONLY ONE__ keyword per line
|
