A monotone data augmentation algorithm for longitudinal data analysis via multivariate skew-t,skew-normal or t distributions

Abstract

The mixed effects model for repeated measures has been widely used for the analysis of longitudinal clinical data collected at a number of fixed time points. We propose a robust extension of the mixed effects model for repeated measures for skewed and heavy-tailed data on basis of the multivariate skew-t distribution, and it includes the multivariate normal, t, and skew-normal distributions as special cases. An efficient Markov chain Monte Carlo algorithm is developed using the monotone data augmentation and parameter expansion techniques. We employ the algorithm to perform controlled pattern imputations for sensitivity analyses of longitudinal clinical trials with nonignorable dropouts. The proposed methods are illustrated by real data analyses. Sample SAS programs for the analyses are provided in the online supplementary material.

Keywords

Block sampling controlled imputations mixed effects model for repeated measures monotone data augmentation penalized complexity prior tipping point analysis

Get full access to this article

View all access options for this article.

References

Siddiqui

Hung

JHM

O'Neill

. MMRM vs. LOCF: A comprehensive comparison based on simulation study and 25 NDA datasets. J Biopharmaceut Stat 2009; 19: 227–246.

Mallinckrodt

Lane

Schnell

, et al. Recommendations for the primary analysis of continuous endpoints in longitudinal clinical trials. Drug Inform J 2008; 42: 303–319.

ICH E9. Statistical principles for clinical trials: ICH harmonized tripartite guideline. Stat Med 1999; 18: 1905– 1942.

CHMP. Guideline on adjustment for baseline covariates in clinical trials. London: CHAMP, 2015.

Laird

Lange

Stram

. Maximum likelihood computations with repeated measures: application of the EM algorithm. J Am Stat Assoc 1987; 82: 97–105.

Tang

. Closed-form REML estimators and sample size determination for mixed effects models for repeated measures under monotone missingness. Stat Med 2017; 36: 2135–2147.

Mehrotra

. Specification of covariance structure in longitudinal data analysis for randomized clinical trials. Stat Med 2010; 29: 474–488.

Gurka

Edwards

Muller

. Avoiding bias in mixed model inference for fixed effects. Stat Med 2011; 30: 2696–2707.

Jacqmin-Gadda

Sibillot

Proust

, et al. Robustness of the linear mixed model to misspecified error distribution. Computat Stat Data Analys 2007; 51: 5142–5154.

10.

Gomez

Schaalje

Fellingham

. Performance of the Kenward-Roger method when the covariance structure is selected using AIC and BIC. Commun Stat-Simul Computat 2005; 34: 377–392.

11.

Little

Yau

. Intent-to-treat analysis for longitudinal studies with drop-outs. Biometrics 1996; 52: 1324–1333.

12.

Carpenter

Roger

Kenward

. Analysis of longitudinal trials with protocol deviation: a framework for relevant, accessible assumptions, and inference via multiple imputation. J Biopharmaceut Stat 2013; 23: 1352–1371.

13.

Mallinckrodt

Roger

Chuang-stein

, et al. Missing data: turning guidance into action. Stat Biopharmaceut Res 2013; 5: 369–382.

14.

Tang

. An efficient multiple imputation algorithm for control-based and delta-adjusted pattern mixture models using SAS. Stat Biopharmaceut Res 2017; 9: 116–125.

15.

Rubin

. Inference and missing data. Biometrika 1976; 63: 581–92.

16.

ICH E9 (R1) addendum on estimands and sensitivity analysis in clinical trials to the guideline on statistical principles for clinical trials. http://www.ema.europa.eu/docs/en_GB/document_library/Scientific_guideline/2017/08/WC500233916.pdf, 2017.

17.

CHMP. EMA guideline on missing data in confirmatory clinical trials (EMA/CPMP/EWP/1776/99). London: CHAMP, 2010.

18.

National Research Council. Panel on handling missing data in clinical trials. Committee on National Statistics, Division of Behavioral and Social Sciences and Education. The prevention and treatment of missing data in clinical trials. Washington, DC: The National Academies Press, 2010.

19.

Little

RJA

. Pattern-mixture models for multivariate incomplete data. J Am Stat Assoc 1993; 88: 125–134.

20.

Tang

. An efficient monotone data augmentation algorithm for multiple imputation in a class of pattern mixture models. J Biopharmaceut Stat 2017; 27: 620–638.

21.

