| mygllm | R Documentation |
Fits a log-linear model for collapsed contingency tables.
mygllm(y, s, X, maxit = 1000, tol = 1e-05, E = rep(1, length(s)))
y |
Vector of observed cell frequencies. |
s |
Scatter matrix. s[i] is the cell in the observed array that corresponds to cell i in the full array. |
X |
Design matrix. |
maxit |
Maximum number of iterations. |
tol |
Convergence parameter. |
E |
Full contingency table. Should be initialized with either ones or a priori estimates. |
This is an implementation and extension of the algorithm published by Haber (1984). It also incorporates ideas of David Duffy (see references).
A priori estimates of the full contingency table can be given as
start values by argument E. This can reduce
execution time significantly.
Estimated full contingency table.
Andreas Borg, Murat Sariyar
Michael Haber, Algorithm AS 207: Fitting a General Log-Linear Model, in: Applied Statistics 33 (1984) No. 3, 358–362.
David Duffy: gllm: Generalised log-linear model. R package version 0.31. https://cran.r-project.org/package=gllm
emWeights, which makes use of log-linear fitting for
weight calculation.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.