Fits log-linear models for incomplete contingency tables, including some latent class models, via EM and Fisher scoring approaches.

1 | ```
gllm(y,s,X,method="hybrid",em.maxit=1,tol=0.00001)
``` |

`y` |
is the observed contingency table. |

`s` |
is a vector of indices, one for each cell of the full (unobserved)
contingency table, representing the appropriate cell of |

`X` |
is the design matrix, or a formula. |

`method` |
chooses the EM, Fisher scoring or a hybrid (EM then scoring) method for fitting the model. |

`em.maxit` |
is the number of EM iterations. |

`tol` |
is the convergence criterion for the LR criterion. |

The generalized log-linear model allows for modelling of incomplete contingency tables, that is tables where one or more dimensions have been collapsed over. These include situations where imprecise measures have been calibrated using a "perfect" gold standard, and the true association between imperfectly measured variables is to be estimated; where data is missing for a subsample of the population; latent variable models where latent variables are "errorless" functions of observed variables - eg ML gene frequency estimation from counts of observed phenotypes; specialised measurement models eg where observed counts are mixtures due to perfect measures and error prone measures; standard latent class analysis; symmetry and quasi-symmetry models for square tables.

The general framework underlying these models is summarised by Espeland
(1986), and Espeland & Hui (1987), and is originally due to Thompson &
Baker (1981). An observed contingency table *y*, which will be
treated as a vector, is modelled as arising from an underlying complete
table *z*, where observed count *y_j* is the sum of a number of
elements of *z*, such that each *z_i* contributes to no more
than one *y_j*. Therefore one can write *y=F'z*, where *F*
is made up of orthogonal columns of ones and zeros.

We then specify a loglinear model for *z*, so that
*log(E(z))=X'b*, where *X* is a design matrix, and *b* a
vector of loglinear parameters. The loglinear model for *z* and
thus *y*, can be fitted using two methods, both of which are
available in `gllm`

. The first was presented as AS207 by Michael
Haber (1984) and combines an iterative proportional fitting algorithm
for *b* and *z*, with an EM fitting for *y*, *z* and
*b*. The second is a Fisher scoring approach, presented in Espeland
(1986).

The `gllm`

function is actually a simple wrapper for `scoregllm()`

.

A list with components:

`iter` |
the number of scoring iterations until convergence |

`deviance` |
the final model deviance (-2 log likelihood) |

`df` |
the model degrees of freedom |

`coefficients` |
the model parameter estimates |

`se` |
the standard errors for the model parameter estimates |

`V` |
the variance-covariance matrix for the model parameter estimates |

`observed.values` |
the observed counts in |

`fitted.values` |
the expected counts under the fitted model |

`residuals` |
Pearsonian residuals under the fitted model |

`full.table` |
the expected counts for the full (unobserved) table. |

Espeland MA (1986). A general class of models for discrete multivariate
data. *Commun. Statist.-Simula* 15:405-424.

Espeland MA, Hui SL (1987). A general approach to analyzing epidemiologic
data that contains misclassification errors. *Biometrics*
43:1001-1012.

Haber M (1984). AS207: Fitting a general log-linear model. *Appl
Statist* 33:358-362.

Thompson R, Baker RJ (1981). Composite link functions in generalized
linear models. *Appl Statist* 30: 125-131.

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 | ```
#
# latent class analysis: two latent classes
#
# Data matrix 2x2x2x2x2 table of responses to five binary items
# (items 11-15 of sections 6-7 of the Law School Admission Test)
#
y<-c( 3, 6, 2, 11, 1, 1, 3, 4,
1, 8, 0, 16, 0, 3, 2, 15,
10, 29, 14, 81, 3, 28, 15, 80,
16, 56, 21, 173, 11, 61, 28, 298)
#
# Scatter matrix: full table is 2x2x2x2x2x2
#
s<- c(1:32,1:32)
#
# Design matrix: x is the latent variable (2 levels),
# a-e are the observed variables
#
x<-as.integer(gl(2,32,64))-1
a<-as.integer(gl(2,16,64))-1
b<-as.integer(gl(2,8 ,64))-1
c<-as.integer(gl(2,4 ,64))-1
d<-as.integer(gl(2,2 ,64))-1
e<-as.integer(gl(2,1 ,64))-1
res1<-gllm(y,s,~x*(a+b+c+d+e),method="em",tol=0.01)
res1
#
# An example of model fitting: gametic association between two diallelic loci
#
# Data matrix
#
y<-c( 187,386,156,
352,310,20,
136,0 ,0)
#
# Scatter matrix
#
s<- c( 1, 2, 2, 3,
4, 5, 5, 6,
4, 5, 5, 6,
7, 8, 8, 9)
#
# Design matrix
#
X<- matrix(c( 1,0,0,0,0,0,1,
1,0,1,0,0,0,0,
1,0,1,0,0,0,0,
1,0,2,0,1,0,0,
1,1,0,0,0,0,0,
1,1,1,0,0,1,0,
1,1,1,0,0,0,1,
1,1,2,0,1,1,1,
1,1,0,0,0,0,0,
1,1,1,0,0,0,1,
1,1,1,0,0,1,0,
1,1,2,0,1,1,1,
1,2,0,1,0,0,0,
1,2,1,1,0,1,1,
1,2,1,1,0,1,1,
1,2,2,1,1,2,2), byrow=TRUE, ncol=7)
colnames(X)<-c("Intercept", "A", "B", "P1", "P2", "Delta", "Epsilon")
res2<-gllm(y,s,X[,c(1:6)],method="hybrid",em.maxit=1,tol=0.00001)
res2
#
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.