Schafer JL. Analysis of incomplete multivariate data. London: Chapman Hall, 1997.

22.

Tang

. Controlled pattern imputation for sensitivity analysis of longitudinal binary and ordinal outcomes with nonignorable dropout. Stat Med 2018; 37: 1467–1481.

23.

Tang

. A monotone data augmentation algorithm for multivariate nonnormal data: with applications to controlled imputations for longitudinal trials. Stat Med 2019; 38: 1715–1733.

24.

Azzalini

. A class of distributions which includes the normal ones. Scand J Stat 1985; 12: 171–178.

25.

Azzalini

Capitanio

. Distributions generated by perturbation of symmetry with emphasis on a multivariate skew t distribution. J Royal Stat Soc, Ser B 2003; 65: 367–389.

26.

Verbeke

Lesaffre

. The effect of misspecifying the random effects distribution in linear mixed models for longitudinal data. Computat Stat Data Analysis 1997; 23: 541–556.

27.

Pinheiro

Liu

. Efficient algorithms for robust estimation in linear mixed-effects models using a multivariate t-distribution. J Computat Graphics Stat 2001; 10: 249–276.

28.

Zhang

Davidian

. Linear mixed models with flexible distributions of random effects for longitudinal data. Biometrics 2001; 57: 795–802.

29.

Arnau

Bendayan

Blanca

, et al. The effect of skewness and kurtosis on the robustness of linear mixed models. Behavior Res Meth 2013; 45: 873–879.

30.

Boos

Hughes-Oliver

. How large does n have to be for z and t intervals?. Am Stat 2000; 54: 121–128.

31.

von Hippel

. Should a normal imputation model be modified to impute skewed variables?. Sociol Meth Res 2013; 42: 105–138.

32.

Azzalini

Dalla Valle

. The multivariate skew-normal distribution. Biometrika 1996; 83: 715–726.

33.

Liu

. Missing data imputation using the multivariate t distribution. J Multivariate Analys 1995; 53: 139–158.

34.

Jara

Quintana

Martin

. Linear mixed models with skew-elliptical distributions: a Bayesian approach. Computat Stat Data Analys 2008; 52: 5033–5045.

35.

Lachos

Dey

Cancho

. Robust linear mixed models with skew-normal independent distributions from a Bayesian perspective. J Stat Plan Inference 2009; 139: 4098–4110.

36.

Arnold

Beaver

. Skewed multivariate models related to hidden truncation and/or selective reporting. Test 2002; 11: 7–54.

37.

Arellano-Valle

Branco

Genton

. A unified view on skewed distributions arising from selections. Canadian J Stat 2006; 34: 581–601.

38.

Liseo

Loperfido

. A note on reference priors for the scalar skew-normal distribution. J Stat Plan Inference 2006; 136: 373–389.

39.

Bayes

Branco

. Bayesian inference for the skewness parameter of the scalar skew-normal distribution. Brazilian J Probabil Stat 2007; 21: 141–163.

40.

Branco

MGG Liseo

. Objective Bayesian analysis of skew-t distributions. Scand J Stat 2013; 40: 63–85.

41.

Pewsey

. Problems of inference for Azzalini's skew-normal distribution. J Appl Stat 2000; 27: 859–870.

42.

Arellano-Valle

Azzalini

. The centred parametrization for the multivariate skew-normal distribution. J Multivariate Analys 2008; 99: 1362–1382.

43.

Liseo

Parisi

. Bayesian inference for the multivariate skew-normal model: a population Monte Carlo approach. Computat Stat Data Analys 2013; 63: 125–138.

44.

Fonseca

TCO

Ferreira

MAR

Migon

. Objective Bayesian analysis for the Student-t regression model. Biometrika 2008; 95: 325–333.

45.

Liu

. Parameter expansion for data augmentation. J Am Stat Assoc 1999; 94: 1264–1274.

46.

Liu

Sabatti

. Generalised Gibbs sampler and multigrid Monte Carlo for Bayesian computation. Biometrika 2000; 87: 353–369.

47.

Lee

McLachlan

. On mixtures of skew normal and skew t-distributions. Adv Data Analysis Classification 2013; 7: 241–266.

48.

Fruhwirth-Schnatter

Pyne

. Bayesian inference for finite mixtures of univariate and multivariate skew-normal and skew-t distributions. Biostatistics 2010; 11: 317–336.

49.

Huang

Wand

