Identifying patterns in treatment response profiles in acute bipolar mania: a cluster analysis approach
© Lipkovich et al; licensee BioMed Central Ltd. 2008
Received: 07 June 2007
Accepted: 29 July 2008
Published: 29 July 2008
Patients with acute mania respond differentially to treatment and, in many cases, fail to obtain or sustain symptom remission. The objective of this exploratory analysis was to characterize response in bipolar disorder by identifying groups of patients with similar manic symptom response profiles.
Patients (n = 222) were selected from a randomized, double-blind study of treatment with olanzapine or divalproex in bipolar I disorder, manic or mixed episode, with or without psychotic features. Hierarchical clustering based on Ward's distance was used to identify groups of patients based on Young-Mania Rating Scale (YMRS) total scores at each of 5 assessments over 7 weeks. Logistic regression was used to identify baseline predictors for clusters of interest.
Four distinct clusters of patients were identified: Cluster 1 (n = 64): patients did not maintain a response (YMRS total scores ≤ 12); Cluster 2 (n = 92): patients responded rapidly (within less than a week) and response was maintained; Cluster 3 (n = 36): patients responded rapidly but relapsed soon afterwards (YMRS ≥ 15); Cluster 4 (n = 30): patients responded slowly (≥ 2 weeks) and response was maintained. Predictive models using baseline variables found YMRS Item 10 (Appearance), and psychosis to be significant predictors for Clusters 1 and 4 vs. Clusters 2 and 3, but none of the baseline characteristics allowed discriminating between Clusters 1 vs. 4. Experiencing a mixed episode at baseline predicted membership in Clusters 2 and 3 vs. Clusters 1 and 4. Treatment with divalproex, larger number of previous manic episodes, lack of disruptive-aggressive behavior, and more prominent depressive symptoms at baseline were predictors for Cluster 3 vs. 2.
Distinct treatment response profiles can be predicted by clinical features at baseline. The presence of these features as potential risk factors for relapse in patients who have responded to treatment should be considered prior to discharge.
The clinical trial cited in this report has not been registered because it was conducted and completed prior to the inception of clinical trial registries.
Manic episodes in bipolar disorder requiring hospitalization contribute to substantial personal and economic burden on patients, their families and society , which could be alleviated with rapid resolution of acute manic symptoms . Unfortunately, patients with acute mania can respond differentially to treatment due to underlying symptomatic dimensions that vary individually, which can result in failure to obtain or sustain symptom remission . Behavioral subtypes in mania have been described  and were found to respond differently to treatment  even though the pattern of manic symptom reduction in response to treatment was not found to differ. However, this observation was based on endpoint measures, and any treatment-associated differences in the pattern of symptom change that occurred earlier during treatment may have been missed . What might be of potential utility to the clinician would be to identify distinct time-course patterns of acute treatment response.
The objective of the present study was to identify groups of patients with similar response profiles and to construct predictive models for these groups using baseline data.
The present study is unique in that patients with acute mania were assessed daily during the first week of treatment which allowed the capture of symptom changes that characterized their longer-term response to treatment.
Study design and patient population
This was a post-hoc analysis of patients (n = 222) with bipolar mania or mixed-episode, with or without psychotic features, from a randomized, double-blind clinical trial comparing olanzapine (n = 125; 5–20 mg/d) and divalproex (n = 123; 500–2500 mg/d) in a 3-week acute phase  followed by a 44-week maintenance phase . Details of this study have been reported elsewhere, and the following is a brief summary. Acute mania was assessed as having Young-Mania Rating Scale (YMRS)  total scores ≥ 20 at baseline. Depressive symptoms were assessed with the Hamilton Rating Scale for Depression (HAMD) . Patients were assessed daily during the first week, weekly for the following 2 weeks, and then bi-weekly for the remainder of the study (a total of 7 weeks for this analysis).
During this study period, remission was defined as having YMRS total scores ≤ 12, and sustained remission was defined as having YMRS total scores ≤ 12 for at least 2 subsequent visits including the last one (discontinuation or at Week 7, whichever occurred first). Relapse (ongoing mania or recurring mania post-baseline after treatment initiation) was defined as having YMRS total scores ≥ 15 at any assessment, and "sustained relapse" (sustained manic episode) was defined as having YMRS total scores ≥ 15 for at least 2 subsequent visits including the last one (discontinuation or at Week 7, whichever occurred first). To facilitate clinical interpretation of the resulting clusters, sustained remission and sustained relapse were evaluated as binary outcomes for each cluster.
The objective of the present analysis was to identify distinct response-to-treatment profiles in acute mania that may be representative of certain types of patients, their disease history, baseline conditions, or treatment regimen. An easy-to-implement way of identifying different patterns of treatment response is by clustering estimated regression coefficients fitted to each individual patient's response profile [9–11]. To achieve this, a method of hierarchical clustering based on Ward's distance was applied to the estimated coefficients of orthogonal polynomial regressions of the third degree that were fitted to YMRS total scores during the first 7 weeks of treatment. To smooth YMRS profiles, splines were applied to the data (prior to fitting orthogonal polynomials) using penalized least-squares estimates from a non-parametric regression (SAS® PROC TPSPLINE). This procedure adapts to the pattern of individual profiles by varying model complexity for individual curves. While attempting to reach good fit to the data, this procedure also avoids excessive roughness or rapid variation [12, 13]. The outlined procedure of approximating data using orthogonal basis functions (here polynomials) and clustering estimated coefficients appears quite naturally in the context of functional data analysis [9–11, 14]. Specifically, when using orthogonal polynomials, the coefficients have the following natural interpretation, independently of other terms included: the coefficient at linear term represents an overall trend in the outcome profile, the coefficient at quadratic term is the rate of change in the expected outcome at time t, and so on .
The number of clusters that best represented patient response to treatment as revealed by hierarchical clustering was selected by comparing the percent variation accounted for by the clusters (R2) achieved at each level of data aggregation. The level(s) that corresponded to subsequent sharp deterioration of R2 were selected for further inspection, and the clusters that were the most interpretable in terms of clinical relevance were selected for final presentation.
One difficulty encountered with this dataset was the large number of patient dropouts. In order to regularize the data prior to fitting the curves, profiles of patients who discontinued during Weeks 2–7 of treatment were completed using the average of 50 imputed values generated with Bayesian regression fitted to the previous outcomes at every time point past the first week (SAS® PROC MI). This imputation procedure is consistent with the likelihood of observed data and is not expected to introduce any bias when evaluating representative curve shapes, as long as data are missing at random (that is, the probability of dropouts does not depend on unobserved outcomes after the observed outcomes have been accounted for).
To identify baseline variables that may have been contributing to membership in a particular cluster, logistic regression (SAS® PROC LOGISTIC) with stepwise variable selection was used. Baseline variables included in the analysis were: HAMD total scores, individual YMRS items scores, gender, ethnicity, age, weight, body mass index (BMI), index episode and/or illness features (mixed, psychotic, rapid cycling), number of previous manic, mixed, or depressive episodes in the previous year/or lifetime, number of hospitalizations, age at onset of illness, and treatment during the study (olanzapine vs. divalproex).
Baseline characteristics of patients who had at least one week of active treatment.
(n = 111)
(n = 111)
YMRS total, mean (SD)
HAMD-21 total, mean (SD)
Gender (% male)
Age, year, mean (SD)
Ethnicity (% white)
Weight, kg, mean (SD)
Mixed episode (%)
Psychotic features (%)
Rapid cycling (%)
Number of manic episodes in previous year, mean (SD)
Cluster 1 (Figure 2A) included patients (n = 64) most of whom (95%) did not attain sustained remission of manic symptoms, and who largely (over 50% of them) had a sustained relapse. In this cluster, mania symptoms gradually improved during the first week (mean YMRS total scores = 19.6 [SD 8.8] at Day 7), and 46% of the patients were discharged from the hospital by the end of the first week. In the following weeks, mania symptoms of patients in Cluster 1 improved slightly, but by the end of Week 7 most of these patients still had not remitted indicating poor response to treatment.
Cluster 2 (Figure 2B) included patients (n = 92) who experienced rapid improvement in their symptoms in the first 2 days, and by Day 7, mean YMRS total score was 10.2 (5.8), and 86% were discharged. Most of the patients in this cluster (66%) reached a sustained remission and continued to improve for the duration of the study, and only 5% had a sustained relapse.
Cluster 3 (Figure 2C) included patients (n = 36), who also experience rapid improvement in their symptoms within the first 2 days of treatment, but their symptoms had worsened slightly by Day 7, even so 84% were discharged during the initial week of treatment. Most (82%) of the patients in Cluster 3 had met relapse criterion at Week 2 (n = 34, mean YMRS total score of 21.8 [8.4]), and their YMRS scores did not improve substantially thereafter: 72% of them had maintained relapse until the end of the 7-week period (or discontinuation). There was only one patient (3%) in this cluster who had a sustained remission.
Cluster 4 (Figure 2D) included patients (n = 30) whose mania symptoms improved gradually during the first week and 66% were discharged. Most (83%) of the patients in this cluster were remitted by Week 2, and 80% maintained their remission for the duration of the 7-week period (or until discontinuation). There was only one patient (3%) who had a sustained relapse.
Baseline characteristics of patients associated with each cluster.
n = 64
n = 92
n = 36
n = 30
Age, y mean (SD)
White Ethnicity, %
Weight, kg mean (SD)
YMRS total, mean (SD)
YMRS-09: Disruptive-aggressive behavior, mean (SD)
YMRS-10: Appearance, mean (SD)†
HAMD-21 total, mean (SD)
Mixed episode, %
Psychotic features, %
Rapid cycling, %
# manic episodes previous year, mean (SD)
Randomized to olanzapine, %
Patients in Clusters 2 and 3 also shared certain clinical characteristics at baseline, but had different outcomes (Figure 1): less than one-third of the patients in these clusters had psychotic symptoms, over 50% were in a mixed episode, and over 60% were rapid cycling. Cluster 3 patients were poor responders with an average of 8 manic episodes in the year preceding study entry, and over two-thirds were randomized to divalproex. Cluster 2 patients, who were responders, had on average 4 manic episodes in the preceding year and over 60% were randomized to olanzapine. Membership in Cluster 3 vs. 2 was significantly predicted by randomization to divalproex (adjusted OR = 5.0, CI: 1.88–13.3, p = .001), experiencing a larger number of manic episodes over the previous 12 months (adjusted OR = 1.11, 95% CI: 1.04–1.18, p = .001), lack of disruptive-aggressive behavior (YMRS Item 9, adjusted OR = 0.57, 85% CI: 0.38–0.85, p = .006) and to some extent by more prominent depressive features at baseline (greater HAMD total score; adjusted OR = 1.10, 95% CI: 1.02–1.18, p = .013).
As we found that patients in Cluster 3 differentiated from those in Cluster 2 by having a substantially larger number of previous manic episodes during the preceding year, one may question whether the study investigators considered the baseline history of manic episodes when making a decision to discharge a patient. From visually inspecting the curves along with associated rates of discharge for each cluster, it seems that the investigators tended to discharge patents with robust improvement during the first week, since Clusters 2 and 3 with rapid early response have the highest proportion of early discharges. To formally evaluate this relationship we fitted logistic regression models with early discharge (yes/no as a response) and YMRS scores during the first 7 days and various baseline and historical characteristics evaluated as potential predictors. The stepwise logistic model revealed that the most important predictors of early discharge were baseline YMRS total score (adjusted odds ratio (OR) = 0.89, p < .001), and change in YMRS from baseline to evaluation at Days 2 (adjusted OR = 0.94, p = .026) and 7 (adjusted OR = 0.93, p = .003). Interestingly, the previous number of manic episodes had a marginally significant but positive effect on early discharge (adjusted OR = 1.08, p = .067), which suggests that while early discharge was driven by mania improvement during the first 7 days, the patient's history of recurrence of mania was not considered in making early discharge decisions. Of course, from our data, we cannot evaluate the possible impact on the outcome, had such information been considered and patients with higher number of manic episodes not been discharged.
In this post-hoc analysis of a double-blind clinical trial  in patients with bipolar mania or mixed-episode, we identified 4 clusters with different response profiles that could be grouped into 2 contrasting sets of patterns. The first set involved contrasting patients who rapidly improved during the first week of treatment, some of whom either relapsed (Cluster 3) or continued improving (Cluster 2), and a second set contrasting patients who showed less rapid response during the first week of treatment, one group of which did not gain remission (Cluster 1) and the other group continued to improve to remission (Cluster 4). The presence of psychotic features, not being in a mixed episode and the YMRS Appearance item were the most significant predictors of slower initial improvement as represented by Clusters 1 and 4 vs. Clusters 2 and 3. Patients in Clusters 1 and 4 differed in their rates of rapid cycling at baseline, 39.1% vs. 56.7%. Interestingly Cluster 4, which had the higher rates for rapid cycling, eventually gained remission. For Cluster 3 vs. 2, the larger number of previous manic episodes and randomization to divalproex treatment were significant predictors of relapse following rapid reduction in symptoms. The similarity between clinical courses for Clusters 1 and 4 and for Clusters 2 and 3 during the first week of treatment and their markedly different response profiles is interesting, and may have implications about the importance of clinical observation early in treatment of patients in order to avert potential relapse.
Treatment of manic episodes may be difficult when they are complicated by the presence of other features such as depression, psychosis, and anxiety. Attempts have been made to characterize subtypes of mania that may be associated with these other features and that respond differently to treatment. In a study in untreated patients with acute mania, Swann et al  utilized factor analysis of behavioral rating scales scores followed by cluster analysis, which yielded 4 mania subtypes they described as anxious-depressive, psychotic, classic euphoria, and irritable. These subtypes were subsequently found to respond differently to treatment with divalproex or lithium . The anxious-depressed subtype had higher distressed appearance scores and was resistant to treatment, the psychotic and classic euphoria subtypes responded to either lithium or divalproex, and the irritable subtype responded better to divalproex than lithium.
In the current study, we used data from a clinical trial focusing primarily on identifying patterns of response that were based on individual mania profiles during the active treatment and not specifically on clustering patients by pre-treatment mania features. Therefore, the clinical characteristics of the clusters presented herein were not as distinct or as definitive as those for the naturalistic subtypes described by Swann et al . However, the response profiles for the clusters in the present study are supported by previous reports. Cluster 4, which was comprised of more psychotic patients, with more severe mania baseline scores, and less severe depressive scores responded to treatment with either divalproex or olanzapine similar to the psychotic mania subtype described by Swann et al . Patients in Cluster 3, who had more severe depressive symptoms at baseline, relapsed within 2 weeks and remained resistant to treatment similar to the anxious-depressed subtype . In addition, Cluster 3 was characterized by having more prior manic episodes, which has been reported by Welge et al  to be a predictor of poor response. Cluster 2, which responded well to treatment with either olanzapine or divalproex, may share features with the classic euphoria subtype of Swann et al . Cluster 1 characteristics and response profiles do not appear to share similarities with any of the naturalistic mania subtypes proposed by Swann et al .
This was a post-hoc analysis of a study not designed to assess the complexities of treatment response patterns and warrants further study. Furthermore, replication of this method in larger populations of patients may be needed to validate clusters identified based on this limited data. A further limitation is that this cluster analysis deterministically identified mania profiles ignoring the inherent uncertainty of class membership. For example, some patients may have had profiles on the borderline of two response patterns. Ignoring such uncertainty may result in overstating the odds ratios from the logistic regressions predicting class membership. Additionally, the relationships found between subscores of YMRS and class membership probabilities may be somewhat overstated because the baseline YMRS total score was included as part of the outcome profile creating an overlap. It is important to emphasize that the findings of this study are exploratory and not confirmatory; therefore interpretation of all hypothesis testing should be done with caution and merely indicate possible existing patterns that need to be validated in different datasets, perhaps using more powerful clustering procedures like latent class finite mixture modeling [13, 16–18]. This type of modeling would evaluate the profiles representing distinct clusters, as well as predictors of class membership, while properly accounting for class uncertainty within a single estimation step. It requires certain parametric assumptions about the curves and underlying mixture distributions, and the results may be highly sensitive to the initial partitioning of observations into clusters. The approach of the present paper to cluster response profiles was exploratory and data-driven, and is similar to that of Tarpley et al. [10, 11]. In addition, our purpose was to search for interesting patterns, including those that were unusual and of relatively low occurrence potentially representing an anomaly, which is more in line with the philosophy of exploratory functional data analysis .
In an attempt to validate our exploratory findings, we reanalyzed our data utilizing a fully parametric maximum-likelihood based approach proposed in Jones et al  (that is very similar to that of Muthen and Shedden ) using the freely available SAS PROC TRAJ software  developed by Jones and colleagues. The clustering structure of the data obtained using the latent class mixtures was very similar to our original findings, and the variables identified as significant predictors of class membership were virtually identical to those we identified using simple logistic regression. It is important to point out, however, that the methods used in the present analysis have not been validated in a separate study and may have limitations that have yet to be realized.
Cluster analysis was useful to identify interesting patterns in treatment response profiles that would be difficult to detect using only pre-specified clinical definitions of relapse or response. Patient history and features of manic episode were important for discriminating between these patterns. In particular, the previous number of manic episodes was found to be a risk factor for relapse in patients who initially responded to treatment and therefore should be considered prior to discharge.
Funding provided by Eli Lilly and Company.
- Olie JP, Levy E: Manic episodes: the direct cost of a three-month period following hospitalization. Eur Psychiatry. 2002, 17: 278-86. 10.1016/S0924-9338(02)00680-6.View ArticlePubMedGoogle Scholar
- Tohen M, Waternaux CM, Tsuang MT: Outcome in mania: a 4-year prospective follow-up of 75 patients utilizing survival analysis. Arch Gen Psychiatry. 1990, 47: 1106-11.View ArticlePubMedGoogle Scholar
- Swann AC, Bowden CL, Calabrese JR, Dilsaver SC, Morris DD: Pattern of response to divalproex, lithium, or placebo in four naturalistic subtypes of mania. Neuropsychopharmacology. 2002, 26: 530-6. 10.1016/S0893-133X(01)00390-6.View ArticlePubMedGoogle Scholar
- Swann AC, Janicak PL, Calabrese JR, Bowden CL, Dilsaver SC, Morris DD, Petty F, Davis JM: Structure of mania: depressive, irritable, and psychotic clusters with different retrospectively-assessed course patterns of illness in randomized clinical trial participants. J Affect Disord. 2001, 67: 123-32. 10.1016/S0165-0327(01)00447-5.View ArticlePubMedGoogle Scholar
- Tohen M, Baker RW, Altshuler LL, Zarate CA, Suppes T, Ketter TA, Milton DR, Risser R, Gilmore JA, Breier A, Tollefson GA: Olanzapine versus divalproex in the treatment of acute mania. Am J Psychiatry. 2002, 159: 1011-17. 10.1176/appi.ajp.159.6.1011.View ArticlePubMedGoogle Scholar
- Tohen M, Ketter TA, Zarate CA, Suppes T, Frye M, Altshuler L, Zajecka J, Schuh LM, Risser R, Brown E, Baker RW: Olanzapine versus divalproex sodium for the treatment of acute mania and maintenance of remission: a 47-week study. Am J Psychiatry. 2003, 160: 1263-71. 10.1176/appi.ajp.160.7.1263.View ArticlePubMedGoogle Scholar
- Young RC, Biggs JT, Ziegler VE, Meyer DA: A rating scale for mania: reliability, validity and sensitivity. Br J Psychiatry. 1978, 133: 429-35.View ArticlePubMedGoogle Scholar
- Hamilton M: A rating scale for depression. J Neurol Neurosurg Psychiatry. 1960, 23: 56-62.View ArticlePubMedPubMed CentralGoogle Scholar
- Tarpey T: Linear transformations and the k-means clustering algorithm: Application to clustering curves. Am Stat. 2007, 61: 34-40. 10.1198/000313007X171016.View ArticlePubMedPubMed CentralGoogle Scholar
- Tarpey T, Petkova E, Ogde RT: Profiling placebo responders by self-consistent partitions of functional data. J Am Stat Assoc. 2003, 98: 850-8. 10.1198/016214503000000783.View ArticleGoogle Scholar
- Tarpey T, Kinateder K: Clustering functional data. J Classification. 2003, 20: 93-114. 10.1007/s00357-003-0007-3.View ArticleGoogle Scholar
- Edwards LJ, Stewart PW, MacDougall JE, Helms RW: A method for fitting regression splines with varying polynomial order in the linear mixed model. Stat Med. 2006, 25: 513-27. 10.1002/sim.2232.View ArticlePubMedGoogle Scholar
- Luan Y, Li H: Clustering of time-course gene expression data using a mixed-effects model with B-splines. Bioinformatics. 2003, 19: 474-82. 10.1093/bioinformatics/btg014.View ArticlePubMedGoogle Scholar
- Ramsay JO, Silverman BW: Functional Data Analysis. 1997, New York Springer, 310-View ArticleGoogle Scholar
- Welge JA, Keck PE, Meinhold JM: Predictors of response to treatment of acute bipolar manic episodes with divalproex sodium or placebo in 2 randomized, controlled, parallel-group trials. J Clin Psychopharmacol. 2004, 24: 607-12. 10.1097/01.jcp.0000145342.38170.7f.View ArticlePubMedGoogle Scholar
- McLachlan G, Peel D: Finite Mixture Models. 2000, John Wiley & SonsView ArticleGoogle Scholar
- Muthen B, Shedden K: Finite mixture modeling with mixture outcomes using the EM algorithm. Biometrics. 1999, 55: 463-9. 10.1111/j.0006-341X.1999.00463.x.View ArticlePubMedGoogle Scholar
- Jones B, Nagin D, Roeder K: A SAS procedure based on mixture models for estimating developmental trajectories. Social Methods Res. 2001, 29: 374-93. 10.1177/0049124101029003005.View ArticleGoogle Scholar
- SAS PROC TRAJ. Accessed January 2008., [http://www.andrew.cmu.edu/user/bjones/index.htm]
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-244X/8/65/prepub