Generalized linear array model

Generalized linear array model

In statistics, the generalized linear array model(GLAM) is used for analyzing the data sets with array structure. It based on the generalized linear model with the regression matrix written as a Kronecker product.

Overview

In the article published in the Journal of the Royal Statistical Society series B, 2006, Currie, Durban and Eilers introduced the generalized linear array model or GLAM. GLAMs provide a structure and a computational procedure for fitting generalized linear models or GLMs whose model matrix can be written as a Kronecker product and whose data can be written as an array. In a large GLM, the GLAM approach gives very substantial savings in both storage and computational time over the usual GLM algorithm.

Suppose the data mathbf Y is arranged in a d-dimensional array with size n_1 imes n_2 imesldots imes n_d; thus,the corresponding data vector mathbf y = extbf{vec}(mathbf Y) has size n_1n_2n_3cdots n_d. Suppose also that the regression matrix mathbf X = mathbf X_dotimesmathbf X_{d-1}otimesldotsotimesmathbf X_1.

The standard analysis of a GLM with data vector mathbf y and regression matrix mathbf X proceeds by repeated evaluation of the scoring algorithm

mathbf X' ilde{mathbf W}_deltamathbf Xhat{oldsymbol heta} = mathbf X' ilde{mathbf W}_delta ilde{mathbf z}

where ilde{oldsymbol heta} represents the approximate solution of oldsymbol heta, and hat{oldsymbol heta} is the improved value of it; mathbf W_delta is the diagonal weight matrix with elements

w_{ii}^{-1} = left(frac{partialeta_i}{partialmu_i} ight)^2 ext{var}(y_i),

and mathbf z = oldsymboleta + mathbf W_delta^{-1}(mathbf y - oldsymbolmu) is the working variable.

Computationally, GLAM provides array algorithms to calculate the linear predictor, oldsymboleta = mathbf X oldsymbol heta and the weighted inner product mathbf X' ilde{mathbf W}_deltamathbf X without evaluation of the model matrix mathbf X .

Example: In 2 dimensions, let mathbf X = mathbf X_2otimesmathbf X_1 then the linear predictor is written mathbf X_1 oldsymbolTheta mathbf X_2' where oldsymbolTheta is the matrix of coefficients; the weighted inner product is obtained from G(mathbf X_1)' mathbf W G(mathbf X_2) and mathbf W is the matrix of weights; here G(mathbf M) is the row tensor function of the r imes c matrix mathbf M given by

G(mathbf M) = (mathbf M otimes mathbf 1') * (mathbf 1' otimes mathbf M) where * means element by element multiplcation and mathbf 1 is a vector of 1's of length c.

These low storage high speed formulae extend to d-dimensions.

Applications: GLAM is designed to be used in d-dimensional smoothing problems where the data are arranged in an array and the smoothing matrix is constructed as a Kronecker product of d one-dimensional smoothing matrices.

References

* I.D Currie, M. Durban and P. H. C. Eilers (2006) Generalized linear array models with applications to multidimensional smoothing,"Journal of Royal Statistical Society - Series B", 68, part 2, 259-280.


Wikimedia Foundation. 2010.

Игры ⚽ Поможем решить контрольную работу

Look at other dictionaries:

  • Array data type — Not to be confused with Array data structure. In computer science, an array type is a data type that is meant to describe a collection of elements (values or variables), each selected by one or more indices that can be computed at run time by the …   Wikipedia

  • Probit model — In statistics, a probit model is a popular specification of a generalized linear model, using the probit link function. A probit regression is the application of this model to a given dataset. Probit models were introduced by Chester Ittner Bliss …   Wikipedia

  • List of statistics topics — Please add any Wikipedia articles related to statistics that are not already on this list.The Related changes link in the margin of this page (below search) leads to a list of the most recent changes to the articles listed below. To see the most… …   Wikipedia

  • List of mathematics articles (G) — NOTOC G G₂ G delta space G networks Gδ set G structure G test G127 G2 manifold G2 structure Gabor atom Gabor filter Gabor transform Gabor Wigner transform Gabow s algorithm Gabriel graph Gabriel s Horn Gain graph Gain group Galerkin method… …   Wikipedia

  • Biological neuron model — A biological neuron model (also known as spiking neuron model) is a mathematical description of the properties of nerve cells, or neurons, that is designed to accurately describe and predict biological processes. This is in contrast to the… …   Wikipedia

  • Parametric array — The parametric array is a nonlinear transduction mechanism that generates narrow, nearly sidelobe free beams of low frequency sound, through the mixing and interaction of high frequency sound waves, effectively overcoming the diffraction limit (a …   Wikipedia

  • Effect size — In statistics, an effect size is a measure of the strength of the relationship between two variables in a statistical population, or a sample based estimate of that quantity. An effect size calculated from data is a descriptive statistic that… …   Wikipedia

  • Info-gap decision theory — is a non probabilistic decision theory that seeks to optimize robustness to failure – or opportuneness for windfall – under severe uncertainty,[1][2] in particular applying sensitivity analysis of the stability radius type[3] to perturbations in… …   Wikipedia

  • Principal component analysis — PCA of a multivariate Gaussian distribution centered at (1,3) with a standard deviation of 3 in roughly the (0.878, 0.478) direction and of 1 in the orthogonal direction. The vectors shown are the eigenvectors of the covariance matrix scaled by… …   Wikipedia

  • Maximum spacing estimation — The maximum spacing method tries to find a distribution function such that the spacings, D(i), are all approximately of the same length. This is done by maximizing their geometric mean. In statistics, maximum spacing estimation (MSE or MSP), or… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”