1 / 32

Lecture 7 Model Checking for Linear Mixed Models for Longitudinal Data

Lecture 7 Model Checking for Linear Mixed Models for Longitudinal Data. Ziad Taib Biostatistics, AZ MV, CTH May 2009. Outline of lecture 7. Introduction to model checking Model checking for the linear model Model checking for the linear mixed models for longitudinal data. Select.

lourdes
Download Presentation

Lecture 7 Model Checking for Linear Mixed Models for Longitudinal Data

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 7Model Checking for Linear Mixed Models for Longitudinal Data Ziad Taib Biostatistics, AZ MV, CTH May 2009 Name, department

  2. Outline of lecture 7 • Introduction to model checking • Model checking for the linear model • Model checking for the linear mixed models for longitudinal data Name, department

  3. Select Summarize Model Class Some Models Conclusions Data Stop 1. Introductio to model checking • The process of statistical analysis might take the form

  4. In the above process, however, even after a careful selection of model class, the data themselves may indicate that the particular model is unsuitable. Thus, it seems to be reasonable to introduce model checking to the original process. The news process of statistical analysis is Select Model Class Some Models Conclusions Summarize Data Stop Model Checking

  5. The inadequacy indicated by model checking could take two forms and is part of the technique of model checking. • The detection of systematic discrepancies. It may be e.g. that the data as a whole show some systematic departure from the fitted model. An example of this type is informal checking using residuals. • The detection of isolated discrepancies. It may be that a few data values are discrepant from the rest. This can be done using measures of leverage or measures of influence

  6. 2. The linear model Model checking for linear models uses mainly the following statistics: • The fitted values: • The mean residual sum of square: • The residual:

  7. Residual checking • Plot residuals against mean

  8. Checks for Isolated Departures from the Model • In the case of the standard linear model, the Cook’s distance can be used to assess the influence of observation i, by considering the parameter estimate without the contribution from the i’th observation:

  9. 3. Model checking in linear mixed models

  10. 3.1 Model selection: likelihood • When choosing between different models we want to be able to decide which model fits our data best. If the models compared are nested within each other it is possible to do a likelihood ratio test where the test statistic has an approximate distribution. The test statistic for the likelihood statistic is, • where DF are the degrees of freedom which is the difference in number of parameters for the models and L1 and L2 are the likelihoods for the first and second model respectively. Name, department

  11. If the two models compared are not nested with each other but contain the same number of parameters they can be compared directly by looking at the log likelihood and the model with the biggest likelihood value wins • If the two models are not nested and contain different number of parameters the likelihood can not be used directly. It is still possible to compare these models with some of the methods described below. • The bigger the likelihood is the better the model fits data and we use this when we compare different models • Since we are interested in getting as simple models as possible we also have to consider the number of parameters in the structures. A model with many parameters usually fits data better than a model with less number of parameters Information. Name, department

  12. 3.2 Model selection: Information criteria • It is possible to compute so called information criteria and there are different ways to do that and here we show two of these, Akaikes information criteria (AIC) and Bayesian information criteria (BIC). The idea with both of these is to punish models with many parameters in some way. We present the information criteria the way they are computed in SAS. • The AIC value is computed as below where q is the number of parameters in the covariance structure. Formulated this way, a smaller value of AIC indicates a better model. • The BIC value is computed using the following formula where q is the number of parameters in the covariance structure and n is the number of effective observations, which means the number of individuals. Like for AIC a smaller value of BIC is better than a larger. Name, department

  13. Model fit • It is possible to define a goodness of fit measure similar to, R, the coefficient of determination often used for linear models. It is called Concordance Correlation coefficient (CCC). Unlike the AIC or the BIC, the CCC does not compare the model at hand to other models, thus it does not require that other models be fitted. • For simple linear regression we have

  14. 3.3 Residuals for linear mixed models • In model selection, we accept the model with the best likelihood value in relation to the number of parameters but we still do not know if the model chosen is a good model or even if the normality assumption we have made is realistic.To check this we can look at two types of plots for our data, • normal plots • residual plots to check • normality of the residuals and the random effects • if the residuals seem to have a constant variance • outliers Name, department

  15. The predicted values and residuals can be computed in many different ways. Some of these are accounted for in the in what following. • Recall that the general linear mixed model is of the form: • Assuming we have ML estimates of the fixed parameters and EB predictions of the random parameters Name, department

  16. We can “estimate” the residuals according to the following three methods: • Each type of residual is useful to evaluate some of the assumptions of the model. Name, department

  17. Marginal residuals • Can be used to assess linearity of the response wrt explanatory variables. A random behviour around zero is a sign of linearity. Conditional residuals • Plots of e/s against Y can be used to assess homogeneity of the variances as well as normality. EBLUP • Plots of bi against subject indices can be used to find outliers. Plot elements in bi b to assess normality and check for outliers. Name, department

  18. 3.4 An example To illustrate the above procedures, we analyze data from a study conducted at the School of Dentistry of the University of So Paulo, Brazil, designed to compare a low cost toothbrush (monoblock) with a conventional toothbrush with respect to the maintenance of the capacity to remove bacterial plaque under daily use. The data in the table correspond to bacterial plaque indices obtained from 32 children aged 4 to 6 before and after tooth brushing in four evaluation sessions. Name, department

  19. Following Singer et al. (2004) who analyze a different data set from the same study, we considered fitting models of the form Name, department

  20. (3.1) (3.2) (3.3) Three possible models i subject d session j type of toothbrush pre post

  21. The model reduction procedure can be based on likelihood ratio tests (LRT) and AIC and BIC: The LRT p-values corresponding to the reduction of (3.1) to (3.2) and of (3.2) to (3.3) were, respectively 0.3420 and 0.1623. The AIC (BIC) for the three models are AIC (BIC) (3.1) 95.0 68.6 (3.2) 102.8 86.7 (3.3) 105.6 92.1 Based on these results, we adopt (3.3) to illustrate the use of the proposed diagnostic procedures.

  22. To check for the linearity of effects, we plot the marginal residuals versus the logarithms of the pretreatment bacterial plaque index in Figure 2. The figure supports the regression model for the transformed response (log of the bacterial plaque index) Figure 2 Name, department

  23. The figure suggests something is wrong with observations #12.2 and #29.4

  24. Name, department

  25. References • Atkinson, C. A. (1985). Plots, transformations, and regression: an introduction to graphical methods of diagnostic regression analysis. Oxford University Press, Oxford. • Cook, R. D. and Weisberg, S. (1982). Residuals and influence regression. Chapman & Hall, New York. • Cox, D. R. and Snell, E. J. (1968). A general definition of residuals (with discussion). Journal Royal Statistical Society B 30, 248–275. • Fei, Y. and Pan, J. (2003). Influence assessments for longitudinal data in linear mixed models. In 18th international workshop on Statistical Modelling. G. Verbeke, G. Molenberghs, M. Aerts and S. Fieuws (eds.). Leuven: Belgium, 143–148. • Grady, J. J. and Helms, R.W. (1995). Model selection techniques for the covariance matrix for incomplete longitudinal data. Statistics in Medicine 14, 1397–1416. Name, department

  26. References 6. Jiang, J. (2001). Goodness-of-fit tests for mixed model diagnostics. The Annals of Statistics 29, 1137–1164. 7. Lange, N. and Ryan, L. (1989). Assessing normality in random effects models. The Annals of Statistics 17, 624– 642. 8. Longford, N. T. (2001). Simulation-based diagnostics in random-coefficient models. Journal of the Royal Statistical Society A 164, 259–273. 9. Nobre, J. S. and Singer, J. M. (2006). Fixed and random effects leverage for influence analysis in linear mixed models. (Submitted; http://www.ime.usp.br/jmsinger). 10. Oman, S. D. (1995). Checking the assumptions in mixed-model analysis of variance: a residual analysis approach. Computational Statistics and Data Analysis 20, 309–330.

  27. References 11. Verbeke, G. and Lesaffre, E. (1997). The effect of misspecifying the random-effects distributions in linear mixed models for longitudinal data. Computational Statistics and Data Analysis 23, 541–556. 12. Waternaux, C., Laird, N. M., and Ware, J. H. (1989). Methods for analysis of longitudinal data: blood-lead concentrations and cognitive development. Journal of the American Statistical Association 84, 33–41. 13. Weiss, R. E. and Lazaro, C. G. (1992). Residual plots for repeated measures. Statistics in Medicine 11, 115–124. 14. Wolfinger, R. (1993). Covariance structure selection in general mixed models. Communications in Statistics-Simulation 22, 1079–1106.

  28. Any Questions? Name, department

More Related