> Hmmm. Thinking some more, I might not have answered your (unstated)
> question. What do your mean by GLM?
>
> => I meant generalized linear model as well. Thanks for the references.
    The first one was mentioned first in my life time after keeping asking
the same question.


> I mean the Generalized Linear Model, not the General Linear Model. The
> Generalized one allows for non-normal responses and different
> mean-variance relationships and is the GLM of Nelder and Wedderburn
> (1972, J. Royal Statistical Society, Series A, 135(3),370-384) and the
> monograph by McCullagh and Wedderburn (1989, Generalized Linear Models,
> Chapman & Hall/CRC). The R function glm() fits these kinds of model.
>
> The General Linear Model (IIRC) was the linking of linear regression and
> anova into a single entity. The R function lm() fits these kinds of
> models.
>
> The linear model is a special case of the Generalized Linear Model when
> the Gaussian error is used with the identity link function. Hence a
> Gaussian GLM (my GLM) with the identity link fitted by glm() will give
> the same results as lm(), but it will do so in a very inefficient
> manner. As this was what your code was doing I suggested using lm()
> instead.
>

=> Thanks for the important suggestion.
    I tried both lm and glm until now but no one commented on the method.

>
> Elaine
>

	[[alternative HTML version deleted]]

