1 / 43

Lecture 7: Parametric Models for Covariance Structure (Examples)

Lecture 7: Parametric Models for Covariance Structure (Examples). 1. Model for the mean. 2. Model for the covariance matrix. 2. Model for the covariance matrix (cont’d). Which model to pick?. Which model to pick? (cont’d). Example: CD4+ Level.

amosr
Download Presentation

Lecture 7: Parametric Models for Covariance Structure (Examples)

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Lecture 7:Parametric Models for Covariance Structure (Examples)

  2. 1. Model for the mean

  3. 2. Model for the covariance matrix

  4. 2. Model for the covariance matrix (cont’d)

  5. Which model to pick?

  6. Which model to pick? (cont’d)

  7. Example: CD4+ Level HIV attacks CD4+ cells, which regulate the body’s immuneresponse against infectious agents We have 2376 values of CD4+ cell numbers plotted against time since seroconversion for 369 infected men enrolled in the MAC Study

  8. Example: CD4+ Level (cont’d) • Goals: • Estimate the average time course of CD4+ cell depletion • Identify factors which predict CD4+ cell changes • Estimate the time course for an individual man taking into account the measurement error in CD4+ cell determinations • Characterize the degree of heterogeneity across men in the rate of progression

  9. Example: CD4+ LevelGoal 1: Estimate average time course of CD4+ cell depletion The model for the covariance matrix is a model of serial correlation.

  10. Example: CD4+ LevelGoal 2: Identify factors predictive of CD4+ cell changes The model for the covariance matrix is still a model of serial correlation, however we have changed the model for the mean.

  11. Example: CD4+ LevelParameter Interpretation

  12. Example: CD4+ LevelGoal 3: Estimate time course for an individual man, accounting for measurement error in CD4+ cell counts This is a model with a random intercept and slope + serial correlation + measurement error.

  13. Example: CD4+ LevelParameter Interpretation

  14. Random Effects Model:Interpretation of coefficients Heterogeneity between subjects atbaseline Heterogeneity between subjects in rate of change

  15. Example: CD4+ LevelGoal 4: Characterize degree of heterogeneity across men in progression rate (From the previous slide)

  16. Example: Protein contents of milk samples

  17. Barley (25 cows) • Mixed (27 cows) • Lupins (27 cows)

  18. Example: Protein contents of milk samples (cont’d)

  19. 0.02

  20. Example: Protein contents of milk samplesModel for the Mean

  21. Example: Protein contents of milk samplesModel for the Covariance Matrix

  22. Table 5.1 (b)

  23. Example: Protein contents of milk samplesDoes diet affect the mean response profile?

  24. Example: Protein contents of milk samplesIs there a rise in the mean response towards the end of the experiment?

  25. Example: Protein contents of milk samples Is there a rise in the mean response towards the end of the experiment? (cont’d)

  26. Example: Protein contents of milk samples Is there a rise in the mean response towards the end of the experiment? (cont’d)

  27. 0.02

  28. Example: Body weight of 26 cows Dataset consists of body weights of 26 cows, measured at 23 unequally-spaced times over a period of about 22 months. Treatments were allocated in a 2x2 factorial design: Control (4) Iron-dosing (4) Infection (9) Iron + Infection (10)

  29. Control N=4 Iron N=4 Iron + Infection N=10 Infection N=9 Log Y: Variance-stabilizing transformation

  30. Example: Body weight of 26 cows (cont’d) • Look at your data • Estimate empirical variogram • What do you see? • Measurement variance small • Substantial between-cow variability • Gaussian correlation model appropriate

  31. Empirical variogram of the OLS residuals from a saturated model for the mean response • Small measurement error • Experimental correlation • Random effects

  32. Example: Body weight of 26 cowsModel for the Mean

  33. Example: Body weight of 26 cowsModel for the Covariance Matrix

  34. Example: Body weight of 26 cowsQuestions Q1: Can we conclude linear (vs. quadratic) growth? …The quadratic curve is appropriate. Q2: Is there a main effect for iron?…NO

  35. Example: Body weight of 26 cowsQuestions (cont’d) Q3: Is there a main effect for infection?…YES Q4: Is there an interaction between iron and infection?…NO

  36. Example: Body weight of 26 cowsQuestions (cont’d) We re-fit the model with only the infection term: • Conclusions: • Highly significant effect of infection • No significant effect of iron • No significant effect of interaction

More Related