Assessing Informative Drop-out in Models for Repeated Binary Data by Lee Shean Er B.Sc, University of Guelph 1998 A THESIS SUBMITTED IN PARTIAL F U L F I L L M E N T OF THE REQUIREMENTS FOR THE D E G R E E OF M a s t e r of Science in THE FACULTY OF GRADUATE STUDIES (Department of Statistics) we accept this thesis as conforming to the required^ standard The University of British Columbia March 2001 © Lee Shean Er, 2001 In presenting this degree at the thesis in University of partial fulfilment of of department this thesis for or by his or scholarly purposes may be her representatives. permission. of The University of British Columbia Vancouver, Canada Date DE-6 (2/88) $ i M l for an advanced Library shall make it agree that permission for extensive It publication of this thesis for financial gain shall not Department requirements British Columbia, I agree that the freely available for reference and study. I further copying the is granted by the understood that head of copying my or be allowed without my written Abstract Drop-outs are a common problem in longitudinal studies. In terms of statistical models for the data, there are three types of drop-out mechanisms: drop-out occurring completely at random (CRD), drop-out occurring at random (RD) and informative drop-out (ID). The drop-out mechanism is classified as C R D if the drop-out mechanism is independent of the measurements; as R D if the drop-out mechanism depends only on the observed but not the unobserved measurements, and as ID if the drop-out mechanism depends on both the observed and unobserved measurements. C R D and R D are referred to as ignorable because the drop-out mechanism can be ignored for the purpose of making inferences about the observed measurements, while ID is non-ignorable. Analyses based on an assumption of ignorable drop-out, when in reality the drop-out mechanism is non-ignorable, can lead to misleading or biased results. Likelihood-based models for continuous and categorical longitudinal data subject to non-ignorable drop-out have been developed. In this thesis, we focus on exploring likelihood-based models for binary longitudinal data subject to informative drop-out. The two modelling approaches considered are a selection model proposed by Baker (1995) and a transition model proposed by L i u et al. (1999). We apply these models to a data set from a multiple sclerosis (MS) clinical trial. The aims of the analyses are to investigate whether there is an indication of informative drop-out in this data, and to assess the sentivity of inferences concerning the treatment effects to the underlying drop-out mechanisms. We do not attempt to provide a definitive ii analyses of the data set, but rather to explore a variety of models which incorporate informative drop-out. iii Contents Abstract ii Contents iv List of Tables viii List of Figures xii Acknowledgements xiii Dedication xiv 1 Introduction 1 1.1 Background of this Thesis 1 1.2 Methods of Analyses 2 1.3 Outline of this Thesis 5 2 Data Description 7 2.1 Description of the Berlex Clinical Trial . . '. 2.1.1 2.2 7 Drop-out Rate in this Clinical Trial Description of the Data 9 . 12 2.2.1 Drop-out Patterns 12 2.2.2 Binary Outcome Variables 14 iv 2.2.3 2.3 Baseline Covariates . 17 Questions of Interest : . . . 3 Classification of Missing Values in Longitudinal Data 4 18 22 3.1 Introduction 22 3.2 Types of Drop-outs 24 Selection Model 4.1 26 Baker's Selection Model for Binary Longitudinal Data with Informative Non-response 4.2 26 Selection Model for Binary Longitudinal Data with Informative Dropout . . . . 30 4.2.1 Outcome Model 30 4.2.2 Drop-out Model 32 4.2.3 Likelihood Function '. . 5 Transition Model 34 36 5.1 Introduction 36 5.2 The Liu et al. Transition Model for Binary Longitudinal Data with Informative Drop-out 37 5.2.1 Outcome Model 38 5.2.2 Drop-out Model 39 5.2.3 Likelihood Function . : . . 40 6 Identifiability in Models for Incomplete Binary Data 41 6.1 Introduction 41 6.2 Discussion in Fitzmaurice et al. (1996) and Glonek (1999) 43 6.2.1 Fitzmaurice et aVs Suggested Procedures 44 6.2.2 Glonek's Necessary and Sufficient Conditions 47 v 6.3 6.4 Discussion of Model Identifiability for Incomplete Binary Responses in Baker (1995) 50 Discussion of Model Identifiability 51 6.4.1 Identifiability of/i ({p,p},y3 | a;; 773) . . . 53 6.4.2 Identifiability of /i2({p}> 2 | x; r} ) . . 56 6.4.3 Identifiability of / i i ( { }iYi I \ Vi) 59 3 v 2 x 7 Application to the Data 7.1 Introduction 7.2 Baker's Selection Model: With Only Treatment Groups and Time as 7.3 7.4 8 61 61 Covariates 62 7.2.1 The Quasi-Newton (QN) Algorithm 66 7.2.2 Results 68 7.2.3 Summary 85 Baker's Selection Model: Extensions of the Drop-out Model 86 7.3.1 Results 87 7.3.2 Summary 95 Baker's Selection Model: Extension of the Outcome M o d e l . . . . . . 96 7.4.1 Results 97 7.4.2 Summary , • • 100 7.5 Overall Summary for Baker's Selection Model 100 7.6 The Liu et al. Transition Model 103 7.6.1 Results 7.6.2 Summary '. 110 Conclusions 8.1 Conclusions 8.2 Further Work 104 113 113 '. Bibliography 120 122 vi Appendix A Proof for Condition (6.4) 127 Appendix B Detailed Results for the Selection Models Described in Section 7.2 129 Appendix C Detailed Results for the Selection Models Described in Section 7.3 138 Appendix D Detailed Results for the Selection Models Described in Section 7.4 143 Appendix E Detailed Results for the Liu et al. Transition Models Described in Section 7.6 151 vn List of Tables 2.1 Cumulative Number of Drop-outs After 1, 2 and 3 Years on Study . 10 2.2 Number of Patients in the 4 Drop-out Cases (x=present,o=absent) . 13 2.3 Summary of the Our Annual Data 14 2.4 Frequency Table of the Exacerbation Counts for Patients with At Least One Outcome : 15 2.5 Number of Female and Male Patients in Each Treatment Group . . . 4.1 A l l Possible Patterns for Incomplete Data for the Case where 3 Observations were Intended for Every Unit: x = observed, o = missing. 7.1 18 28 Drop-out Models under Different Drop-out Mechanisms: 1/ denotes inclusion of a parameter and rn denotes parameters which are restricted to be equal 64 7.2 Negative Log-likelihood Values for Five Outcome Model Specifications 69 7.3 Non-response Probability for the Third Response Using Model ID1 with 7703 -> - 0 0 and 7703 + 7723 = A 7.4 77 Results for Model ID1 Evaluated on the Boundary: 7703 - » — 00 and V03 + ma = A 7.5 78 Results for Model ID2 Evaluated on the Boundary: 7703 —> - 0 0 , 7 7 0 2 —> - 0 0 , 7703 + 772 = A i and 7702 + m = A viii 2 79 7.6 Results for Model ID4 Evaluated on the Boundary: r/03 —> - 0 0 and »703 + V23 = A 7.7 79 Results for Model ID5 Evaluated on the Boundary: - 0 0 , T703 -h T72 = Ai 7702 + and m = A 7703 —> - 0 0 , 7 7 0 2 —• 80 2 7.8 Negative Log-likelihood Values for Models in Table 7.1 82 7.9 Non-response Probability for the Second and Third Responses . . . . 89 7.10 Results for Model T R T * L U R Evaluated on the Boundary: - 0 0 , r/02 -> - 0 0 , 7702 + % = Ai 7703 + and 773 = A 7703 - » 90 2 7.11 Results for Model T R T + L O R + L U R Evaluated on the Boundary: 7703 ->• - 0 0 , 7702 - > - 0 0 , 7702 + 774 = Ai and 7703 + 774 = A 91 2 7.12 Results for Model L O R * L U R Evaluated on the Boundary: W i t h 7703 - > -00,7702 Vi + V3 = A -> - 0 0 , m -°°> V02 + m = 7703 + Ai, m =. A 92 3 7.13 Results for Model T R T + L U R Evaluated on the Boundary: - 0 0 , 7702 - > and 2 - 0 0 , 7702 + 774 = Ai 7703 + and 774 = A 7703 ->• 93 2 7.14 The L R T Statistics in the Forward Stepwise Procedure 98 7.15 Data sets used for assessing the sensitivity of the results when considering log(BOD) in addition to treatment group and gender as a covariate 99 7.16 The Observed and Expected Cell Counts for Baker's Selection Model with Drop-Out Model ID5 ("*" denotes missing) 102 7.17 Results for L i u Transition Model with Drop-out Model ID1 Evaluated on the Boundary: 7722 = A 7703 ->• -00,7702 -> - 0 0 , 7703 + 7723 = A i and 7702 + 106 2 7.18 Results for L i u Transition Model with Drop-out Model ID2 Evaluated on the Boundary: 7703,7702 - > - 0 0 , Ai = 7703 + f]i and A 2 = 7702 + V2 • 106 7.19 Results for L i u Transition Model with Drop-out Model ID5 Evaluated on the Boundary: 7703,7702 - > - 0 0 , ix Ai = 7703 + 772 and A 2 = 7702 + ?72 • 107 7.20 Goodness-of-fit Statistics for L i u Transition Model with Drop-out Models ID1, ID2, ID3 and ID5 107 7.21 The Observed and Expected Cell Counts for the L i u Transition Model with Drop-Out Model ID5 ("*" denotes missing) Ill 8.1 Estimated Chance of Exacerbations Based on Baker's Selection Model 115 8.2 Estimated Chances of Exacerbations Based on the L i u et al. Transition Model 8.3 116 Estimated P r ( Y / = l,Y * = 1) and P r ^ * = 1,Y * = 1,Y * = 1) by t 2 3 Treatment Groups 118 B.l Results for Model ID1 130 B.2 Results for Model ID2 131 B.3 Results for Model ID3 132 B.4 Results for Model ID4 133 B.5 Results for Model ID5 134 B.6 Results for Model ID6 B.7 Results for Model RD1 B.8 Results for Model RD2 136 B.9 Results for Model RD3 136 , 135 ." 135 B.10 Results for Model C R D 1 137 B. l l Results for Model C R D 2 137 C. l Results for Drop-out Model: T R T * L U R 139 C.2 Results for Drop-out Model: T R T + L O R + L U R 140 C.3 Results for Drop-out Model: L O R * L U R 141 C. 4 Results for Drop-out Model: T R T + L U R 142 D. l Results for Case 1 in Table 7.14 Evaluated at the Boundary: 7703 —> - 0 0 , 7702 - 0 0 , 7703 + 772 = Ax, 7702 + 772 = A x 2 144 D.2 Results for Case 2 in Table 7.14 Evaluated at the Boundary: -co, 77 2 -> - 0 0 , 770.3 + 772 = A i , 770.2 + *72 = A 0 145 2 D.3 Results for Case 3 in Table 7.14 Evaluated at the Boundary: -00, 7702 - > - 0 0 , 7703 + 772 = A i , 7702 + 772 = A 7702 -> - 0 0 , 7703 + 772 = A i , 7702 + 772 = A 146 7702 - 0 0 , 7703 + 772 = A i , 7702 + 772 = A 7703 —> 147 2 D.5 Results for Case 5 in Table 7.14 Evaluated at the Boundary: -00, 7703 —>• 2 D.4 Results for Case 4 in Table 7.14 Evaluated at the Boundary: -00, 7703 -> 7703 -> 148 2 D.6 Results for Case 5 in Table 7.14 Evaluated at the Boundary (364 patients): 7703 -> - 0 0 , 770:2 ->• - 0 0 , 7703 + 772 = Ai, 7702 + 772 = A 2 . . . 149 D. 7 Results for Model ID2 Evaluated at the Boundary (364 patients): 7703 - » • - 0 0 , 7702 -> - 0 0 , 7703 + 772 = A i , 7702 + 772 = A 2 150 E. l Results for Liu Transition Model with Drop-out Model ID1 152 E.2 Results for Liu Transition Model with Drop-out Model ID2 152 E.3 Results for Liu Transition Model with Drop-out Model ID3 153 E.4 Results for Liu Transition Model with Drop-out Model ID5 153 E.5 Results for Liu Transition Model with Random Drop-out (RD) . . . 154 E.6 Results for Liu Transition Model with Drop-out Completely At Random (CRD) 154 xi List of Figures 2.1 Histogram for Length on Study (3-month bins) 2.2 Kaplan-Meier Survival Curves for Time on Study: Over 3-year Treatment Period 2.3 11 Proportion of Patients Experiencing Exacerbations Over Time by Treatment Group Based on Dichotomous Annual Data 2.4 10 16 Boxplots of Age, Duration of MS and EDSS at Baseline by Treatment Group 19 2.5 Histogram of B O D and Boxplots of B O D by Treatment Group ( n = 364) 20 7.1 A Two-Dimensional Profile Log-likelihood Surface for Model ID1 . . 73 7.2 A Two-Dimensional Profile Log-likelihood Surface for Model ID1 with Boundary Constraint 7703 —> — 00 and 7703 + 7723 = A xii 76 Acknowledgements First and foremost, I would like to thank my supervisor, Professor John Petkau, for his patience, constant support and invaluable guidance throughout the development of this manuscript. I could not have chosen a better person to work with. I am also grateful to Professor Harry Joe for agreeing to review my work on such short notice. A huge thanks also to Professor Nancy Heckman for her constant encouragement, especially during my first semester in UBC. I would also like to express my gratitude to Ms. Christine Graham for her help. In addition, I would like to thank the entire UBC Statistics Department for making my stay enjoyable. Most importantly, I am indebted to Ryan Woods, as without his love, encouragement and belief in me, I would not have completed my Master's degree. L E E SHEAN E R The University of British Columbia March 2001 xiii my beloved late mother, Goh Seok l m xiv Chapter 1 Introduction 1.1 Background of this Thesis The defining characteristic of a longitudinal study is a sample design which specifies repeated observations on the same individual (or experimental unit). However, failure to obtain a full set of observations on a given individual (or unit), resulting in incomplete data and/or unbalanced designs, is a common problem in longitudinal studies. The form of missingness in longitudinal studies is typically drop-outs, in which sequences of measurements on some individuals terminate prematurely. This drop-out phenomenon is reflected in a data set collected over a 3-year period in a multicenter multiple sclerosis (MS) clinical trial sponsored by Berlex Laboratories of Richmond, California. The work presented in this thesis is motivated by this data set. A detailed description of the clinical trial and the data set to be analyzed can be found in Chapter 2. Multiple sclerosis is a serious disease of the central nervous system (CNS), the nerves that comprise the brain and spinal cord. The term "multiple sclerosis" refers to multiple areas of patchy scarring, or plaques, that result from the destruction of myelin. Myelin is a white substance which forms a sheath around the spinal cord. When the myelin sheath is destroyed, signals transmitted throughout the CNS are 1 disrupted which leads to the occurrence of an acute attack, or exacerbation. During these exacerbations patients can suffer from a variety of symptons such as blurred vision, a sensation of numbness or loss of control of the movements in parts of the body. To date, the cause of MS is unknown and no cure exists. A number of treatments examined over the past decade have reduced rates of acute exacerbations and slowed progression in disability. In fact, the Berlex trial, where the treatment investigated was Interferon /3-lb, was the first to demonstrate beneficial effects of a treatment for MS patients. Patients withdrew from the Berlex trial due to reasons such as lack of efficacy, toxicities in excess of prespecified toxicity levels, or other side effects to the treatment. Nevertheless the overall drop-out rate did not exceed that anticipated at the trial's inception. The intent-to-treat analyses of the trial data were per- formed under the assumption that the drop-out occurred completely at random, as is customary in clinical trials. Methods have been developed for explicitly modelling non-response (not only restricted to drop-outs) under the more general assumption that the non-response may not have occurred completely at random. Our main objective is to investigate the sensitivity of the conclusions concerning the treat. ment effects to different assumptions about the nature of the drop-out mechanisms. Diggle and Kenward's (1994) classification of drop-out mechanisms, modified from Rubin (1976) and Little and Rubin (1987), is described in Chapter 3. Comparison of different models also allows us to study the nature of the drop-out mechanism in this data set. 1.2 Methods of Analyses Likelihood-based methods are commonly used for incomplete data, including for the analysis of longitudinal data with drop-outs. According to the Diggle and Kenward (1994) terminology, drop-out mechanisms can be classified as completely random drop-out (CRD), random drop-out (RD) or informative drop-out (ID). For a 2 C R D mechanism, drop-out is independent of the outcome (or measurement) process; for a R D mechanism, drop-out is independent of the unobserved outcomes, but depends on the observed outcomes. For an ID mechanism, drop-out depends on both the observed and unobserved outcomes. Likelihood-based methods yield valid results in the presence of C R D or R D provided the model used for the measurement process is valid, and the observed information matrix is used rather than the expected information matrix. If, however, the drop-out mechanism is ID, modelling the drop-out process is necessary to permit valid inferences; see Laird (1988). Modelling different drop-out mechanisms can provide insight into the nature of the withdrawal process. It can also be used to investigate the sensitiv- ity of inferences to the underlying assumptions. In the past decades, researchers have proposed a number of methods for quantitative longitudinal data (usually normally distributed) and categorical longitudinal data subject to non-random dropout. Laird (1988) provided an excellent discussion of how the drop-out process can affect the inferences about both continuous and categorical measurement processes. For continuous longitudinal data, Wu and Carroll (1988) considered ID in a random effects model, with the data for each experimental unit following a linear time trend whose intercept and slope vary between individuals according to a bivariate Gaussian distribution. Their likelihood-based method permits the comparison of non-ID and ID drop-out mechanisms. Schluchter (1992) outlined a new approach based on a log-normal survival model when the primary outcome is the rate of change in a continuous variable subject to informative censoring. More recently, Diggle and Kenward (1994) proposed a general model-based approach for analyzing continuous longitudinal data that combines a multivariate linear model for the response with a logistic regression model for the drop-out process. This is the first paper to develop a modelling strategy that explicitly accommodates C R D and R D as special cases within an ID model. The issue of how to deal with ID in categorical longitudinal data is not yet 3 resolved. Further, potential technical difficulties may arise in the likelihood-based methods for correlated categorical data due to the discreteness of the responses. Baker and Laird (1988) developed a log-linear model for categorical response subject to non-ignorable non-response in a sample survey setting and drew attention to the existence of boundary solutions. A number of authors have focused their attention on the multivariate binary data case to better understand some of the potential difficulties for correlated categorical data. Both Baker (1995) and Fitzmaurice et al. (1996) used a multivariate binary model where the marginal probabilities for the responses are specified as logistic regressions. However, these authors modelled the associations among the responses differently. These models for the outcomes were combined with logistic models for ignorable and non-ignorable drop-out mechanisms to analyze multivariate binary data. Both papers also highlighted the issue of identifiability of these models. Baker (1995) provided outlines of the proof of model identifiability for certain models he considered. Fitzmaurice et al. (1996) gave some suggestions on how to examine the identifiability of non-ignorable dropout models. More recently, Ten Have et al. (1998) presented mixed effects models for longitudinal binary responses with informative drop-out analogous to the Wu and Carroll (1988) models for longitudinal continuous data. Liu et al. (1999) adapted the method proposed by Diggle and Kenward (1994) for the analysis of a binary longitudinal outcome. Most of the likelihood-based models mentioned are formulated within the selection modelling framework (Little and Rubin, 1987). A selection model factors the joint distribution of the measurement and response processes into the marginal measurement distribution and the response distribution, conditional on the measurements. Molenberghs et al. (1999) discuss the strengths and limitations of selection models for non-random missingness in the categorical data setting. There are other ways to specify the joint distribution. For categorical responses, a loglinear approach incorporates the measurement and response processes into a sin- 4 gle log-linear model. A time-ordered approach factors the joint probability into a product of conditional probabilities ordered in time and a pattern-mixture model (Little, 1993) factors the joint distribution into the marginal response distribution and the measurement distribution, given the response distribution. The latter two approaches can be applied to both continuous and categorical repeated measurements data. Ekholm (1998) re-analyzed the children's obesity data set considered by Baker (1995) using a pattern-mixture model. Michiels et al. (1999) studied similarities and differences of modelling incomplete data within the selection and pattern-mixture settings assuming a missing at random mechanism. Pseudo-likelihood and non-parametric approaches have also been proposed for carrying out analyses under different types of drop-out mechanisms. Using Dale's model (1986) for ordinal categorical longitudinal data, Kenward et al. (1994) demonstrated that, in the presence of R D , the generalized estimating equations (GEEs) approach proposed by Liang and Zeger (1986) may give misleading results. Robins et al. (1995) showed that appropriately weighted G E E s overcome this problem, but not in the presence of ID. More recently, Sun and Song (2000) proposed a non-parametric approach for analyzing the data from a clinical trial of adult schizophrenics with informative censoring. 1.3 Outline of this Thesis Our attention will be on multivariate binary data. This special form of the data allows us to focus on the aforementioned issues that arise mainly in correlated categorical data. We are also interested in studying the nature of the drop-out process in our data. For this purpose, we choose to work with models within the selection modelling framework. The remainder of this thesis is outlined as follows: Chapter 2 presents the description of the Berlex trial and the binary responses which will comprise the data set to be analyzed. Chapters 4 and 5 discuss Baker's selection model and the L i u t 5 transition model respectively. These models can be used to examine various types of drop-out mechanisms in our data. The definitions of the drop-out mechanisms is provided in Chapter 3. Non-ignorable (or informative) drop-out (or non-response) models are generally harder to implement due to potential analytical problems such as model identifiability issues. Chapter 6 focuses on the issue of identifiability in models for incomplete binary responses. Proofs of identifiability of some of our models are also included in the chapter. The detailed results of our analyses are reported in Chapter 7. We conclude the thesis with some general discussion in Chapter 8; this includes comments on the two models and suggestions of other possible methods for analyzing the data. 6 Chapter 2 Data Description 2.1 Description of the Berlex Clinical Trial The Berlex clinical trial was a phase III trial of the effect of Interferon /3-lb on relapsing-remitting multiple sclerosis (MS) patients. The primary outcome measure was the rate of exacerbations. This was a multicenter, randomized, double-blind, placebo-controlled trial with three parallel treatment groups. The study was originally planned with a 2-year treatment period; the trial was later extended to 3 years (because by the end of the second year, some patients had been on the study for almost three years due to different starting dates). The study was carried out in a double-blind fashion for the full three years. The data from the first 2 years of the study established that the Interferon /3-lb treatment groups had decreased exacerbation rates and increased proportions of patients remaining exacerbation-free. These beneficial results were also found in the 3 year data. This was the first trial to unequivocally identify an effective treatment for relapsing-remitting MS. Interferon /3-lb has emerged as a therapeutic option in MS and has been hailed as a major advance in the management of this disorder. This trial consisted of 372 patients from 11 centers in the United States and Canada on three parallel treatment groups: placebo (PL), low dose (LD) and high 7 dose (HD). The dosage for L D and H D were 1.6 and 8.0 million international units (MIU) respectively. A l l patients were between the ages of 18 and 50 years, had been diagnosed with MS at least 1 year prior to entry to the study, had Krutzke Expanded Disability Status Scale (EDSS) scores of 5.5 or less, and had experienced at least 2 exacerbations in the previous 2 years. Moreover, all had been clinically stable for at least 30 days prior to entry and had received no medications to speed up the recovery from relapse such as A C T H (adrenocorticotrophic hormone) or prednisone during this period. Patients were randomized to the three treatment groups within each center and divided almost evenly within each center. A l l patients were blinded to the treatment assignments. O f these 372 patients, 123 received P L , 125 received L D and 124 received H D of Interferon /3-lb by injection every other day. Two neurologists were appointed at each center: one who performed the periodic examinations was not aware of the drug side effects, and another who knew about the side effects and injection reactions was responsible for reviewing laboratory findings for toxicity and for overall patient care. Patients were scheduled to be evaluated every 12 weeks except for the first few months of the study, where evaluations were more frequent. In addition, visits were made when symptoms occurred suggesting the possibility of an MS exacerbation. A Scripps Neurological Rating Scale (NRS) score and a Kurtzke E D S S score were determined in each evaluation. For all patients in the study, the beginning and end dates of all exacerbations as well as the EDSS scores obtained at each visit were recorded. Besides these clinical outcomes, each patient also had a baseline cranial magnetic resonance imaging (MRI) and this was repeated annually. The patients at one of the centers (the University of British Columbia) had cranial MRIs repeated at 6-week intervals for the first 2 years. 8 2.1.1 Drop-out Rate in this Clinical Trial Since some beneficial results of Interferon /3-lb were found after 3 years of study, patients who remained in the study were offered the high dose treatment for another 2 years. Thus, the entire study continued for over 5 years, but many patients dropped out during this period. Figure 2.1 shows a roughly constant rate of dropout during the first three and a half years, except for a large number of drop-outs at the end of the second year (the original intended end of the study). The plot also indicates the drop-out rate increased dramatically after the end of the 3-year treatment period. Because of the potential difficulty in interpreting the 5 year data (e.g. how should patients who switched from one treatment to another be treated in the analysis and how should the results obtained be interpreted), we employ the 3-year treatment period data to perform a variety of analyses in this thesis. We can represent the information shown in Figure 2.1 in another fashion. Figure 2.2 displays Kaplan-Meier survival curves describing the proportion of patients remaining on study by treatment group; the dash, solid and dotted vertical lines indicate the end of the 1-year, 2-year and 3-year periods respectively. The most drop-outs over the 3-year period occurred in the low dose group (approximately 40%). Roughly 20% of patients withdrew from the trial during the first 2 years in all three groups. A number of patients in each treatment arm dropped out around the end of the 2-year period, but the proportion remaining for most of the third year of the study is roughly 70% in both the P L and H D groups and roughly 60% in the L D group. Table 2.1 summarizes the numbers of patients who dropped-out by the end of the first, second and third year of the study. More details concerning the clinical trial can be found in the published reports of the I F N B Multiple Sclerosis Study Group [27, 35, 36]. 9 Figure 2.1: Histogram for Length on Study (3-month bins) End of 2-year Period End of 3-year Period CD <= •5 o- 500 1000 1500 Length on Study in Months Table 2.1: Cumulative Number of Drop-outs After 1, 2 and 3 Years on Study After 1-Year After 2-Year After 3-Year Group Number Proportion Number Proportion Number Proportion PL 13 11% 27 22% 41 33% LD 11 9% 30 24% 49 39% HD 17 14% 29 23% 35 28% 10 Figure 2.2: Kaplan-Meier Survival Curves for Time on Study: Over 3-year Treatment Period 2.2 Description of the D a t a The main objective of this thesis is to explore models for longitudinal binary responses incorporating different types of drop-out mechanisms in the context of the Berlex trial. We consider the exacerbation variable as the response variable of interest in our analysis. We choose to represent these data in binary form on an annual basis (whether exacerbations occurred in each 1-year interval) to allow a specific focus on models for binary responses. In other words, the data for each patient will be represented by three binary responses indicating whether they experienced any exacerbations during the 1-year intervals. One of the reasons for proceeding in this way, as opposed to refining the time intervals to 6-month intervals say, is to reduce the number of possible different derived sequences of the binary responses as well as the number of drop-out patterns. This allows a focus on the key ideas for modelling such data. This will become clearer in later chapters. The rest of this section is structured as follows. In the next subsection, we describe the annual drop-out patterns. We then discuss how the binary responses are derived. We conclude the section with a brief description of the baseline covariates to be included in our analyses. 2.2.1 Drop-out Patterns The data described in the previous section involve a total of 372 patients. Each patient's termination date from the study was recorded. To derive our annual data set, the data on patients who dropped out are handled as follows: • Scenario 1 If the patient's termination date was prior to 365 days on study, then we will treat these patients as if they dropped out at the beginning of the study. In other words, these patients have no outcomes in our annual data set. 12 Table 2.2: Number of Patients in the-4 Drop-out Cases (x=present,o=absent) 0 Year 1 Year 2 Year 3 Number of Patients 1 o o o 41 2 X 0 o X X o X X X 45 39 247 3 4 • Scenario 2 If the patient's termination date was after 365 days but prior to 730 days, then we will treat these patients as if they dropped out at the end of the first year of the study. These patients have one outcome in our annual data set. • Scenario 3 If the patient's termination date was after 730 days and prior to 1095 days, then we will treat these patients as if they dropped out at the end of the second year of the study. That is, these patients are missing only the third year outcome in our annual data set. • Scenario 4 If the patient's termination date exceeded 1095 days, then we will treat these patients as if they completed the 3-year study and thus all three annual outcomes were observed. Table 2.2 summarizes the total number of patients according to the four scenarios of available annual outcomes over the 3 year period ("x" denotes present and "o" denotes absent). Table 2.3 displays the breakdown of the 372 patients in our annual data set by treatment groups and gender according to the total number of patients entering the study, and dropping-out at the beginning, the end of the first year and the end of the second year of the study. Patients were quite evenly distributed across the 3 treatment arms at the beginning of the study, as were the patients who dropped 13 Table 2.3: Summary of the Our Annual Data By Treatment Groups By Gender Number of Patients PL LD HD Males Females Entering the Study 123 124 113 259 Drop-out At Beginning Drop-out At E n d of Year 1 Drop-out At E n d of Year 2 13 14 14 125 11 19 19 17 12 6 16 9 8 25 36 31 out at the beginning of the study and at the end of year 1. However, fewer patients in the H D group dropped-out at the end of year 2 than in the P L and L D groups. In summary, the L D group has the highest drop-out rate, followed by the P L group, and both rates increase slightly over time. As expected, the drop-out rate in the H D group is the lowest and it decreases over time. Table 2.3 also shows the drop-out rates for females and males are fairly consistent over time, although the drop-out rates for females are a bit higher than for males. In the next two sections, we provide a more detailed description of the binary outcome variable and baseline covariates of interest. A l l the corresponding descriptive statistics presented are based on our annual data. 2.2.2 Binary Outcome Variables As mentioned earlier in the chapter, the start and end dates of any exacerbations patients experienced during the study were recorded. For our purposes, we do not use the end dates even though they could contain valuable information. A l l exacerbations are attributed to the annual period in which they began. Recall we divided the time period of the study into three 1-year intervals. Since these intervals are quite wide, some patients experienced multiple exacerbations within these intervals. The number of exacerbations experienced by patients within these annual intervals ranges from 0 to 6; the frequency of these counts by yearly interval is summarized in Table 2.4. Most patients experienced either no exacerbations or a small number of 14 Table 2.4: Frequency Table of the Exacerbation Counts for Patients with At Least One Outcome Number of Exacerbations 0 1 2 3 4 5 6 Interval Group Year 1 ALL PL LD HD 121 32 39 ALL PL LD 122 Year 2 Year 3 50 Number of Patients 103 36 35 32 63 20 26 27 14 7 13 8 4 2 0 1 2 0 2 331 110 114 17 6 1 1 0 107 86 28 28 53 16 22 15 8 2 4 2 1 0 1 30 15 5 0 5 3 1 1 286 96 95 95 12 4 0 0 1 0 0 0 0 0 0 0 1 HD 39 39 44 ALL 122 74 38 PL LD HD 37 36 49 23 26 18 8 12 25 5 3 2 0 0 247 82 76 89 exacerbations; only a few patients had 4 or more exacerbations within a year. Based on this information and for simplicity of analysis, it seems reasonable to dichotomize these data as no exacerbation or at least 1 exacerbation experienced. Clearly there is some loss of information associated with dichotomizing these data. One way to retain the information is to treat the counts of the total number of exacerbations as if they are Poisson random variables and perform analyses based on the counts. However, we will not explore such analyses in this thesis. Figure 2.3 shows the proportion of patients experiencing exacerbations over time by treatment group based on these dichotomized annual data. In general, the proportion of patients experiencing exacerbations decreased over the 1-year periods in all groups. Further, the H D group has the lowest proportions among the 3 treatment arms throughout the study. The proportion of patients experiencing exacerbations is slightly higher in the P L group than in the L D group. This plot also suggests a dose-response relationship in these data. 15 Figure 2.3: Proportion of Patients Experiencing Exacerbations Over Time by Treatment Group Based on Dichotomous Annual Data Year 16 2.2.3 Baseline Covariates We are primarily interested in the assessment of the treatment effects on the binary outcome variables described in the previous section, but patterns in the data over time and the effects of several baseline covariates are also of interest. The baseline covariates we considered are: • gender, • age, • duration of MS, • Kurtzke Expanded Disability Status Score (EDSS), and • burden of disease (BOD). In general, more females than males suffer from M S . This phenomenon is reflected in this trial; as shown in Table 2.5, the female-to-male ratios are roughly 2.5, 2.1, and 2.3 in the P L , L D and H D groups respectively. Figure 2.4 shows the boxplots of age, duration of MS and EDSS at baseline by treatment group. The ages range between 18 and 50 years. The median age at baseline in the H D group is slightly smaller than in the other groups, but the distribution of the ages is quite similar for the 3 treatment groups. The boxplots also indicate that about 50% of the patients had ages between 30 and 40 years in each treatment group. The duration of MS ranges between 1 and 31 years and the median is slightly higher in the H D group. The boxplots for baseline EDSS indicate a fairly balanced distribution across the three groups, with scores ranging from 0 to 5.5 in each group. There are two distinct forms of magnetic resonance imaging (MRI) scans of interest in MS studies: Tl-weighted scans and T2-weighted scans. A Tl-weighted scan uses a small injection of the chemical gadolinium into the patient's bloodstream. The presence of gadolinium will enhance the appearance of active lesions (areas of inflammation on the blood/brain barrier) on the brain stem, and facilitate their 17 Table 2.5: Number of Female and Male Patients in Each Treatment Group Female Male Total Treatment Group PL L D HD 88 85 86. 35 40 38 123 125 124 detection. A T2-weighted scan provides clearer definition of the actual size and shape of each lesion without any gadolinium injection into the bloodstream, which usually blurs the border of the lesions. The MRI measure of interest in this thesis, known as burden of disease (BOD), is a measure of the total volume of all lesions on the T2-weighted scan. In our data set, there are 8 patients who did not have a BOD measurement at baseline; 3 from the P L group, 4 from the L D group, and 1 from the HD group. Excluding these 8 patients, the histogram of BOD at baseline and the boxplots of BOD at baseline by treatment group are shown in Figure 2.5. The distribution of BOD is highly skewed to the right. There are only 2 patients who did not have any lesions at baseline (BOD — 0), but there are 5 patients with BOD greater than 10,000 (mm ): 2 belong to the P L group and 3 belong to the LD 2 group. This is also reflected in the boxplots in Figure 2.5. Excluding the 3 patients in the LD group who had the largest BOD readings, the general distribution of the BOD measurements is quite similar in each treatment arm. 2.3 Questions of Interest Having introduced the annual data set to be analyzed, we now describe the study questions we plan to address in this thesis. Recall that the main focus of this thesis is to explore models for analyzing repeated binary data incorporating different drop-out mechanisms. Although the drop-out rate in our annual data is moderate, we would like to investigate the 18 Figure 2.4: Boxplots of Age, Duration of MS and EDSS at Baseline by Treatment Group Age at BaseBne by Treatment Duration of MS by Treatment 5 19 EDSS at Baseline by Treatment Figure 2.5: Histogram of B O D and Boxplots of B O D by Treatment Group (n = 364) Histogram of B O D at Baseline o co o CD o cvj 5 10 BOD (/1000) B O D at Baseline by Treatment 8 o Q O m o , _ m 20 15 most appropriate form of model for the drop-out process; in particular, to explore whether there is an indication of informative drop-out. It is also of interest to assess the sensitivity of inferences concerning the treatment effects (primarily) to the form of the models for the drop-out mechanism, and to explore the importance of baseline covariates for our annual data. Chapter 3 provides a discussion of different drop-out mechanisms. We describe general methodology for analyzing incomplete binary data in Chapters 4 and 5 . Chapter 6 sheds some light on potential identifiability problems in such models. Chapter 7 contains all the results from the analyses we performed and we conclude the thesis with some discussion. 21 Chapter 3 Classification of Missing Values in Longitudinal Data 3.1 Introduction Longitudinal studies are usually characterized by collecting a set of measurements on an individual unit at prespecified points in time; in many cases (typically in clinical trials), the set of prespecified points in time are the same for all units. Missing values arise whenever one or more of the intended measurements from units within the study are incomplete. Such missing data are a common problem in longitudinal studies, particularly when the experimental units are human subjects and collecting data involves a visit to a hospital or clinic, or the time between intended measurements is lengthy. It is important to distinguish between unbalanced data and missing values. Unbalanced data result when the set of times of intended measurements is not common to all units; for example, if one chose in advance to take measurements every half hour on one-half of the subjects and every hour on the other half. Such unbalanced data could also be described as incomplete but there are no missing values from the viewpoint of the design of data collection. Missing data also arise in 22 unbalanced data; however, there are deeper conceptual issues as to why the values are missing, and more specifically whether the missingness is related to the questions posed by the study. Little and Rubin (1987) have provided a useful classification of missing value mechanisms. Let Y* denote the complete set of measurements for one unit which would have been obtained if there were no missing values. Partition this set into Y* = (Y(°),Y( )) with Y(°) denoting the measurements actually obtained and m Y( ) the measurements which would have been available if they had not been m missing, for whatever reason or cause. Let R denote a set of indicator random variables, denoting which elements of Y* fall into Y(°) and which into Y^™). We can 1 then specify a probability model for the missing value mechanism as the probability distribution of R conditional on Y* = (Y^°), Y^" )). In the terminology used by 1 Little and Rubin, the missing value mechanism is classified as: 1. completely random if R is independent of both Y(°) and Y^ ^; 2. random if R is independent of Y^ ); m m 3. informative if R is dependent on Y ^ \ m We will abuse the notation / to denote a probability density (or mass) function throughout this thesis; the function being referred to will be clear from the context. For likelihood-based inference, the important distinction is between and random informative missing values. To see this, f{y^°\ y^ \ r), the joint probability m density function (pdf) of /(y (o) ,y (Y<°>, Y( ),R), can be factored as m (m) ,r) /(y (o) ,y (m) )/(r |y ( o ) ,y ( m ) ). (3.1) For a likelihood-based analysis, we need the joint pdf of the observed random variables, (Y(°),R), which can be obtained by integrating (3.1) over all possible values for the unobserved random variables (3.2) 23 If the missing value mechanism is random, f(r \ y^,y^) is independent of y( ) and (3.2) becomes m /(r I y )/(y ). (o) (o) (3.3) Taking logarithms in (3.3), the log-likelihood function is L = log/(r|y( ))+log/(yW), 0 (3.4) which can be maximized by separate maximization of the two terms on the righthand side provided the parameters appearing in /(r | y(°)) and in f{y^) are disjoint. Since the first term contains no information about the distribution of Y ^ , we can ignore it for the purpose of making inferences about Y(°). Because of the above result, both completely random and random missing value mechanisms are sometimes referred to as ignorable. On the other hand, informative missing value mechanisms are referred to as non-ignorable because such a missing value mechanism cannot be ignored when making inferences about Y(°). 3.2 Types of Drop-outs We have distinguished between unbalanced data and missing values. Now let us focus on different types of missing values. Missing values can occur either intermittently or as drop-outs. Suppose we intend to obtain a sequence of n measurements, say Y i , Y , . . . , Y , on a particular unit. We say that missing values occur as drop2 n outs if whenever measurement Yj is missing, so are the measurements, Y^ for all k > j; otherwise the missing values are intermittent. In this thesis, we are particularly interested in studying drop-out mechanisms. Drop-outs are a common phenomenon in longitudinal studies. They typically arise not as a result of censoring applied to the measurements on the experimental unit, but because some units prematurely terminate their participation in the study. 24 A unit's withdrawal may be for reasons directly or indirectly connected to the measurement process. Thus, it is of interest to investigate whether the drop-out-process is related to the measurement process. Following the Little and Rubin (1987) discussion of the classification of missing value mechanisms, Diggle and Kenward (1994) modified the above definitions slightly to describe drop-out processes as: (a) Completely Random Drop-out (CRD): if the drop-out mechanism is indepen- dent of the measurement process; (b) Random Drop-out (RD): if the drop-out mechanism is independent of the unobserved measurements, but depends on the observed measurements; (c) Informative Drop-out (ID): if the drop-out mechanism depends on both the observed and unobserved measurements. Both C R D and R D are referred to as ignorable drop-outs, while ID is referred to as non-ignorable drop-out. In next two chapters, we give an overview of the selection modelling approach for longitudinal binary data subject to non-ignorable non-response. The basic idea of a selection model is to factor the joint distribution of the measurement variables and the non-response indicator variables, / ( Y * , R ) , into / ( R | Y*)/(Y*), where /(Y*) is known as the outcome model and / ( R | Y*) is known as the drop-out model. The only distinction between the next two chapters is in the model for the outcome (or measurement) process. 25 Chapter 4 Selection Model 4.1 Baker's Selection Model for Binary Longitudinal Data with Informative Non-response Diggle and Kenward (1994) provided a general methodology for dealing with continuous responses subject to informative, or non-ignorable, drop-outs in a longitudinal study. Baker (1995) provided a discussion of a related model that accounts for non-ignorable non-response. The methodology is connected to that presented by Diggle and JKenward, however Baker's model is for repeated binary data and the non-response is allowed to occur in various patterns, not only as drop-outs. For simplicity, we limit our discussion to repeated binary data collected at 3 time points, as this coincides with the structure of our data set. Our model is a simplified version of Baker's model as we are only interested in monotonic nonresponse patterns, i.e. drop-outs. We first introduce the concepts of incomplete (observed) and complete data. Let t index the time points where measurements are intended to be taken. In this particular context, t represents the three successive 1-year period measurements that are to be taken, coded as t = 1,2,3. Let X denote a vector of covariates at time t t, and denote X = (X'^X^X^). The vector of random variables for the complete .26 i i I (possibly unobserved) data is (^T> *2*> > -^1) R3, X ) , where Y * is the binary outcome variable at time t which takes on values 0 or 1, and t Rt is an indicator variable of non-response at time t with sample space {a,p} where "a" denotes absent and "p" denotes present. The vector of random variables for the incomplete (observed) data is (Yi,Y y ,X), 2 ) 3 where Y has sample space {0,1, a}. The complete and incomplete random variables t are related as follows: ' Y * i£Rt = p, t a if Rt = a. There are several approaches for modelling the joint distribution of the complete data, Pr(Y * = y{,Y2* = y *,Y * = y ^ 1 2 3 - ruR2 = r , i ? = r | X ) . Baker 2 3 3 chose to use a selection model in which the joint distribution is factored into the probability of the outcomes multiplied by the probability of the non-response indicators, given the outcomes; that is, Pr(i2i = n , R 2 = r 2 > R 3 = r | Y{ = y{,Y * = y* ,Y * = y* ,X) 2 3 2 3 3 x Pr(Y *=yI,Y *= /2,Y * = yS | X ) . 1 2 l 3 Now, denote the outcome model as P r ( Y ; = y{,Y2* = y*2,Y3* = y*3 | X ) = /*(y*,y *,y* | x;0), 2 (4.1) where 6 is a vector of parameters. Also, denote the non-response model as Pr(i2i =ri,R2 = r , R = r | Y{ = y{,Y * = y* ,Y * = y* ,X) 2 3 3 2 2 3 = q(ri,r ,rz\yl,y ',yZ,x;r)), 2 27 2 3 (4.2) Table 4.1: A l l Possible Patterns for Incomplete Data for the Case where 3 Observations were Intended for Every Unit: x = observed, o = missing. Pattern yi V2 V3 1 2 3 4 x X X X X o X o X o X X 5 6 7 8 X 0 o o X o o o X o o o where r\ is a vector of parameters. This construction assumes that the parameters of the outcome and non-response models are distinct (Diggle and Kenward, 1994). This relates back to the idea of ignorable and non-ignorable drop-out mechanisms discussed in Chapter 3 (see p. 18). Under R D and C R D , inferences based on the observed data are valid even though the' drop-out mechanism is ignored, but this is not true for an informative drop-out mechanism. Table 4.1 displays all the possible realizations of the incomplete data in this particular scenario, where "x" denotes the measurement is observed and "o" denotes the measurement is unobserved. Using the outcome and non-response models, we can write down the probability of these 8 realizations of incomplete data as follows: f(yl,yhyt I = f*{y{^yl I «;0)<?(P,P,P | y^y^yh^v), (4.3) i f(Vi,V2,d\x;0,ri) = ^2[f*{y* y2,yt \ x;6) x q{p,p,a \ yl,y^,y^x;r}) u , 3/3=0- (4.4) l / ( y i » o , i / 3 | as;0,T|) = ^2[f*{y{,y2,yt \ x-,6) x q{p,&,p \ yl,y^,y^,x;rj)j, 2/2=0 (4.5) 28 f(a,y ,yl \ x;0,r}) 2 = ^ \ x;G) x <j(a,p,p | 1/1,1/2.3/3, a:; 17) , [/"'(yl,y ,y 2 3 (4.6) 1 f{y{,a,a | X;0,TJ) l = [/*(2/i,y2,2/3 I s.fl) ^ <z(p,a,a | y*,y ,yZ,x;r)) x , 2 y ' = 0 2/*=0 (4-7) l f{a,y ,a \ x; 9,r)) 2 = l ^ ^ [/*(yt>y2^3 I 2/1=0 =o ?(a,P,a | y*, y ,y* , x; rj) , x 2 3 (4.8) l /(a,a,2/3 | x;0,7/) = l ^ ^ [/*(yi%2/2,2/3 I 2/1=0 2/5=0 a;;0) x g(a,a,p | y*,^,y^,a;;T7) , (4.9) l f(a,a,a \ x;6,n) l l ^ ^ [/*(yI,y2>2/3 | x ; 0 ) x g(a,a,a | y i ^ ^ a : ; * / ) ] yj=0 2/^=0 y|=0 = (4.10) Baker specified the outcome model, / * (yj, j / ^ 2/3 I 'i x m terms of a marginal model which models the marginal probabilities as functions of covariates, and an association model which models the temporal associations, using the idea suggested by Ekholm (1991, 1992). The non-response model, 17(7-1,r2,r3 | y\,y ,7/3,x;77), was 2 modelled by employing a general time-order causal model. With the assumption that drop-out does not depend on future events, the time-ordered causal model for three time points has the form q{n,r ,r 2 3 I yl,y%,y$,x;ri) = Pr(# = r \R =r R 3 3 l x Pr(i? = r 2 u 2 = r , Y{ = yj, Y * = y* ,Y * = y*, x) 2 2 2 2 3 I Ri = n,Y{ = yl,Y * = y* ,x) 2 x P r ( i ? ! = n | y * = j/J,x). 1 2 (4.11) To complete the specification of the non-response model, Baker modelled each of these conditional probabilities as a logistic regression. Under non-ignorable 29 non-response, these logistic regressions involve the unobserved outcomes as well as the observed outcomes and the covariates. 4.2 Selection Model for Binary Longitudinal Data with Informative Drop-out Our main focus is to explore models for incomplete binary responses subject to informative drop-out for our annual data set as described in Section 2.2. The approach sketched in the previous section can be modified to serve our purpose. We consider the case where drop-out occurs either at the first, second or third time point. Our data are then limited to 4 of the 8 possible patterns listed in Table 4.1: patterns 1, 2, 5 and 8. Patterns 2, 5 and 8 form a monotone pattern of non-response, and are also known as drop-outs. Following Baker, the probabilities for these 4 incomplete data patterns are given by equations (4.3), (4.4), (4.7) and (4.10). In other words, our model is a simplified version of Baker's more general selection model. In the next few subsections, we specify particular forms for the outcome model, f*(yl,V2,V3 \ x;0), and the drop-out model, q(ri,r ,r 2 3 \ yl,y^,y^,x;ri). The likelihood function is then assembled according to these models. 4.2.1 Outcome Model Baker's outcome model y ,Vz I x;0) is specified in terms of two models: a 2 marginal model (model for the univariate marginal probability) and an association model (model for the multivariate probability). There are several approaches to constructing marginal and association models with binary longitudinal data. Baker used the parameterization introduced by Ekholm (1991, 1992) which expresses f*{y*, 2/2'^3 I ;9) x a s a linear combination of marginal and association models. Let 0 = {/3, a}, where /3 = {PiiP^Pz} 30 and o = {012,013,023,0123} are vectors of parameters associated with the marginal and association models, respectively. We model the logit of the marginal probability, Pr(Y" * = 1 | x) for t t = 1,2,3, as a linear function of the covariates. More precisely, if gt{x;/3 ) = t Pv(Y * = 1 | x), then the marginal model is given by logit{gt(x;/3)} t = X'/3 . We t t denote the association model as g t(x;a t) = Pr(Y * = l , Y * = 1 | x), for {s,t} = s s s t {1,2}, {1,3}, {2,3} and gm(x; a ) = Pr(Y{ = 1,Y * = 1,Y * = 1 | x), where l23 logit {0 st (x; a )} = X' a st t 2 3 and \og\t{g {x\ 0 1 2 3 ) } = - X st X23 1 2 3 The probabilities ai 32 for the different possible outcomes can then be expressed as follows: x;0) — ru,i,o r(i,o,i x-0) = 5i2(x;ai ) -5i23(x;ai23), x;0) = 5i3(x; a y ) - 5i23(x; a'123), ra,o,o x;0) = gi(x;fi ) - g\2{x;cti ) - 9iz(x;aw) + gi 3(x;a ), r (0,1,1 x;0) = g 3(x; a ) - gx (x; a ), no, 1,0 x-0) = 92{x; 0 ) - gn(x; r (0,0,1 x;0) = P3(a:;/3 ) r x;0) = l-9i(x;Pi) r (1,1,1 (0,0,0 ffi23(x;ai23), 2 x 2 2 23 23 2 3 2 123 i23 a i 2 ) - g 3(x; a 2 - gu(x;ai ) -g 3(x;a ) 3 2 2 3 ) + 5i23(x; 0123), + gi (x; cx ), 23 23 U3 -g (x;8 ) -g (x;(3 ) + gi (x;a ) 2 2 3 3 2 12 + g\z{x; an) + g {x; a ) - gmix; a ). 23 23 i23 The above probabilities must sum to 1. Further, each of these probabilities must be bounded between 0 and 1, so that there are many constraints on the parameters. Note that the parameters in the marginal and association models can be interpreted as various types of odds ratios. More detailed interpretation of some of these parameters are provided in Chapter 7 where we discuss the results of the application of this model to our annual data set. But the parameters in the association models do not necessarily have direct interpretation relating to the strength of dependence among the responses. In other words, the magnitude of the parameter estimates may not explicitly reflect whether the responses are positively or negatively associated. Evaluating the correlations among the responses based on this 31 model are straightforward, although somewhat tedious. 4.2.2 Drop-out M o d e l We now consider the model for the drop-out process. We adopt Baker's idea for modelling the drop-out process as presented in (4.11). Let r _ i = {ri,r , • • • ,r -i} t 2 t denote the previous pattern of non-response indicators up to time t — 1 and y£ = {y\,y , • • • iVt} denote the outcomes up to and including the outcome at time t. Let 2 rj denote vectors of parameters associated with drop-out at time t, where t = 1,2,3. t Further, denote M r t - i , y t I 5»7t) = Pr(J?t = aj a I r -i,yt,aj) t (4.12) and model logit{/jt(r _i,y£ | x;rf )} as a linear function of y£ and x. The dropt t out process is ignorable if ht(rt-i,yl,x;r) ) t depends only on observed outcomes and covariates. More specifically, if ht(rt_i,yt | x;r) ) depends only on covariates, t the drop-out is completely random; that is, the drop-out mechanism is referred to as CRD. If, on the other hand, /i((-rt-i,y£ | x;rj ) depends on the observed t outcomes, and perhaps covariates, but not on the unobserved outcomes, the dropout mechanism is referred to as random drop-out (RD). The drop-out is informative if h (r -i, yl | x;r) ) depends on the unobserved outcomes, and perhaps the observed t t t outcomes and covariates as well. Various authors such as Baker and Laird (1988), Fitzmaurice, Laird and Zahner (1996), and Glonek (1999), have drawn attention to the issue of identifiability for non-ignorable non-response models. If there are more independent parameters than available degrees of freedom, a model is clearly not identifiable. But with nonignorable non-response, even some models with fewer independent parameters than available degrees of freedom are not identifiable (Baker and Laird, 1988). Baker (1995) established sufficient conditions for certain non-ignorable non-response models for three repeated binary outcomes to be identifiable. In particular, for r t _ i equal to {p, p}, {p, a}, {a, p}, or {p}, he considered logistic regressions in which 32 the dependence on outcomes is limited to two predictors: y|,OR' * n e ^ * observed as response (LOR), and y* , the last unobserved response (LUR). Models that inLUR clude the predictor y* LUR are non-ignorable. The values of y* QR ^ VLUR depend a n c L on the previous patterns of non-response. For rt_i equal to {a, a}, {a}, and { }, ^t(rt-i Vt I x;t] ) depends only on covariates. More details on Baker's sufficient r t conditions is provided in Chapter 6. For our case, r _ i can take on six patterns: {p, p}, {p, a}, {a, a}, {p}, {a}, t and { }. Since the only type of missing responses in our data set corresponds to drop-outs, we only need to model those /i (rt_i, y£ | x; r) ) where r t - i equals {p, p}, t t {p}, or { }. More precisely, for the cases of drop-outs: • consider r : 2 (i) when r = {p, p} => y* 2 (ii) when r = {p, a} 2 = y%, y* LOR LUR = y* , as in Baker (1995); 3 Pr(_R = a | r ,y3,a;) = 1; 2 3 (iii) when r = {a, a} => Pr(i?3 = a | r , y | , x ) = 1. 2 2 • consider r i : (i) when n = {p} => y* = yj, y* LOR = y%, as in Baker (1995); LUR (ii) when r i = {a} => Pr(i? = a | r i , y 2 , x ) = 1. 2 • when ro = { }, Baker (1995) suggested hi({ },y\ | a:;?7i) should depend only on the covariates, not on the observed and unobserved outcomes. As in Baker (1995), we allow the models for /it(r _i,y£ | x;rj ) when r _ i t t t equals {p,p} or {p} to be nested within one of the following: 1. Covariates (COV) * L U R [= COV + L U R + COV x LUR]: logitf/i^rt-i, yj | x; r) )] t = rf*- + x' rf^y + rf^ 1 + 33 x ' Vcov*LUR x VLURI R y* LUR ( -!3) 4 2. C O V + L O R + L U R : logit[fc (r _i,yt t t | x; rj )} t r/ *-» + x' rj* \, + r = 0 + LUR VLORVLOR (4.14) VLURI N 3. L O R * L U R [= L O R + L U R + L O R x LUR]: logit[Mr -i,yt* t I x; r) )] = rf*- + rf^ y* 1 t + ^ LOR* LOR LUR VLOR + rf ^ y* L R LUR (4.15) V*LUR- The drop-out model considered in Diggle and Kenward (1994) is a special case of the model C O V + L O R + L U R . They assumed the drop-out mechanism only depended on L O R and L U R , and that the effects of these two predictors were the same across different drop-out occasions. Note that the same covariates can appear in both the drop-out and the outcome models. 4.2.3 Likelihood Function We assemble these models into an explicit expression for the logarithm of the likelihood. Let n yy yit 2! 3iX y\ at time denote the total number of subjects with outcome 1, y 2 at time 2, j/3 at time 3 and categorical covariate at level x. Further, denote V = {^h,^,^}- We can then express the log-likelihood as L(0,r)) — VJ L (0, x r/), X where log f{a,a,a \ x,0,ri) 1 1 1 + E E •yl,y^,a,x log f{y*i,y* ,a I x;6,n) 2 2/1=0^=0 1 1 1 +£ £ £ log f{yl,y* ,y*3 I x;0,-n), 2 »I =ow;=o 2/5=0 34 (4.16) where the four functions/(yj, j/5 > J/3 I and f(a,a,a f{yt,yh a I x]0,v), f{y{,a,a\ x;6,r)), \ x;6,rj), are specified in (4.3), (4.4), (4.7) and (4.10) respectively. We obtain the maximum likelihood estimates (MLEs) of the parameters, 9 and 77, by minimizing the negative log-likelihood using a quasi-newton minimization routine [26]. 35 Chapter 5 Transition M o d e l 5.1 Introduction In this chapter, we model the outcome (or measurement) process using a transition model coupled with several models for the drop-out process as described in Chapter 4. The idea of using a transition model to describe the outcome process is motivated by Liu, Waternaux and Petkova (1999), who investigated the effect of human immunodeficiency virus (HIV) status on neurological impairment on a cohort of HIV positive and negative gay men. These subjects were followed for 5 years and assessed every 6 months. The primary outcome is the presence or absence of neurological impairment which varies over time. Predictors of outcome include fixed and time-varying covariates, such as age at baseline, HIV status, disease progression and time of assessment: Nearly half of the subjects dropped out before the end of the study for reasons that might have been related to the missing neurological data. Liu et al. (1999) adapted the likelihood-based approach proposed by Diggle and Kenward (1994) for the analysis of a Gaussian longitudinal outcome with informative drop-out to analyze these binary longitudinal responses. More precisely, they assumed a first-order Markov chain transition model for the binary longitudinal responses combined with different logit models for the occurrence of drop-out. 36 Transition models are often used for equally-spaced longitudinal data when the interest is in prediction (Diggle and Kenward, 1994), and Liu et al. (1999) proposed such a model for the outcome process as their interest was in predicting neurological impairment. Our data set consists of yearly observations on the presence or absence of exacerbations in MS patients. According to Liu et al. (1999), "In biomedical research, sequences of measurements are often fairly short and, in many cases, a first-order transition model is reasonable". Thus in this thesis, we embrace their idea of modelling the repeated binary responses withfirst-ordertransition models. In the next section, we give an overview of the Liu et al. transition model for the outcome process and propose to combine this with Baker's ideas for modelling the drop-out process. We then briefly present the general expression of the loglikelihood under these models to conclude the chapter. 5.2 The Liu et al. Transition Model for Binary Longitudinal Data with Informative Drop-out In this section, we illustrate the general approach of afirst-ordertransition model. To keep the discussion simple and consistent with Chapter 4, we assume each subject is followed at three equally-spaced time points. As in Chapter 4, (Yf, Y , Y *, R±, R , R , X ) , is the vector of random vari2 3 2 3 ables for the complete data, and (Yi, Y2, Y 3 , X ) is the corresponding set of random variables for the incomplete data. The relationship between them is: Yt = Y * if t Rt = p and Yt = a if Rt — a. The joint distribution for the complete data is factored as Pr(Y * = ylY * = y* ,Y * = y' ,R x 2 2 3 3 = Pv(Ri = r R u 1 =r R u = r ,R 2 2 3 = r ,R 2 2 - r | 3 = r \ X = x) 3 yl,y ,y$,x) 3 2 xPr(Y *=yI,Y *=i/ *,Y *=y5 |x), 1 2 37 2 3 where Pr(Y * = y\,Y 1 2 Pr(i?i = n,R 2 = y ,Y * — y \ x) is known as the outcome model, and 2 = r ,R 2 a 3 = r | y\,y ,yl,x) 3 is the drop-out model. 2 3 As for Baker's selection model, the basic idea is to construct a model for both the outcome and drop-out processes. These models then specify a model for the incomplete data. The log-likelihood function is expressed in terms of these models and a maximum likelihood is employed to estimate the model parameters. The only difference from the previous chapter is that here we model the outcome process with transition models. 5.2.1 Outcome Model Denote Ht = {y*,... , y^-i} as the responses up to but not including time t. The joint distribution of the equally-spaced outcome variables given the covariates, i.e. Pr(Y * = y\,..., Y * = yl \ x), can be decomposed as 1 t Pr(y *=yI,...,y * = i/?|a ) 1 t J = Pv(Y * = y* \ H ,x) t t t xPv(Y U^yU\H ^,x) t t x...xPr(y '=yJ|a')1 (5-1) A transition model of order q > 0 postulates that the conditional distribution of y\, given the history Ht, depends only on the observations: y j _ , . . . , yl-i- A firstg order (q — 1) transition model for the case of three repeated responses is of form Pr(y -=yJ,y *=.y2,y *=y3*|x) 1 2 3 = Pr(y * - y* | y* , x) 3 3 x Pr(Y 2 *=y 2 * x PviY^yHx). 2 \y\,x) (5.2) Liu et al. (1999) proposed using a first-order transition model for specifying the joint distribution of an equally-spaced binary outcome process. They employed a specific model for the conditional probabilities of the binary elements in the complete outcome vector y* in which the conditional probabilities are assumed to depend only 38 on the covariates observed at the immediately previous time point, denoted x t _ i . The form of the model is logit{Pr(y; = 1 \y* _ x)} = t lt A-xPi+hvU, (5-3) where /3i and /3 are parameters to be estimated. The parameter /3 represents the 2 2 log odds ratio for presence at time t given presence at time t — 1, against presence at time t given absence at time t — 1. At first glance, the assumed form of dependence on the covariates seems a bit peculiar since the covariates measured at time t should have a stronger influence on the response y\ than the covariates measured at time t — 1. However, they noted that in most biomedical studies, there will be no information available after a subject drops out of the study; that is, if yf is not observed, then x would not be t observed either. They chose to overcome this (potential) data limitation problem by the aforementioned approach. In summary, the outcome model is specified in terms of conditional distribution of yl with the assumption that it depends only on y^_ and x t - i . The structure x of the associations among the responses is more restricted than in Baker's selection model in that this model assumes the association between Y{ and Y * to be the same 2 as the association between Y * a n 2 5.2.2 d Y *. 3 Drop-out Model Similarly, the drop-out model Pr(Ri — r\, R = r , R3 = r$ | y\, y^, 2/3, x) is specified 2 2 in terms of Rt given ( r _ i , H , y* , x) for t = 2,3 as in Chapter 4. Liu et al. (1999) t t t modelled these conditional probabilities as: \ogit{h (r -i,H ,yl t t t \ x;r} )} t = ry + mVt-i + Wit0 (-) 5 This is a special case of (4.13) and (4.14) in which the drop-out mechanism is assumed to be independent of the covariates and the parameters in the logistic regressions are the same regardless of r t - i . In their data set, the first observation 39 4 was always observed, so they did not need to consider drop-out models for the case where r t - i equals {a} or { }. But for our purposes, the drop-out probability with rt_i equals {a} is always 1 . For the case where rt_i equals { }, we model ht({ })£/* | x,r/ ) t according to Baker's (1995) suggestion; that is, this probability should depend only on the baseline covariates, not on the outcome measurements. 5.2.3 Likelihood Function The general expression of the log-likelihood is the same as (4.15). Similarly, the four models for the incomplete data, /(y^y^yt), /(2/*>y2' )> /(yi> > ) a a a a n d f(a,a,a) have the forms (4.3), (4.4), (4.7) and (4.10) respectively, which are specified in terms of the models described in the previous two subsections. L i u et al. (1999) used the S-PLUS function ms to obtain the maximum likelihood estimates ( M L E ) for the parameters in their problem. As in Chapter 4, we use a quasi-newton minimization routine to obtain the M L E s for /3 and rj. In the next chapter, we discuss potential identifiability problems in nonignorable non-response models for incomplete binary data before proceeding to use the Baker's selection model and the Liu et al. transition model to analyze our annual data in Chapter 7. '40 Chapter 6 Identifiability in Models for Incomplete Binary Data 6.1 Introduction Analyses based on an assumption of ignorable non-response when the non-response mechanism is informative (or non-ignorable) can lead to misleading or biased results. Thus in the past decade, various authors have developed models for continuous and categorical response data subject to non-ignorable non-response. In particular, likelihood-based analyses have been widely employed since there is a choice of whether or not to introdure an explicit model for the non-response mechanism. Little and Rubin (1987) noted that, by incorporating a model for non-response in a likelihood-based approach, valid inferences can be obtained when the non-response mechanism is non-ignorable provided the non-response model correctly represents the non-response mechanism. Most of these papers have emphasized the formulation and implementation of those models. However, it has been observed that such models present certain analytical difficulties. In particular, it can happen that the parameters of the non-ignorable models are not identifiable or the maximum likelihood solutions can lie on the boundary of the parameter space. 41 Baker and Laird (1988) drew attention to the issue of boundary solutions to the maximum likelihood equations in a non-longitudinal setting. They illustrated this issue with the pre-election data from four successive Roper polls carried out to predict the proportion of voters preferring Truman in the 1948 presidential election. The four variables used in their analyses were time of survey (XT = July, August, September, October), economic class of voter (XE = A , B, C , D), voter preference (Y = Truman, Dewey, other), and expression of preference (R = yes, no). They employed two different log-linear models to describe the related regressions: marginal outcome model for the XTXEY the margin (a 4 x 4 x 3 array) which describes the regression of Y on XT and XE, and the non-response model for the full contingency table XTXEYR, which describes the regression of R on XT, XE, and Y. For this framework, they showed that with non-ignorable non-response models, overparameterized and saturated models may not yield a perfect fit and the likelihood equations can be satisfied by boundary values even when all observed counts are strictly positive. As discussed in Chapter 4, Baker (1995) used a selection model to analyze data from the Muscatine Risk Factor Study to investigate the effects of gender and age on obesity in schoolchildren who ranged between ages of 5 and 13 years. In these data, each child was intended to have three binary responses at 2-year intervals indicating whether or not they were obese at that point in time. However, there was a substantial amount of non-response due to no consent from the parents or the child not being in school on the day of the examination. In this special setting, Baker (1995) obtained sufficient conditions for non-ignorable non-response models to be identifiable. Following Baker's ideas, we establish sufficient conditions for non-ignorable drop-out models to be identifiable in the last section of this chapter. For the context of models for incomplete multivariate binary data, Fitzmaurice et al. (1996) suggested some simple procedures for examining local and global identifiability in models with non-ignorable non-response. A summary of this por- 42 tion of that paper is given in the following section. More recently, Glonek (1999) formulated the specific application considered in Section 3 of Fitzmaurice et al. (1996) in a somewhat more general fashion to discuss the identifiability issue for models for incomplete binary data. He derived necessary and sufficient conditions for certain simple non-ignorable non-response models (including some of the models considered by Fitzmaurice et al. for their application) to be identifiable. His results show that these models are identifiable except at a set of special parameter values where the conditions fail to hold. The consideration of model identifiability is an issue that should be resolved prior to estimation, because it does not make sense to attempt interpretation of an estimate of a parameter that is not statistically identifiable. In Section 6.2, we describe the procedures suggested by Fitzmaurice et al. (1996) for checking model identifiability. We also describe the necessary and sufficient conditions obtained by Glonek (1999) and the implications of these results for Fitzmaurice et a/.'s suggested approaches to examining the identifiability of non-ignorable non-response models. Baker's (1995) development of sufficient conditions for the identifiability of certain non-ignorable non-response models for the case where the data consist of three repeated binary responses with all possible patterns of non-response is briefly summarized in Section 6.3. We conclude this chapter by applying Baker's ideas to the special situation of interest here of models corresponding to monotone non-response patterns. 6.2 Discussion in Fitzmaurice et al. (1996) and Glonek (1999) Fitzmaurice et al. (1996) proposed a likelihood-based regression model for analyzing incomplete multivariate binary responses based on the multivariate binary model proposed by Fitzmaurice and Laird (1993). The latter model is extended to 43 accommodate incomplete data by assuming a logistic model for the non-response mechanism which depends on covariates and on both the observed and unobserved responses. This idea is motivated by Diggle and Kenward (1994) and Molenberghs, Kenward, and Lesaffre (1997). Throughout Fitzmaurice et al. (1996), monotone non-response is assumed. Various authors have pointed out that the identifiability is an important yet unresolved issue in non-ignorable non-response models. As Fitzmaurice et al. (1996) stated, "So far, no general and practically useful necessary and sufficient conditions for identifiability are available". Fitzmaurice et al. (1996) suggested some simple procedures for examining the identifiability status of non-ignorable models for the case of discrete response variables; these are described in the next subsection. The following subsection describes Glonek's results and the implications of those results for the procedures suggested by Fitzmaurice et al. (1996). 6.2.1 Fitzmaurice et aZ.'s Suggested Procedures Fitzmaurice et al. (1996) indicate what they mean by a non-identifiable model. Consider a non-ignorable model with parameters (0, 77), where 9 and 77 are the vectors of parameters associated with the outcome model and the non-response model respectively. If it is the case that there are distinct parameter vectors (0o,^7o) 7^ ( ^ l i ^ i ) such that / ( y , r i I 0 ,Vo) o i O = \ f{yoi,ri 0 ,rj ) 1 1 for all y j (the vector of observed responses for the i-th subject) and r-j (the vector of Q response indicators for the i-th subject), then L(0Q,r) ) = L(0i,77 ) and the model Q 1 is not statistically identifiable. Showing algebraically that all of the parameters in non-ignorable models are identifiable is not trivial (Fitzmaurice et al., 1996). If there are more parameters to be estimated than available degrees of freedom in the data, the model is clearly not identifiable. But having no more parameters to '44 be estimated than the available degrees of freedom is not sufficient to guarantee identifiability for non-ignorable non-response models (Baker and Laird, 1988). Fitzmaurice et al. (1996) suggested some simple procedures for examining the identifiability of non-ignorable non-response models. Since local identifiability (the model is identifiable in a subspace of the entire parameter space) is a necessary condition for a model to be globally identifiable (the model is identifiable throughout the entire parameter space), a first step is to examine the local identifiability status of the model by checking that the Fisher information matrix is nonsingular. Rothenberg (1971) has shown that, subject to certain regularity conditions, if the Fisher information matrix is nonsingular, then the model is locally identifiable. This idea of using the Fisher information matrix to determine the identifiability status of a model was described in the context of latent class models by Goodman (1974). • Checking for Local Identifiability Fitzmaurice e£ al. (1996) suggested selecting a reasonable set of parameter values for (0, rj) and evaluating the Fisher information matrix at this particular set of parameter values. This can be accomplished by taking the expectation of the outer-product of the score equations, summing all the possible realizations weighted by their respective probabilities. In other words, for each possible realization of ( Y j , R j , X ; ) , calculate the sample covariance matrix of the scores and weight these contributions by their respective joint probabilities. By summing over all possible realizations, the Fisher information matrix is obtained. The information matrix can then be checked to see whether it is nonsingular at this set of parameter values. • Checking for Global Identifiability Having established local identifiability, Fitzmaurice et al. (1996) recommend assessing global identifiability with the following procedure: 1. Select a set of reasonable values for the parameters (0, rj) (e.g. the estimated values) and use them to generate an artificial sample comprising one observa45 tion for each possible realization of (Yj, R;, X,). 2. Solve for (9, r)) from the likelihood equations obtained by weighting the contribution for each possible realization by its respective probability. If the resulting estimate (9, r)) does not equal (8, TJ), then the model is not globally identifiable, and those parameters that give a different value are not statistical identifiable. If the estimate (9, r)) equals (9, 77) for a whole grid of reasonable values for (9, 77), then the model is most likely identifiable (Fitzmaurice et al., 1996). Fitzmaurice et al. (1996) provide a simple example intended to show that the model identifiability problem exists even when the number of parameters is no more than the available degrees of freedom from the data. For the i-th patient there are two binary responses, Yn and Yj , and a dichotomous covariate, Xi. Yn is always 2 observed but Yj is subject to missingness. Thus for each value of Xi, there are 6 2 possible outcomes for (YJI, Yj ): (0,0), (0,1), (0, a), (1,0), (1,1), (l,a). Consequently, 2 the observed data have 10 degrees of freedom, 5 for each of the two possible values of ^ . The outcome model they considered is not fully saturated and they also considered several non-ignorable non-response models. More specifically, the outcome model consists of two parts: a marginal model (for the means of the responses) and an association model. The (unrestricted) marginal model is parametrized as lo&t{E(Yij)} = Boj + ByXi-,. for j = 1,2, but the association between Yn and Yj is assumed to be constant 2 across Xi, i.e. the conditional log odds ratios are assumed to be constant across Xi. Thus, the outcome model involves 5 parameters. This outcome model is coupled with 8 non-ignorable non-response models having at most 5 parameters. With B42 46 denoting the response indicator for Yi , these models are: 2 1. logit{Pr(ifc = P)> = 2. logit{Pr(i? = P)> = r)o + riiXi + r\ Y 3. logit{Pr(Zfo = P)} = Vo + ViYn + r) Yi2 logit{Pr(i2 2 = P)> = r/o + T/iXi + 7/2Y + r) Xi x Y 5. logit{Pr(# = P)> = 7/0 + 6. l o g i t f P r ^ = P)} = Vo + m u + r] Yi2 + m n 7. l o g i t f P r ^ = P)} = 7/o + 7/iXj + 8. l o g i t f P r ^ = P)} = 7/o + rjiXi + r\ Yix + 7/3^2 + m i x Y 2 4. i i2 I i2 2 i2 2 i2 T/I Xi + TfcYn + 3 i2 nY 3 Y i2 Y 2 x i2 Y n Yn + r) Y + r) Xi x Yn 2 3 i2 A x 2 i2 Based on the use of their suggested procedures, Fitzmaurice et al. (1996) claimed that only three of these eight non-response models (Models 1, 2, and 4) are statistically identifiable. However, they do not indicate how they selected reasonable sets of values for the parameters and how many sets they checked to reach their conclusions. 6.2.2 Glonek's Necessary and Sufficient Conditions Glonek (1999) attacks the model identifiability problem from a different point of view. He formulated the problem considered in Fitzmaurice et al. (1996) in a more general fashion to address the issue of identifiability. Two binary responses, Y\ and I 2 , and a categorical covariate X with I levels are considered. Only Y is subject to 2 non-response and R2 is the response indicator for Y (R2 = p if Y2 is observed and 2 R 2 = a otherwise). The outcome model is denoted as n ijk = Pr(Yi =j,Y = k | X = i) 2 for j, k = 0,1. The non-response model is denoted as p ijk = Pv(R = p\Y =j,Y 2 l 47 2 = k,X = i). Thus, the observations corresponding to the i-th level of the covariate X are multinomial across six cells with probabilities 8ijk = KijkPijk = Pr(Yi = j, Y = k, both responses observed | X = i) Qij* = 7Tijo(l - 2 Pijo) + 7Tiii(l - Piji) = Pr(Yi = j , Y2 unobserved | X — i). The simple example used by Fitzmaurice et al. (1996) to illustrate their suggested procedures for checking model identifiability is of this form. As described in the previous subsection, for the case of a binary covariate (1 = 2), they considered a restricted model for 7Tjj involving no three-factor interaction and eight different fc models for pij^. Combined with an unrestricted model for 7 ^ , Glonek (1999) considered homogeneous non-response models of two forms: Pijk = Pjk, (6.1) and Pijk = Pik- > (6.2) In the first of these models, the probability of response is independent of the covariate, while in the second, the probability of response does not depend on the first response variable. Non-response models 1, 3 and 6 of Fitzmaurice et al. (1996) are of the first form, whereas models 1, 2, and 4 are of the second form; models 5, 7, and 8 are of more general forms. For the case 1 = 2 with non-response model (6.1), Glonek showed that the condition P r ( F = l\ 2 Y =j,X 1 = l) ± Pr(Y = 1 | Y = j,X 2 x = 2) (6.3) for j = 0,1, is necessary and sufficient for the parameters of the model to be identified. The condition (6.3) would generally be satisfied, even under the restriction of no three-factor interaction incorporated into the Fitzmaurice et al. (1996) outcome 48 model. However, the restriction does not imply the condition (6.3); the condition could fail to hold for specific values of the parameters. Hence, their outcome model combined with any of their non-ignorable non-response models 1, 3, and 6 is identifiable except at those special values of the parameters where (6.3) fails to hold. Similarly for the case 1 = 2 with non-response model (6.2), a necessary and sufficient condition for the parameters of the model to be identified is P r ( y = 1 | Yi = 0,X 2 for i = 1,2. = i) ^ P r ( y = 1 | Y = 1,X = i) 2 x (6.4) The proof is provided in Appendix A . Again, even under the restric- tion of no three-factor interaction in the outcome model, the condition (6.4) would generally be satisfied. But the restriction does not imply the condition. Hence, the Fitzmaurice et al. (1996) outcome model combined with any of their non-ignorable non-response models 1, 2 and 4 is identifiable except at those special values of the parameters where (6.4) fails to hold. Contrary to the conclusions of Fitzmaurice et al. (1996), Glonek was able to establish that with these homogeneous non-response models, the Fitzmaurice et al. models 1, 2, 3, 4 and 6 for this simple example are identifiable except at a set of special values of the parameters. (He did not address the issue for the Fitzmaurice et al. non-response models 5, 7 and 8.) Thus, Glonek established that the identifiability status of these models depends on the particular values of the parameters. Glonek also provided a simple example with a non-homogeneous non-response model where this phenomenon occurs. This is problematic for inference since it may happen for a particular set of data that the maximum likelihood estimates are welldefined in the sense that the parameters are identified while, in fact, the true values of the parameters that generated the data are not. In such cases, it is clear that local calculations performed at the M L E will not bring to light this underlying non-identifiability. This phenomenon is different from the structural type of nonidentifiability that would lead to rank deficiency in the Fisher information matrix, as considered by Fitzmaurice et al. (1996). 49 Hence, the procedures suggested by Fitzmaurice et al. (1996) are not adequate to resolve the issue of identifiability. Our annual data setting is slightly different from the problem Glonek considered. We have three binary responses, Y"i, Y , Y , and all are subject to non-response. 2 3 The derivation of the necessary and sufficient conditions for the identifiability of nonignorable non-response models following Glonek's ideas appears to be much more complicated in our setting. However, we were able to establish sufficient conditions for certain nonignorable models to be identified in our setting, by following the ideas illustrated in Baker (1995). We briefly describe Baker's ideas in the next section and conclude this chapter with a description of the sufficient conditions we established. 6.3 Discussion of Model Identifiability for Incomplete Binary Responses in Baker (1995) In Chapter 4, we described Baker's selection model for three repeated binary responses. He pointed out that all models with ignorable non-response are identifiable, but identifiability becomes a concern with non-ignorable non-response models. To restrict his models to those that are identifiable, he introduced two predictors for the non-response model: ULORI the last observed response (LOR), and y* uR, the L last unobserved response (LUR). Non-response models that include the predictor V*LUR a r e n o -ig o ble. n n ra Recall that Baker modelled Pr( R = ri,R J 1 2 = r ,R 2 = r \ y^y^y^x) 3 3 in terms of conditional probabilities assuming the non-response does not depend on future events; that is, Pv(Ri=ri,R = r ,R 2 P(R 2 3 = r \yl,y2,yl,x) 3 3 = r \R =r R 3 1 u x P(R = r \R,= 2 2 = = r , Y{ = y\,Y * = y *, Y * = vl,x) 2 2 ,Y{ n 2 2 3 = y{,Y * = y* ,x) 2 x P ( i 2 i = n \Y? = y\,x). 2 • 50 (6.5) Each of these conditional probabilities is modelled as a logistic regression that depends on y* R and y*iuR- The values of these predictors are determined by the LO previous observation pattern, r<_i = {ri, r , . . . , rt-\}. He claimed that the non2 ignorable non-response models are identifiable if the following conditions are satisfied: A. When rt-i equals {a, a}, {a}, or { }, the corresponding conditional non- response probabilities should depend only on covariates. B. When rt-i equals {p, p}, {p, a}, {a, p}, or {p}, the non-response models should be nested within one of the following three types: (a) COV * LUR; (b) COV + LOR + LUR; (c) LOR * LUR. Baker allowed the model parameters to differ for each of the previous observation patterns. Some of the details of the verification of identifiability are presented in the appendix of his paper. Our situation is slightly different from that Baker considered. He had 7 non-response history patterns to consider, i.e. {p, p}, {p, a}, {a, p}, {a, a}, {p}, {a}, and { }. Since the non-response in our data set is monotonic, we need to consider only three different non-response history patterns: {p,p}, {p}, and { }. In the following section, we present verifications of the identifiability of the non-ignorable non-response models considered in our context. 6.4 Discussion of Model Identifiability Our data set is a special case of Baker's general data structure as we have only monotone non-responses, i.e. drop-outs. In particular, we have four monotone nonresponse patterns to consider: {p, p, p}, {p, p, a}, {p, a, a}, and {a, a, a}. Recall that 51 77 = {T7 6 = {/3,a} and 5l T J , 2 AS TJ }. 3 in Chapter 4, we model the incomplete data in terms of the product of the outcome model, Pv(Y{ = yl,Y * = y* ,Y * = y* | X ) = /*(y?,y *,y * | x;9), 2 2 3 2 3 and the drop-out model, Pr(i?! = n , R = r , i i = r I Yj* = y*, Y * = y^,Y * = y *,X) 2 2 2 2 2 = 3 3 <?(ri,r ,r I t/i,y ,2/ ,x;r7), 2 3 2 3 where tj(ri, r , r | y*, y , y , x; 77) is specified as in (6.5). Recall also that y£ denotes 2 3 2 3 the outcomes up to and including occasion t and r t - i denotes the non-response history prior to time t. The three conditional non-response probabilities are denoted as follows: Pr(i? = a I r = {p,p},y3,x) 3 = 2 Pr(i2 = a | r i = {p},y^,x) = 2 Pr(i?i = a I r = { },y^x) 3 h ({p}, y 2 = 0 M{p,p},y;s I x\rj ) 2 \ x; TJ ) 2 M { },yi | a;;»h). Consequently, the drop-out models for the four monotone non-response patterns are ?(P,P,P I y*3,x;rj) = [1 - M { p , p } , y I x;rj )][l - h {{p},y* \ 3 3 2 X;TJ )] 2 2 x[i-M{},yil*;»h)] 9(p,P,a j y*2,x;r)) = M { p , p } , y I x;rj )[l - h {{p},y* \ x;r} )} 3 3 2 2 x [ l - / i i ( { },yi ?(p,a,a I yi.xjT/) = h ({p},y \ x;r) )[l - hi{{ },y\ tj(a,a,a I *;?/) = h {{ },y{ \ x;77^. 2 x 2 2 52 | x;^)] 2 For the case of categorical covariates, the kernel of the log-likelihood function is L{0,r]) = 2~2 x(0,v), where L X 1 L (0,v) x 1 1 =S E E n w!,y5.»5.x g{/*(l/i.y2.y3 l o ~M{P>P}>y5 I x;ri )] I 3 yl=0y*=0y*=0 • x [1-M{p},y2 1 1 \x\V2)][l-hi{{},yl + XI !/J.y5.^ g{Xl f*(y*>y2,y% n lo yJ=0j/J=0 i + 51 i i 3 }.yi I ^m)]} /*(j/i,l/2.S/3 I x , t f ) / i 2 ( { p } , y 2 I x;ri ) 2 !/2=0 J/*=0 x [ 1 } i + n , , log[ a 51 yr>w s{H l o 2/i* =0 0 \ x;r] ) Iz,0)M{P»P},ys j/*=0 x [i - M { p } , y 2 1 a;;»?2)][i - ^ i ( { n Ix;^)]} 1 , y l |ar;»7i)]} l l £ ^ /*(yi,!/ ,!/3 « J = 0 y j = 0 j/*=0 aiX 2 I x,0)hi({ },y\ | X;T7I)}- (6-6) In the following subsections, we discuss the identifiability of the drop-out models by verifying whether the conditional probabilities /i ({p, p},V3 | x;r) ) 3 ^2({p},y2 I 'i"2) x a n 3 } d h\({ },y{ I z;»7i) are identifiable under the conditions de- scribed in the previous section. That is, ^ ({p,p},y3 | x;rj ) could depend on y , 3 3 3 y\ and x, while /i2({p}>y I 'i t2) could depend on y*, yj*. as l l as a;. However, x r w e 2 h\({ },y\ I a;;»7i) is allowed to depend only on x. 6.4.1 Identifiability o f h ({p, p}, y£ \ x; rj ) 3 3 The contribution of h ({p,p},y$ \ x;rj ) to L (0,rf) is given by 3 3 x 53 1 1 1 S E E io {r (vh2/2)2/3 x [i - i + I z,0)[l ~ M{P>Pj>y | X ; T 7 ) ] g yi=o y^=o y*=0 M{ },y21 P 3 3 ^;r? )][i - M { },y*i I s;»h)]} 2 i 5Z %i.j/2^ g{Z^ /*(yi)y2>2/ I x,0)M{p,p},ys I lo J/i*=Oj/*=0 3 x ^ 3 ) y|=0 x [1 - ^ ( { p } , y I s;rj )][l - M { },yj | x;^)]} 2 a 2 To simplify the notation, substitute i for y*, j for y n j / * ,2/2-2/3",* = ^sijfc a n d ny* y* a,x t t x j k d & f°r 2/ a n 3 d denote We further define Wxij- 1/2,1/31 x,o)[i - /i3({p,p},y31 Pxijk = , . = a n 2 (6.7) «;»?3)] x [I - M { p } , y 2 1 «;»72)] [i - M { }.yi I *;*h)] M{P,P>,y I x;T73) 1 - /i3({p,p},y I x;r7 )' _ . 3 3 3 where the notation reflects that /i ({p,p},y | x;rj ) 3 3 does not depend upon y\. 3 Then (6.7) can be re-expressed as: i l l X Yl i ]L i mxi k ^ i Pxijk+x ] L i=0 j'=0 A;=0 l lo s{Xl piijfc^jfc}- fc=0 j=0 j=0 This is identical to the log-likelihood for a contingency table {m ijk} with a supplex mentary margin {w ij } corresponding to cases where k was not observed. Therefore, x the expected cell counts for m ijk and w ij are jJ. ijk l x x X = Pxijk( x+++ + x++) m w a n d 2~2 fJ>xijk<t*xjk, respectively. We address the identifiability of hs({p,p},y 3 | x;rj ) 3 for the three specific forms of non-response models introduced by Baker (1995). In each case, ci saturated outcome model is assumed. • COV * LUR This model has 4>xjk = <l>xk, implying two distinct parameters for each level of x. 54 A perfect fit requires p ijk xijk = m x xij w = 1 1 d w ij = VJ ii ijk fixk- Hence, we require fe=o x X X) xijk^xk- Thus, for each level of x, we have four equations in the two m k=o unknowns, <f> o level. Hence, • a n x and <§>\- The parameters are overdetermined even if x has only one x hs({p,p},Vg | x;r}3) is identifiable under this specification. COV + LOR + LUR In this model, we can represent <j> jk — 4>x <f>jfik-If we denote X 4>\\\ = <t>, <f>uo = 4> <t>K, <t>\oi = 4> <f>J, </»ioo = <t> <t>J <t>K, then if x has only 2 levels, we can write (f>2l\ = <f> 4>X, <f>210 = <f> <$>X <t>K, ^201 = </> <f>X <1>J, <p200 = <t> <t>X <t>J 4> 1 A perfect fit requires w ij = VJ m ijk f>x<f>j l>k- For x = 1 (level 1), we have the < x ( x following equations: win = m i n i <p + m uw tono = mnoi $ <pj + <f> 4>K mnoo 4>K = mion <f> wioo = miooi <j> 4>j +TOIOOO<i> <t>J <t>K IOIOI + mioio <j> <t>j (6.8) <f> (6-9) (6.10) 4>K ( -ll) 6 We can solve the two linear equations (6.8) and (6.10) for the two unknowns <j> and (f>K- Substituting these solutions into (6.9) and (6.11) yields two equations for (f>j and thus (j>j is overdetermined. The equation: 10211 then yields a value for = <f>x- "12111 fi <j>X + ™2110 4> 4>X <j>Ki Indeed, each of the w ij equations yields an equation for 2 <PxIf x has more than 2 levels, we would write <f> as fif for the first level, fix as (j>2forthe second level, etc. In other words, there is one parameter for each 55 level of x and these parameters can all be identified. Thus, this specification for ^3({PiP};y*t • I\}) x r 3 is identifiable. LOR* LUR This model has 4>xjk = fijk as there is no dependence on the covariate x. Thus, there are only four parameters for all the levels of x. As before, a perfect fit requires w ij x l = YJ mxijkfijk which represents four linear equations in the same four fc=o unknowns for each level of x. Hence, /i ({p, p } , y 3 | x;r] ) 3 3 is also identifiable under this parameterization. In summary, /i3({p,p},y3 | x;rj ) can be identified if its form is one of the 3 three types considered above. 6.4.2 I z;^) Identifiability of h ({p},y2 2 h,2{{p},y \ x\rf ) The verification of the identifiability of ^3({p,P},y denote 3 I s;^)- v i =n* , x y l n 2 2 is similar to that for addition to the notation from the previous subsection, and >a>a x I x]r} ) i-M{p},ySI*;i )' M{p},y*i l x i ] 2 2 The contribution of /i2({p}>yj! i l I ^i^) l ^2j=Q ^2 ^2 i=0 J mxi k § i + l0 k=0 t ° Lx(0,ri) P x i 3 in (6.6) can be expressed as: l l Yl i &{X^ Pxijkfixjk } ol Wxi k i=0 j=0 1 1 k=0 1 + ^W jlog|^^/J j I x i=0 j j f c [ l + txjkhxij}- (6.12) j=0 k=0 This is identical to the log-likelihood function for a contingency table two supplementary margins, namely {w ij} x (where {m ijk} k was not observed) x with and {vxi} (where neither of j and k were observed). Therefore, the expected cell counts for l ™ ijk, xij and v w x l xi are n k xij = p ijk{m x x+++ l . X] Z) A xijfc(l i + (pxjkhxij, respectively. j=0k=0 56 +w x++ + v ), VJ Pxijkfixjk and x+ • COV * L U R This model has (f> j = <p x k and jxij 1 1 xk l Wxij = 2~2 fc=o Uxijk<i>xk and v = xi j=o k=0 = 2~2 = wXIJ j j. x A perfect fit requires Miijfc(l + <f>xk)lxj- aj xi = m j, k xi k Hence, we require (6.13) y^fnxijk^xk, fc=o and l vi X l ^2J2 xijk(l m = + (6.14) 4>xk)lxj- j=0 jfc=0 For a fixed level of x, (6.13) represents four linear equations in the two unknowns, (j) i, CJ) Q and x x indicating these are overdetermined. With solutions for (6.14) represents two linear equations in the two unknowns, j o x c^ x0 and j i. x and <j> i, x Thus, /i2({p},y2 I x;r7 ) is identifiable under this model. 2 • COV -I- LOR + LUR In this model, we can represent for (f> jk X <f> j = (fi 4>j (t> x k x and k j ij = j 7^ jj. x x The equations are identical to the earlier case for this model and so are identifiable provided the covariate takes on at least two levels. It remains to show that the parameters j ij x can also be identified. The equations for v { are x 1 Vxi = 1 ^2 ^2 rn k(l xij + <j> jk)lx x 7i Ij j=0 k=0 1 = ^T,lx li Ij M , xij (6.15) 3=0 1 M j = ^ m ij (l + </>Xjk) is treated as known since solutions for the </>'s X k fc=o exist. Suppose x has 2 levels. Using the same representation for jxij as was used where xi -57 for (fixjk earlier, these equations become no = M o 7 7/ 7 J + Afioi 7 7/, m = M V20 = M200 7 ix 7/ U + M201 7 7x 7/, v\ = M \o 7 ix U + Mn 2 (6-16) W no 7 jj + M 2 m 7, 2 (6.17) (6-18) 7 ix- (6.19) Taking the ratio of (6.18) to (6.16) to eliminate 7 7/ and of (6.19) to (6.17) to eliminate 7 leads to two equations in 7^ and 77 from which 7x is easily eliminated. This leads to a quadratic equation in 7^; that is, A jj + B A = M100M210 - VOR-M110.M200 B = (M101M210 + M211M100) - VOR(MUIM2OO C = M101M211 - VOR = no/vn 7—• + C = 0, where + M201M110) V RM IM O 20 U1 W20/W21 A perfect fit requires real roots, or B 2 — 4 A C > 0. Thus, /i2({p} y2 5 I 'i l ) x r 1S 2 identifiable under this model provided the covariate takes on at least two levels and the equation B — 4 A C > 0 is satisfied. 2 • LOR *LUR This model has 4> jk — 4>jk and 7 ^ = jij. x Thus, there are 4 distinct parameters of each type for all the levels of x. These 8 parameters can be identified from the equations for a perfect fit: 1 y^'m i <f> , = Wxij x jk jk (6.20) k=0 and 1 Y2 j=0 1 XI xijk{ik=0 m 58 + <t>3khij- (6-21) For each x, (6.20) corresponds to 4 linear equations in the same 4 unknowns as | x;rj ). in the verification for hs({p,p},y3 Substituting these solutions for the 3 </>s' into (6.21) leads to 2 linear equations in the same 4 unknowns for each x. The 4 jij parameters are determined as long as x has 2 or more levels. Hence, ^2({p},y2 I 'i l2) x 6.4.3 r 1S identifiable provided the covariate x has 2 or more levels. Identifiability of hi({ },y{ \ x^i) In addition to the notation from the previous subsection, denote z = n , x aA ayX and i - M { }>yi I x;f?i) The contribution of h\({ },yj | x\r) ) to L (8,r)) can then be expressed as: x i i x i X 5Z ^2 i log p xijk m x i j k ^ = 0 ^ = 0 ^ = 0 i l XX + ^glX^'* ^^} 1 fc=o 2^=0 ^ = 0 i j=0 i 1 1 i fc=0 1 +^xlog| X) X 5Z / ^j'fc( + ^ f c ) ( + TxijMx}9 1 1 (6.23) 3/J=0 2 / ; = 0 j / * = 0 A perfect fit requires l Wxij (6-24) y^,m ijk<l>xjk, = X fe=o vi = X l l Y2^2™> ijk(l + ^xjkhxij j=o k=0 x (6.25) and z l l l = E E E i*( «j=o»;=oy5=o m i i x 1 + M( 1 + 7«i)^-. (6-26) • COV * LUR This implies (f> jk = </> k and j j x x xi = ^ j, while the equations (6.19) for w ij and x 59 x (6.20) for v { are the same as before. The argument in the previous subsections x shows that the <f> k and 7 j are identified. With these solutions, (6.21) becomes a x X single equation in one unknown, namely 8 . Thus, 5 is also identified. In other X words, h\({ },y\ | z;T7i) X is identifiable. • COV + LOR + LUR In this model, we can represent <f> jk = <f>x <f>j <j>k and j ij — 7x 7 i 7j- The argument x x for the identifiability of the <> / and 7 parameters is identical to that in the previous subsection. Additionally, we have a 6 parameter for each level of x in (6.21). In X other words, there exists a solution for S provided the solutions for the <f> and 7 x parameters exist. Hence, h\({ },y\ \ is identifiable. • LOR * LUR This model implies 4> jk = <f>jk (4 parameters for all levels of x), "f ij = lij (4 x x parameters for all levels of x) and 6 = 5 (1 parameter for all levels of x). The X argument for the identifiability of the (f> and 7 parameters is again identical to that in the previous subsection. The additional parameter, 6, can be determined from (6.21) provided solutions exist for the <> / and 7 parameters. Hence, hi({ }, | x; r^) is identifiable. Thus, we have shown that, when coupled with a saturated outcome model, the parameters in the drop-out models of the three forms suggested by Baker (1995) are identifiable. Notice, that we only consider the case where the covariates are categorical. In the next chapter, we analyze our annual data set with the models mentioned in the previous chapters. 60 Chapter 7 Application to the Data 7.1 Introduction In this chapter, we implement the selection model approach for our annual MS data as described in Chapter 2. Recall our study questions of interest are: • to investigate the most appropriate form of drop-out model for our annual data (in particular, to explore whether the data provide evidence of informative drop-out); • to assess the sensitivity of inferences concerning the treatment effects (and other covariate effects) to the form of drop-out model employed; • to explore the influence of baseline covariates. Recall that the basic idea of a selection model is to factor the joint distribution for the response variables (Y) and the indicator variables corresponding to whether or not the response variables are observed (R) as follows: /(Y,R) = . / ( R | Y)/(Y). (7.1) Thus, the selection model approach involves the specification of a model for the outcomes, / ( Y ) , and for the drop-out pattern conditional on the outcomes, / ( R | Y ) . 61 The outline of this chapter is as follows: Section 7.2 considers a simple structure for Baker's selection model where only treatment group and time are included as covariates in the outcome model. This outcome model is coupled with a LOR-f-LUR type of drop-out model. In Section 7.3, we consider three more general model specifications for the drop-out process in conjunction with the same outcome model: COV * LUR, COV + LOR + LUR, and LOR * LUR. We extend this simple model by incorporating other baseline covariates described in Section 2.2.3 into the outcome model in Section 7.4. The latter two sections can be viewed as further explorations of Baker's selection model. We conclude the chapter with a brief discussion of the use of the Liu et al. transition model for the outcome model. 7.2 Baker's Selection Model: With Only Treatment Groups and Time as Covariates As described in Chapter 4, Baker (1995) suggested specifying the outcome model in terms of marginal and association models. The drop-out process is modelled using a time-dependent causal model assuming the non-response does not depend on future events. • R e p e a t e d B i n a r y Outcomes w i t h Informative D r o p - o u t • o Outcome Model The outcome model / * (y^y^Vz I x;0) is expressed in terms of marginal and as- sociation models. As is apparent from Figure 2.3, the proportion of patients with exacerbations seems to vary across the treatment groups and with time, so the marginal model employed is logit{«7t(as;/3)} = fa + faLD + faHD + fat, (7.2) where t — 1,2,3, and LD and HD are indicator variables to represent the treatment groups. For patients in the LD group, LD = 1 and HD = 0. Similarly, LD = 0 and 62 HD = 1 if patients belong to the HD group. For patients in the P L group, both LD and HD take on value 0. We propose modelling the 2-way and 3-way associations with different intercept parameters to describe different degrees of association. We further assume the association among the responses is related to the treatment arms. For simplicity, these treatment effects are taken to be the same for all associations. • Models for 2-way Association: logit{g {x; a )} st = a + ct\LD + a HD st st (7.3) 2 where st = {12,13,23}. • Model for 3-way Association: logit { # 1 2 3 {x; 0 1 2 3 ) } = « i 2 3 + OL\LD + a HD. (7.4) 2 Both the marginal and association models remain the same throughout the analyses in this section regardless of the assumption on the drop-out mechanism. The adequacy of this non-saturated outcome model for our data has been confirmed by comparing it to various more general models. This information is presented in the next subsection. o Drop-out Model We model the drop-out process using time-dependent causal models assuming the non-response does not depend on future events. We allow different regression parameters for the logistic regressions specifying the different conditional probabilities of absence, /it(rt_i,y£ | x,n ); see (4.12). To simplify the notation, we introduce t two subscripts for these regression parameters: logit{/i (r = {p,p},y | 3 2 3 aj,T} )} 3 = logit{/i (ri = {p},y | x, r» )} = 2 2 2 63 r/03 + VnV + mzVz 2 r/ + rjnvt + V22V 02 2 (7.5) Table 7.1: Drop-out Models under Different Drop-out Mechanisms: J denotes inclusion of a parameter and Vi denotes parameters which are restricted to be equal Drop-out Mechanism ID RD CRD Parameter Model 1 2 3 4 5 6 1 2 3 1 2 V03 V V vo V V no V V 110 V23 V m m v V02 Vl2 V V 7 V2 V2 Vi Vi Vo Vol V V V V2 V2 V Vo - V2 V2 Vo - V - - V2 V2 Vo - V V - Vo Vi Vi - V - - - m Vi V V V22 - V - - Vo - - Vo . - V V V V Vo Vo where the first subscript indexes the specific parameter in the model, while the second subscript indexes the year the drop-out occurred. According to Baker (1995), if the conditional non-response probability in the first year, Pr(i?i = a | y*,x) = hi(ri = { },y* | x,r]i), depends only on the covariates, then the non-ignorable nonresponse models under consideration will be identifiable. In our case, the model for M r i = { }>2/i I iV\) x becomes: logit{/»i(n = { },y{\ x, )} Vl = r/oi. (7.6) These drop-out models belong to Baker's LOR + LUR class of models. For simplicity, we have taken the drop-out mechanism to be independent of the available covariates. We relax this assumption in Section 7.3. To explore the adequacy of simpler models, we consider five other model specifications which are obtained by letting certain parameters be equal or be equal to zero. The ID models to be considered are summarized in the first six rows of Table 7.1. 64 • Repeated Binary Outcomes with Ignorable Drop-out • To investigate the types of drop-out in our annual data, we also fit the data to models under ignorable drop-out assumptions, i.e. with RD and CRD models. o Random Drop-out We consider the three RD models summarized in Table 7.1. Modifying an ID model by setting the parameters associated with the unobserved response to zero leads to an RD model. For instance, RD1 (Model 1 under RD) is obtained by setting ??23 = ?722 = 0 in ID1 (Model 1 under ID). RD2 and RD3 are similarly obtained from ID2 and ID3. o Completely Random Drop-out The two CRD models considered in Table 7.1 are obtained by simplifying the RD models. Under CRD, the drop-out mechanism is independent of the measurement process. Thus, CRD1 (obtained by setting 7713 = 7712 = 0 in RD1, or 771 = 0 in RD2) and CRD2 (obtained by setting 771 — 0 in RD3) each consist of only intercept parameters. For both the RD and CRD cases, we also have the opportunity to examine the sensitivity of the covariate effects (treatment and time) under different forms of the RD or CRD models. These outcome and drop-out models can be assembled into explicit expressions for the logarithm of the likelihood (see (4.16)). The maximun likelihood estimates of the parameters in these models are obtained by minimizing the negative log-likelihood function using a quasi-newton (QN) minimization procedure. This procedure is briefly described in the following subsection. The corresponding results are summarized in next subsection. 65 7.2.1 The Quasi-Newton (QN) Algorithm The QN algorithm used to maximize the log-likelihood is a variable metric algorithm. All variable metric methods seek to minimize a certain function S(6) (in our case, S(6) is the negative log-likelihood function) of p parameters by means of a sequence of basic iterative steps (7.7) 6' = 0- kBg where g is the gradient of the function S, B is a matrix defining a transformation of the gradient and k is a step length. Consider the set of nonlinear equations formed by the gradient at a minimum 9(0') = 0. (7.8) As in the one-dimensional root-finding problem, one can use a linear approximation from the current 0, that is 9(O')^g(0)+H(O)(O'-O) (7.9) where H(0) is the Hessian matrix (the matrix of second derivatives of the function S). For convex functions, H will be positive definite. From (7.8), (7.9) becomes 0' ttO-H^WgiO) (7.10) which is Newton's method for a function of p parameters. This is equivalent to (7.7) with B - H _ 1 and k = 1. Newton's method is generally preferable if second derivatives can be analytically computed. But the implementation of Newton's method may induce errors when closed form expressions for the second derivatives do not exist as it involves composing subroutines for evaluating p first derivatives, p second derivatives and 2 a matrix inversion. For these reasons, Newton's method does not recommend itself for some problems. 66 If H _ 1 could be approximated directly from the first derivative information available at each step of the iteration, this would save a great deal of work in computing both the matrix H and its inverse. This is precisely the role of the matrix B in the iteration defined by (7.7). The transformed gradients in the matrix B are used to generate linearly independent search directions; equivalently, these search directions are conjugate to each other with respect to H. Further, the step parameter k is rarely fixed; its value is usually determined by some form of a linear search. In particular, the role of A; is to allow a search for values of 6' at which the function value is reduced, i.e. S(6') < S(6). Since the second derivatives required in Newton's method are approximated in the iteration (7.7), this algorithm is known as a guasi-Newton method. We employ the QN algorithm suggested in Nash (1979). It involves specific choices of the formula for updating the matrix B and of the linear search procedure for obtaining the updated values of 6'. A n 'acceptable point' search procedure suggested by Fletcher (1970) and a matrix-updating formula for B due to Broyden (1970a, 1970b), Fletcher (1970) and Shanno (1970) are employed. Generally speaking, the algorithm first goes through a linear search to find one value for 6 which gives a smaller function value than that at the previous value for 6. The approximation to the Hessian matrix is then updated accordingly. The algorithm stops when all the parameter values on consecutive iterations are sufficiently close. For our purposes, the absolute difference between the parameter values of consecutive iterations must be smaller than 1 0 . A detailed outline of this algorithm can be -7 found in Chapter 15 of Nash (1979). Note that in this version of the QN algorithm, the matrix B is initialized as a unit matrix. This simple choice nevertheless has the advantage of generating the steepest descent direction (Nash, 1970). To ensure rounding errors which occur in updating the matrix B and forming the search directions, t, through the equation t = 0' - 0 = -kBg 67 have not accidentally given a direction in which the function S cannot be reduced, a reset of B to a unit matrix is suggested in any of the following cases: (i) t 9 > 0; that is, the direction of the search is 'uphill'; T (ii) 6' = 6; that is, no change is made in the parameters by the linear search along t; (iii) t {g(0') — g(Q)} < 0; that is, an updating contrary to the objective of the T method to reduce S along t (t g(6') is expected to be greater (less negative) T than t g(6)), indicating a danger that matrix B may no longer be positive T definite. If either (i) or (ii) occurs during the first step after B has been set to the unit matrix, the algorithm is taken to have converged. All results described in this thesis are obtained using this Q N algorithm implemented in C. The results for the models described in the beginning of this section are discussed in the next subsection. 7.2.2 • Results Adequacy of the Outcome M o d e l To verify the adequacy of our reduced (non-saturated) outcome model, we consider four more general outcome model specifications. These outcome models are 1. Saturated: a saturated marginal model (9 distinct parameters) and a saturated association model of the same form as (7.3) and (7.4) but with regression parameters that differ for each of the 2-way and 3-way association models (12 distinct parameters); 2. Semi-saturated I: a saturated marginal model and a reduced association model with common treatment effects in the 2-way associations (8 distinct parameters); 68 Table 7.2: Negative Log-likelihood Values for Five Outcome Model Specifications Outcome Model Saturated Semi-saturated I Semi-saturated II Semi-saturated III Reduced Negative Log-likelihood 928.923 930.450 931.680 930.304 933.407 Number of Parameters 28 24 22 23 17 3. Semi-saturated II: a saturated marginal model and a reduced association model with common treatment effects for all associations (6 distinct parameters). Note that this reduced association model is exactly (7.3) and (7.4); 4. Semi-saturated III: a reduced marginal model assuming linearity in time (4 parameters) and a saturated association model (12 distinct parameters). Note that this reduced marginal model is exactly (7.2). The negative log-likelihood values presented in Table 7.2 correspond to these outcome models coupled with the drop-out model (7.5). The likelihood ratio test (LRT) indicates the reduction from the fully saturated outcome model to semi-saturated I is reasonable (LR statistic = 3.05 on degrees of freedom (df) = 4; p-value = 0.55). To examine whether the treatment effects in the association model can be taken to be common across all associations, we compare semi-saturated I to semi-saturated II. The L R statistic of 2.46 (df = 2; p-value = 0.29), indicates the reduction is permissible. The result based on a direct comparison between the saturated and semi-saturated II models also agrees (LR statistic = 5.51 on df = 6; p-value = 0.48). This indicates that an association model with common treatment effects for all associations is reasonable for our data set. The further reduction to our reduced outcome model is also allowed (LR statistic = 3.45 on df = 5; p-value = 0.63). 69 As our primary focus is on the marginal model, a more interesting comparison is between the semi-saturated III and saturated outcome models. In the context of a saturated association model, this provides an assessment of whether the reduced marginal model (7.2), which incorporates additive treatment effects and a linear pattern over time for the log odds of having exacerbations, is reasonable. The LRT allows this reduction (p-value = 0.74). As should be expected from the earlier comparisons, the semi-saturated III model can be further reduced to our non-saturated model (p-value = 0.40). Both sequences of model reductions lead to the same conclusion: the reduction to the model presented in the beginning of this section is permitted. This reduced model also provides an adequate fit to our data. The usual goodness-of-fit statistics based on the 15 different possible patterns of binary responses for each treatment arm lead to G = 24.65 and X 2 2 = 22.80 on 25 degrees of freedom (p- values = 0.48 and 0.59 respectively). Thus, we can proceed confidently with further work using this reduced model as a starting point in the investigations. • Informative Drop-out (ID) The detailed results corresponding to the six ID models described in Table 7.1 can be found in Appendix B: Tables B . l to B.6. These tables include the sets of starting values (SV) used, and the maximum likelihood estimates for the parameters (Est), the corresponding standard errors (SE), and the negative log-likelihood computed at the M L E which are all provided as part of the output from the Q N minimization procedure. The number of iterations needed to achieve convergence is also cited in the tables. In each of these tables, regardless of the starting values in the Q N procedure, the corresponding negative log-likelihoods computed at the parameter estimates (at convergence) are the same (at least up to the 4 significant decimal digits displayed). However, in Tables B . l , B.2, B.4 and B,5, not all the reported MLEs are the same 70 (see especially for parameters 7703 and 7/23 in Tables B . l and B.4, and parameters 7703, 772 and 7702 in Tables B.2 and B.5). Also, in these four tables, the SEs for the estimates vary quite a bit across different sets of starting values. This phenomenon might be due to how the Hessian matrix is approximated in the minimization procedure. As mentioned earlier, the Hessian matrix is approximated based on the search directions for the parameter estimates obtained in each successive iteration. To illustrate, consider starting value Sets #1 and #4 in Table B . l . For Set #1, the estimated Hessian matrix was reset to a unit matrix at the a positive definite matrix. 56th iteration due to it not being The final SEs as displayed thus depend on both the parameter estimates at convergence and the corresponding search directions at the subsequent iterations, i.e. the 57th iteration until convergence was achieved (at the 71st iteration). The estimated Hessian matrix for Set #4, however, was reset to a unit matrix three times during the process of minimization (at the 4th, 9th, and 75th iterations), with convergence established at the 91st iteration. Since the process of minimization for the two sets was quite different, this might be the reason why the estimated SEs differ considerably from one set of starting values to another. The substantially different values of the estimates obtained with different sets of starting values for some of the parameters in models ID1, ID2, ID4 and ID5 indicates a more fundamental difficulty. example. Consider the results for model ID1, for Table B . l shows the parameter 7703 is always estimated as being large negative, while 7723 is always estimated as large positive. Furthermore, for all four sets of starting values, the sum of these two parameter estimates equals a constant value, —1.548. This suggests the maximum likelihood estimates for this data set satisfy the constraint 7703 + 1723 = —1.548, with the M L E occurring on the boundary of the parameter space (7703 = — 00 or 7723 = 0 0 ) . Recall that the non-response probability for the third observation is modelled as a logistic regression on last observed outcome (y^) and last unobserved outcome (7/3); see (7.5). When 7703 = —00, 7703 + *723 = —1-548 and 7713 is finite, the probability that the third observation 71 is missing is estimated to be zero if the history is either {y = 0, y | = 0} or {y = 2 1,2/3 = 2 0}' but non-zero for the remaining two histories. The same phenomenon is observed for model ID4 in Table B.4, but with 7703 +1723 = —1.165. This phenomenon is also apparent for models ID2 (Table B.2) and ID5 (Table B.5), but manifests itself in a slightly different fashion. Here, the parameters parameter 772 7703 772 and 7/02 are always estimated as being large negative, while the is always estimated as large positive. However, the sum of always equals a constant, and the sum of and 7702 772 7703 and equals another constant. The pair of constants differ from model ID2 to model ID5. Thus under models ID2 and ID5, the probabilities for the second and third observations to be missing are estimated to be zero when the past observations are either {y* = 0, y = 0} or 2 {y{ = 1, y =0}, and when the history is either {y = 0, y | = 0} or {y = 1, y =0}, 2 2 2 3 respectively. In the next few paragraphs, we discuss the issue of boundary solutions for model ID1 in greater detail. The corresponding discussion for models ID2, ID4, and ID5 is omitted as the details are essentially identical to model ID1. But the results for these three models evaluated at the boundary solutions are also presented. • Discussion of Boundary Solutions Consider model ID1. The estimates obtained for 7703 and vary across different starting values, but in each case 7723 displayed in Table B . l 7703 + 7723 = —1.548. Further the negative log-likelihood remains the same up to the four decimal digits displayed. We believe that the M L E is located on the boundary of the parameter space. To confirm this conjecture, we first use a graphical visualization of the negative loglikelihood function incorporating the special feature (e.g. 7703 is estimated with large negative value, while 7723 is estimated with large positive values, and the sum of the two is always the same) observed in Table B . l . 72 Figure 7.1: A Two-Dimensional Profile Log-likelihood Surface for Model ID1 Figure 7.1 is a graphical representation of the profile log-likelihood surface for the parameters 7703 and 7723 in model ID1. This three-dimensional plot is produced by maximizing the log-likelihood over all parameters except values of 7703 and 7723, 7703 and 7723. For fixed we apply the QN minimization procedure to the negative log- likelihood function. This log-likelihood value is then plotted against these values for 7703 and 7723 using the S-PLUS function "persp". We chose the values for 7703 and 7723 to be a sequence of numbers between —.20 and 20 with increment size of 0.5. This yields a 81 by 81 grid of log-likelihood values. Notice that there seems to be a steady, but very shallow, decrease in this surface along a line (where the grid where 7703 and 7723 7703 + 7723 = —1.548) in take on values ranging from —20.0 to —0.5, and from 0.0 to 20.0, respectively. This seems to agree with the results presented in Table B . l . We also computed the log-likelihood on the boundary of the parameter space to check that the log-likelihood values obtained in Table B . l are what one would obtain at the suggested point on the boundary. Because the parameter estimates appear to satisfy the constraint in terms of 7703 and 7723 = ~V03 7703 + V23 = —1.548, it is useful to re-parameterize + A , where A is a finite-valued parameter. As 7703 approaches —00, the log-likelihood is a function of the remaining parameters and A . For the probability of non-response, Pr(i?3 = a | {p, p}, y , x), we substitute the 3 values presented in Table 7.3 to obtain the reduced log-likelihood function. Applying the QN minimization routine to this reduced negative log-likelihood function yields the results summarized in Table 7.4. The estimates for the model parameters are essentially the same as those presented in Table B . l and the loglikelihood value also agrees. Thus, both Figure 7.1 and this computation of the log-likelihood at the indicated boundary point seem to support our conjecture that the parameter estimates for model ID1 occur on the boundary of the parameter space. The same values of the estimates reported in Table 7.4 were obtained with different choices of starting values and these minimizations required many fewer iterations than those presented in Table B . l . Further, the estimated Hessian matrix 74 was never reset to a unit matrix during these minimizations. Notice that the standard errors for 7/02 and 7722 in Table 7.4 are relatively large. One might suspect this reflects a potential boundary solution phenomenon for the reduced log-likelihood even though these estimates did not vary with the sets of starting values chosen (see also Table B.l). Perhaps these large standard errors are simply indicating that our data set does not contain sufficient information to obtain precise estimates for these parameters. We explored this further graphically. Figure 7.2 shows the profile log-likelihood surface for the parameters 7722 of the reduced model ID1. The values for 7702 and 7722 7702 a n d were chosen to be a sequence of numbers between —20 and 20 with increment size of 0.5. The plot is not very informative in terms of revealing the existence of optimal solutions. The rotating option in "persp" allowed us to view Figure 7.2 from different directions and convinced us of the existence of optimal solutions in the interior of the parameter space for this reduced log-likelihood function. For further assurance, we also calculated the log-likelihood values at various points in the neighbourhood of the suggested estimates for 7702 and 7722; these values are all larger than 933.407. Thus we are certain that this situation does not indicate a boundary solution, but simply indicates a lack of information in the data to precisely estimate these parameters. One can easily show, in a similar fashion, that the parameter estimates for models ID2, ID4 and ID5 also occur on the boundary of the parameter space. The corresponding results for these three models computed at the suggested boundary points are presented in Tables 7.5, 7.6, and 7.7. Note that the parameter estimates in the outcome models for ID2 and ID5 are the same. With the imposed boundary constraints, the log-likelihood functions can be expressed as the sum of a function of the parameters in the outcome model and a function of the parameters in the drop-out model. Hence, the parameters in the outcome and drop-out models can be maximized separately. Compared to the minimizations summarized in Tables B.2, B.4 and B.5, the convergence for these three cases is achieved with many fewer 75 Figure 7.2: A Two-Dimensional Profile Log-likelihood Surface for Model ID1 with Boundary Constraint 7703 —> — 0 0 and 7703 + 7723 = A 76 Table 7.3: Non-response Probability for the Third Response Using Model ID1 with 7703 ->• -oo and 7703 + 7723 = A y*2 y*s logit{Pr(i? = a | {p,p},y ,a:)} 3 3 0 0 0 1 A 1 0 —00 1 1 —00 TJ13 + A iterations. Further, the estimated Hessian matrices were never reset to a unit matrix during the course of minimization. As expected, the standard errors for 7702 and 7722 in Table 7.6 behave similarly as in Table 7.4. This is again verified (by the same approach) not to reflect a boundary solution. On the other hand, the standard errors for all the estimates in Tables 7.5 and 7.7 look quite reasonable. This feature of boundary solutions does not appear in models ID3 and ID6. For both models, the solutions obtained by the QN minimization are located in the interior of the parameter space. Different sets of starting values lead to the same parameter estimates and similar standard errors for the estimates, as shown in Tables B.3 and B.6. Even though the Hessian matrix was never reset to unity during the minimization process, the small discrepancy in the estimated SEs is expected due to the way the Hessian matrix is approximated. For these two models, the convergence is achieved between 17 and 21 iterations, which is much faster than for the models where the solutions are located on the boundary of the parameter space. This concludes the discussion concerning the existence of boundary solutions. • Results for the I D Models Now we examine if the treatment effects are sensitive to the form of the informative drop-out model based on the results presented in Tables 7.4, 7.5, B.3, 7.6, 7.7 and B.6. Our primary focus is on the treatment effects in the marginal model for the exacerbation rates even though treatment effects are also incorporated in 77 Table 7.4: Results for Model ID1 Evaluated on the Boundary: V03 + V23 = A Parameter Estimate A> 0i (LD) 02 (HD) 0.876 -0.028 -0.489 -0.122 0.206 0.200 0.195 0.074 -0.020 0.031 -0.136 -0.534 -0.113 -0.657 0.170 0.168 0.183 0.187 0.213 0.221 A 0.558 -1.548 0.409 0.347 ?702 »7l2 -3.360 0.140 2.218 0.417 1.860 2.615 ??01 -2.089 0.167 03 (time) «12 «13 «23 «123 ai a 2 "13 Neg. Loglik r SE 933.407 (# Iter = 25) 78 7703 —• - 0 0 and Table 7.5: Results for Model ID2 Evaluated on the Boundary: - 0 0 , 7703 + 772 = Ai and 7702 + 772 = A Parameter Po (LD) fa (HD) fa (time) Pi "12 "13 "23 "123 Q!2 »7i Ai A »7oi Neg. Loglik 2 2 7703 fa fa (LD) fa (HD) fa (time) "12 "13 "23 "123 "1 " 2 A V02 V22 V01 Neg. Loglik 7702 -> • Estimate SE 0.886 0.204 -0.017 0.195 -0.484 0.194 -0.118 0.074 -0.004 0.163 -0.010 0.161 0.173 -0.111 -0.511 0.177 -0.103 0.208 -0.649 0.217 0.286 0.275 -1.356 0.264 -1.499 0.258 -2.089 0.164 933.922 (# Iter = 20) Table 7.6: Results for Model ID4 Evaluated on the Boundary: »?03 + 7723 = A Parameter -> - 0 0 , Estimate SE 0.880 0.189 -0.024 0.190 -0.487 0.187 -0.120 0.071 -0.013 0.145 -0.022 0.137 -0.126 0.151 -0.524 0.152 -0.109 0.202 -0.654 0.214 -1.165 0.181 -3.819 3.002 2.464 3.217 -2.089 0.165 934.432 (# Iter = 27) 79 7703 —• —00 and Table 7.7: Results for Model ID5 Evaluated on the Boundary: - 0 0 , 7703 + 772 = A i and 7702 + 772 = Parameter A> Pi (LD) (HD) 03 (time) 02 «12 "13 "23 "123 ai 0:2 Ai A Vol Neg. Loglik 2 A 7703 —> — 0 0 , 7/02 —> 2 Estimate SE 0.202 0.886 -0.017 0.198 -0.484 0.192 -0.118 0.073 -0.004 0.162 -0.010 0.160 -0.111 0.172 -0.511 0.176 -0.103 0.211 -0.649 0.217 -1.165 0.182 -1.293 0.168 -2.089 0.165 934.473 (# Iter = 21) the association model. The structure of the ID drop-out model does not change the conclusions about the treatment effects in the marginal model. All six models conclude that the exacerbation rates in the LD and P L groups at any given time are not significantly different (approximate two-sided p-value > 0.62 based on 0i in each case). On the other hand, the exacerbation rate in the HD group is estimated to be significantly lower than in the P L group at all time points (two-sided p-value < 0.02 based on 02 in each case). The odds of experiencing exacerbations in the P L group are roughly 1.6 times higher than in the HD group. There is a weak suggestion of a linear decrease with time in the log odds of experiencing exacerbations under models ID1, ID2, ID4 and ID5 (two-sided p-value « 0.10 in each model), but the estimates of 03 in both ID3 and ID6 provide a strong indication of a linear decrease over time (two-sided p-values < 0.008). The conclusions regarding the treatment effects in the association model are similar. All six models indicate that the odds of having exacerbations at two 80 occasions or at all three occasions in the study are not significantly different between the L D and P L groups (two-sided p-values > 0.38 based on a\). But the models suggest that the odds in the HD group are significantly smaller than in the P L group (two-sided p-values < 0.004). Under models ID1, ID2, ID4 and ID5, the estimates of the intercept parameters, c*i2 and ai3, are fairly similar while 0:23 is slightly more negative. As would be expected, the estimate for the intercept in the 3-way association model is most negative. The situation is similar for models ID3 and ID6, although the estimates are slightly more negative. Note that the estimates for a i 2 , " 1 3 and « 2 3 are not very different, suggesting a possibility of a common intercept parameter for all the 2-way association models. However, the reduction to a model with the same intercept parameter for all 2- and 3-way association models may not seem reasonable since the estimate for a\23 is always quite different from the others. Further, we could explore explicitly whether the responses are positively or negatively associated by comparing the joint probabilities of the responses with those obtained under the independence assumption. If the joint probabilities are larger than the product of the marginal probabilities, then there is some positive dependence among the responses; otherwise, the responses are negatively correlated. See Chapter 8 for more details. We now consider selecting a parsimonious ID model to describe our data. Table 7.8 summarizes the negative log-likelihood and available degrees of freedom for all models listed in Table 7.1. Based on the LRT, the reduction from model ID1 to ID2 is permissible (p-value = 0.60), indicating the dependence on the previous and current observations is similar at time points 2 and 3. Using model ID2 as the base model and comparing to model ID3 examines whether the odds of dropping out (for the same history) change over time; that is, the hypothesis is 7703 = 7702 = 7701 = ??o- But the LRT statistic indicates this reduction is not reasonable (p-value = 0.03). Note that one can also assess the reduction from model ID1 directly to ID3, although this assessment is not as sensitive as the comparison between models ID2 and ID3. 81 Table 7.8: Negative Log-likelihood Values for Models in Table 7.1 Drop-out Mechanism ID RD CRD Model 1 2 3 4 5 6 1 2 3 1 2 Negative Log-likelihood 933.407 933.922 937.349 934.432 934.473 938.464 936.833 937.250 937.457 940.422 941.040 Degrees of Freedom (df) 25 27 29 27 28 30 27 28 30 29 31 The associated p-value is 0.096, indicating only fairly weak evidence against reducing from model ID1 to ID3. Thus, based on the more sensitive assessment, we conclude that model ID2 is the simplest permissible ID model among these three. To consider further model reductions, we next compare model ID2 to ID5. The LRT statistic suggests this reduction is reasonable (p-value = 0.29). The overall reduction from model ID1 to ID5 also agrees (p-value = 0.54). In model ID5, the drop-out probabilities do not depend on the last observed response, only on the last unobserved response. The further reduction from model ID5 to ID6 is not allowed (p-value = 0.02). We conclude that model ID5 is the simplest of these six informative dropout models that can be used to describe our annual data set. The two reduced models, ID2 and ID5, both fit the data adequately. For model ID2, G = 25.94 and 2 X = 23.81 on 27 degrees of freedom (p-values = 0.52 and 0.64 respectively). For 2 model ID5, G = 26.53 and X 2 2 = 24.09 on 28 degrees of freedom (p-value = 0.54 and 0.68 respectively). Note that all parameter estimates in the outcome model are the same for drop-out models ID2 and ID5. This phenomenon is induced by the 82 imposed boundary constraints mentioned earlier which allow separate maximizations for the parameters in the outcome and drop-out models. • Ignorable Drop-out Under the assumption of ignorable drop-out (either RD or CRD), the maximum likelihood estimates obtained by the QN minimization are in the interior of the parameter space. The results are summarized in Tables B.7 to B . l l . As expected, the parameter estimates in the measurement process are the same in all the RD and CRD models. Hence, the conclusions about the treatment effects in the marginal model for the exacerbation rates do not differ across the different specifications of these drop-out models. Only the HD group has a different effect on the exacerbation rates compared to the P L group (two-sided p-value « 0.01 based on /3 ); the odds of having exacerbations in the P L group are about 1.6 times 2 the odds in the HD group. There is a strong indication of a linear decrease over time in the log odds of having exacerbations (two-sided p-value « 0.001 based on ft). The treatment effects express themselves similarly in the association model. There are no apparent differences between the LD and P L groups in the odds of having exacerbations at two and three occasions (two-sided p-value « 0.40 based on di), but the HD and PL groups differ (two-sided p-value ~ 0.004 based on d ). The 2 intercept parameter estimates are quite similar, although slightly more negative, to those obtained under models ID3 and ID6. Again, the estimated values for a i , 2 ai3 and a 3 are reasonably similar, and the estimate for ai 3 is somewhat more 2 2 negative. This indicates a model which assumes a common intercept parameter for all the 2-way association models and a separate intercept parameter for the 3-way association may be reasonable for our data. We next consider selecting a simpler model among the three RD models. Based on the LRT, the model reduction from RD1 to RD2 is permissible (p-value — 83 0.36). One can also reduce model RD2 to RD3 (p-value = 0.81). The LRT statistic comparing model RD1 to RD3 also indicates the reduction to model RD3 is reasonable (p-value = 0.74). Thus, model RD3 is the simplest permissble model under the RD assumption. Similarly, if a CRD mechanism is assumed, model CRD2 can be used instead of CRD1 to describe our annual data (p-value = 0.54). • Types of Drop-out in the D a t a In the earlier part of this section, we determined that reductions from model ID1 to models ID2 and ID5 are permissible, with model ID5 being the simplest possible model among the six ID models considered. These three models can be used to examine whether the drop-out mechanisms in our data is ID, RD or CRD according to the classification by Little and Rubin (1987). To assess whether the drop-out occurred at random (RD), we can compare model ID1 to RD1. This comparison examines 7723 = 7722 = 0. The L R statistic of 6.85 (df = 2; p-value = 0.03) provides evidence against this reduction. As already established, it is reasonable to have common regression parameters describing dropout at the different time points (reduce from ID1 to ID2). Hence, the comparison between model ID2 and RD2 should provide a more sensitive assessment of our question. In this case, we investigate whether 772 — 0 and the result agrees with the previous assessment. (LR statistic = 6.66, df = 1; p-value = 0.01). The less sensitive comparison of model ID1 to RD2 also sugguests one should not reduce to the simpler model (LR statistic = 7.69, df = 3; p-value = 0.05). Thus, the data indicate that the drop-out did not occur at random. As reduction to an RD model is not allowed, presumably reduction to a CRD model will also not be allowed. For the sake of completeness, we perform various assessments to examine this. Model CRD1 can be compared to model ID1, ID2 and ID5 to examine the dependence between the drop-out and the outcome processes. The L R test comparing models ID1 and CRD1 clearly indicates the reduction is 84 not permissible (LR statistic = 14.03, df = 4; p-value = 0.007). The L R statistics for examining the reduction from model ID2 and ID5 to CRD1 are 13.00 (df = 2; p-value = 0.002) and 11.90 (df = 1; p-value < 0.001), respectively. As expected, the comparison to ID5 provides the strongest evidence. Thus, the data provide strong evidence against the hypothesis that the drop-out process is independent of the outcome process. According to these comparisons, one cannot reduce from the ID models to any of these RD and CRD models. We can thus confidently conclude that the drop-out process in our data is informative. 7.2.3 Summary We fitted six ID models and the maximum likelihood solutions for four of these models lie on the boundary of the parameter space. This phenomenon does not occur in the case where the drop-out mechanism is assumed to be ignorable. Based on L R tests, we conclude that the drop-out mechanism in our data is informative and model ID5 is determined to be the simplest possible model for our data. The treatment effects appear in both the marginal and association models. However, we focus primarily on the treatment effects in the marginal model. Under model ID5, the HD group has a lower rate of exacerbations compared to the P L group. The odds ratio of having exacerbations in the HD group relative to the P L group is estimated to be 0.62 and the corresponding approximate 95% confidence interval (CI) is (0.42, 0.90). The indication of a linear decrease in the odds of having exacerbations over time is quite weak; the approximate 95% CI for fa is (—0.26,0.03). The treatment effects in the association model convey a similar story: the odds of experiencing exacerbations at two occasions and at all three occasions in the L D group are not significantly different from the P L group, but these odds are clearly lower in the HD group. Interestingly, these conclusions are not very sensitive to the underlying drop- 85 out mechanisms for this data set. In particular, the parameter estimates (and standard errors) in the outcome model obtained with the ID assumption are fairly similar to those obtained with the ignorable drop-out assumptions. 7.3 Baker's Selection Model: Extensions of the Dropout Model In this section, we are interested in investigating the impact of different specifictions of the drop-out model on inferences concerning the treatment effects. The outcome model remains the same as in the previous section, and is coupled with the dropout models considered in Baker (1995); that is, COV + LOR + LUR, COV * L U R and L O R * LUR. Since the only covariates to be used are the treatment groups indicators, we replace COV with TRT throughout this section. We have established that models ID1, ID2 and ID5 can be used to describe our annual data but no reduction to the RD and CRD models is allowed. Model ID1 is of form LOR + LUR, with different parameters associated with each time of occurrence of the drop-outs. Models ID2 is obtained from model ID1 by assuming the regression parameters to be common at each time of occurrence of the dropouts, while model ID5 corresponds to the further assumption that the drop-out probabilities do not depend on LOR. In this section, we retain the feature of common regression parameters in all drop-out models considered. The three non-nested ID models considered for M r t - i , y ? I ',Vt) x a r e : 1. T R T * LUR: For t = 2,3 ( r _ i equal to {p,p} or {p}), we have t logit[ftt(rt-i,y? I x;ri )] t = not + m^D + n HD 2 + r LDy* +r HDy* , ]4 86 t ]5 t + r/ y * 3 t (7.11) and for t = 1 ( r _ i equal to { }), the model is t logit[/»i({ }>yi I w,Vi)] = Voi+mLD + r) HD; 2 (7.12) 2. T R T + LOR + LUR: For t = 2,3, the model is logit[/it(rt_i, y* | x; r) )] = r) + rj\LD + rj HD t t 0t 2 + mvi-i + ViVt, (7.13) and for t = 1, we have logit[/ii({ },yl | a;;»h)] = Vol + ViLD + mHD; (7.14) 3. LOR * LUR: For t = 2,3, the model is l o g i t [ M r - i , y*t I t »7t)] = Vot + mVt-i + VMt + mVt-iVt> (- ) 7 15 and for t = 1, we have logit[/ii({ },y\ | x;T7j] = 7701. (7.16) One can view these models as expansions of models ID2 and ID5. More specifically, all three drop-out models are expansions of model ID5. Further, models TRT * LUR and TRT + LOR + LUR can also be considered as expansions of model ID2. Hence we can compare these models to models ID2 or ID5 for examining the improvement of the fit with these more general models. The results are presented in next subsection. 7.3.1 Results Tables C . l to C.3 in Appendix C display detailed summaries of the results corresponding to the three extended drop-out models. For each drop-out model, we 87 report the starting values used to obtain the parameter estimates, the estimated standard errors, negative log-likelihood values, and the number of iterations required to achieve convergence. The phenomenon observed in models ID2 and ID5 can also be seen in these drop-out models. For the drop-out model TRT * LUR (see Table C . l ) , the same parameter estimates are obtained regardless of the starting values used except for the intercept parameters, associated with LUR values, and 773 (773). Parameters 7703 and 7703 7702 and 77025 and the parameter are estimated as large negative is estimated as large positive. Further the estimates of 7703 and 773 al- ways sum to -1.226, and 7702+773 = -1.350. Similarly, for model T R T + L O R + L U R (see Table C.2), the intercept parameters, 7703 and 7702, are estimated as being large negative, and the estimated value for 774 (the regression parameter corresponding to LUR) is large positive, but 7703 + 774 = —1.430 and 7702 + 774 = —1.573. The situation for model LOR * LUR is more complicated. Here we have the same phenomenon described for both models TRT * LUR and TRT + LOR + LUR, but the estimates of 771 (the parameter corresponding to LOR) and 773 (the parameter associated with the interaction term, LOR x LUR) also appear to satisfy the constraint, 771 + 773 = 0.286. The parameter estimates obtained from the fourth set of the starting values, in particular, indicate that the maximum likelihood solution corresponds to 771 —> — 00 with 771 + 773 = 0.286. To make comparison to models ID2 or ID5, we need to verify that the maximum likelihood solutions for these extended models occur at the suggested points on the boundary of the parameter space. Re-parameterizing in a similar fashion as previously, the conditional drop-out probabilities at years 2 and 3 can be expressed as in Table 7.9. We then substitute these expressions into the log-likelihood functions for the three models. To obtain the MLEs, we minimize the negative log-likelihood functions using the QN procedure. The results are reported in Tables 7.10 to 7.12. The minimizations reported in Tables C . l to C.3 required a large number of iterations for convergence and, in each case, the estimated Hessian matrix was 88 Table 7.9: Non-response Probability for the Second and Third Responses Model: T R T * L U R With r/02 + r/3 = A i , rj03 + m = A LOR LUR logit{Pr(/? = a | {p},y* ,x)} 0/1 0 —oo 0/1 1 2 logit{Pr(ii3 = a | {p,p},Y3,a;)}. 2 A i + (m + r) )LD + (r/ + m)HD 4 2 2 —oo A 2 + (r?i + 7] )LD + (r) + r) )HD 2 A b Model: T R T + L O R + L U R With 7/02 + 774 = A i , 7703 + 774 = A LOR LUR 0/1 0 0 1 1 1 logit{Pr(i? = a | {p},y* ,x)} 2 2 logit{Pr(i? = a | {p,p},yg,x)} 2 3 —00 • —00 A i + 771 + V2HD A i +rnLD + r) HD + T] 2 A + 771 L i ? + r) HD A + 771 L D + r) HD + 773 2 3 2 2 2 Model: L O R * L U R With 7702 + T? = A i , 7703 + 772 = A , 771 +773 = A 2 2 3 LOR LUR logit{Pr(i? = a | {p},yS,aj)} logit{Pr(i? = a | ( p , p } , y 3 , x ) } 0/1 0 —00 —00 0 1 Ai 1 1 Ai + A 2 3 A + A 2 A 3 2 3 reset to a unit matrix in the course of the computations. These features were not found for the minimizations reported in Tables 7.10 to 7.12. In particular, the number of iterations needed in Tables 7.10 to 7.12 is, on average, only one-third the number required in Tables C . l to C.3. Furthermore, the estimated Hessian matrix in Tables 7.10, 7.11 and 7.12 was never reset to unity throughout the minimization process. The parameter estimates and the log-likelihood values in the corresponding tables in these two sets are identical to the number of digits displayed, but the loglikelihood is always slightly larger at the boundary point than at the interior points located by the original minimizations. Hence, we have shown that the maximum likelihood solutions for these extended models are indeed located at the suggested points on the boundary. 89 Table 7.10: Results for Model TRT * L U R Evaluated on the Boundary: 7702 -> - 0 0 , 7702 + 773 = A i and 7703 +773 = Parameter /V 01 (LD) P2 (HD) 03 (time) "12 "13 "23 "123 "1 "2 V01 Vi(LD) V2(HD) Vi(LD x LUR) Vs(HD x LUR) Ai A Neg. Loglik 2 A 2 Estimate SE 0.204 0.886 -0.017 0.201 -0.484 0.198 -0.118 0.075 -0.004 0.165 -0.010 0.165 -0.111 0.180 0.182 -0.511 -0.103 0.213 -0.649 0.221 -2.136 0.293 -0.203 0.433 0.296 0.394 0.571 0.521 -0.620 0.518 -1.350 0.244 -1.226 0.251 931.223 (# Iter = 23) 90 7703 -> - 0 0 , Table 7.11: Results for Model TRT + L O R + L U R Evaluated on the Boundary: 7703 - > - 0 0 , 7702 ->• - 0 0 , 7702 + 774 = Parameter A) Pi (LD) P2 (HD) flz (time) "12 "13 "23 "123 "1 "2 7701 m(LD) m{HD) m(LOR) Ai A Neg. Loglik 2 Ai and 7703 + 774 = A 2 Estimate SE 0.886 0.198 -0.017 0.195 -0.484 0.190 -0.118 0.074 -0.004 0.157 -0.010 0.155 -0.111 0.169 -0.511 0.173 -0.103 0.207 -0.649 0.215 -2.156 0.223 -0.209 0.238 -0.023 0.249 0.290 0.202 -1.573 0.283 -1.430 0.242 933.350 (# Iter = 25) 91 Table 7.12: Results for Model L O R * L U R Evaluated on the Boundary: With 7/03 -> -oo, 7/02 -> -oo, T/I -» -oo, 7/ 2 + T/ = A i , 7/ 3 + T/ = A and 7/1 + 7/3 = A 0 Parameter Pi (LD) P2 (HD) 03 (time) "12 "13 "23 "123 "1 "2 7/01 Al A A Neg. Loglik 2 3 2 0 2 2 3 Estimate SE 0.886 0.206 -0.017 0.196 -0.484 0.194 -0.118 0.074 -0.004 0.163 -0.010 0.160 -0.111 0.172 -0.511 0.177 ^0.103 0.208 -0.649 0.217 -2.089 0.167 -1.499 0.265 -1.356 0.265 0.286 0.277 933.922 (# Iter = 21) There is an interesting point to note before moving on to the comparisons between these models and the models described in the previous section. Tables 7.10 to 7.12 (see also Tables C . l to C.3) display identical estimates for all the parameters in the outcome model. In fact, these parameter estimates are identical to those reported in Tables 7.5 and 7.7 (see also Tables B.2 and B.5) for models ID2 and ID5, respectively. The explanation for this is simple: for these drop-out models, the conditional probabilities that the second and third observations are missing are estimated to be zero when LUR = 0 (for both values of LOR). This simplifies the log-likelihood functions and allows the parameters in the outcome model and in the drop-out model to be maximized separately. As the 5 models share the same specification for the outcome process, it is then no surprise that the estimates of the parameters in the outcome model are identical even though the model specifications for the drop-out process differ. 92 Table 7.13: Results for Model TRT + L U R Evaluated on the Boundary: 7702 -> - 0 0 , 7702 + 774 = A i and 7703 + 774 = Parameter Po Pi {LD) P2 (HD) Pz (time) "12 «13 «23 0!123 Oil Ci2 VOl Vi(LD) V2(HD) Ai A Neg. Loglik 2 A 2 Estimate SE 0.886 0.207 -0.017 0.193 -0.484 0.194 -0.118 0.076 -0.004 0.163 -0.010 0.159 -0.111 0.171 -0.511 0.175 -0.103 0.206 -0.649 0.218 -2.136 0.214 0.191 0.232 -0.051 0.248 -1.349 0.211 -1.222 0.226 933.910 (# Iter = 27) 93 7703 ->• - 0 0 , To examine if one of these more complicated models should be employed for the drop-out process, we compare models TRT * LUR, TRT + L O R + LUR and LOR * LUR to models ID2 and ID5. By comparing model TRT * LUR to model ID5, we are examining whether the additional treatment effects (771,772) interaction between the treatment effects and the last unobserved response and the (774,775) provide a significant improvement on the fit of model ID5. The L R statistic (6.50 on df = 4; p-value = 0.16) indicates that there is not strong evidence that we should employ model TRT * LUR instead of model ID5. Table 7.10 suggests the two interaction terms contribute the major improvement in expanding the model from ID5 to TRT * LUR. Further, the comparisons of model ID5 with models TRT + LOR + LUR and LOR * LUR seem to agree with this observation (LR statistics = 2.25 and 1.10, df = 3 and 2; p-values = 0.52 and 0.58, respectively). That is, neither the terms LOR and TRT nor the terms L O R and LOR x LUR contribute significant improvement to the fit of model ID5. Thus model TRT + L U R (obtained by setting 774 = 775 = 0 in model TRT * LUR) is an interesting intermediate model between models ID5 and TRT * LUR. The detailed results for model TRT + LUR are provided in Table C.4, while the maximum likelihood estimates evaluated on the suggested point on the boundary of the parameter space is presented in Table 7.13. Comparing model TRT * LUR to TRT -I- L U R examines the contribution of the interaction terms, TRT x LUR. The corresponding LRT statistic is 5.37 on 2 degrees of freedom (p-value = 0.07), indicating fairly weak evidence against the hypothesis that the interaction terms are negligible. The cautious approach in this situation might be to retain the more general model, i.e. TRT * LUR, rather than reducing to the simpler TRT + LUR. But the evidence is not compelling, so we choose to reduce to the simpler TRT + LUR as the drop-out model. We then further examine whether the reduction from model TRT + LUR to ID5 is reasonable. Not surprisingly, in view of the earlier comparisons of model TRT * L U R to ID5, the LRT shows that the data provide no evidence to conclude 94 that the additional TRT covariates improve the fit of model ID5 (p-value = 0.57). We have already identified that models TRT + LOR + LUR and LOR * LUR do not improve the fit of model ID5. We can also examine whether these extended drop-out models provide improvements to model ID2 (LOR 4- LUR). The LR statistics are 1.15 and 0.00 (due to possible round-off error) on 2 and 1 degrees of freedom, respectively, indicating insufficient evidence to conclude that these extended drop-out models improve upon the fit of ID2 to our data set. Thus, neither the addition of TRT nor of LOR x LUR, provides a meaningful improvement in fit to ID2 (LOR + LUR). Hence, the simpler models ID2 or ID5 can be used to describe the drop-out process in our annual data set. 7.3.2 Summary We explored various ways of modelling the drop-out process in our data. More specifically, the three models considered can be viewed as extensions of ID2 and ID5, two of the permissible drop-out models described in the previous section. We introduce treatment effects and interactions terms into the drop-out model with a view to examining whether there is any impact on the conclusions about the treatment effects. Because some of the conditional drop-out probabilities at years 2 and 3 are estimated to be zero for each of these three drop-out models (see Table 7.9), the estimates of the parameters in the outcome model from these three dfop-out model specifications are identical to those obtained under models ID2 and ID5. It is also of interest to investigate whether a more general model specification for the drop-out process improves the fit. The results indicate that the simpler drop-out models ID2 or ID5 are adequate for our annual MS data. Thus, models ID2 and ID5 would be used throughout the next section. 95 7.4 Baker's Selection Model: Extension of the Outcome Model In this section, we explore extensions of the outcome model considered in the two previous sections based on including other baseline covariates such as gender, age, duration of MS, EDSS and BOD, in addition to the treatment arms and time. The main purpose of this section is to investigate whether or not inclusion of other baseline covariates in the model has any impact on the conclusions about the treatment effects identified in Section 7.2. For simplicity, we only consider the five baseline covariates described in Section 2.2.3, and these are introduced only into the marginal model for the exacerbation rates. The structure of the associations among the measurements is assumed to remain as previously described. This is thought reasonable as our primary interest focuses on the impact of additional covariates on the conclusions about the treatment effects in the marginal model for the exacerbation rates. The baseline covariates are included one at a time into the marginal component of the outcome model. The forward stepwise procedure for inclusion of the baseline covariates in addition to the treatment and time effects is carried out in the following fashion: (1) . Consider each covariate for inclusion in the marginal model and examine if it has a significant effect; (2) . If any covariates have significant effects, include the most significant covariate in the marginal model and repeat (1). Stop when no remaining covariates are found to be significant; (3) . If no covariates have significant effects, terminate the procedure. Even though EDSS score is an ordinal variable, for simplicity, we treat it as a continuous variable in our analysis. The BOD at baseline is skewed to the right .96 as is evident in Figure 2.5. Further, this covariate has a much larger scale than the other covariates. To avoid potential difficulties these features could induce in the estimating procedure, we use a logarithm transformation of the baseline BOD. Baseline BOD and its logarithm are highly associated (the correlation between them is roughly 0.7 based on the 362 patients who had baseline BOD greater than zero). Among the ID models considered with the original form of the outcome model, we found that the reduced models ID2 and ID5 were adequate. The extensions considered in Section 7.3 did not improve the fit significantly, so these same drop-out models will be considered here. The inclusion of additional covariates in the outcome model contemplated here could improve the overall fit, in which case it would again be of interest to examine whether the drop-out process is ID, RD or CRD. As noted in Section 7.2, model ID2 is more suitable for this purpose. Hence, model ID2 is used to describe the drop-out process throughout this section. 7.4.1 Results The results of the forward stepwise procedure to examine the role of each baseline covariate are summarized in Table 7.14. These log-likelihood values correspond to maximum likelihood estimates on the boundary of the parameter space as in the earlier fitting with models ID2 and ID5. Detailed summaries for the several cases reported in Table 7.14 appear in Tables D . l to D.5 of Appendix D. The minimization process for obtaining the estimates reported in Tables D . l to D.5 are similar to those described earlier. These maximum likelihood estimates were, on average, obtained at the 24th iteration and the Hessian matrix was never reset to a unit matrix in any of the minimizations. The first baseline covariate in addition to the treatment group included in the model is gender of the patients (Gender). The LRT indicates gender is not an important covariate when estimating the exacerbation rate. This agrees with Wald test (see Table D . l : z-score = 1.17, p-value = 0.24). The effects of baseline EDSS 97 Table 7.14: The LRT Statistics in the Forward Stepwise Procedure Case 1 2 3 4 5 Neg. Loglik for Model with LD + HD + time: 933.922 Additional COV Neg. Loglik LRT Gender 933.244 1.357 EDSS 933.901 0.043 Dur 0.154 933.768 Age 933.354 1.137 log(BOD) 933.088 1.668 log(BOD) 1.414 933.215 log(BOD) 933.083 1.677 log(BOD) 933.211 1.421 p-value 0.24 0.84 0.69 0.29 0.20 0.23 0.20 0.23 Comment Based Based Based Based on on on on Imputed Imputed Imputed Imputed Set Set Set Set 1 2 3 4 (EDSS), duration of MS at baseline (Dur), and age at baseline (Age), are similarly not significant; see Table 7.14. As mentioned before, there are 8 patients with missing BOD at baseline. In addition, 2 patients did not have any lesions at baseline, i.e. their baseline BOD value is zero. This creates a minor difficulty for converting baseline BOD to the log scale. Since the smallest non-zero baseline BOD value is 9, we impute a value between 0 and 9 for these 2 patients and perform a sensitivity analysis to determine whether the specific value chosen has any impact on the conclusion of our analysis. The arbitrary values chosen are 1.0 and 4.5. For the 8 patients who did not have any reading on BOD at baseline, one way to impute values for them is with the expectation-maximization (EM) algorithm, utilizing the other baseline covariates. For our purposes, it is sufficient to use the following values to fill in the 8 missing values and perform a sensitivity assessment: • the average of the log of the baseline BOD from 362 patients (excluding the 10 patients mentioned earlier), i.e. 7.085 (BOD = 1194.516); • the average of the log of the BOD from 364 patients (2 patients with zero 98 Table 7.15: Data sets used for assessing the sensitivity of the results when considering log(BOD) in addition to treatment group and gender as a covariate Data Set Imputed Set Imputed Set Imputed Set Imputed Set 1 2 3 4 In terms of BOD The 8 Patients The 2 Patients 1194.516 1.0 1194.516 4.5 1148.905 1.0 1158.439 4.5 • baseline BOD imputed to have a value 1.0), i.e. 7.047 (BOD = 1148.905); • the average of the log of the BOD from 364 patients data (2 patients with zero baseline BOD imputed to have a value 4.5), i.e. 7.055 (BOD = 1158.439). The four different combinations of values for imputing the 8 missing values and the 2 zero baseline BOD values are listed in Table 7.15. All four imputed data sets lead to a similar conclusion: log(BOD) is not a statistically important factor; see Table D.5 for the detailed results. Since the other baseline covariates are demonstrated to be not important for estimating the rate of exacerbations, we can also perform an alternative assessment for the significance of log(-BO-D). In particular, the 8 patients with missing baseline BOD are withheld from the analysis and the 2 patients with zero baseline BOD are imputed to have values of 1.0 and 4.5. The results evaluated on the boundary of the parameter space are displayed, in Table D.6. To perform a LRT, we re-fit model ID2 with this reduced data set; see Table D.7. The conclusion from this assessment remains the same as in the previous analyses. The L R statistics corresponding to the data sets with zero baseline BOD imputed as 1.0 and 4.5 are 1.83 and 1.56 on 1 degree of freedom (p-values = 0.17 and 0.21, respectively). The Wald-test for 04 also leads to the same conclusion (z-scores — 1.30 and 1.20, with p-values = 0.19 and 0.23, respectively). As expected, the parameter estimates associated with the drop-out process 99 are identical in Tables D . l to D.5. The reason is exactly as in the previous section. Because the conditional drop-out probabilities at years 2 and 3 are estimated to be zero, the log-likelihood functions in all five cases can be expressed as the sum of a function of the parameters for the outcome model and a function of the parameters for the drop-out model. Hence, the MLEs for the parameters in the two processes can be obtained separately. Since we employ the ID2 drop-out model in all five cases, the parameter estimates are expected to be identical. 7.4.2 Summary In the previous sections, the outcome model includes only the treatment groups and time as covariates. Here we consider also including the five baseline covariates, gender of the patients, EDSS, duration of disease, age and BOD, into the marginal model for the exacerbation rates. Model ID2 is used to described the drop-out process throughout the section. We found that none of these five baseline covariates contribute significantly to the fit in estimating the exacerbation rates. 7.5 Overall Summary for Baker's Selection Model We have used Baker's selection modelling approach to address various questions, and we provided a brief summary of our findings at the end of Sections 7.2, 7.3 and 7.4. In this section, we briefly describe what we have learned about the data according to the results obtained with the simplest acceptable model. In Section 7.2, we first determined that the non-saturated outcome model described in (7.2) — (7.4) is sufficient for our data by comparing it to various more general outcome models. This outcome model was then used throughout the section, coupled with drop-out models of the type L O R + LUR, to address questions of interest. We discovered that the maximum likelihood solutions for four (models ID1, ID2, ID4 and ID5) of the six informative drop-out models are located on the boundary of the parameter space. This results in identical parameter estimates for 100 the outcome model associated with drop-out models ID2 and ID5 for our data set. This boundary phenomenon does not arise in any of the ignorable drop-out models, i.e. the RD and CRD models. Based on likelihood ratio tests, we concluded the drop-out mechanism in our data set is informative. Models ID1, ID2 and ID5 are permissible and adequate models for modelling the drop-out process in our data. Model ID5, the simplest permissible informative drop-out model, indicates that the drop-out process in our data depends on the outcome process only through the last unobserved measurement (LUR). In Section 7.3, we explored several drop-out models that can be viewed as generalizations of models ID2 and ID5. In particular, we allowed the drop-out process to depend on the treatment groups. We found that these general drop-out models do not provide significant improvement to the fit of models ID2 or ID5. Thus, our drop-out process can be described by the simpler models ID2 and ID5. In Section 7.4, we addressed the question of the significance of other baseline covariates such as gender, EDSS, duration of MS, age and BOD in estimating the rate of exacerbations. These covariates were considered for inclusion only in the marginal component of the outcome model. Based on the forward stepwise procedure, none of these covariates were found to contribute significantly to the fit. Consequently, the simplest Baker's selection model consists of an outcome model composed of (7.2) — (7.4), and the drop-out process described by model ID5; see Table 7.7. This model fits the data quite adequately (p-value > 0.54). The observed and expected counts for the 15 (observation patterns) by 3 (treatment groups) contingency table are presented in Table 7.16. None of the expected cell counts are zero even though this model estimates some of the conditional probabilities of drop-out to be zero. The discrepencies between the observed and the expected counts are generally small, indicating the data are well-described by the model. Thus, we make inferences based on our data using this model in Chapter 8. 101 Table 7.16: The Observed and Expected Cell Counts for Baker's Selection Model with Drop-Out Model ID5 ("*" denotes missing) Pattern PL LD HD (0,0,0) 14 (13.5) 9 (9.1) 15 (15.3) (0,0,1) 3 (3.0) 5 (5.0) 11 (7.7) (0,1,0) 6 (5.2) 7 (7.4) 12 (10.3) (0,1,1) (1,0,0) 5 (6.4) 7 (6.3) 7 (5.3) 9 9 (13.9) 12 6 (8.6) (1,1,0) 8 (10.7) 10 11 (9.5) (10.2) 9 (1,0,1) (6.7) (10.2) (10.7) 13 (9.0) (1,1,1) (0,0,*) 25 (24.5) 18 (23.4) 16 (15.8) 0 2 (0.9) 1 1 (2.4) (2.0) (3.2) 3 5 0 2 (1.6) 1 (1.6) (2.0) (3.2) (1,1,*) (0, *, *) 11 (7.7) 10 (7.3) 3 (4.9) 2 (3.7) 7 (4.3) 4 (4.7) (1,*,*) (*,*, *) 12 (11.8) 12 (11.3) 8 (8.1) 13 (13.6) 11 (13.8) 17 (13.7) (0,1,*) (1,0,*) (2.7) Goodness-of-fit Tests G' = 26.53 on 28 degrees of freedom; p-value = 0.54 X> = 24.09 on 28 degrees of freedom; p-value = 0.68 1 102 7.6 The Liu et al. Transition Model In this section, we apply the Liu transition model to our annual data. Recall that Liu et al. (1999) employ a first-order transition model to model the outcome process. Further they assume that each of the conditional probabilities, Pr(Y * = y% | Y' !_ = yjL^x), does not depend on the covariates measured at time t t 1 t which seems somewhat unusual (see Chapter 5 for details). In our case, we can proceed with their idea without making such an assumption about the dependence on the covariates measured at time t as we consider only covariates measured at baseline. For the drop-out process, we consider three models: ID1, ID2 and ID3 as described in Table 7.1. Based on the LRT, we can select the simplest permissible model among the three. The basic idea of these models is similar to those considered in Liu et al. (1999) in the sense that the drop-out probabilities are assumed to depend only on the response observed prior the drop-out (LOR) and the response which would be observed if drop-out had not occurred (LUR). But in their data set, the first observation is always observed. Thus their models are slightly different than ours as they do not need a model for the case where the response pattern r _ i t is equal to { }. • Repeated Binary Outcomes with Informative Drop-out • o Outcome Model A first-order transition model is assumed for the binary longitudinal data. This means that the current measurement, yt, is related only to the previous measurement, y £ _ i , for t = 2,3, as well as to the baseline covariates of interest. Here only the treatment assignment and time are considered in the analysis since the results in the previous section indicate that gender of the patients, baseline EDSS, age at baseline, duration of MS, and baseline BOD were not important covariates in estimating the 103 rates of exacerbation. Thus, the outcome model employed can be expressed as: logit{Pr(Y7* = 1 | y / _ ! = yU,*t)} = A) + PiLD + faHD + 0 t + fay* _ 3 t x (7.17) o Drop-out Model Models similar to ID3 and ID6 from Table 7.1 were considered in Liu et al. (1999). Here we propose to model the drop-out process using models ID1, ID2, ID3 and ID5. We choose to focus on these three ID models out of the six listed in Table 7.1 because they allow straightforward investigation for the form of the drop-out mechanisms according to the terminology by Little and Rubin (1987). Furthermore, it will be interesting to determine if this leads to the same choice of the ID models for the drop-out process, namely ID2 and ID5, as the Baker selection model approach. • Repeated Binary Outcomes with Ignorable Drop-out • To investigate the impact of different drop-out mechanisms on the treatment effects, we also consider drop-out models assuming the drop-out occurred at random (RD) and completely at random (CRD). The RD and CRD models are the same in Table 7.1. ~ Likelihood ratio tests can be performed to examine the type of drop-out in our annual data based on these models. The results for the parameter estimates under different drop-out mechanisms are presented in the subsequent subsection. We conclude the section with a brief summary. 7.6.1 • Results Informative Drop-out The maximum likelihood solutions for models ID1 and ID2 lie on the boundary of the parameter space, while those for model ID3 exist in the interior. The detailed 104 results for these models are summarized in Tables E . l to E.6 of Appendix E. The boundary solutions for models ID1 and ID2 occur in a similar fashion as in Baker's selection model (see Tables E . l and E.2). We present the MLEs computed on the boundary for drop-out models ID1 and ID2 in Tables 7.17 and 7.18, respectively. These reported estimates are obtained with many fewer iterations than those in Tables E . l and E.2. Moreover, the estimated Hessian matrix in both cases was never reset to unity throughout the minimization process. Notice that the MLEs for the parameters in the outcome model are identical for drop-out models ID1 and ID2. This is again because some of the conditional probabilities of drop-out at years 2 and 3 are estimated to be zero and hence the parameters in the outcome and drop-out models can be estimated separately. The G and X 2 2 goodness-of-fit statistics shown in Table 7.20 provide some evidence of lack-of-fit in each case. Although the evidence is not compelling, the fit of these models for our data are somewhat questionable; perhaps a more complicated association structure or a more general drop-out model should be employed. However, our objective is not to perform a definitive analysis on our annual data, but rather to explore different approaches for modelling incomplete longitudinal binary data with informative drop-outs. Hence, despite their somewhat questionable fit, we do not elaborate on these models but rather go on to consider the best choices within this collection of models. All three ID models lead to similar conclusions about the treatment effects. In particular, the chance that an exacerbation would be experienced, given the past history (whether or not an exacerbation occurred at the previous time point), is not significantly different between the LD and P L groups (all p-values > 0.47). Nevertheless, the LD effect is estimated to be much stronger in model ID3 than in models ID1 and ID2. A l l three models conclude that the HD group has a lower chance than the P L group to experience an exacerbation, given the past history (p-values < 0.01). Further, there is a strong suggestion of a linear decrease over 105 Table 7.17: Results for Liu Transition Model with Drop-out Model ID1 Evaluated on the Boundary: 7703 -oo, 7702 Parameter Po Pi (LD) P2 (HD) •03 (time) Pi V01 V13 V12 Ax A Neg. Loglik 2 - 0 0 , 7703 + 7723 = A i and 7702 + 7722 = A 2 Estimate SE 1.007 0.206 -0.040 0.168 -0.462 0.167 -0.324 0.095 0.692 0.161 -2.089 0.165 0.558 0.413 0.374 0.048 -1.548 0.348 0.314 -1.327 942.259 (# Iter = 16) Table 7.18: Results for Liu Transition Model with Drop-out Model ID2 Evaluated on the Boundary: 7703,7702 - » - 0 0 , A i = 7703 + 772 and A = 7702 + 772 2 Parameter PQ (LD) P2 (HD) 03 (time) PI Pi V01 Vi Ai A Neg. Loglik 2 Estimate SE 1.007 0.199 -0.040 0.167 -0.462 0.167 -0.324 0.094 0.692 0.161 -2.089 0.169 • 0.286 0.296 -1.356 0.279 -1.499 0.273 942.687 (# Iter = 15) 106 Table 7.19: Results for Liu Transition Model with Drop-out Model ID5 Evaluated on the Boundary: 7703, V02 - » - 0 0 , A i = 7703 + V2 and A = 7702 + 772 2 Parameter A) 0i (LD) 0 (HD) 03 (time) 04 2 A A Neg. Loglik x 2 Estimate SE 1.007 0.206 -0.040 0.169 -0.462 0.168 -0.324 0.095 0.692 0.161 -2.089 0.166 -1.356 0.185 -1.499 0.168 943.239 (# Iter = 14) Table 7.20: Goodness-of-fit Statistics for Liu Transition Model with Drop-out Models ID1, ID2, ID3 and ID5 Model ID1 ID2 ID3 ID5 Degrees of Freedom 30 32 34 33 G' p-value 42.56 0.06 42.77 . 0.10 48.09 0.06 46.10 0.06 2 107 X 40.91 41.15 46.06 45.76 2 p-value 0.09 0.13 0.08 0.07 time in the log odds of having exacerbations given the past history (p-value < 0.001 based on fa in each model). The association parameter fa is highly significant (all p-values < 0.001). Under models ID1 and ID2, the odds of having an exacerbation given there was an exacerbation at the previous visit are 2.0 times the odds of having an exacerbation given there was no exacerbation at the previous visit; the corresponding approximate 95% CI for the odds ratio is (1.46,2.74). Under model ID3, the odds ratio is estimated as 1.8 and the approximate 95% CI is (1.27,2.48). The LR statistic for the reduction from model ID1 to model ID2 is 0.86 on 2 degrees of freedom (p-value = 0.65) and hence is permissible. However, we cannot further reduce model ID2 to model ID3 (LR statistic = 6.43, df = 2; p-value = 0.04). Thus, the simplest ID model among these three is ID2, which is the same conclusion obtained with Baker's selection model. Recall that with Baker's selection model, drop-out model ID5 is a reasonable reduction of model ID2. Thus, it is of interest to perform this assessment with the Liu transition model. The parameter estimates obtained from the Q N minimization with drop-out model ID5 are summarized in Table E.4. The results indicate a similar feature of boundary solutions as in model ID2. Table 7.7 presents the maximum likelihood estimates obtained at the suggested boundary points for model ID5. The LR statistic indicates that the term corresponding to the last observed response included in ID2 does not provide an important improvement to the fit (LR statistic = 1.10, df = 1; p-value = 0.29). Further, while the goodness-of-fit of model ID5 is slightly less satisfactory than for ID2 (see Table 7.20), the evidence against the adequacy of model ID5 is not overly compelling. These conclusions are qualitatively similar to those obtained with Baker's selection model. • Ignorable Drop-out The results for the RD and CRD models are displayed in Tables E.5 and E.6, respectively. As expected, the parameter estimates for the outcome model are identical 108 under both drop-out mechanisms. All parameter estimates are located in the interior of the parameter space. Under the assumption that the drop-out process is ignorable, the Wald tests suggest that the chance a patient would have an exacerbation given the past history is similar in the LD and P L group (p-value 0.50). But the risk is significantly lower in the HD group than in the P L group (p-value « 0.01). As in the ID case, the suggestion of a linear decrease over time in the log odds of having exacerbations given the past history is quite strong (z-score « —4.4 based on fa; p-value < 0.001). The odds of having an exacerbation given there was an exacerbation in the previous period are about 1.8 times the odds of having an exacerbation given there was no exacerbation in the previous period; the approximate 95% CI for the odds ratio is (1.32, 2.50). We perform LRTs for selecting the simplest RD and CRD models. The reduction from model RD1 to RD2 is permissible (p-value = 0.36), but the further reduction from model RD2 to RD3 is not allowed (p-value = 0.007). Under the CRD assumption, CRD1 is identified as the simplest possible model, as the reduction from CRD1 to CRD2 is not permissible (p-value = 0.004). These choices differ from those for Baker's selection model; see Section 7.2. • Types of Drop-out We established that models ID1 and ID2 are reasonable for describing our data. To investigate the types of drop-out in our annual data set, we compare these models with some RD and CRD models. For assessing if the drop-out mechanism is of type RD, model ID1 can be compared to model RD1 and similarly, model ID2 can be compared to model RD2. The reduction from ID1 to RD1 is permissible (LR statistic = 3.68, df = 2; pvalue = 0.16). However, the more sensitive assessment comparing model RD2 to ID2 (since the reduction from ID1 to ID2 is reasonable) provides a less definite 109 conclusion; the L R statistic equals 3.66 on 1 degree of freedom (p-value = 0.06). With a 5% level of significance, we would not reject the hypothesis that 772 = 0, but with only a slightly larger acceptable type I error, we would reject the hypothesis. Thus further investigation is required. The L R test indicates one cannot reduce from model RD2 to C R D 1 (LR statistic = 6.14, df = 1; p-value = 0.01). The reduction from model ID2 to C R D 1 is also not permitted (LR statistic = 9.80, df = 2; p-value = 0.007). Thus we need to make a decision based on the comparison between model ID2 and RD2. In such an ambiguous situation, one would usually prefer not to reduce from ID2 to RD2 because the simpler model may be more susceptible to potential bias in the results. As mentioned earlier, model ID2 can be further reduced to ID5. The comparison of model ID5 to C R D 1 confirms that model C R D 1 is not appropriate for our data (LR statistic = 8.70, df = 1; p-value = 0.003). Thus we conclude that the drop-out process in our data appears to be informative. 7.6.2 Summary We considered a first-order transition model for modelling the outcome process, coupled with the same drop-out models considered in Section 7.2. Based on the likelihood ratio tests, it appears that the drop-out process in our model cannot be ignored. Model ID5 is identified as the simplest drop-out model that is acceptable for our data. Based on this model, we computed the expected cell counts for the 15 (observation patterns) by 3 (treatment arms) contingency table; see Table 7.21. Despite some of conditional drop-out probabilities being estimated as zero, the expected counts are all nonzero. Notice that the differences between the observed and expected counts in some cells are quite large. For instance, the differences in cells (0,0,0) and (1,1,0) for the P L group and in cell (1,1,1) for the L D group are larger than 5.0 in magnitude. This is also reflected in the values of G 2 110 and X , 2 both J Table 7.21: The Observed and Expected Cell Counts for the Liu Transition Model with Drop-Out Model ID5 ("*" denotes missing) Pattern (0,0,0) (0,0,1) (0,1,0) (0,1,1) (1,0,0) (1,0,1) (1,1,0) (1,1,1) (0,0,*) (0,1,*) (1,0,*) (1,1,*) (o, *, *) (1,*,*) (*, *, *) PL 14 3 6 5 9 12 8 25 0 2 1 11 2 12 13 (7.4) (6,1) (5.8) (9.5) (9.3) (7.6) (14.4) (23.6) (1-6) (2.4) (2.0) (6.1) (3.9) (9.8) (13.6) LD 9 5 7 7 9 10 11 18 1 3 5 10 7 12 11 (8.1) (6.4) (6.1) (9.6) (9.7) (7.7) (14.6) (23.1) (1.6) (2.5) (2.0) (6.0) (4.1) (9.8) (13.8) 15 11 12 7 9 6 13 16 1 0 2 3 4 8 17 HD (15.6) (8.1) (8.3) (8.6) (13.2) (6.9) (14.0) (14.5) (2.1) (2.2) (1.8) (3.7) (4.3) (7.2) (13.7) Goodness-of-fit Tests & = 46.10 on 33 degrees of freedom; p-value = 0.06 X = 45.76 on 33 degrees of freedom; p-value = 0.07 2 111 indicating a potential lack-of-fit of this model. One could explore more complicated drop-out models or association structures to improve the fit of the model, but such extensions are not our main interest. Thus, we go on to make inferences based on our data using this model in the concluding chapter. 112 Chapter 8 Conclusions 8.1 Conclusions The main focus of this thesis has been on exploring likelihood-based methods for analyzing longitudinal binary responses under informative (or non-ignorable) dropout. The two modelling approaches, considered were Baker's selection model and the Liu et al. transition model. Both models belong to a general class of models known as selection models. A selection model factors the joint distribution for the response variables (Y) and the indicator variables denoting whether the response variables were observed (R) as /(Y,R) = /(R|Y)/(Y), (8.1) where / ( R | Y ) is the model for the drop-out process and / ( Y ) corresponds to the model for the measurement (or outcome) process. The main difference between Baker's selection model and the Liu transition model resides in the model specification for the measurement process. Baker's selection model uses a parameterization proposed by Ekholm (1991, 1992) to accomodate longitudinal binary measurements. That is, the outcome model is expressed in terms of a model for the (univariate marginal) probabilities of the responses and 113 an association model for the temporal associations among the responses. The L i u transition model, however, employs a first-order Markov chain transition model for the measurement process. The conditional distribution of response at time t (yt) given the history of the responses up to time t — 1 is assumed to depend only on the response at the previous time point (yt-i)- These outcome models are coupled with a drop-out model specified as a time-ordered causal model incorporating the assumption that the drop-out does not depend on future events. Given that the two approaches model the outcome process differently, this raises the question of the advantages and disadvantages of the two approaches. If the objective of the study is to study the effects of covariates on the marginal probabilities of the responses, marginal models provide a direct answer to this question. However, transition models should be used when the interest is in prediction (Diggle et al., 1994). Baker's selection model incorporates a more general structure for the strength of association among the responses than the L i u transition model. The structure for the associations among the responses in the L i u transition model is completely specified in terms of a single lagged effect. (Additional lagged effects could be added to the model but the nature of the association structure is limited by this parameterization.) For Baker's selection model, the expression for the outcome model for a sequence with more than three responses becomes more complicated, and the number of parameters increases rapidly. This is particularly so for the association model if no assumptions are made regarding the nature of the association structure among the responses. Unlike Baker's selection model, the number of parameters in the L i u transition outcome model need not change with the length of the response sequence. Both models were applied to our annual version of the Berlex exacerbation data described in Chapter 2 to examine the sensitivity of the estimated effects of Interferon /3-lb on the exacerbation rates in relapsing-remitting MS patients to various assumed forms for the drop-out mechanisms. More fundamentally, we were 114 Table 8.1: Estimated Chance of Exacerbations Based on Baker's Selection Model Treatment Group Year 1 Year 2 Year 3 PL 0.68 0.66 0.63 LD 0.68 0.65 0.63 HD 0.57 0.54 0.51 interested in studying the nature of the drop-out process in this clinical trial. Using Baker's selection modelling approach, we verified that the relationships expressed in (7.2) — (7.4) are sufficient for describing the outcome process in our data. This outcome model coupled with drop-out model ID5 is determined to be the most parsimonious yet adequate model among other more general models considered. In other words, the drop-out process in our data is informative and it depends on the last unobserved response, but not on the last observed response. Based on this model, we conclude that the low dose effect is not significant. The odds of having exacerbations in the L D group are reduced only by 1.7% relative to the odds of having exacerbations in the P L group. The corresponding approximate 95% CI for the precent reduction in the odds is (—44.9%, 33.3%). The high dose effect, however, is evidently different from the placebo effect. The odds of having exacerbations in the H D group are roughly 38.4% lower than the odds in the P L group (95% CI: 10.1%, 57.7%). Under the model assumption that the log odds of having exacerbations changes linearly over time, the odds are estimated to decrease by 11.1% per year in each group. The approximate 95% CI for the relative reduction in odds over time is (—2.6%, 23.0%), indicating the reduction is not statistically significant. The estimated chances of having exacerbations at each occasion presented in Table 8.1 also reflect these conclusions. The chances of experiencing exacerbations are almost the same in the L D and P L groups, but are much smaller in the H D group. In each group, these chances decrease only slightly over time. As for the association models, the L D and P L groups seem to have similar 115 Table 8.2: Estimated Chances of Exacerbations Based on the Liu et al. Transition Model Exacerbation Experienced in Previous Period Year 3 Treatment Group Year 1 Year 2 PL 0.74 0.67 0.80 LD 0.79 0.73 0.67 HD 0.64 0.71 0.57 No Exacerbation Experienced in Previous Period Year 3 Treatment Group Year 1 Year 2 PL 0.66 0.59 0.51 LD 0.66 0.58 . 0.50 HD 0.56 • 0.47 0.39 chances of having exacerbations at exactly two or all three time points during the study, but these chances are lower in the HD group. The odds ratios for the LD and P L groups are estimated as 0.90, reflecting a 9.8% reduction in the odds in the LD group. The corresponding approximate 95% CI is (—36.4%,40.3%), implying the LD effect is not statistically significant. On the other hand, the odds in the HD group are only about half the odds in the P L group. The approximate 95% CI for the decrease in the odds in the HD group is (20.1%, 65.8%). The estimates for the intercept parameters, a\2, and 0:23, are all quite small. This suggests a possible reduction to a model with all the 2-way associations in each treatment group being the same, i.e. au = 0:13 = 0 2 3 . On the other hand, a separate intercept parameter for the 3-way association appears to be useful as di23 is considerably larger in magnitude. Notice that, the estimated joint probabilities of the responses obtained from our model are slightly larger than those obtained under the independence assumption, indicating that there is some positive dependence among the responses; see Table 8.3. With the Liu et al. transition approach, the simplest acceptable drop-out model is also identified to be ID5, again indicating the drop-out mechanism in our 116 data is informative. Even though the outcome model, and hence the parameters being estimated, are different than in Baker's selection model, the conclusion regarding the treatment effects remain quite similar. For fixed t and previous response yjLj, the odds of having exacerbations are reduced by 3.9% (95% CI: —33.8%, 31.0%) in the LD group and by 37.0% (95% CI: 12.5%, 54.6%) in the HD group relative to the PL group. This indicates that only the high dosage of Interferon /3-lb effectively reduces the odds of experiencing exacerbations in MS patients. Similarly, the parameters fa and ^4 can also be interpretated as log odds ratios. In particular, exp(fa) represents the ratio of the odds of having exacerbations at time t + 1 as relative to time t for a patient with the same history at times t — 1 and t {yt-i = yt)• This odds ratio is estimated as 0.72 with approximate 95% CI (0.60, 0.87). The odds of having exacerbations given exacerbations in the previous period are 2.00 (= ex.p(fa)) times the odds given no exacerbations in the previous period; the corresponding 95% CI for the odds ratio is (1.46, 2.74). The estimated chances of experiencing exacerbations given the previous history presented in Table 8.2 also indicate similar conclusions regarding the treatment effects: the risks are much smaller in the HD group than in the LD and P L groups. Given that exacerbations were observed in the previous period (i.e. y^_ = 1), the x relative differences in the chances between the HD and PL groups are 11%, 14% and 15% at years 1, 2, and 3, respectively. For the case where no exacerbations were detected in the previous period (i.e. yl_\ = 0), the relative differences are slightly larger: 15%, 20% and 24% at years 1, 2, and 3, correspondingly. Table 8.3 displays the values of Pi(Y* = l,Y * = 1) where {s,t} = {1,2}, t {1,3}, {2,3} and Pr(Yj* = 1, Y * = 1, Y * = 1) obtained from Baker's selection model 2 z and the Liu et al. transition model. The estimates are generally similar for the two approaches except for the estimated probability of exacerbations at visits 1 and 3 and at all three visits. The differences are more substantial for the former estimated probabilities; the magnitudes of the (absolute) differences are 0.08, 0.06, 0.12 in the 117 Table 8.3: Estimated Pr(Y * = l,Y * = 1) and Pr(Y * = 1,Y * = 1,Y * = 1) by Treatment Groups s t x Baker's Selection Model PL LD Pr(Y * = l,y * = l) 0.50 0.47 Pr(Y * = l,Y * = l) 0.50 0.47 Pr(Y * = l,Y * = l) 0.47 0.45 Pr(Y * = l,Y * = l,Y * = l) 0.38 0.35 Assuming Independent Responses PL LD Pr(Y * = l,Y * = l) 0.45 0.44 Pr(Y * = l,Y * = l) 0.43 0.43 Pr(Y * = l,Y * = l) 0.42 0.41 Pr(Y * = l,Y * = l,y * = l) 0.28 0.28 2 3 HD 0.34 0.34 0.32 0.24 3 HD 0.31 0.29 0.28 0.15 Liu et al. Transition Model PL LD Pr(Y * = l,Y * = l) 0.49 0.48 0.42 0.41 Pr(Y *.= l,Y * = l) Pr(Y * = l,Y * = l) 0.47 0.45 Pr(Y * = l,Y * = l,Y * = l) 0.33 0.32 HD 0.36 0.28 0.32 0.20 1 2 1 3 2 3 1 2 1 2 1 3 2 1 3 2 1 2 1 3 2 1 3 2 3 3 PL, LD and HD groups respectively. In the intent-to-treat analyses reported in [35] (which assumed the drop-out occurred completely at random), the exacerbation rate was defined as the number of exacerbations experienced in one year. This is different from the exacerbation rate referred to throughout this thesis (the chance of having one or more exacerbations in a year). Nevertheless, it is of interest to compare the two sets of estimated treatment effects in terms of the relative change in the exacerbation rates. From the intent-to-treat analyses, the exacerbation rates in the PL, L D and HD group were 1.21, 1.05 and 0.84, respectively. Thus, the rates were 13% and 31% lower for the LD and HD groups relative to the P L group. Based on Baker's selection 118 model, the odds of having exacerbations are reduced by 1.7% and 38.4% in the LD and HD groups, respectively. Similarly, they are reduced by 3.9% and 37.0% under the Liu et al. transition model. The relative changes for the low dose effect are quite different between our approaches and the intent-to-treat analyses, but the variation is not as large for the high dose effect. Even though the magnitudes of the relative changes are quite different, the results convey a similar conclusion; that is, the effect of the high dosage of Interferon /3-lb is much more evident than that of the low dosage. We also found that there is a weak positive association over time in the presence/absence of exacerbations, and that the influence of the association is present over more than 1 time period. In the previous chapter, we provided the results from goodness-of-fit tests for both Baker's selection model and the Liu et al. transition model. The tests provided no evidence to suggest any lack-of-fit of Baker's selection model for our data. However, the adequacy of the Liu et al. transition model (p-values = 0.06 and 0.07 for G and X 2 2 respectively) is questionable. The discrepency between some of the observed and expected counts obtained from the Liu et al. model is quite large (see Table 7.21). This seems to suggest the restrictive assumption on the form of the associations among the responses in the Liu transition model may not be adequate for our data; that is, a higher-ordered transition model could possible be used instead. Alternatively, this may suggest a more general model for the drop-out process should be employed. Between Baker's selection model and the Liu transition model, Baker's selection model seems much more satisfactory as it fits the data quite well (see Table 7.20). In summary, analyses based on an assumption of ignorable non-response when the non-response mechanism is informative could lead to misleading results. By incorporating a non-response model in a likelihood-based approach, valid inferences can be obtained when the non-response mechanism is non-ignorable provided the non-response model correctly describes the non-response mechanism (Little and 119 Rubin, 1987). However, this approach is not without analytical difficulties. The parameters of the non-ignorable models may not be identifiable or the solutions to the likelihood equations (which may not be the maximum) may lie on the boundary of the parameter space. In Chapter 6, we showed that, with a saturated outcome model, the informative models of types C O V * L U R , C O V + L O R + L U R and L O R * L U R where C O V represents categorical covariates, are identifiable. In the course of our analyses in Chapter 7, we demonstrated that the maximum likelihood solutions for some of our non-ignorable models were located on the boundary of the parameter space. This boundary phenomenon did not occur in any of the ignorable non-response models considered. 8.2 • Further Work Other Approaches of Interest There are approaches other than selection models that can be used for analyzing incomplete data. In particular, the pattern-mixture modelling framework proposed by Little (1993) has become an area of active research. The pattern-mixture approach specifies the joint distribution of the measurement and response processes in terms of the marginal distribution of the responses multiplied by the distribution of measurements, conditional on the response patterns. Pattern-mixture models are natural when the interest is in population strata defined by missing data patterns, but these models are typically underidentified (Little, 1993). Thus the models require restrictions or prior information to identify the parameters. Unlike selection models, with the pattern-mixture approach one can avoid specifying the form of the missing data mechanism as it is incorporated indirectly via parameter restrictions (Little, 1993). This is a possible attractive feature over the selection model ap- proach as the latter is vulnerable to misspecification of the form of the missing-data mechanism. Further, pattern-mixture models are closer to the form of the data and sometimes simpler to fit. Thus, it would be of interest to re-analyze our annual data 120 with this approach and compare the results to those reported here. • Generalizations of the Data We chose to express the exacerbation data in terms of annual binary outcome variables. One could perform similar analyses on the binary data with more refined time intervals; for instance, semi-annual intervals. This semi-annual data may contain more information and may provide more precise estimates for the parameters. As mentioned at the outset, there is a loss of information associated with dichotomizing the data. To retain all the information, one could analyze the count data presented in Table 2.4 treating these as realizations of Poisson random variables [18, 19]. One could also use this approach with finer time-intervals, semi-annual intervals say. The conclusions obtained from these annual and semi-annual count data might be more informative than those based on the dichotomized data. 121 Bibliography [1] Baker, S.G. and Laird, N . M . (1988). Regression analysis for categorical variable with outcome subject to nonignorable nonresponse. Journal Statistical Association of the American 83, 62-69. [2] Baker, S.G. (1995). Marginal regression for repeated binary data with outcome subject to non-ignorable non-response. Biometrics 51, 1042-1052. [3] Broyden, C . G . (1970a). The convergence of a class of double-rank minimization algorithms, pt 1. Journal of the Institute of Mathematics and Its Applications 6, 76-90. [4] Broyden, C . G . (1970b). The convergence of a class of double-rank minimization algorithms, pt 2. Journal of the Institute of Mathematics and Its Applications 6, 222-231. [5] Dale, J . (1986). Global cross-ratio models for bivariate discrete ordered responses. Biometrics 42, 909-917. [6] Diggle, P.J. and Kenward, M . G . (1994). Informative drop-out in longitudinal data analysis. Applied Statistics 43, 49-93. [7] Ekholm, A . (1991). Fitting regression models to a multivariate binary response. In: A Spectrum nomics, of Statistical and Population Thought: Essays Genetics in Honour 122 in Statistical of Johan Fellman, Theory, Eco- G . Rosenqvist, K . Juselius, K . Nordstrom, J . Palmgren (eds), 19-32. Helsinki: Swedish School of Economics and Business Administration. [8] Ekholm, A . (1992). Discussion of: Multivariate regression analysis for categorical data by K . Liang, S.L. Zeger, and B. Qaqish. Journal Statistical Association of the American 81, 354-365. [9] Ekholm, A . (1998). The muscatine children's obesity data reanalysed using pattern mixture models. Applied Statistics 47, 251-263. [10] Fitzmaurice, G . M . and Laird, N . M . (1993). A likelihood-based method for analysing longitudinal binary responses. Biomeirika 80, 141-151. [11] Fitzmaurice, G . M . , Laird, N . M . and Zahner, E . P . (1996). Multivariate logistic models for incomplete binary responses. Journal Association of the American Statistical 91, 99-108. [12] Fletcher, R. (1970). A new approach to variable metric algorithms. The puter Journal Com- 13, 317-322. [13] Glonek, G . F . V . (1999). O n identifiability in models for incomplete binary data. Statistics & Probability Letters 41, 191-197. [14] Goodman, L . A . (1974). Exploratory latent structure analysis using both identifiable and unidentifiable models. Biometrika 61, 215-231. [15] Kenward, M . G . , Lesaffre, E . and Molenberghs, G . (1994). A n application of maximum likelihood and generalized estimating equations to the analysis of ordinal data from a longitudinal study with cases missing at random. Biometrics 50, 945-953. [16] Laird, N . M . (1988). Missing data in longitudinal studies. Statistics in 7, 305-315. 123 Medicine [17] Liang, K . Y . and Zeger, S.L. (1986). Longitudinal data analysis using generalized linear models. Biometrika 73, 13-22. [18] Lindsey, J . K . (1997). Applying Generalized Linear Models. Springer-Verlag, New York. [19] Lindsey, J . K . (1999). Models for Repeated Measurements. Oxford University Press, New York. [20] Little, R . J . A . and Rubin, D . B . (1987). Statistical Analysis with Missing Data. John Wiley, New York. [21] Little, R . J . A . (1993). Pattern-mixture models for multivariate incomplete data. Journal of the American Statistical Association 88, 125-134. [22] L i u , X . , Waternaux, C . and Petkova, E . (1999). Influence of human immunodeficiency virus infection on neurological impairment: an analysis of longtudinal binary data with informative drop-out. Applied Statistics 48, 103-115. [23] Michiels, B., Molenberghs, G . and Lipsitz, S.R. (1999). Selection models and pattern-mixture models for incomplete data with covariates. Biometrics 55, 978-983. [24] Molenberghs, G . , Kenward, M . G . and Lesaffre, E . (1997). The analysis of longitudinal ordinal data with non-random dropout* Biometrika [25] Molenberghs, G . , Goetghebeur, 84, 33-44. E . J . T . , Lipsitz, S.R. and Kenward, M . G . (1999). Nonrandom missingness in categorical data: strengths and limitations. The American Statistician 53, 110-118. [26] Nash, J . C . (1979). Compact Numerical and Function Minimisation. Methods for Computers: Adam Hilger L t d , Bristol. 124 Linear Algebra [27] Paty, D . W . , L i , D . K . B . , The U B C M S / M R I Study Group, and The IFNB Multiple Sclerosis Study Group (1993). Interferon B-lb is effective in relapsingremitting multiple sclerosis: II. M R I analysis results of a multicenter, randomized, double-blind, placebo-controlled trial. Neurology 4 3 , 662-668. [28] Robins, J . M . and Rotnitzky, A . (1995) Semiparametric efficiency in multivariate regression models with missing data. Journal Association of the American Statistical 9 0 , 122-129. [29] Rothenberg, T . J . (1971). Identification in parametric models. Econometrica 39, 577-591. [30] Rubin, D . B . (1976). Inference and missing data. Biometrika 6 3 , 581-592. [31] Schluchter, M . D . (1992). Methods for the analysis of informatively censored longitudinal data. Statistics in Medicine 11, 1861-1870. [32] Shanno, D . F . (1970). Conditioning of quasi-Newton methods for function minimization. Mathematics 24, 647-656. of Computation [33] Sun, W . and Song, P. (2000). Statistical analysis of repeated with informative cersoring times. Statistics in Medicine. measurements To appear. [34] Ten Have, T . R . , Kunselman, A . R . , Pulkstenis, E . P . and Landis, J.R. (1998). Mixed effects logistic regression models for longitudinal binary response data with informative drop-out. Biometrics 54, 367-383. [35] The I F N B Multiple Sclerosis Study Group (1993). Interferon B-lb is effective in relapsing-remitting multiple sclerosis: I. Clinical results of a multicenter, randomized, double-blind, placebo-controlled trial. Neurology 4 3 , 655-661. [36] The I F N B Multiple Sclerosis Study Group (1995). Interferon B-ib in the treatment of multiple sclerosis: final outcome of the randomized controlled trial. Neurology 4 5 , 1277-1285. 125 [37] Wu, M.C. and Carroll, R.J. (1988). Estimation and comparison of changes in the presence of informative right censoring by modelling the censoring process. Biometrics 44, 175-188. 126 Appendix A Proof for Condition (6.4) As in Section 6.2.2, there are two binary responses, Y\ and Y , with only Y subject 2 to non-response. j, k = The outcome model is Pr(Yi 0 , 1 . The non-response model, PT(R = 2 assumed to be homogeneous in that is, Y\\ j,Y = p | Y\ = pi . pij k = 2 = k 2 k j, Y2 \ X = = k, X i) = = i) for iVijk, = pij , k is Thus, the joint probabilities for the observed data are Pr(Yi = j , Y 2 = Pr(Yi = k, both observed | X = i) = 9 unobserved | X = i) = 0 j, Y2 ijk ijt = = itijkpik ir (l ij0 - p ) i0 + 7^1(1 - p ), ix and the marginal probabilities for Y\ are ITij. Let (f>ik = 1/pik = 7Tjj0 + Ttijl = + dijQ + 9ij\. and assume 1 = 2. The 0,^ must satisfy the following system of equations: V #100 #101 0 0 #110 #111 0 0 0 0 #200 #201 0 0 #210 #211 \ / ^ \ 127 010 \ f 7T10- 011 TTll- 020 7T20- 021 / \ 7T21- \ (A.l) / Given the multinomial probabilities 6, there is a unique solution for the fak provided the coefficient matrix is non-singular; that is, provided the determinant of the coefficient matrix does not equal to 0. The determinant of the coefficient matrix, ( 0 m 0 i o o — 0 i o i 0 n o ) ( 0 2 i i 0 2 o o — ^201^210)1 will be non-zero provided #1110100 - 0 0ioi#iio (A.2) and #211#200 - #2010210 7^ 0- (A-3) To satisfy (A.2), we require 01110100 01010110 TI'lOlPllTrilO/'lO TTlllPllTTlooPlO 7Tlll(7riO- - TTlOl) * 7!"lll/7I"ll- Pr(Y = l 1 y 2 1 = 7!"l0l(7Tll- - TTlll) TTlOl/^lO- l , X = l) P r ( y = 1 1 Y = 0,X = 1) 2 X Similarly, to satisfy (A.3) requires P r ( y = 1 I Y = 1,X = 2) ^ P r ( y = 1\Y 2 X 2 = 0,X = 2). 1 Thus the necessary and sufficient condition for the coefficient matrix to be nonsingular is P r ( y = 11 y 2 = i,x = i) ?Pr(y 2 = I\ Y =O,X 1 = i) for i = 1,2. Thus, the fak are identifiable unless this condition fails to hold. Note that, in contrast to the argument leading to condition (6.3), the argument leading to this condition remains the same if the number of levels of the categorical covariate X is greater than 2 (J > 2). 128 Appendix B Detailed Results for the Selection Models Described Section 7.2 129 Table B . l : Results for Model ID1 Parameter Po Pi P2 03 "12 "13 "23 "123 "1 "2 »703 Vl3 V23 V02 Vl2 V22 Voi Neg. Loglik Parameter 00 0i 02 03 "12 "13 "23 "123 "1 "2 V03 Vl3 V23 V02 Vl2 V22 V01 Neg. Loglik Set 1 SV Estimate (SE) 0.82 0.876 (0.919) 0.00 -0.028 (1.036) 0.00 -0.489 (0.896) -0.26 -0.122 (0.568) -0.60 -0.020 (0.579) -0.63 -0.031 (0.840) -0.77 -0.136 (0.959) -0.534 (0.702) -1.15 0.00 -0.113 (1.188) 0.00 -0.657 (0.938) -1.95 -14.421 (1.025) 0.00 0.558 (1.003) 0.00 12.874 (1.015) -1.95 -3.360 (1.042) 0.00 0.140 (1.001) 0.00 1.860 (1.030) -1.95 -2.089 (1.057) 933.407 ( # Iter = 71) Set 3 SV Estimate (SE) 0.90 0.876 (0.799) -0.03 -0.028 (0.266) 0.00 -0.489 (0.286) -0.12 -0.122 (0.394) -0.02 -0.020 (0.287) -0.04 -0.031 (0.25.1) -0.15 -0.136 (0.434) -0.50 -0.534 (0.292) 0.00 -0.113 (0.287) 0.00 -0.657 (0.310) -2.00 -15.400 (0.799) 0.56 0.558 (0.865) 0.00 13.853 (0.794) -2.40 -3.360 (0.806) 0.15 0.140 (0.783) 2.00 1.860 (0.813) 0.00 -2.089 (0.266) 933.407 ( # Iter = 59) 130 Set 2 Estimate (SE) 0.876 (0.816) -0.028 (0.584) -0.489 (0.357) -0.122 (0.388) -0.020 (0.404) -0.031 (0.378) -0.136 (0.486) -0.534 (0.446) -0.113 (0.656) -0.657 (0.391) -15.848 (0.784) 0.558 (0.769) 14.301 (0.775) -3.360 (0.690) 0.140 (0.540) 1.860 (0.787) -2.089 (0.563) ( # Iter = 70) Set 4 SV Estimate (SE) 0.876 0.876 (0.330) -0.028 -0.028 (0.338) -0.489 -0.489 (0.341) -0.122 -0.122 (0.079) -0.020 (0.302) -0.020 -0.031 -0.031 (0.311) -0.136 -0.136 (0.321) -0.534 -0.534 (0.334) -0.113 -0.113 (0.355) -0.657 -0.657 (0.371) -20.000 -15.171 (0.737) 0.558 0.558 (0.503) 1.000 13.624 (0.739) -3.360 -3.360 (0.800) 0.140 0.140 (0.828) 1.860 1.860 (0.813) -2.089 -2.089 (0.171) ' 933.407 ( # Iter = 91). SV 0.90 -0.02 -0.50 -0.12 -0.02 -0.03 -0.14 -0.50 -0.11 -0.66 -1.00 0.50 1.00 -2.00 0.14 2.00 -2.00 933.407 Table B.2: Results for Model ID2 Parameter ft 0i 02 ft C*12 ai3 "23 am ai 02 »703 »?1 V02 Vol Neg. Loglik Parameter ft ft ft . ft "12 "13 "23 «123 Oil Q-2 V03 Vl m V02 V01 Neg. Loglik Set 1 SV Estimate (SE) 0.82 0.886 (0.553) 0.00 - 0,017 (0.538) 0.00 - 0.484 (0.538) -0.26 - 0.118 (0.248) -0.60 - 0.004 (0.341) -0.63 - 0.010 (0.325) -0.77 - 0.111 (0.373) -1.15 - 0.511 (0.368) 0.00 - 0.103 (0.592) 0.00 - 0.649 (0.647) -1.95 -14.421 (0.888) 0.00 0.286 (0.815) 0.00 13.065 (0.708) -1.95 -14.564 (0.889) -1.95 - 2.089 (0.990) 933.922 (# Iter = 67) Set 3 SV Estimate (SE) 0.90 0.886 (0.788) -0.02 - 0.017 (0.734) -0.50 - 0.484 (0.751) -0.12 - 0.118 (0.324) 0.00 - 0.004 (0.562) -0.01 - 0.010 (0.694) -0.11 - 0.111 (0.593) -0.50 - 0.511 (0.599) -0.10 - 0.103 (0.792) -0.60 - 0.649 (0.847) -6.00 -14.384 (1.033) -0.30 0.286 (0.957) 6.00 13.028 (0.990) -4.00 -14.527 (0.904) -2.00 - 2.089 (0.958) 933.922 (# Iter = 56) 131 Set 2 Estimate (SE) 0.886 (0.738) - 0.017 (0.603) - 0.484 (0.700) - 0.118 (0.379) - 0.004 (0.487) - 0.010 (0.683) - 0.111 (0.848) - 0.511 (0.607) - 0.103 (0.713) - 0.649 (0.818) -14.760 (1.075) 0.286 (0.844) 13.404 (0.875) -14.903 (0.982) - 2.089 (0.995) (# Iter = 64) Set 4 SV Estimate (SE) 0.886 0.886 (0.530) -0.017 -0.017 (0.389) -0.484 -0.484 (0.377) -0.118 -0.118 (0.286) -0.004 -0.004 (0.370) -0.010 -0.010 (0.420) -0.111 -0.111 (0.546) -0.511 -0.511 (0.489) -0.103 -0.103 (0.441) -0.649 -0.649 (0.426) -14.384 -13.732 (0.616) 0.000 0.286 (0.318) 0.000 12.376 (0.584) 0.000 -13.875 (0.642) 0.000 -2.089 (0.838) 933.922 (# Iter = 72) SV 0.88 -0.02 -0.48 -0.12 0.00 -0.01 -0.11 -0.51 -0.10 -0.65 -1.95 0.00 0.00 -1.95 -2.00 933.922 Table B.3: Results for Model ID3 Parameter SV 00 0i 02 03 0.88 -0.02 -0.48 -0.12 0.00 -0.01 -0.11 -0.51 -0.10 -0.60 -1.95 0.00 0.00 «12 «13 «23 "123 a.\ vo Vi V2 Neg. Loglik Set 1 Estimate (SE) 0.986 -0.097 -0.475 -0.230 -0.082 -0.189 -0.345 -0.706 -0.191 -0.648 -2.195 0.416 0.222 (0.210) (0.199) (0.196) (0.083) (0.169) (0.177) (0.195) (0.200) (0.219) (0.226) (0.159) (0.286) (0.449) 937.349 (# Iter = 20) 132 SV 1.00 -0.10 -0.50 -0.20 -0.08 -0.20 -0.30 -0.70 -0.20 -0.60 -2.00 0.40 0.20 Set 2 Estimate (SE) 0.986 -0.097 -0.475 -0.230 -0.082 -0.189 -0.345 -0.706 -0.191 -0.648 -2.195 0.416 0.222 (0.210) (0.200) (0.195) (0.083) (0.171) (0.178) (0.198) (0.202) (0.219) (0.224) (0.161) (0.290) (0.459) 937.349 (# Iter = 21) Table B.4: Results for Model ID4 Parameter ft ft ft ft "12 "13 «23 ttl23 ai «2 V03 V23 V02 V22 VOl Neg. Loglik Parameter ft ft ft ft "12 «13 "23 «123 «i «2 »703 ?723 %2 %2 »7oi Neg. Loglik Set 1 SV Estimate (SE) SV 0.82 0.880 (0.731) 0.88 0.00 -0.024 (0.252) 0.00 0.00 -0.487 (0.241) -0.50 -0.12 -0.26 -0.120 (0.386) -0.013 (0.216) 0.00 -0.60 -0.022 .(0.225)" -0.63 -0.02 -0.77 -0.126 (0.422) -0.13 -1.15 -0.524 (0.257) -0.52 0.00 -0.109 (0.298) -0.10 0.00 -0.654 (0.274) -0.65 -1.95 -14.818 (0.730) -4.00 0.00 13.654 (0.751) 2.00 -1.95 -3.819 (0.728) -3.80 2.464 (0.761) 0.00 0.00 -1.95 -2.089 (0.209) -2.08 934.432 (# Iter = 67) 934.432 Set 3 SV Estimate (SE) SV 0.880 0.880 (0.880) 0.880 -0.024 -0.024 (0.720) -0.024 -0.486 -0.487 (0.703) -0.487 -0.120 -0.120 (0.434) -0.120 -0.010 -0.013 (0.438) -0.013 -0.022 -0.022 (0.771) -0.022 -0.125 -0.126 -0.126 (0.845) -0.524 -0.524 -0.524 (0.641) -0.109 -0.109 (0,776) -0.109 -0.654 (0.820) -0.654 -0.650 0.000 -15.432. (0.986) -15.432 0.000 14.268 (0.988) 14.268 -3.820 -3:819 (0.725) -3.819 2.464 (0.766) 2.464 2.460 -2.090 -2.089 (1.002) -2.089 934.432 (# Iter = 61) 934.432 133 Set 2 Estimate (SE) 0.880 (0.217) -0.024 (0.197) -0.487 (0.202) -0.120 (0.075) -0.013 (0.174) -0.022 (0.168) -0.126 (0.179) -0.524 (0.189) -0.109 (0.208) -0.654 (0.230) -16.284 (1.015) 15.119 (1.008) -3.819 (2.606) 2.464 (2.754) -2.089 (0.142) (# Iter = 63) Set 4 Estimate (SE) 0.880 (0.170) -0.024 (0.165) -0.487 (0.171) -0.120 (0.074) -0.013 (0.118) -0.022 (0.121) -0.126 (0.139) -0.524 (0.131) -0.109 (0.174) -0.654 (0.195) -15.432 (1.363) 14.268 (1.368) -3.819 (1.973) 2.464 (2.093) -2.089 (0.156) (# Iter = 24) Table B.5: Results for Model ID5 Parameter So Si S2 S3 "12 "13 "23 "123 "1 "2 V03 m V02 V01 Neg. Loglik Parameter •Bo Si S2 S3 "12 "13 "23 "123 "1 "2 Voz m V02 V01 Neg. Loglik Set 1 SV Estimate (SE) 0.82 0.886 (0.926) 0.00 -0.017 (0.898) -0.484 (0.874) 0.00 -0.26 -0.118 (0.543) -0.60 -0.004 (0.971) -0.63 -0.010 (0.965) -0.77 -0.111 (0.990) -0.511 (0.891) -1.15 0.00 -0.103 (0.901) 0.00 -0.649 (0.940) -1.95 -13.737 (1.000) 12.573 (1.002) 0.00 -1.95 -13.866 (1.000) -2.089 (1.000) -1.95 934.473 (# Iter = 55) • Set 3 SV Estimate (SE) 0.886 0.886 (0.213) -0.017 -0.017 (0.206) -0.484 (0.192) -0.484 -0.118 (0.074) -0.118 -0.004 -0.004 (0.171) -0.010 -0.010 (0.167) -0.111 -0.111 (0.178) -0.511 -0.511 (0.185) -0.103 -0.103 (0.220) -0.649 (0.214) -0.649 -15.608 -15.608 (0.582) 14.443 14.443 (0.568) -15.737 -15.737 (0.569) -2.089 -2.089 (0.166) 934.473 (# Iter = 20) 134 Set 2 SV Estimate (SE) 0.886 (0.169) 0.890 -0.017 -0.017 (0.178) -0.484 (0.187) -0.480 ^0.120 -0.118 (0.070) -0.004 (0.134) 0.000 -0.010 -0.010 (0.140) -0.110 -0.111 (0.153) -0.510 -0.511 (0.151) -0.100 -0.103 (0.195) -0.649 (0.204) -0.650 -4.000 -15.608 (301.091) 14.443 (301.092) 0.000 -2.000 -15.737 (301.091) -2.089 (0.166) -2.080 934.473 (# Iter = 60) Set 4 SV Estimate (SE) 0.900 0.886 (0.400) -0.017 -0.017 (0.573) -0.484 (0.482) -0.480 -0.118 (0.214) -0.120 -0.004 (0.305) 0.000 -0.010 -0.010 (0.306) -0.110 -0.111 (0.456) -0.510 -0.511 (0.460) -0.100 -0.103 (0.619) -0.649 (0.534) -0.650 0.000 -13.737 (0.942) 12.573 (0.587). 0.000 0.000 -13.866 (0.896) -2.089 (0.577) -2.000 934.473 (# Iter = 71) Table B.6: Results for Model ID6 Parameter So Si ' Si s 3 • "12 "13 "23 "123 Oil "2 no. m Neg. Loglik Set 1 SV Estimate (SE) 0.90 0.962 (0.209) -0.10 -0.080 (0.197) -0.50 -0.483 (0.198) -0.20 -0.201 (0.077) -0.05 -0.057 (0.171) -0.10 -0.137 (0.172) -0.20 -0.279 (0.187) -0.50 -0.646 (0.193) -0.10 -0.172 (0.213) -0.60 -0.655 (0.223) -1.95 -2.206 (0.151) 0.661 (0.264) 0.00 938.464 (# Iter = 19) Set 2 SV Estimate (SE) 1.00 0.962 (0.208) -0.08 -0.080 (0.198) -0.50 -0.483 (0.195) -0.20 -0.201 (0.078) -0.06 -0.057 (0.168) -0.10 -0.137 (0.170) -0.30 -0.279 (0.186) -0.60 -0.646 (0.190) -0.20 -0.172 (0.216) -0.70 -0.655 (0.224) -2.00 -2.206 (0.161) 0.70 0.661 (0.269) 938.464 (# Iter = 17) Table B.7: Results for Model RD1 Parameter So Si s 2 S3 "12 "13 "23 "123 "1 "2 V03 ni3 no2 nn noi Neg. Loglik Set 1 SV Estimate (SE) 0.82 0.999 (0.211) 0.00 -0.106 (0.196) 0.00 -0.470 (0.196) -0.26 -0.246 (0.080) -0.60 -0.097 (0.166) -0.63 -0.219 (0.169) -0.77 -0.384 (0.182) -1.15 -0.742 (0.188) 0.00 -0.201 (0.216) 0.00 -0.643 (0.229) -1.95 -2.416 (0.335) 0.00 0.878 (0.396) -1.95 -2.117 (0.300) 0.00 0.401 (0.360) -1.95 -2.089 (0.165) 936.833 (# Iter = 25) 135 Set 2 SV Estimate (SE) 1.00 0.999 (0.216) -0.12 -0.106 (0.194) -0.47 -0.470 (0.194) -0.20 -0.246 (0.082) -0.10 -0.097 (0.173) -0.22 -0.219 (0.172) -0.38 -0.384 (0.188) -0.74 -0.742 (0.195) -0.20 -0.201 (0.217) -0.64 -0.643 (0.227) -2.41 -2.416 (0.336) 0.90 0.878 (0.396) -2.11 -2.117 (0.293) 0.40 0.401 (0.350) -2.08 -2.089 (0.164) 936.833 (# Iter = 23) Table B.8: Results for Model RD2 Parameter ft ft 02 ft "12 «13 «23 "123 ai OL2 V03 m V02 Vol Neg. Loglik Set 1 SV Estimate (SE) 0.82 0.999 (0.204) 0.00 -0.106 (0.193) 0.00 -0.470 (0.187) -0.26 -0.246 (0.080) -0.60 -0.097 (0.166) -0.63 -0.219 (0.162) -0.77 -0.384 (0.174) -1.15 -0.742 (0.183) 0.00 -0.201 (0.217) 0.00 -0.643 (0.229) -1.95 -2.239 (0.258) 0.00 0.625 (0.261) -1.95 -2.278 (0.250) -2.08 -2.089 (0.169) 937.250 (# Iter = 26) Set 2 SV Estimate (SE) 0.999 0.999 (0.208) -0.106 -0.106 (0.197) -0.470 -0.470 (0.196) -0.246 -0.246 (0.077) -0.097 -0.097 (0.168) -0.219 -0.219 (0.170) -0.384 -0.384 (0.186) -0.742 -0.742 (0.193) -0.201 -0.201 (.0.217) -0.643 -0.643 (0.227) -2.239 -2.239 (0.247) 0.625 0.625 (0.264) -2.278 -2.278 (0.253) -2.089 -2.089 (0.167) 937.250 (# Iter = 17) Table B.9: Results for Model RD3 Parameter ft ft ft ft ai2 "13 "23 "123 ai Vo Vi Neg. Loglik Set 1 SV Estimate (SE) 0.82 0.999 (0.206) 0.00 -0.106 (0.197) 0.00 -0.470 (0.179) -0.26 -0.246 (0.078) -0.60 -0.097 (0.167) -0.63 -0.219 (0.168) -0.77 -0.384 (0.180) -1.15 -0.742 (0.187) 0.00 -0.201 (0.227) 0.00 -0.643 (0.219) -1.95 -2.153 (0.133) 0.00 0.518 (0.196) 937.457 (# Iter = 22) 136 Set 2 SV Estimate (SE) 0.999 0.999 (0.208) -0.106 -0.106 (0.195) -0.470 -0.470 (0.196) -0.246 -0.246 (0.073) -0.097 -0.097 (0.168) -0.219 -0.219 (0.163) -0.384 -0.384 (0.181) -0.742 -0.742 (0.191) -0.201 -0.201 (0.217) -0.643 -0.643 (0.225) 0.000 -2.153 (0.121) 0.000 0.518 (0.188) 937.457 (# Iter = 24) Table B.10: Results for Model CRD1 Parameter Bo Bi B Bs 2 "12 "13 "23 "123 ai a 2 V03 V02 Vol Neg. Loglik Set 1 SV Estimate (SE) 0.82 0.999 (0.203) 0.00 -0.106 (0.173) 0.00 -0.470 (0.189) -0.26 -0.246 (0.078) -0.60 -0.097 (0.172) -0.63 -0.219 (0.169) -0.77 -0.384 (0.186) -1.15 -0.742 (0.192) 0.00 -0.201 (0.196) 0.00 -0.643 (0.231) 0.00 -1.846 (0.161) 0.00 -1.849 (0.152) 0.00 -2.089 (0.154) 940.322 (# Iter = 29) Set 2 SV Estimate (SE) 0.999 (0.212) 1.00 -0.10 -0.106 (0.199) -0.40 -0.470 (0.197) -0.20 -0.246 (0.080) 0.00 -0.097 (0.169) -0.20 -0.219 (0.171) -0.40 -0.384 (0.186) -0.70 -0.742 (0.193) -0.20 -0.201 (0.221) -0.60 -0.643 (0.229) -2.00 -1.846 (0.170) -2.00 -1.849 (0.159) -2.10 -2.089 (0.165) 940.322 (# Iter = 21) Table B . l l : Results for Model CRD2 Parameter Bo Bi B B 2 3 «12 . «13 «23 "123 ai • a Vo Neg. Loglik 2 Set 1 SV Estimate (SE) 0.82 0.999 (0.281) 0.00 -0.106 (0.211) 0.00 -0.470 (0.228) -0.26 -0.246 (0.078) -0.60 -0.097 (0.232) -0.63 -0.219 (0.221). -0.77 -0.384 (0.227) -1.15 -0.742 (0.250) 0.00 -0.201 (0.228) 0.00 -0.643 (0.239) -1.95 -1.933 (0.106) 941.040 (# Iter = 21) 137 Set 2 SV Estimate (SE) 1.00 0.999 (0.210) -0.10 -0.106 (0.192) -0.40 -0.470 (0.193) -0.20 -0.246 (0.079) 0.00 -0.097 (0.169) -0.20 -0.219 (0.172) -0.40 -0.384 (0.185) -0.70 -0.742 (0.190) -0.20 -0.201 (0.217) -0.60 -0.643 (0.224) -2.00 -1.933 (0.095) 941.040 (# Iter = 19) Appendix C Detailed Results for the Selection Models Described Section 7.3 138 Table C . l : Results for Drop-out Model: TRT * L U R Parameter ft Pi(LD) fo{HD) 03 (time) "12 "13 "23 "123 "i "2 V03 V02 Voi Vi(LD) V2(HD) Vs(LUR) r] (LD x LUR) m(HD x LUR) Neg. Loglik 4 Set 1 SV Estimate (SE) 0.80 .0.886 (1.178) -0.10 -0.017 (0.779) -0.50 -0.484 (0.876) -0.20 -0.118 (0.422) -0.004 (1.143) -0.08 -0.20 -0.010 (1.025) -0.30 -0.111 (0.877) -0.70 -0.511 (0.617) -0.20 -0.103 (1.096) -0.70 -0.649 (0.885) -1.95 -13.608 (1.053) -1.95 -13.732 (1.069) -1.95 -2.136 (1.235) 0.00 -0.203 (1.232) 0.00 0.296 (1.044) 12.382 (1.122) 0.00 0.00 0.571 (1.001) 0.00 -0.620 (1.101) 931.223 (# Iter = 60) 139 Set 2 Estimate (SE) 0.90 0.886 (0.926) -0.02 -0.017 (0.988) -0.484 (1.095) -0.50 -0.10 -0.118 (0.911) -0.004 (1.365) 0.00 -0.01 -0.010 (1.198) -0.10 -0.111 (1.265) -0.50 -0.511 (1.193) -0.10 -0.103 (0.987) -0.60 -0.649 (1.143) -1.00 -15.118 (1.047) -1.00 -15.242. (1.044) -2.10 -2.136 (1.026) -0.20 -0.203 (1.020) 0.30 0.296 (1.133) 13.892 (1.348) 1.00 0.60 0.571 (1.025) -0.60 -0.620 (1.100) 931.223 (# Iter = 65) . sv Table C.2: Results for Drop-out Model: T R T + L O R + L U R Parameter SV Bo Bi(LD) B {HD) 8z(time) 0.80 -0.10 -0.50 -0.20 «12 -0.08 -0.20 -0.30 -0.70 -0.20 -0.70 2 "13 «23 t*123 Oil Oi2 V03 V02 Vol Vi(LD) (HD) m V3(LOR) rj^(LUR) Neg. Loglik Set 1 Estimate (SE) 0.886 (0.602) -0.017 (0.687) -0.484 (0.650) -0.118 (0.267) -0.004 (0.561) -0.010 (0.899) -0.111 (0.715) -0.511 (0.572) -0.103 (0.732) -0.649 (0.774) SV 0.90 -0.02 -0.50 -0.10 0.00 -0.01 -0.10 -0.50 -0.10 -0.60 Set 2 Estimate (SE) 0.886 (1.280) -0.017 (1.112) -0.484 (1.012) -0.118 (0.683) -0.004 (1.067) -0.010 (1.112) -0.111 (1.089) -0.511 (1.104) -0.103 (0.978) -0.649 (1.188) -1.95 -1.95 -13.920 (1.010) -14.063 (0.998) -2.00 -1.00 -14.006 (1.019) -14.149 (1.368) -1.95 -2.10 0.00 0.00 -2.156 (0.923) 0.209 (0.934) -0.023 (0.974) -2.156 (1.102) -0.209 (1.287) 0.00 0.00 0.290 (0.873) 12.490 (0.777) 933.350 (# Iter = 65) 140 0.20 -0.20 0.30 1.00 -0.023 (1.381)' 0.290 (1.145) 12.576 (2.942) 933.350 (# Iter = 60) Table C.3: Results for Drop-out Model: LOR * LUR Set 1 SV Estimate (SE) 0.886 (0.664) 0.80 ft -0.10 Pi(LD) ' -0.017 (0.411) -0.50 -0.484 (0.639) f3 (time) -0.20 -0.118 (0.361) -0.08 -0.004 (0.448) "12 -0.20 -0.010 (0.539) "13 -0.30 -0.111 (0.711) "23 -0.70 -0.511 (0.643) "123 -0.20 -0.103 (0.456) "1 -0.70 -0.649 (0.751) "2 -1.95 -12.674 (0.953) V03 -1.95 -12.817 (0.936) V02 -1.95 -2.089 (0.979) VOI Vi(LOR) -0.10 -0.711 (0.850) V2(LUR) 0.00 11.318 (0.785) 0.20 m{LOR x LUR). 0.996 (0.852) Neg. Loglik 933.922 (# Iter = 61) Set 2 SV Estimate (SE) 0.90 0.886 (0.650) -0.02 -0.017 (0.733) -0.50 -0.484 (0.800) -0.10 -0.118 (0.188) 0.00 -0.004 (0.487) -0.01 -0.010 (0.766) -0.10 -0.111 (0.586) -0.50 -0.511 (0.610) -0.10 -0.103 (0.788) -0.60 -0.649 (0.847) -1.00 -14.340 (0.992) -1.00 -14.483 (1.001) -2.10 -2.089 (0.998) 0.10 0.152 (1.002) 1.00 12.984 (0.995) 0.20 0.134 (0.998) 933.922 (# Iter = 66) Set 3 SV Estimate (SE) 0.80 0.886 (0.908) -0.10 -0.017 (0.882) -0.50 -0.484 (0.963) -0.20 -0.118 (0.593) -0.08 -0.004 (0.956) -0.20 -0.010 (0.945) -0.30 -0.111 (0.887) -0.70 -0.511 (0.884) -0.20 -0.103 (0.918) -0.70 -0.649 (0.985) -1.95 -14.338 (1.002) -1.95 -14.481 (1.004) -1.95 -2.089 (1.002) -0.50 -0.496 (1.001) 12.982 (1.009) -0.10 0.10 0.782 (1.001) 933.922 (# Iter = 60), Set 4 SV Estimate (SE) 0.90 0.886 (0.704) -0.02 -0.017 (0.722) -0.484 (0.628) -0.50 -0.10 -0.118 (0.345) -0.004 (0.710) 0.00 -0.01 -0.010 (0.806) -0.10 -0.111 (0.812) -0.50 -0.511 (0.705) -0.10 -0.103 (0.787) -0.60 -0.649 (0.793) -4.00 -13.092 (0.922) -3.00 -13.235 (1.002) -2.10 -2.089 (0.992) 0.10 -3.387 (0.916) 2.00 11.736 (0.861) 3.672 (0.929) -0.20 933.922 (# Iter = 59) Parameter 3 Parameter ft Pi(LD) (3 (HD) 03 (time) 2 "12 . "13 "23 "123 "1 "2 V03 V02 VOI Vi(LOR) V2(LUR) m(LOR x LUR) Neg. Loglik 141 Table C.4: Results for Drop-out Model: T R T + L U R Parameter SV Set 1 Estimate (SE) SV Set 2 Estimate (SE) Bo Bi(LD) 8 (HD) Bz(time\ 0.80 -0.10 -0.50 -0.20 0.886 -0.017 -0.484 -0.118 (0.202) (0.195) (0.192) (0.074) 0.80 -0.02 -0.50 -0.10 0.886 -0.017 -0.484 -0.118 (0.927) (0.995) (0.995) (0.600) Ct\2 -0.08 -0.20 -0.30 -0.70 -0.20 -0.70 -0.004 -0.010 -0.111 -0.511 -0.103 -0.649 (0.157) (0.156) (0.168) (0.171) (0.207) (0.215) 0.00 -0.01 -0.10 -0.50 -0.10 -0.60 -0.004 -0.010 -0.111 -0.511 -0.103 -0.649 (0.988) (0.981) (0.981) (0.976) (1.000) (0.997) -1.95 -1.95 -1.95 0.00 0.00 0.00 -2.140 -15.364 -15.236 0.191 -0.051 14.014 (0.151) (0.862) (0.807) (0.205) (0.219) (0.876) -2.10 -2.00 -3.00 -2.140 -13.697 -13.569 0.191 -0.051 12.347 (1.001) (1.001) (1.001) (1.002) (1.002) (1.021) 2 «13 C*23 "123 oti a.2 Vol V02 V03 m(LD) m{HD) r)i(LUR) Neg. Loglik 933.910 (# Iter = 67) 142 0.20 -0.05 1.00 933.910 (# Iter = 54) Appendix D Detailed Results for the Selection Models Described Section 7.4 143 Table D . l : Results for Case 1 in Table 7.14 Evaluated at the Boundary: 7/03 -> —00, 7702 -> - 0 0 , 7 / 0 3 + m = A i , 7/02 + m = Parameter & /?3 (time) 04 (gender) «12 "13 «23 "123 «1 «2 r/i(LOi?) Ai • A *7oi Neg. Loglik 2 A 2 SE Estimate 0.861 0.206 -0.007 0.196 -0.495 0.193 -0.122 0.073 0.052 0.045 -0.002 0.164 -0.011 0.159 0.172 -0.111 -0.511 0.176 -0.094 0.208 -0.661 0.218 0.286 0.275 -1.356 0.264 -1.499 0.262 -2.089 0.167 933.244 (# Iter = 24) 144 Table D.2: Results for Case 2 in Table 7.14 Evaluated at the Boundary: 7702 -> - 0 0 , 7/03 + V2 = A i , 7702 + 7/2 = Parameter Bo Bi (LD) B (HD) 2 63 (time) At (EDSS) "12 «13 «23 «123 Oil Oi2 m(LOR) Ai A Voi Neg. Loglik 2 A 2 Estimate SE 0.894 0.205 -0.016 0.193 -0.483 0.190 -0.117 0.073 -0.004 0.017 -0.004 0.157 -0.011 0.155 -0.111 0:168 -0.511 0.170 0.207 -0.101 -0.649 0.215 0.286 0.251 -1.356 0.209 -1.499 0.237 0.164 -2.089 933.901 (# Iter = 24) 145 7703 - » — 0 0 , Table D.3: Results for Case 3 in Table 7.14 Evaluated at the Boundary: 7703 —> — 0 0 , 7702 -> - 0 0 , 7703 + 772 = A i , 7702 + 772 = A Parameter 2 Estimate SE Po Pi (LD) P2 (HD) Ps (time) Pi (duration) 0.877 -0.025 -0.484 -0.119 0.002 0.205 0.196 0.194 0.073 0.003 "12 0.000 -0.008 -0.108 0.163 0.161 0.173 0.177 "13 "23 "123 "1 "2 Vi(LOR) Ai A 2 Voi Neg. Loglik -0.508 -0.110 -0.651 0.208 0.216 0.286 0.279 -1.356 -1.499 0.263 0.265 -2.089 0.165 933.768 (# Iter = 24) 146 Table D.4: Results for Case 4 in Table 7.14 Evaluated at the Boundary: 7703 —> — 0 0 , 7/02 -> - c o , 7703 + 7/2 = A i , 7702 + % = A Parameter A> 2 Estimate SE 0.745 -0.025 -0.486 -0.118 0.004 0.241 -0.001 -0.005 -0.110 -0.507 -0.108 -0.654 0.162 0.162 0.174 0.286 -1.356 -1.499 0.241 Ai A V01 -2.089 0.138 (LD) (HD) /?3 (time) Bi (age) 01 02 "12 "13 "23 • "123 "1 "2 m(LOR) 2 Neg. Loglik 0.195 0.185 0.073 0.004 0.179 0.208 0.209 0.247 0.232 933.354 (# Iter = 26) 147 Table D.5: Results for Case 5 in Table 7.14 Evaluated at the Boundary: 7702 -> - c o , 7703 + 772 = A i , 7702 + m = A Parameter Po Pi(LD) P (HD) 2 Pz(time) p (log(BOD)) 4 "12 "13 "23 "123 "1 "2 Vi(LOR) Ai A »7oi Neg. Loglik 2 Parameter Po Pi(LD) P (HD) 2 Pz(time) P (log(BOD)) 4 "12 "13 "23 "123 "1 "2 Vi(LOR) Ai A ' V01 Neg. Loglik 2 7703 —> — 0 0 , 2 Imputed Set 1 Estimate SE 0.718 0.253 0.008 0.193 -0.474 0.190 -0.108 0.075 0.020 0.016 -0.022 0.164 -0.031 0.161 -0.131 0.172 -0.530 0.178 -0.071 0.207 -0.624 0.216 0.286 0.274 -1.356 0.264 -1.499 0.263 -2.089 0.164 933.088 (# Iter = 23) Imputed Set 2 Estimate SE 0.726 0.251 0.006 0.195 0.192 -0.475 -0.109 0.073 0.019 0.016 -0.021 0.162 -0.029 0.159 -0.130 0.171 -0.529 0.175 -0.074 0.210 -0.626 0.218 0.286 0.277 -1.356 0.266 -1.499 0.261 -2.089 0.166 933.215 (# Iter = 23) Imputed Set 3 Imputed Set 4 Estimate SE Estimate SE 0.717 0.251 0.725 0.249 0.009 0.190 0.006 0.195 -0.474 0.189 -0.475 0.189 -0.108 0.075 -0.109 0.073 0.020 0.016 0.019 0.016 -0.022 -0.021 0.160 0.160 -0.031 0.158 -0.029 0.155 -0.131 • 0.172 -0.129 0.167 -0.530 0.173 -0.529 0.171 -0.071 0.204 -0.073 0.210 -0^624 0.215 -0.625 0.214 0.286 0.278 0.286 0.283 -1.356 0.266 -1.356 0.268 0.264 -1.499 -1.499 0.267 -2.089 0.166 -2.089 0.166 933.083 (# Iter = 23) 933.211 (# Iter = 23) 148 Table D.6: Results for Case 5 in Table 7.14 Evaluated at the Boundary (364 patients): 7/03 -> " O O , 7/02 -> - O O , 7/ + 7/ = A , 7/ + 7/ = A 03 2 X Imputed with 1.0 Parameter Estimate SE 0.686 0.248 00 0i{LD) 0.063 0.193 B (HD) -0.446 0.196 8% (time) 0.074 -0.108 fa(log(BOD)) 0.020 0.016 -0.061 0.166 "12. 0.162 -0.059 "13 -0.161 0.175 "23 -0.567 0.180 "123 -0.022 0.208 "1 -0.587 0.220 "2 0.274 0.301 ni(LOii) -1.342 0.264 Ai -1.491 0.261 A -2.149 0.168 V01 Neg. Loglik 914.912 (# Iter = 24) 2 2 149 02 2 2 Imputed with 4.5 Estimate SE 0.694 0.249 0.061 0.193 0.192 -0.447 0.074 -0.108 0.020 0.016 -0.060 0.160 -0.058 0.159 0.172 -0.160 -0.565 0.176 -0.025 0.206 0.212 -0.589 0.301 0.276 -1.342 0.264 -1.491 0.263 -2.149 0.171 915.047 (# Iter = 26) Table D.7: Results for Model ID2 Evaluated at the Boundary (364 patients): - 0 0 , 7702 -> - 0 0 , 7703 + 772 = A i , 7702 + 772 = Parameter 00 Pi(LD) 03 (time) "12 "13 "23 "123 Oil "2 m(LOR) Ai A »7oi Neg. Loglik 2 A 2 Estimate SE 0.859 0.206 0.035 . 0.193 -0.457 0.201 0.074 -0.118 -0.041 0.171 -0.037 0.161 0.174 -0.140 -0.546 0.179 -0.056 0.207 -0.614 0.224 0.301 0.276 -1.342 0.250 -1.491 0.266 -2.149 0.170 915.825 (# Iter = 21) 150 7703 ->• Appendix E Detailed Results for the Liu e t a l . Transition Models Described in Section 7.6 151 Table E . l : Results for Liu Transition Model with Drop-out Model ID1 Parameter ft 0i ft ft ft Voz Vl3 V23 V02 Vl2 ' V22 Vol Neg. Loglik SV Estimate (SE) 0.89 1.007 (0.558) -0.12 -0.040 (0.876) -0.50 -0.462 (1.336) -0.42 -0.324 (0.369) 0.90 0.692 (0.707) -1.95 -12.083 (2.727) 0.00 0.558 (1.013) 0.00 10.535 (0.888) -1.95 -18.645 (1.144) 0.00 0.048 (2.246) 0.00 17.318 (0.913) -1.95 -2.089 (0.828) 942.259 (# Iter = 137) SV Estimate (SE) 1.00 1.007 (0.447) -0.04 -0.040 (0.475) -0.50 -0.462 (0.372) -0.30 -0.324 (0.404) 0.70 0.692 (0.361) -2.00 -10.451 (1.197) 0.60 0.558 (0.767) 2.00 8.903 (0.970) -1.00 -22.807 (0.746) 0.05 0.048 (0.646) 2.00 21.480 (0.746) -2.00 -2.089 (0.363) 942.259 (# Iter = 97) Table E.2: Results for Liu Transition Model with Drop-out Model ID2 Parameter ft ft ft ft ft V03 V02 Vol Vi V2 Neg. Loglik SV Estimate (SE) 0.89 1.007 (0.817) -0.12 -0.040 (0.748) -0.50 -0.462 (0.604) -0.42 -0.324 (0.483) 0.90 0.692 (0.631) -1.95 -14.182 (0.991) -1.95 -14.324 (1.388) -1.95 -2.089 (0.924) 0.00 0.286 (0.710) 0.00 12.826 (0.812) 942.687 (# Iter = 58) 152 SV Estimate (SE) 1.00 1.007 (0.873) -0.04 -0.040 (0.710) -0.50 -0.462 (0.625) -0.30 -0.324 (0.481) 0.70 0.692 (0.697) -2.00 -14.175 (0.908) -1.00 -14.318 (0.909) -2.00 -2.089 (0.886). 0.30 0.286 (0.415) 3.00 12.819 (0.594) 942.687-(# Iter = 61) Table E.3: Results for Liu Transition Model with Drop-out Model ID3 Parameter A> 01 SV 0.89 -0.12 03 -0.50 -0.42 04 0.90 Vo -1.95 0.00 0.00 02 Vi V2 Neg. Loglik Estimate (SE) SV Estimate (SE) 1.123 (0.211) 1.10 -0.13 -0.40 -0.128 (0.176) -0.437 (0.173) -0.50 0.60 -0.443 (0.103) 0.573 (0.177) -2.00 0.50 -0.30 -2.023 (0.161) 0.542 (0.313) -0.262 (0.609) -0.128 -0.437 -0.443 0.573 (0.178) (0.174) (0.102) (0.172) -2.023 (0.163) 0.542 (0.311) -0.262 (0.610) 939.471 (# Iter = 17) 1.123 (0.212) 939.471 (# Iter = 15) Table E.4: Results for Liu Transition Model with Drop-out Model ID5 Parameter SV Estimate (SE) SV Estimate (SE) 00 0.89 1.007 (0.931) 1.10 01 -0.12 -0.03 1.007 (0.929) -0.040 (0.994) 02 03 -0.50 -0.42 -0.040 (1.026) -0.462 (0.984) -0.40 -0.30 -0.462 (0.992) -0.324 (0.485) 04 0.90 0.60 0.692 (0.971) -0.324 (0.463) 0.692 (0.975) V03 -1.95 • V02 -14.704 (0.863) -14.832 (0.972) -1.00 -1.95 Voi -1.95 0.00 -2.089 (1.028) 13.539 (0.742) -1.00 -13.967 (0.970) -14.096 (0.894) -2.089 (0.992) 0.00 12.803 (0.570) V2 Neg. Loglik 943.239 (# Iter = 55) 153 -0.00 943.239 (# Iter = 57) Table E.5: Results for Liu Transition Model with Random Drop-out (RD) Parameter Po '0i 02 03 04 V03 ??02 Voi Vo RD1 Estimate 1.113 -0.118 -0.445 -0.431 0.596 -2.416 -2.117 -2.089 SE 0.206 0.170 0.165 0.096 0.160 0.327 0.297 0.108 - Vl3 V12 0.878 0.401 Vi - Neg. Loglik RD2 Estimate 1.113 -0.118 -0.445 -0.431 0.596 -2.239 -2.278 -2.089 - - 0.367 0.337 - 944.101 (# Iter = 21) RD3 Estimate 1.113 -0.118 -0.445 -0.431 0.596 SE 0.209 0.172 0.172 0.099 0.164 0.248 0.261 0.163 - 0.625 0.261 944.518 (# Iter = 17) - -2.068 • SE 0.212 0.175 0.173 0.099 0.169 - 0.132 - 0.432 0.195 939.578 (# Iter = 14) Table E.6: Results for Liu Transition Model with Drop-out Completely At Random (CRD) Parameter 00 0i 02 03 04 V03 V02 Vol Vo Neg. Loglik CRD1 Estimate SE 0.209 1.113 0.174 -0.118 -0.445 0.172 • -0.431 0.099 0.596 0.169 -1.846 0.172 -1.849 0.160 -2.089 0.166 - - 947.589 (# Iter = 13) 154 CRD2 Estimate SE 0.209 1.113 0.172 -0.118 -0.445 0.172 -0.431 0.099 0.164 0.596 - - - - -1.880 0.096 942.074 (# Iter = 11)
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Assessing informative drop-out in models for repeated...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Assessing informative drop-out in models for repeated binary data Er, Lee Shean 2001
pdf
Page Metadata
Item Metadata
Title | Assessing informative drop-out in models for repeated binary data |
Creator |
Er, Lee Shean |
Date Issued | 2001 |
Description | Drop-outs are a common problem in longitudinal studies. In terms of statistical models for the data, there are three types of drop-out mechanisms: drop-out occurring completely at random (CRD), drop-out occurring at random (RD) and informative drop-out (ID). The drop-out mechanism is classified as CRD if the drop-out mechanism is independent of the measurements; as RD if the drop-out mechanism depends only on the observed but not the unobserved measurements, and as ID if the drop-out mechanism depends on both the observed and unobserved measurements. CRD and RD are referred to as ignorable because the drop-out mechanism can be ignored for the purpose of making inferences about the observed measurements, while ID is non-ignorable. Analyses based on an assumption of ignorable drop-out, when in reality the drop-out mechanism is non-ignorable, can lead to misleading or biased results. Likelihood-based models for continuous and categorical longitudinal data subject to non-ignorable drop-out have been developed. In this thesis, we focus on exploring likelihood-based models for binary longitudinal data subject to informative drop-out. The two modelling approaches considered are a selection model proposed by Baker (1995) and a transition model proposed by Liu et al. (1999). We apply these models to a data set from a multiple sclerosis (MS) clinical trial. The aims of the analyses are to investigate whether there is an indication of informative drop-out in this data, and to assess the sentivity of inferences concerning the treatment effects to the underlying drop-out mechanisms. We do not attempt to provide a definitive analyses of the data set, but rather to explore a variety of models which incorporate informative drop-out. |
Extent | 6198152 bytes |
Genre |
Thesis/Dissertation |
Type |
Text |
File Format | application/pdf |
Language | eng |
Date Available | 2009-07-27 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0089811 |
URI | http://hdl.handle.net/2429/11271 |
Degree |
Master of Science - MSc |
Program |
Statistics |
Affiliation |
Science, Faculty of Statistics, Department of |
Degree Grantor | University of British Columbia |
Graduation Date | 2001-05 |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-ubc_2001-0024.pdf [ 5.91MB ]
- Metadata
- JSON: 831-1.0089811.json
- JSON-LD: 831-1.0089811-ld.json
- RDF/XML (Pretty): 831-1.0089811-rdf.xml
- RDF/JSON: 831-1.0089811-rdf.json
- Turtle: 831-1.0089811-turtle.txt
- N-Triples: 831-1.0089811-rdf-ntriples.txt
- Original Record: 831-1.0089811-source.json
- Full Text
- 831-1.0089811-fulltext.txt
- Citation
- 831-1.0089811.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0089811/manifest