Topic 30: Random Effects Modeling

Transcription

1 Topic 30: Random Effects Modeling

2 Outline One-way random effects model Data Model Inference

3 Data for one-way random effects model Y, the response variable Factor with levels i = 1 to r Y ij is the j th observation in cell i, j = 1 to n i Almost identical model structure to earlier one-way fixed effect ANOVA. Difference in level of inference

4 Level of Inference In one-way fixed effects ANOVA, interest was in comparing the factor level means In random effects scenario, interest is in the population of factor level means, not just the means of the r study levels Need to make assumptions about population distribution Will take random draw from the population of factor levels for use in study

5 KNNL p 1036 KNNL Example Y is the rating of a job applicant Factor A represents five different personnel interviewers (officers), r=5 n=4 different applicants were interviewed by each officer The interviewers were selected at random and the applicants were randomly assigned to interviewers

6 Read and check the data data a1; infile 'c:\...\ch5ta01.dat'; input rating officer; proc print data=a1; run;

7 The data Obs rating officer

8 The data Obs rating officer

9 Plot the data title1 'Plot of the data'; symbol1 v=circle i=none c=black; proc gplot data=a1; plot rating*officer/frame; run;

10

11 Find and plot the means proc means data=a1; output out=a mean=avrate; var rating; by officer; title1 'Plot of the means'; symbol1 v=circle i=join c=black; proc gplot data=a; plot avrate*officer/frame; run;

12

13 Random effects model Y ij = μ i + ε ij the μ i are iid N(μ, σ μ ) the ε ij are iid N(0, σ ) μ i and ε ij are independent Y ij ~ N(μ, σ μ + σ ) Two sources of variation Key difference Observations with the same i are not independent, have covariance σ μ

14 Random effects model This model is also called Model II ANOVA A variance components model The components of variance are σ μ and σ The models that we previously studied are called fixed effects models

15 Random factor effects model Y ij = μ + τ i + ε ij τ i ~ N(0, σ μ ) ***** ε ij ~ N(0, σ ) Y ij ~ N(μ, σ μ + σ )

16 Parameters There are three parameters in these models μ σ μ σ The cell means (or factor levels) are random variables, not parameters Inference focuses on these variances

17 Primary Hypothesis Want to know if H 0 : σ μ = 0 This implies all µ i in model are equal but also all µ i not selected for analysis are also equal. Thus scope is broader than fixed effects case Need the factor levels of the study to be representative of the population

18 Alternative Hypothesis We are sometimes interested in estimating σ μ /(σ μ + σ ) This is the same as σ μ /σ Y In some applications it is called the intraclass correlation coefficient It is the correlation between two observations with the same i Also percent of total variation of Y

19 ANOVA table The terms and layout of the ANOVA table are the same as what we used for the fixed effects model The expected mean squares (EMS) are different because of the additional random effects but F test statistics are the same Be wary that hypotheses being tested are different

20 EMS and Parameter Estimates E(MSA) = σ + nσ μ E(MSE) = σ We use MSE to estimate σ Can use (MSA MSE)/n to estimate σ μ Question: Why might it we want an alternative estimate for σ μ?

21 Main Hypotheses H 0 : σ μ = 0 H 1 : σ μ 0 Test statistic is F = MSA/MSE with r-1 and n T -r degrees of freedom, reject when F is large, report the P-value

22 Run proc glm proc glm data=a1; class officer; model rating=officer; random officer/test; run;

23 Model and error output Source DF MS F P Model Error Total 19

24 Random statement output Source Type III Expected MS officer Var(Error) + 4 Var(officer)

25 Proc varcomp proc varcomp data=a1; class officer; model rating=officer; run;

26 Output MIVQUE(0) Estimates Variance Component rating Var(officer) Var(Error) Other methods are available for estimation, mivque is the default

27 Proc mixed proc mixed data=a1 cl; class officer; model rating=; random officer/vcorr; run;

28 Output Covariance Parameter Estimates Cov Parm Est Lower Upper officer Residual /( )=.53

29 Output from vcorr Row Col1 Col Col3 Col

30 Other topics Estimate and CI for μ, p1038 Standard error involves a combination of two variances Use MSA instead of MSE r-1 df Estimate and CI for σ μ /(σ μ + σ ), p1040 CIs for σ μ and σ, p Available using Proc Mixed

31 Applications In the KNNL example we would like σ μ /(σ μ + σ ) to be small, indicating that the variance due to interviewer is small relative to the variance due to applicants In many other examples we would like this quantity to be large, e.g., Are partners more likely to be similar in sociability?

32 Last slide Start reading KNNL Chapter 5 We used program topic30.sas to generate the output for today