790 likes | 852 Views
Hypothesis testing for the GLM. The General Linear Hypothesis. Testing the General Linear Hypotheses The General Linear Hypothesis. H 0 : h 11 b 1 + h 12 b 2 + h 13 b 3 +... + h 1 p b p = h 1 h 21 b 1 + h 22 b 2 + h 23 b 3 +... + h 2 p b p = h 2 ...
E N D
Hypothesis testing for the GLM The General Linear Hypothesis
Testing the General Linear Hypotheses The General Linear Hypothesis H0: h11b1 + h12b2 + h13b3 +... + h1pbp = h1 h21b1 + h22b2 + h23b3 +... + h2pbp = h2 ... hq1b1 + hq2b2 + hq3b3 +... + hqpbp = hq where h11, h12, h13, ... , hqp and h1, h2, h3, ... , hq are known coefficients. In matrix notation
Examples • H0: b1 = 0 • H0: b1 = 0, b2 = 0, b3 = 0 • H0: b1 = b2
Examples • H0: b1 = b2 , b3 = b4 • H0: b1 = 1/2(b2 + b3) • H0: b1 = 1/2(b2 + b3), b3 = 1/3(b4 + b5 + b6)
TheLikelihood Ratio Test The joint density of is: The likelihood function The log-likelihood function
Defn (Likelihood Ratio Test of sizea) Rejects H0:q w against the alternative hypothesis H1:q w. when and K is chosen so that
Note The Lagrange multiplier technique will be used for this purpose We will maximize.
or finally or
Thus the equations for are Now or and
Thus Note Now
The Likelihood Ratio Test of sizea) Rejects H0:q w against the alternative hypothesis H1:q w. when and K is chosen so that
Thus The LR tests rejects H0:q w when
Theorem If is true then has an F distribution with n1 = q d.f. in the numerator and n2 = n – p d.f. in the denominator.
Proof Recall Thus and Hence Also recall
Finally Thus is independent of has an F distribution with n1 = q d.f. in the numerator and n2 = n – p d.f. in the denominator.
Special Case Testing vs Thus the test would reject when
Testing if a sub vector is zero (an alternative formulation)
An Alternative form of the F statistic Exercise: Show
Example: General Linear Model with intercept Assume we have collected data on Y, X1, X2, X3 and Y = b0 + b1X1 + b1X2 + b1X3 + e The Data
Confidence intervals, Prediction intervals, Confidence Regions General Linear Model
i.e. is a (1 – a)100 % confidence interval for Special cases
i.e. are (1 – a)100 % confidence interval for s2 and s.