. Simple marginally noninformative prior distributions for covariance matrices. Bayesian Analysis 2013; 8: 439–452.

50.

Gelman

. Prior distributions for variance parameters in hierarchical models. Bayesian Analysis 2006; 1: 515–534.

51.

Tang

. An efficient monotone data augmentation algorithm for Bayesian analysis of incomplete longitudinal data. Stat Probabil Lett 2015; 104: 146–52.

52.

Simpson

Rue

Riebler

, et al. Penalising model component complexity: a principled, practical approach to constructing priors. Stat Sci 2017; 32: 1–28.

53.

Daniels

Kass

. Nonconjugate bayesian estimation of covariance matrices and its use in hierarchical models. JASA 1999; 94: 1254–1263.

54.

Permutt

. Sensitivity analysis for missing data in regulatory submission. Stat Med 2016; 35: 876–879.

55.

Koch

. Comparison between two controlled multiple imputation methods for sensitivity analyses of time-to-event data with possibly informative censoring. Stat Biopharmaceut Res 2015; 7: 199–213.

56.

Tang

. On the multiple imputation variance estimator for control-based and delta-adjusted pattern mixture models. Biometrics 2017; 73: 1379–1387.

57.

Tang

. Controlled pattern imputation for sensitivity analysis of longitudinal binary and ordinal outcomes with nonignorable dropout. Stat Med 2018; 37: 1467–1481.

58.

Tang

. Algorithms for imputing partially observed recurrent events with applications to multiple imputation in pattern mixture models. J Biopharmaceut Stat 2018; 28: 518–533.

59.

Barnard

Rubin

. Small-sample degrees of freedom with multiple imputation. Biometrika 1999; 86: 948–955.

60.

Spiegelhalter

Best

Carlin

, et al. Bayesian measures of model complexity and fit. J Royal Stat Soc, Ser B 2002; 64: 583–639.

61.

Ghidey

Lesaffre

Eilers

. Smooth random effects distribution in a linear mixed model. Biometrics 2004; 60: 945–953.

62.

Cabral

CRB

Lachos

Madruga

. Bayesian analysis of skew-normal independent linear mixed models with heterogeneity in the random-effects population. J Stat Plan Inference 2012; 142: 181–200.

63.

Lachos

Ghosh

Arellano-Valle

. Likelihood based inference for skew normal independent linear mixed models. Statistica Sinica 2010; 20: 303–322.

64.

Chib

Greenberg

. Understanding the metropolis-hastings algorithm. Am Stat 1995; 49: 327–335.

65.

Senn

. Testing for baseline balance in clinical trials. Stat Med 1994; 13: 1715–1726.

66.

. On efficiency of constrained longitudinal data analysis versus longitudinal analysis of covariance. Biometrics 2010; 66: 891–896.

67.

. An analytic method for the placebo-based pattern-mixture model. Stat Med 2014; 33: 1134–1145.

68.

Kleinman

Ibrahim

. A semiparametric bayesian approach to the random effects model. Biometrics 1998; 54: 921–938.

69.

Vasconcellos

KLP

Da Silva

. Corrected estimates for student t regression models with unknown degrees of freedom. J Stat Computat Simulat 2005; 52: 409–423.

70.

Lange

Little

RJA

Taylor

JMG

. Robust statistical modeling using the t-distribution. J Am Stat Assoc 1989; 84: 881–896.

71.

Lukacs

Burnham

Anderson

. Model selection bias and Freedman's paradox. Ann Inst Stat Math 2009; 62: 117–125.

72.

Mundry

Nunn

. Stepwise model fitting and statistical inference: turning noise into signal pollution. Am Naturalist 2009; 173: 119–123.

73.

Wang

Luo

. Bayesian multivariate augmented beta rectangular regression models for patient-reported outcomes and survival data. Stat Meth Med Res 2017; 26: 117–125.

74.

Glass

Peckham

Sanders

. Consequences of failure to meet assumptions underlying the fixed effects analyses of variance and covariance. Rev Educ Res 1972; 42: 237–288.

75.

Mehrotra

Liu

, et al. Analysis of longitudinal clinical trials with missing data using multiple imputation in conjunction with robust regression. Biometrics 2012; 68: 1250–1259.

76.

Kotz

Nadarajah

. Multivariate t-distributions and their applications, Cambridge: Cambridge University Press, 2004.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.02 MB

0.00 MB

0.23 MB