1. General Linear Model With correlated error terms S=s2V ≠ s2I

2. The General Linear ModelS≠s2I

6. Summary

7. Example Simple Linear Model where variance is proportional to X2.

11. Testing and Confidence Intervals The Model: can be converted to the model

12. Thus Simultaneous Confidence Intervals (using model (2))

13. Example: Simple Linear Model with no intercept The model

14. Thus

15. Also Special cases

16. General Linear Model Case 2: Sunknown

17. The General Linear ModelSunknown

18. The General Linear ModelSunknown Call this the Ordinary Least Squares (OLS) estimator of Note: Thus the Ordinary Least Squares (OLS) estimator of is always unbiased.

19. This is the Optimal (UMVU) estimator of Note: This is also an unbiased estimator of The Optimal (UMVU) estimator of requires the knowledge of Sin order to calculate it.

20. Theorem: Equivalence of OLS estimator with UMVU estimator

21. Proof

22. Application:Consider the general linear model with intercept In this case the error terms are equally correlated. Also in this case the OLS estimators are equivalent to the UMVU estimators

23. Proof

24. Design Matrix, X, not of full rank

25. The General Linear Model

26. If the rank of X is equal to p then the columns of X are linearly independent and there is a unique way of representing If the rank of X is strictly less than p then there is no unique way of representing

27. Comment: Usually the situation where the rank of X, r < p, arises in the following instances. • The design of the study (the choice of the values of X1, X2, …, Xp) was not careful enough to ensure that X had full rank. • Observations were missing causing the model to be altered Elements of are deleted along with corresponding rows of X, reducing the number of linear independent rows from p to r. • The model was defined in such a way that: • mi = b1xi1 + b2xi2 +… +bpxip • is not uniquely determined by b1, b2,… ,bp.

28. Two Basic approaches: • Impose p – r linear restrictions on the parameters • This allows us to reduce the number of parameters to r. • will have a unique representation if the p – r restrictions are added. • This technique is usually used with ANOVA, MANOVA, ANACOVA models. • Live with the singularity. • Restrict our attention to linear combinations of the parameters that have unique estimators. • The two approaches are essentially the same (lead to the same conclusions).

29. Recall: Linear Equations theory Consider the system of linear equations M(A), the linear space spanned by the columns ofA

30. Then the general solution to the system of linear equations is

31. Maximum Likelihood Estimation leads to the system of linear equations • p equations in p unknowns • called the Normal equations

32. Theorem The Normal equations are consistent. Proof It can be shown that M(XX)  M(X) M(X) M(XX) Theorem The general solution to the Normal equations is

33. Theorem is the same for all solutions of the Normal equations Proof: the general solution to the Normal equations is Since M(XX)  M(X) there exists a p × n matrix L such that X = XXL or X = LXX

35. Definition: (Estimability) The linear function of the parameter vector, is called estimable if there exists a vector such that Example The simple linear model

36. Thus is the only estimable function of b0, b1.

37. Theorem: The following conditions are equivalent 2. For some solution, , of the Normal equation, , is a linear (in ) unbiased estimate of M(XX) M(X)

38. Proof: Assume Then there exists a vector such that M(X) Thus Thus 1. implies 5. (as well as 4.) Now assume 4.

39. Thus 4. implies 3. Thus 4. implies 2. and 1.

40. Suppose we have k normal populations Example: One-way ANOVA (Analysis of Variance) Let yi1, yi2, … , yindenote a sample of n from Let eij= yij - (m + ai), then ei1, ei2, … , eindenotes a sample of n from distribution. where e11, e12, … , eknare kn independent observations from N(0,s2) distribution.

41. Let Matrix Notation

42. Let Then the model is

43. M(X) = then linear space spanned by the vectors

44. Thus the estimable parameters are of the form: The common approach is to add the restriction This reduces the number of parameters to k, and converts the model to full rank.

45. Properties of estimable functions: All linear functions are estimable Proof If rank(X) = p then M(X) = Ep =p-dimensional Euclidean space (which contains all p-dimansional vectors) is estimable if Proof (unique for all solutions of the normal equations) Hence is estimable.

46. If and are estimable then Proof since and are estimable then