loglin(table, margin, start = rep(1, length(table)), fit = FALSE, eps = 0.1, iter = 20, param = FALSE, print = TRUE)
table
|
a contingency table to be fit, typically the output from
table .
|
margin
|
a list of vectors with the marginal totals to be fit.
(Hierarchical) log-linear models can be specified in term of these
marginal totals which give the ``maximal'' factor subsets contained
in the model. For example, in a three-factor model,
The names of factors (i.e., |
start
|
a starting estimate for the fitted table. This optional
argument is important for incomplete tables with structural zeros
in table which should be preserved in the fit. In this
case, the corresponding entries in start should be zero and
the others can be taken as one.
|
fit
| a logical indicating whether the fitted values should be returned. |
eps
| maximum deviation allowed between observed and fitted margins. |
iter
| maximum number of iterations. |
param
| a logical indicating whether the parameter values should be returned. |
print
|
a logical. If TRUE , the number of iterations and
the final deviation are printed.
|
loglin
is used to fit log-linear models to multidimensional
contingency tables by Iterative Proportional Fitting.iter
iterations are performed, convergence is taken to occur when the
maximum deviation between observed and fitted margins is less than
eps
. All internal computations are done in double precision;
there is no limit on the number of factors (the dimension of the
table) in the model.
Assuming that there are no structural zeros, both the Likelihood
Ratio Test and Pearson test statistics have an asymptotic chisquare
distribution with df
degrees of freedom.
Package `MASS' contains loglm
, a front-end to loglin
which allows the log-linear model to be specified and fitted in a
formula-based manner similar to that of other fitting functions such
as lm
or glm
.
lrt
| the Likelihood Ratio Test statistic. |
pearson
| the Pearson test statistic (X-squared). |
df
| the degrees of freedom for the fitted model. There is no adjustment for structural zeros. |
margin
|
list of the margins that were fit. Basically the same
as the input margin , but with numbers replaced by names
where possible.
|
fit
|
An array like table containing the fitted values.
Only returned if fit is TRUE .
|
param
|
A list containing the estimated parameters of the
model. The ``standard'' constraints of zero marginal sums
(e.g., zero row and column sums for a two factor parameter) are
employed. Only returned if param is TRUE .
|
Alan Agresti (1990). Categorical data analysis. New York: Wiley.
table
## Currently no appropriate data sets are available.