 Methodology
 Open Access
 Published:
Realising stratified psychiatry using multidimensional signatures and trajectories
Journal of Translational Medicine volume 15, Article number: 15 (2017)
Abstract
Background
Stratified or personalised medicine targets treatments for groups of individuals with a disorder based on individual heterogeneity and shared factors that influence the likelihood of response. Psychiatry has traditionally defined diagnoses by constellations of cooccurring signs and symptoms that are assigned a categorical label (e.g. schizophrenia). Trial methodology in psychiatry has evaluated interventions targeted at these categorical entities, with diagnoses being equated to disorders. Recent insights into both the nosology and neurobiology of psychiatric disorder reveal that traditional categorical diagnoses cannot be equated with disorders. We argue that current quantitative methodology (1) inherits these categorical assumptions, (2) allows only for the discovery of average treatment response, (3) relies on composite outcome measures and (4) sacrifices valuable predictive information for stratified and personalised treatment in psychiatry.
Methods and findings
To achieve a truly ‘stratified psychiatry’ we propose and then operationalise two necessary steps: first, a formal multidimensional representation of disorder definition and clinical state, and second, the similar redefinition of outcomes as multidimensional constructs that can expose within and betweenpatient differences in response. We use the categorical diagnosis of schizophrenia—conceptualised as a label for heterogeneous disorders—as a means of introducing operational definitions of stratified psychiatry using principles from multivariate analysis. We demonstrate this framework by application to the Clinical Antipsychotic Trials of Intervention Effectiveness dataset, showing heterogeneity in both patient clinical states and their trajectories after treatment that are lost in the traditional categorical approach with composite outcomes. We then systematically review a decade of registered clinical trials for cognitive deficits in schizophrenia highlighting existing assumptions of categorical diagnoses and aggregate outcomes while identifying a small number of trials that could be reanalysed using our proposal.
Conclusion
We describe quantitative methods for the development of a multidimensional model of clinical state, disorders and trajectories which practically realises stratified psychiatry. We highlight the potential for recovering existing trial data, the implications for stratified psychiatry in trial design and clinical treatment and finally, describe different kinds of probabilistic reasoning tools necessary to implement stratification.
Background
There is a growing consensus that psychiatric disorders—defined as syndromes in categorical diagnostic systems such as the International Classification of Diseases [1] and the Diagnostic and Statistical Manual [2]—are heterogeneous in terms of aetiology, presenting psychopathology, and treatment response. Simultaneously, it is now recognized that the search for the aetiology and treatment of psychiatric disorders is not wellserved when traditional categorical definitions of syndromes in the ICD10 and DSM5 are effectively equated with disorders. Instead, the emerging paradigm of stratified psychiatry [3] emphasises multifactorial predictors or moderators of disorders grounded in underlying neurobiology, gene/environment interactions and intermediate endophenotypes such that the final phenotypic expression does not necessarily align with the “classical” disorder specifications of e.g. schizophrenia, bipolar affective disorder and mood disorders.
This has led to the development of the Research Domain operational Criteria (RDoC) [3–8] and the roadmap for mental health research (ROAMER) [9]. These proposals reconstruct psychiatric disorders on the basis of biological mechanism and endophenotypes that describe (1) the aetiology of the disorder, (2) help identify predictors and biomarkers for the disease and/or subtype the disease and (3) variation in response to treatment. This may help resolve the apparent dilemma experienced in everyday clinical practice, where two patients respond differently to the same intervention—one patient’s symptoms and signs improve substantially but another’s remain stubbornly unresponsive. Through the lens of stratified psychiatry, these two patients share some features, but may not necessarily have the same disorder despite a common categorical label of, for example, schizophrenia.
Recent studies examining illness features of mood disorders have shown that combinations of clinical variables predict response to selective serotonin reuptake inhibitors (SSRIs) [10–12]. In substance misuse disorders, cocaine dependence is predicted by a combination of parameters in neurocognitive measures of impulsecontrol [13]. In schizophrenia, antipsychotic treatment response and pathogenesis are predicted by overlapping sets of genes [14].
To constrain our scope, we necessarily focus on the single group of schizophreniform disorders. The proposals that follow are, however, equally applicable to other psychiatric disorders but the grouping of phenotypes will differ. For example, in terms of the DSM5, depressive disorders may be too coarsegrained, but persistent depressive disorder (dysthymia) and major depressive disorder may be appropriate because they share syndromic features. However, depressive episodes with psychotic or catatonic features might be better dealt with separately because of the qualitatively different presentation of these patients.
Given our scope, and to make our examples concrete, we further focus on the neurocognitive symptoms of psychotic disorders that are clinically significant for prognosis and quality of life [15, 16] and have proven notoriously difficult to treat with either cognitive or pharmacological interventions, despite over a decade of clinical trials. Beyond the obvious explanation that some of the compounds or interventions trialled do not affect the relevant neurophysiological substrates of these features of the illness, we propose that the failure stems from methodological problems inherited from broader considerations about disorder definition and the nature of treatment response in individuals recruited into these trials.
To proceed, we first consider how disorder definition can be reconceptualised using clinical signatures and a multidimensional understanding of signs and symptoms of illness which serves recruitment of patients into studies, in other words, patient stratification. Clinical signatures represent an individual’s disease state at a given time, defined along different quantified dimensions (loosely, ‘axes’ in a multidimensional space) of symptoms, signs or quantitative/qualitative measures (such as the individual’s loading or expression for a biomarker) that we define and operationalise from first principles using concepts found in multivariate statistics and pattern recognition. In the context of both clinical trials and interventions more generally, response is then naturally defined as a trajectory between clinical signatures before and after interventions for the same patient. In contrast, in the traditional analysis of trial data, response is defined a priori as an aggregate univariate outcome for a group of patients defined by their shared categorical disorder. We propose that with the richer multidimensional information contained in signatures, treatment response can be framed as modelling trajectories, retaining valuable information that predicts or stratifies how individual patients might respond. We then show how this proposal can be made compatible with traditional analyses by defining outcomes as a function of trajectories, rather than starting with univariate aggregate outcomes.
These principles are then applied in the context of a systematic review of a decade of trials in the treatment of neurocognitive features of schizophreniform disorders, evaluating their suitability or “readiness” for a stratified approach. We conclude with proposals for the future design and analysis of patientlevel data for trials and clinical practice.
Step one: multidimensional definition of disorder
Moving beyond categorical diagnoses
In psychiatry, disorders have historically been defined by necessity as syndromes with operationalised thresholds for symptoms and signs to justify a categorical diagnosis. Modern classification systems such as ICD and DSM follow from the traditions of nineteenth century medicine [17]. Psychiatric disorders appear as categories in the ICD system from version 6, published in 1948 [18] with the first edition of the DSM arriving in 1952. The empirical fieldtrialled approach to definitions reflected in DSMIII (1980) resulted in operational criteria, which are reflected, though less prescriptively, in ICD from version 9 (in use from 1979 to 1994). The development of the ICD follows its origins as a statistical method for indexing causes of mortality and morbidity [18] whereas the DSMIV and later, DSMV, continue to focus on operational criteria for making diagnoses. The ICD in its current revision (version 10) continues to be hierarchical, whereas the DSMV is more dimensional and organised around current understanding of phenotypes. Because of it’s operational, criteriabased approach, the DSM is often used for clinical and research work while outside Europe, the ICD10 is used primarily as a hospital and healthcareprovider coding system, rather than a diagnostic system per se. In either case, clinicians use categorical diagnoses as a technical shorthand, to provide treatments to patients who fit the syndromic description (i.e. given by the DSM criteria or described by the ICD narrative) and who are presumed to have the nominal disorder (e.g. schizophrenia, major depressive disorder). Medications are given regulatory approval if they can be shown to treat the categorical disorder [19, 20], on the basis of clinical trials, in which patients are assigned to treatment based on the same categorical disorders. We will argue that this fails patients given our contemporary understanding of the aetiology of psychiatric disorders, how we design trials and decide on outcomes.
In the midtwentieth century, advances in psychophysics, neuroscience and the emergence of an overarching cognitive science reframed observable behaviour in the context of its proximal causes i.e. in broadly descending levels of abstraction as cognitive, neurophysiological, molecular and genetic systems. Within a classification system, the definition of categorical diagnoses is enforced by mutual exclusivity such that these categorical labels represent a ‘disorder’ [21] and the assignment of more than one of these represents comorbidity. Clinical management is then necessarily ‘mapped’ onto these categories [22] to produce clinical guidelines for treatment. This applies even when multiaxial versions of classification systems are considered.
However, equating disorders to their syndromes is demonstrably artificial. In the past decade, research on schizophrenia has benefitted from insights into genetics of specific features of the disorder; for example, variation in DTNBP1 with severity of general cognitive performance [23], polymorphism of the COMT gene in working memory [24, 25] and executive function that are further implicated in response to antipsychotic medication [24, 26, 27]. Similarly, genomic studies [28, 29] and phenotype clustering [28, 30–32] have shown that the traditional view of diagnostic categories has less utility than symptom and signspecific definition of disorder and illness at the individual level.
Moreover, categorical classification cannot readily accommodate the observation of shared symptoms in the face of divergence in neurodevelopment. Diagnostic categories may exhibit marked differences in neurodevelopment [33–35] but also overlap in terms of shared symptoms: psychotic features are seen in borderline personality disorder [36–39], and bipolar disorder and schizophrenia both exhibit similarities in nonverbal communication [40], affective symptoms [41], cognitive deficits [42, 43], genetic risk factors [44–46] and a broader transdiagnostic ‘psychosis’ phenotype [47]. The observation that signs or symptoms are rarely exclusive is an inevitable stumbling block for categorical classification in psychiatry, because they are a distinct manifestation of neurobiological dysfunction whose idiosyncratic expression in an individual patient is shaped by complex environmental factors in the history of that individual’s illness. Any given sign or symptom is likely shared across categorical disorder boundaries and is rarely observed in only a single disorder. Therefore, psychiatric disorders as entities can only be investigated if signs and symptoms are considered as arising in “bottom up” fashion from dysfunction in neurophysiology, genetics and environment—a principle embodied in both RDoC and ROAMER. The implication for patients is then apparent: the unique cluster of signs and symptoms, their disorder signature, may arise from a number of interacting dysfunctional cognitive systems that further map to a number of underlying neurobiological deficits. Thus, in a group of patients, the meaningful observable variables are the state of their signature components rather than their a priori diagnoses according to a categorical system.
Dimensional definitions of disorder
In recent years, dimensional definitions of disorders have been advocated [41, 48–50] such that disease differentiation, between health and pathology as well as between different disorders with some shared symptoms and signs [50], proceeds by comparing a given patient state (signature) to prototypes along one or more dimensions/scales. Some have noted that this approach is complicated by the absence of a “mathematical, precise resolution of what constitutes ‘sufficiently similar’ patients” [51], which in fact succinctly captures the focus of the current proposal, as it is central to and necessary for the concept of stratified psychiatry [3].
Any proposal for a signature and prototype based approach needs to be capable of modelling differences and similarities between patients. For ease of exposition, but without loss of generality, consider Fig. 1a, where two perpendicular ‘axes’ form a plane (a two dimensional Cartesian space) with each axis representing a monotonically increasing scale from lowtohigh symptom ‘load’ (for example, in schizophreniform illness, the axes could represent psychotic positive and negative symptom load respectively). A patient’s signature at a given time is defined as the point on the plane given by its twodimensional coordinates. Figure 1a shows two patients with very different signatures at time 1, but signatures at time 2 which are close in the plane. Under this model, the cause of these different signatures and trajectories could represent (1) different interventions, (2) different responses to the same intervention or (3) patients with different disorders (but with the same categorical disorder label). We note that although we use symptoms as examples, these axes could as easily represent any quantitative state ‘marker’; for example, a genotypic riskprofile score [29], functional neuroimaging activation of cortical areas [52] or indices of neurocognitive performance [53].
When we populate the space of signatures with multiple patients, clustering of individuals in certain locations may be observed. Figure 1b expands on this showing a two dimensional example of the positive and negative symptoms scores of 100 simulated patients. The patients’ signatures are shown as two tentative clusters (black and grey points). Patients’ similarity to each other (illustrated as dotted lines) is inversely proportional to the Euclidean (‘straight line’) distance between them. This intuition of defining signatures and spaces derives from the formal definition of Hilbert and innerproduct spaces that generalise Cartesian coordinates [54, 55] such that ideas of similarity and distance can be naturally extended to describe arbitrarily highdimensional spaces. Each dimension can be concretely envisaged as an axis orthogonal (perpendicular) to the others and represents any measurement of clinical state. It will later be demonstrated how the 30dimensional ‘space’ of signatures can be formed by translating the Positive and Negative Symptom Scale (PANSS) items into patient signatures.
Operationalising signatures, prototypes and stratification
In Fig. 1b it can be seen that the two groups of patients have a prototype, shown as coloured diamonds that are the centroids (the multivariate statistical means, or first moments) of the two distributions and labelled Class A and B. Note that these prototypes need not represent categorical disorders in the ICD10/DSM5 sense, but, instead, meaningful groups of patients with sufficiently similar signatures. Stratification then becomes the task of assigning membership of a given patient (signature) to one of the two tentative groups (classes), where ‘membership’ can be either a discrete or continuous property.
Stratification has been implemented traditionally by inclusion and exclusion criteria that attempt to control patient heterogeneity by defining the categorical diagnosis of interest, and then specifying confounds or disease features that would alter the interpretation of the trial results. The process of specifying these criteria to arrive at a prototypical patient suitable for a given trial may serve for diseases characterised by few and narrowly defined criteria and low comorbidity, which dictates by necessity low heterogeneity among the individuals that make up that population. This however, is not the case in psychiatry, psychosis being a prime example, which demonstrates high heterogeneity, evident in well documented differences in clinical presentation as well as drug response. Traditional approaches to trial design and statistical data treatment are not equipped to deal with heterogeneity (i.e. nonnormality, outliers); they use inclusion/exclusion criteria to overcome it, in the process sacrificing valuable data from those nonprototypical individuals that fail at the cutoff, which leads to trial outcomes and drugs that work for some but not all. In “Experiment one: patient heterogeneity and stratification” section, we illustrate how datadriven approaches are ideally suited to expose and capitalise on this heterogeneity, by taking into account the entire sample, and rather than binary inclusion/exclusion, treat each individual on the basis of a similarity score to an estimated prototype. This approach maximises the utility of each patient enrolled in a trial, or registered on a database, flexibly adapting to subtly shifting prototypes which may vary between datasets and importantly, represents a powerful tool with which to use heterogeneity rather than seek to limit it.
In the discrete case, ‘hardclassifying’ to stratify patients involves assigning their signature to one of the mutually exclusive discrete classes (A or B) according to a decision rule; for example, in Fig. 1b, we compute the distance of any signature to the prototypes (centroids) of Class A and Class B, then assign the signature to the class with smallest distance. Alternatively, and especially relevant to our proposal, is that membership to a class can be a continuous ‘soft’ property. In this case, assignment is probabilistic, i.e. a signature belongs to class A or B with probabilities 0.7 and 0.3 respectively, where the likelihood of belonging to A or B is some function of the distance. Both hard and soft membership can be made more flexible by formulating the problem in Bayesian decision theory (see [56] Chapter 2 for a thorough treatment) where the probability of class membership for a given signature can be biased towards or away from prototypes by defining more sophisticated products of a likelihood (e.g. the distance measure in this example) and priors over classes (e.g. ‘weighting’ the class membership by prevalence of certain prototypes).
Learning models of signatures
Thus far, we have shown that multidimensional signatures can form a space, how a metric can be defined that exposes similarity or difference between patients and the uses of these definitions for ‘hard’ and ‘soft’ stratification. In any statistical model, relationships between variables are most often acquired, learned or estimated from the data. We consider two situations; first, a collection of signatures, but no a priori information about their relationship and second, a set of signatures associated with other known variables. The first scenario represents a situation where we have signatures and seek a datadriven model by inferring prototypes supported by the data. For example, given a large number of patients whose PANSS scores (a 30item instrument) are used to populate a 30dimensional space, we might seek a set of prototypes that represent clusters of patients that are sufficiently similar such that stratified assignment to intervention occurs by similarity in psychopathology. Of note, this represents exploration and discovery of candidate prototypes and the model is constructed by ‘unsupervised’ learning. In Fig. 1b, stratification of a new, previously unknown patient would be some function of the distance to discovered prototypes (solid green and red lines). This straightforward idea is commonplace in the literature on statistical pattern classification and machine learning—see [57] chapter 13 and 14. The second scenario represents a different situation where we have existing information about the patients, for example, where we have patient signatures before an intervention and an indicator of change following treatment (note that we do not necessarily restrict ourselves to a dichotomised or even univariate definition of response or failure, as we will show later when we consider trajectories). In this case, using existing knowledge of treatment response can help validate an existing stratification protocol or to derive one informed by known outcomes (i.e. we know which patients responded, but not why given their signatures). Thus, the aim is to model the ‘mapping’ from patient signatures before intervention, to treatment response after intervention. In this context, estimating the model proceeds by ‘supervised’ learning and results in a discriminant function such as the solid blue line shown in Fig. 1b where in contrast to the unsupervised case, assignment of signatures to class A and B was known in advance. Previously unseen patients are then prospectively stratified based on which ‘side’ of the discriminant line they fall.
It is worth noting that these concepts of representation, similarity/distance and learning/model estimation have a long tradition in statistics, pattern recognition and machine learning, where they appear in the context of class assignment based on prototypes by learning vector quantisation and selforganising maps [58], binary or multinomial classification using linear and quadratic discriminant analysis [59], generalised linear models such as logistic regression [60], support vector machines [61] and Gaussian process regression and classification [62].
There are important considerations in choosing any datadriven unsupervised clustering or supervised classification algorithm (when a priori assumptions about prototypes exist). For unsupervised clustering methods, the choice of similarity metric between signatures as well as the measure of cluster optimality will dictate the discovered prototypes and need to be justified based on the data and their intended use. Rather than use a supervised approach reflecting a priori information about presumed prototypes, in Experiments 1 and 2, we use unsupervised clustering to allow the data to support candidate prototypes that reveal heterogeneity. The algorithm we chose [63] embodies the aforementioned similarity metric (see “Operationalising signatures, prototypes and stratification” section) and implements an optimality measure based on a further measure that constrains the number of discovered prototypes, such that a candidate prototype differentiates itself from others while binding together similar signatures.
There have been attempts to adopt multidimensional principles to the problem of (1) signature definition (2) treatment outcomes and (3) diagnostic prediction. However, none has applied a multidimensional definition consistently across all levels from disorder definition through to disease state and outcome. For example, while multivariate signatures based on symptomatology were used in studies predicting outcome for firstepisode psychosis [64] and symptom severity and persistence in posttraumatic stress disorder (PTSD) [65], only univariate and dichotomised outcomes were addressed and categorical diagnosis was presumed. When neuroimaging data were used as the signature, machine learning classifiers were used to predict categorical diagnosis [66], transition from an atrisk state to a dichotomised ‘psychosis versus health’ outcome [67] and dichotomised clozapine response [68, 69], or univariate aggregate predicted univariate global assessment of function (GAF) [70]. Only one study [71] used machine learning to predict multiple outcomes in major depressive disorder although again, these were dichotomised and did not model trajectories as multidimensional constructs. We argue that stratified psychiatry should avoid categorical diagnoses and univariate treatment outcomes.
Step two: beyond univariate outcomes
Analyses of clinical trial data typically attempt to find a parsimonious group of independent variables to predict a univariate primary outcome, which is, most often, a continuous scalar variable representing aggregate scores on a clinical instrument or a dichotomised outcome based on a cutoff point applied to some aggregate score. For example, response to an intervention may be defined as at least 50% reduction in the overall/aggregated symptom score. Ideally, there will be predictors that exhibit low collinearity, enabling independent effects on the outcome/dependent variable to be modelled using well understood statistical methods (e.g. logistic regression or survival analysis for dichotomised outcomes, and linear/generalised linear models for continuous outcomes). The appeal is obvious—without some clear, univariate primary endpoint (dependent variable), there would appear to be no tractable way of analysing the data.
However, evidencebased treatment protocols derived from randomised controlled trials (RCT), including metaanalysed evidence and naturalistic effectiveness trials such as CATIE, do not translate reliably to clinical practice [72] and leads to publication bias [73]. This is a direct consequence of categorical disorder definition: patients rarely conform to the strict inclusion/exclusion criteria of RCTs, they often exhibit comorbidity, in itself a consequence of categorical disorder definitions as their symptoms do not conform to the demarcated boundaries of their assumed categorical diagnosis [74] and finally, analysis emphasises average response in the presumed homogenous patient group with the same categorical disorder. In recognition of this gap between the evidence base and clinical practice, researchers have begun exploring individualised/stratified medicine. Recent translational research has focused on biomarkers [75] and signatures [76, 77], where multidimensional clinical metrics can be used to define trajectory of illness [78], response to treatment and relapse [79]. In other words, multivariate outcomes derived from multidimensional signatures and their trajectories.
This problem of average response is shared with pain medicine, where the classification of pain syndromes is a function of diffuse pathophysiological mechanisms (which are also difficult to measure ‘objectively’) and individual response to pain medications varies substantially also as a function of psychological and social factors. A study of 200 fibromyalgia patients described treatment response to pregabalin 450 mg for 14 weeks as a bimodal distribution [80]. Given the categorical diagnosis of fibromyalgia, the notion of average response provides no useful estimate of the effect of the drug because patients tended to either respond well, or not at all.
Trial design stands to benefit from ‘responder analysis’ [81–83] which enrols patients sequentially into a trial decisiontree; failing on one “arm” of the tree for whatever reason, e.g. intolerable sideeffects or failure of analgesia enrols them into an alternative treatment arm. A similar approach was used in the design of the CATIE trial [84, 85] and there is now a literature on adaptive trial design [86–88]. This should come as no surprise given the preceding discussion. Assigning patients to treatments by their categorical diagnosis, rather than by signatures or prototypes leads to loss of valuable information about individuals, and avoidable sample attrition. The definition of aggregate, univariate outcomes to measure average treatment response of the group defined by their assumed categorical diagnosis further obscures clinically meaningful treatment response. To paraphrase Moore et al. [80]—we should “expect [analgesic] failure but pursue success” by accommodating response profiles of individual patients rather than focusing on average (unimodal) response.
Multivariate trajectories preserve response information
To illustrate the impact of (1) preserving multidimensional disease signatures and (2) considering multivariate trajectories and outcomes, consistent with the two outlined themes, we present the following simulated example of 100 patients and a hypothetical intervention which yields an 80% improvement for positive symptoms alone, but where only 50% of the patients respond. In assessing treatment response, we are generally interested in changes in signs and symptoms, so we retain the example of using the PANSS scale. However, changes in other relevant measures of disease states could also be used, for example, change in functional neuroimaging markers. If the usual approach of defining an aggregate, univariate outcome is adopted, such as the sum of positive and negative PANSS scores, we arrive at the distributions in Fig. 2a. Note how at time 1 (prior to the intervention, light grey) there is a clear mode to the distribution. At time 2 (following the intervention) there is a wider variance in the summed score, but the average response does not differ significantly from time 1; in other words, this appears to be a failed trial. However, in Fig. 2b which preserves the twodimensional patient signatures in terms of positive and negative symptom load, there is a clear clustering of signatures, highlighting a group of responders (light grey) from nonresponders (black). Even in a simple two dimensional signature space, a dramatic effect is revealed which is otherwise overlooked using a univariate measure of average response.
This simple demonstration illustrates the challenges faced by standard approaches to analysing treatment response based on linear (and generalised linear) models which by definition model a univariate random variable Y as the mean value of Y given the multiple predictors X [60]. If Y represents the ‘collapsed’ outcome—i.e. summed change in positive and negative scores or some dichotomised version of this, e.g. >30% change—it is unsurprising that the different response profiles which are preserved and clearly visible in the multidimensional signatures (Fig. 2b) are inaccessible and essentially lost from the analysis. The standard way of attempting to recover such effects is to use secondary outcomes, which in this case would be a change in the set of positive symptoms. The traditional treatment of secondary outcomes in the univariate framework is to correct for multiple comparisons to avoid type I errors (false positives). However, this correction may be more or less stringent, depending on the number of a priori hypothesised outcome variables, and ultimately one can never be certain that the outcome measures that survive in a particular analysis will be replicated in different datasets or future trials. The advantage of addressing changes in individual symptoms within a multivariate disease trajectory framework is that rather than assigning different status to primary and secondary outcomes, and hence variable statistical fate following correction, response is treated multidimensionally from the outset, so that individual symptoms are treated equally, as is their relation to their associated multidimensional disease signature.
Operationalising trajectories
Without loss of generality we will restrict examples to two time points which can be taken to be before (time 1) and after the intervention (time 2). A signature represents a patient’s state at a given time, so a trajectory is a sequence of such signatures over time. A geometric interpretation of a single patient’s trajectory is the ‘line’ (a vector) in a multidimensional space connecting at least two signatures (as shown in Fig. 1a). This captures and describes change in a single patient, but provides no information about patterns, regularities or structure across many patients. Since there are potentially an infinite number of trajectories, it is necessary to find structure that enables inference over a tractable number, analogous to formally defining events and probability spaces—for discussion, see [89]. We have already shown how clusters can be learned, their prototypes defined and also, that signatures can be hard or softassigned to these clusters using simple metrics. We can therefore use the finite number of prototypes to assign patients to their nearest cluster at different time points. Trajectories are then modelled as a sequence of ‘movements’ (in multidimensional space) between different clusters at each time point.
Anticipating the experiments later in this paper, assume that at time 1, our collection of signatures supports a number of clusters for which we define prototypes—e.g. as the centroid of the cluster—and assign alphabetic class labels A, B, C and so on. Applying the procedure to the same patient signatures at time 2 (postintervention) yields a different number and structure of clusters which, to distinguish from time 1, we label numerically as classes 1, 2, 3 and so on. Note that these label assignments are purely categorical, and do not imply ordering, weighting or ranking. The probability that a given patient is a member of each of the classes is a function of how similar the patient’s signature is to the corresponding prototypes A, B, C at time 1. Using standard probability theory, a given patient’s trajectory is the likelihood of belonging to one of the classes (1, 2, 3) at time 2, given the probability that they belonged to each of the classes A, B, C, at time 1. Such a concept is captured naturally in the theory of Bayesian networks and graphical models [90] which in this example is simply the conditional probability of classes at time 2 on classes at time 1. This provides a tractable, countablyfinite interpretation of trajectories that are supported by the data.
Trajectories and signature spaces allow flexibility in defining outcomes
Thus far, we have used the term trajectory to refer to movements in spaces of signatures representing response to an intervention and argued that this preserves valuable information (cf. Fig. 2a, b). We have intentionally avoided the term ‘outcome’ because this implies a univariate aggregate measure that is most often dichotomised to represent success or failure of treatment.
However, in our signature and trajectory model, response can acquire the meaning of an outcome once a set of conditions over the dimensions of the signature space is applied. For example, the bottom half of the space in Fig. 1b represents signatures where there is low positive, but variable negative symptom load. Therefore, an outcome in our model is specified as a region in the multidimensional signature space, whose semantics can be described as ‘clinically significant change in positive, but not necessarily negative symptoms’. This can be achieved by using the discriminant line (see Fig. 1b), such that any signature, or trajectory end point, that falls in the region below the discriminant line is deemed a treatment ‘success’ in terms of specific change in positive, but not negative, symptom dimensions. Alternatively, the prototypes for classes A and B can be used to define regions based on ‘soft’ assignment within a radius around each prototype, with class B similarly capturing the meaning of treatment success. There is a clear advantage to defining an outcome as a region demarcated by a discriminant line (or radius around a prototype) comprising essentially multiple thresholds within a multidimensional signature space, instead of applying one arbitrary threshold onto a univariate, aggregate measure, e.g. total PANSS score: by modelling highdimensional trajectories which preserve useful information about change, and define what is clinically meaningful multidimensionally, it becomes possible to account for heterogeneity in treatment response over individuals who respond in different ways, but respond nonetheless. In the context of a clinical trial, the inherently richer, multidimensional specification of what constitutes response, or success, effectively reduces the possibility of a false negative result, which can result from defining success according to a single point on a collapsed univariate continuum as illustrated in Fig. 2.
Summary

1.
Patients present with heterogeneous constellations of symptoms that do not respect disorder categories, even within a classical categorical diagnosis. This heterogeneity carries meaning and clinical utility. Recruitment into trials and interventions would be better defined as targeting dimensions of disorder signs and symptoms rather than categories of disorder. Variation and similarity between patients can be defined as multidimensional signature spaces.

2.
The categorical assumption leads to defining outcomes by collapsing multidimensional clinical state to univariate scalar or dichotomised variables. In doing so, information is lost in favour of simplicity of analysis in trials. By representing clinical states as multidimensional spaces of signatures, data driven techniques can identify prototypes that identify relevant structure in these spaces. Trajectories defined as movement between classes and their prototypes at different times define how patients respond to interventions.

3.
Aggregate or ‘collapsed’ measures of clinical state and outcomes, as well as categorical disorder definition enable measurement of only the mean/average response to an intervention. Therefore, the reasons why an intervention works for some individuals but not others are obscured by solely examining average response for a group. By defining multidimensional signature spaces and utilising the prototypes therein, we preserve information and can flexibly define clinically meaningful response by specifying conditions on multiple dimensions (i.e. by defining a region rather than a single threshold on a collapsed, univariate aggregate measure).
Heterogeneity, stratification and trajectories of patients in the CATIE trial
To illustrate the potential of our multivariate framework we apply it to a real, largescale naturalistic RCT of the efficacy of different antipsychotics. The CATIE trial [85] recruited patients on the basis of a categorical diagnosis of schizophrenia and assigned people to parallel arms for comparing the efficacy of a number of antipsychotic medications. If the originally assigned medication arm failed to treat the patient or side effects were intolerable, they were switched to another. For our experiments, following [91, 92] we extracted patientlevel data with baseline PANSS scores as well as repeated neurocognitive performance measures at baseline and at 2month followup (although the CATIE protocol did not specify PANSS reevaluation at 2 months). This yielded 750 complete datasets. In Experiment 1, we use the baseline PANSS data to illustrate principles of multivariate signatures (and prototypes) yielding structure that could be used for prospective patient stratification. In Experiment 2, we apply our method to the neurocognitive data at times 1 and 2 to illustrate signatures, prototypes and trajectories.
Experiment one: patient heterogeneity and stratification
Our hypothesis is that there is heterogeneity and structure in the 750 CATIE patients which is overshadowed by the assumptions of a categorical diagnosis of schizophrenia. By comparing our multivariate framework to univariate measures of clinical state, we show that datadriven stratification (e.g. a set of prototypes) is possible within this categorically defined group of patients.
Methods
The PANSS naturally forms a 30dimensional space, with variates (items) measuring 7 positive, 7 negative and 16 general symptoms. First, the 750 patients were represented in the most reductive way, as one might to measure clinical state suitable for a univariate analysis, by forming the individual univariate distributions of total (summed) positive, negative and general symptoms. The multivariate approach to defining the signature space was then applied. As it is impossible to visualise a 30dimensional space, the total positive, negative and general scores were used to form bivariate (i.e. twodimensional) representations of all patients, with signatures determined by combinations of positive and negative, positive and general and negative and general symptom scores, analogous to the approach in Fig. 2a, b. This representation preserves more information than the univariate approach, although there is still some loss of information as PANSS signatures form a 30dimensional space.
Then, in the full 30dimensional space of PANSS signatures (which cannot be plotted), we perform datadriven unsupervised clustering using Rodriguez and Laio’s algorithm [63] which finds clusters according to two criteria (using the same Euclidean metric discussed earlier) and has been shown to be robust across a number of low and high dimensional clustering problems. First, a quantity rho representing the local density for each patient’s signature is defined as the number of neighbouring signatures inside a specified radius. Then a quantity, delta, is defined for each signature as the minimum distance to any other signature with a higher local density rho. Prototypes are then identified as the individual signatures where delta is anomalously large, in other words, indicative of good separation from the nearest locallydense clusters. We then assign each of the 750 patients to their nearest prototype in 30dimensional space, illustrating the signatureprototype structure by colour coding each patient by their cluster membership in the bivariate representation.
Results
Figure 3 (top row) shows the univariate (single dimensional) representation for summed symptom domains (positive, negative and general). Each shows an approximately normal distribution but no clear heterogeneity beyond that captured by the mean/mode and variance.
The middle row of Fig. 3 shows bivariate plots of signatures formed by combinations of symptom domains. Note how each combination of variables shows an approximately bivariate normal distribution, but where the centroid and covariance are quite different.
Applying the clustering algorithm [63] identified 8 clusters in the full 30dimensional space of PANSS signatures. Figure 3 (bottom row) shows the bivariate plots from the middle row, but with each signature coloured according to the cluster assigned in the full 30dimensional space. It can be seen that univariate (Fig. 3, top) representations conceal heterogeneity that is visible in higher dimensional representations; as dimensionality of the representation increases, structure and heterogeneity are revealed within the signature space. In contrast to the single univariate distribution, 8 clusters emerged in the 30dimensional space, representing an opportunity for stratification of patients by psychopathology.
Experiment two: trajectories in neurocognitive performance
Similarly to the PANSS signatures, our hypothesis is that there is heterogeneity and structure in patients’ neurocognitive performance at baseline and follow up that can be exploited to define trajectories.
Methods
Neurocognitive performance in the CATIE trial was defined via Zscores that summarise performance on cognitive tasks over five domains; verbal working memory, vigilance, speed of processing, reasoning, and memory [91, 92]. This forms a 5dimensional signature space.
We proceeded as for Experiment 1; first, the univariate distributions of each domain were computed separately, then the 5dimensional space was formed, which can be visualised using bivariate plots of combinations of cognitive domain performance (verbal × vigilance, verbal × speed, vigilance × speed, vigilance × reasoning and so on). Prototypes were derived by applying the clustering algorithm in the 5dimensional space, assigning each patient to the nearest cluster by distance to prototype, and finally, colour coding patients in the bivariate plots according to the assigned cluster. To illustrate trajectories, we apply the clustering algorithm separately at baseline (time 1) and twomonth follow up (time 2). To show how trajectories can be formalised probabilistically, we then compute the probabilities of belonging to each cluster at time 2, given the probability of this patient being in each cluster at time 1; a simple discrete Bayesian network was modelled using the gRain package [93]. We then illustrate how a ‘test’ patient’s trajectory can be found by first computing the patient’s distance to prototypes at time 1. Then, we compute the conditional probabilities of cluster membership at time 2 given their original cluster.
Results
As the number of bivariate combinations of each of the 5 cognitive domains is large (requiring 10 separate plots), in Fig. 4, we display verbal working memory as an example known to be impaired in patients with schizophrenia and firstdegree relatives [94]. Following the approach in Experiment 1, the univariate plot of verbal working memory and corresponding bivariate plots by vigilance, processing speed, reasoning and memory are shown. The top row shows time 1 (baseline) and the bottom row time 2 (followup at 2 months of treatment). If a univariate measure alone is defined for verbal working memory (column A) and scores at time 1 and 2 (top, bottom row) are directly compared, there is little discernible change.
The top row of Fig. 4b–e reveals 3 clusters within the patient group at time 1 supported by the data in the 5dimensional space. Note that these are not completely contained in any part of the two dimensional representations, illustrating how information is lost as dimensionality is reduced. At time 2 (Fig. 4, bottom row, panels b–e) application of the algorithm separately on this space revealed 5 clusters, suggesting a higher degree of heterogeneity after treatment. In other words, patient trajectories can be described as ‘starting’ in one of 3 clusters and diverging to ‘arrive’ into one of a separate set of 5 clusters after treatment.
To quantify trajectories over the entire cohort of 750 patients, we computed for each patient the conditional probability of belonging to each of the five clusters at time 2, given the probability that the patient was in one of the three clusters at time 1. In Fig. 5, we label the three clusters at time 1 as a, b and c to emphasise that they are different to the clusters obtained at time 2 which are labelled numerically, 1–5. Figure 5a shows the overarching structure of trajectories for all 750 patients between clusters at time 1 and time 2. The line weight is proportional to the probabilities of ‘arriving’ in clusters 1 through 5 at two month followup, given the probability of starting at baseline in clusters A, B or C (Fig. 5b).
The same model can be used to predict a patient’s trajectory which is crucial to the idea of stratification solely on the basis of their signature at time 1. Figure 5c shows the location of a ‘case’ patient (black circle), in relation to cluster A, B and C prototypes at time 1. Figure 5d displays the likelihood at time 1 of this test patient belonging to clusters A, B or C as a probability distribution, representing ‘soft’ assignment and stratification. By querying the Bayesian network model with these likelihoods [95, 96], it is possible to extract the conditional probabilities of arriving in any of the 5 clusters at time 2. For this patient, Fig. 5f contains these corresponding predictions: clusters 2 and then 4 are the most likely ‘end points’ of the trajectory for this patient. The actual end point of this patient’s trajectory is close to the prototypes for the mostlikely predicted clusters 2 and 4 (Fig. 5e).
Systematic review of trials for neurocognitive symptoms of schizophrenia
To demonstrate the scope of our proposal and how it might influence trial design, we conducted a systematic review of the Clinical Trials registry to locate trials on cognition in schizophrenia over the period 1st January 2004 through 1st September 2015. The aim of this review was not to assess the quality of trials as reported, for example whether significant effects were obtained or whether these were adequately powered. Instead, we sought to systematically interrogate the extent to which they are characterised by assumptions of categorical and univariate assumptions relating to diagnosis and outcome, which as we have argued are incongruent with a stratified psychiatry. Using the set of criteria outlined below, which were designed in line with our proposed stratification framework, we sought to assess broadly the extent to which this recent body of work actually mirrors, conceptually and methodologically, the growing consensus for the need to stratify.
Method and data extraction
Registered trials recruit for one or more categorical diagnoses, therefore, these definitions were necessarily used to mine the literature. Using the Clinical Trials registry (https://clinicaltrials.gov/), our search terms were as follows: trials registered between 01/01/2004 and 01/09/2015 limited to conditions defined as “schizophrenia” AND (“cognition” OR “cognitive”). Interventions were kept broad, and included trials with titles including “drugs” OR “behavioural” OR “dietary” OR “device” OR “transcranial” OR “fMRI” OR “cognitive” in order to capture the largest range of interventions. Primary or secondary outcome measures were included.
Our initial search yielded 114 trials that met the screening criteria. We then excluded 89 records where results were not accessible or were unreported. Of the remaining 25 registered trials, 2 further studies were excluded because there was no reporting of primary or secondary outcomes that included measures specific to cognition (see Preferred Reporting Items for Systematic Reviews and MetaAnalyses (PRISMA) diagram in Fig. 6) suggesting less relevance.
The design and reported analyses of the remaining 23 studies were examined. For each trial, the full text and results tables were examined by two authors (DWJ and JP). For 8/23 trials, results had also been published in journals and these articles were reviewed and tabulated alongside information from the trials database.
Each trial was assessed according to the following criteria:

1.
Diagnostic categories (e.g. schizophrenia, schizoaffective) and the criteria used, which, for all studies, was either DSMIV, DSMIVTR or not specified.

2.
Within the diagnostic inclusion criteria, which specific symptoms/features (for example, as used in the research domain operational criteria), domains (e.g. positive or negative symptoms) or smaller groups of symptoms were specified as inclusion/exclusion criteria. This represented recognition of the heterogeneity of the categorically labelled disorder.

3.
Which patientspecific variables were included in the design—for example, whether symptom signatures/profiles measured by a specific instrument scale/subscale were considered along with endophenotypes or biomarkers such as genomic data, functional magnetic resonance imaging (fMRI), electroencephalography (EEG) features, hypothesised as mediators of response.

4.
Primary and secondary outcomes and their respective measurement type: if composite/aggregate measures where used such as total or mean scores, or domain specific scores on specific subscales/instruments. For example, if a cognitive battery such as Measurement and Treatment Research to Improve Cognition in Schizophrenia (MATRICS) or Brief Assessment of Cognition in Schizophrenia (BACS) was used as primary outcome, were individual subtests measured and reported, or only the composite “overall” score. The measure of change was also examined in terms of hours, days, weeks, months or years and whether statistically significant results were reported.

5.
Whether the primary and secondary outcomes where tested relative to patientspecific variables—for example, change on immediate recall tests with respect to a given biomarker or accounting for symptom signatures or responder status.
Trials accumulated points depending on the quality of specification of the aforementioned criteria, so that a score was assigned for the following:

1.
If the study included specific symptoms or groups of these (Criterion 2) or patient specific variables (Criterion 3), it scored 1 point. If the study included both Criteria 2 and 3, it scored 2 points.

2.
If the study considered primary or secondary outcomes and change measures, where it would be possible to explore individual response, trajectories and differences (e.g. if domain specific subscales rather than only composite/total scores—see Criterion 4) the study scored a further 2 points. If only total or composite scores were reported, only a single point was scored.

3.
A further 1 point was scored for Criterion 5, where outcomes were analysed accounting for specified subjectspecific variables. For example, if an attempt was made to address individual patientlevel factors in mediating response or change in an outcome.
The range of scores was therefore 0–5. If there was insufficient reported information (either in the trials database, or in accompanying published articles), the above rubric was conservatively applied favouring a lower score. If the design of the study recorded an intention to explore features described in the criteria (e.g. to explore genetic predictors, subgroups in psychotic disorders), but it was not mentioned further (e.g. in reported results) then the rubric was applied generously, favouring the higher score. This allowed for studies that were completed and closed, but where all results had not been analysed and published on clinicaltrials.gov.
As an example, if a study recruited subjects with schizophrenia or schizoaffective disorder with specific positive features and cognitive symptoms (such as thought disorder and impaired verbal working memory), with a hypothesised biomarker such as regional fMRI activity change in then by Criteria 2 and 3, 2 points were scored. If the same study records the a specific medication as the intervention, with primary outcome being change in the biomarker and individual symptom domains, then by Criterion 4, a further 2 points are accumulated. If secondary outcomes were changes in specific symptom profiles (e.g. thought disorder) mediated or grouped by the study biomarker, a further point would be scored by Criterion 5, resulting in a total score of 5. This would describe a study where the proposals of this paper are accounted for or there is potential to reanalyse the data using the same principles. Lower scores indicated that the studies were more vulnerable to the problems described in the themes articulated in this paper.
Results
Table 1 shows the complete set of 23 included trials. Of these, 19 were efficacy or augmentation trials, most with a novel application of a drug: 13 of these were repurposing trials, where a drug with an existing license for a different indication was trialled for cognition in schizophrenia, 5 tested new compounds and 1 was an extension study from adult to a paediatric indication. In total, 7 of the included trials have been published (citations are given in Table 1).
Scoring
By applying the scoring criteria to the 18 studies 1 study NCT00611806 [97] scored 5, the highest score, indicating highest compatibility with the criteria; 5 studies scored 4 indicating high compliance and the potential to use the methods proposed in this paper; 2 studies scored 3 indicating that the trial would be difficult to translate in our proposed framework. Scores of 2 and 1, indicating poor compatibility were achieved by 9 and 1 studies respectively, that is, half of the total reviewed (see Table 1).
We consider that those trials scoring 4 (NCT00548327, NCT00963924, NCT00435370, NCT00506077, NCT00455702) could be used in the framework we describe by reanalysis of the data.
Those scoring 3 or less would require substantial revision of design to use the framework proposed here, either because of design or data recording issues such as not including specific symptoms, or the use of only composite scores to measure patient signatures.
Of note, of the 6 trials scoring 4 or 5, five were published studies, suggesting that peer review processes are sensitive to and aligned with the principles discussed in this paper.
Factors affecting response
Of the 18 studies, 10 studies considered specific disorder features; of these, 6 defined the features using “threshold” definitions—e.g. cognitive impairment or PANSS exceeding a value—and the remainder used qualitative descriptors such as “acute exacerbation” and “depressive subtype”. Ten studies defined individualspecific descriptors which might predict or mediate response; three of these studies used biomarkers (EEG, fMRI or genomics)—the others defined clinical state measures (e.g. negative symptoms, clinically stable on named antipsychotic, illness duration and cognitive performance). Four studies both defined patientspecific descriptors and analysed the primary outcome with respect to these. Only one study considered all of these factors NCT00611806 [97].
Response/outcome measurement
For primary outcomes defined on clinical state (e.g. by PANSS scales), 4/18 studies considered domain scores—for example where individual components or subscales were used instead of mean, total or composite scores. Perhaps unsurprisingly, secondary outcomes reflected more attempts to understand multidimensional (rather than univariate) measurement of clinical state change: 10 studies included “domain scores” or some subscalebased test of the hypotheses, but only three actually analysed how the stated patientspecific factors would impact on the domains.
Discussion and conclusions
In this paper, we have described how relying on categorical disorder definition leads to untenable assumptions about homogeneity in patient populations, impacts on the assessment of treatment and finally, coerces measurement of outcome to group averages sacrificing valuable response heterogeneity. We have shown how to operationalise proposals from the literature on dimensional definition of disorders consistent with proposals such as RDoC. To this end, we have introduced concepts of signatures, prototypes and trajectories, and operationalised these from firstprinciples using concepts from multivariate statistics, pattern recognition and probabilistic reasoning. Using well established data from the CATIE trial as an example, we have shown it is possible to apply these concepts for stratification, to define multidimensional trajectories and derive outcomes.
By systematically reviewing just over a decade of clinical trials registry data, we show that a majority of studies are contaminated or confounded by assumptions that are misaligned with principles of stratified psychiatry. This could explain the limited success in finding interventions that work for patients in the domain of neurocognitive symptoms in psychosis, which have been identified to be significant predictors of quality of life and functional outcome. Only a third of the studies reviewed scored 4 or 5 on our criteria, suggesting compatibility with the principles outlined in our proposal. We now consider specific implications of our proposal, highlighting limitations and future research directions.
Clinical symptomatology and neurocognitive data have been used throughout as examples in aid of developing our arguments, primarily to emphasise commonality between clinical diagnoses, outcomes and stratification. Our approach generalises to other kinds of data such as neuroimaging and genomics (“Dimensional definitions of disorder” and “Learning models of signatures” sections), as both symptomatology and biomarker data are amenable to such treatments and similarly subject to constraints arising from data complexity arguments that exceed the scope of this paper. While these are not trivial problems, they are now tractable with contemporary computational resources and recent advances in algorithm design. Genomics research [98, 99] has seen the application of pattern classification and regression techniques more generally [100, 101], with debate focusing on feature selection [102, 103], sparse sampling [57], the ‘curse’ of dimensionality [104] and asymptotic classification performance as the dimensionality of data increases [105].
Implications for patient recruitment and stratified trials
We now consider the practical application of our proposals for trial design. The first and most obvious implication is to recruit patients for constellations of symptoms/signs in alignment with candidate biological mechanisms from the RDoC. As our systematic review demonstrates, almost all trials recruit for a diagnosis in expectation of a response that follows a unimodal average for the diagnostic group but few explore or even define outcome in terms of features of the illness.
Traditionally, if the aim is to study and treat disorganisation in psychosis [106, 107], a disorder (schizophrenia) and a cutoff measure over some constellation of symptoms (e.g. threshold over an aggregate of PANSS or SANS scores) might be defined, in the face of the known difficulties in robustly defining presumed syndrome or subtype of schizophrenia [108] or neurocognitive correlates [109]. Alternatively, in our framework, we would accept patients who display signs and symptoms clinically consistent with disorganisation syndrome regardless of a diagnosis of schizophrenia or schizoaffective disorder for instance. Similarly, in reexamining existing trial data or using Nof1 trial databases, we would ignore categorical diagnosis. Next, each patient’s signature is defined by their neurocognitive performance and signs/symptoms profile (e.g. PANSS). Prototypes can be defined (1) according to clinical judgement (on groupings of patients’ signatures that clinicians agree are the most representative) or (2) using datadriven approaches (e.g. Experiment One). In both cases, visualisation and exploration of multidimensional signatures is required. Then, a ‘radius’ around these prototypes defines how any individual is ‘typical’ of, or similar to, the prototypes, enabling continuous inclusion based on distance from the prototype or if one requires ‘hard’ inclusion/exclusion, discriminant rules can be defined on the basis of prototypes. More patients can be included but each patient is assigned a ‘weight’ of class membership (with classes being defined by prototypes or discriminant surfaces); increasing study power, but at the potential risk of treatment exposure. Alternatively as in [3], one could datamine existing populations of patients with schizophreniform illness for appropriate candidate prototypes, mirroring Experiment One.
The same approach applies to parallelarm or switchover designs [81, 82], including naturalistic designs such as CATIE. Assume we wish to compare the effectiveness of two treatments, and hypothesise that there may be benefit for patients with a disorganised syndrome and less so for patients with pronounced positive symptoms. We define a treatment arm for patients with a disorganised syndrome whose prototype is distinct from those with more pronounced positive symptoms. Patients can then be assigned based on similarity to prototypes for the treatment arms, or randomly; but their multidimensional signature is retained.
Second, we turn to trajectories. In contrast to traditionally defined univariate aggregate outcomes, the multidimensional approach necessitates specification of changes along each axis of the signature space to define outcomes a priori, e.g. for pretrial registration. Signature spaces can be parcellated by discriminant surfaces or by ‘soft’ membership to clusters that divide the space into regions. These regions can be learned from preexisting data, or specified a priori. Treatment response is embodied in a trajectory where the pretreatment signature moves to a location after treatment that lies in some prespecified region—i.e. one ‘side’ of a discriminant surface or within a radius of a prototype. Defining improvement over the constellation of disorganisation signs and symptoms will result in a region that is distinct from improvement in positive symptoms (which is lost if aggregate univariate outcomes are used). If assignment to treatment arms is random, treatment effectiveness can be crudely estimated as the number of patients who, after treatment, transitioned to each region given their similarity to their pretreatment prototypes after treatment.
In terms of stratifying patients prospectively, recent proposals [110] consider using individual patients as Nof1 trials, thereby rendering routine clinical practice a source of data to maximise data availability and a basis for making inferences. Nof1 trials require formalisation of data collection in routine clinical practice and our framework is compatible with collecting and measuring response in a format useful for prospective stratification. A candidate patient case can be assigned to interventions based on prototypes acquired from historical information. As above, similarity to the prototype of interest (a cluster discovered in previous cases) is propagated forward into the analysis of trajectory and response. If two patients respond similarly, that is, their signatures after intervention are similar, but a third patient does not as their signature lies at a greater distance from the other two, then the respective signature distances preceding the intervention may represent a meaningful explanatory covariate in parsing why the third patient did not benefit; or indeed, predicting that another patient, at a similar distance to the third patient, will not respond to the intervention.
This model is not only useful for data with inherent heterogeneity. Even when the space of patient signatures suggests only a single group, recruitment and stratification is no worse off than when recruiting for categorical disorder. However, signatures and trajectories naturally translate into regionbased outcomes which as discussed, render trials less prone to false negatives due to the flexibility afforded by their multidimensional definition. Thus, even within a relatively homogenous group of patients, trajectories offer a rich source of information on those who benefitted and the extent to which they did so.
Implications for the nosology of psychiatric disorders
We have argued that recruiting individuals for a categorical diagnosis leads to incorrect assumptions about measuring how they respond to treatment. Primarily, we have proposed methods for stratifying patients, rather than redefining categories of disorder, just as RDoC emphasises disease signatures underwritten by biological processes. Specifically, [3] describes how a large collection of people with different a priori categorical disorders are aggregated together because of shared disease features; for example, major depressive disorder, dysthymia (mild depressive disorder) and depressivephase bipolar disorder all share affective symptoms and signs. Then, large collections of endophenotype and biomarker data for this population are used to derive datadriven clusters for disorder classification (i.e. exposing common and distinct biological processes underpinning signs and symptoms).
Implications for analysis of trials
We have argued that simple aggregation and dichotomisation of a multidimensional measure of clinical state (e.g. after a trial intervention) to a univariate scalar outcome discards and obscures relationships that may exist in the original multidimensional representation of clinical state as would any data dimensionalityreducing method [111]. Yet this represents a costly tradeoff for ease of analysis and compatibility with established and wellunderstood statistical techniques; one is forced to derive a univariate outcome in order to utilise multivariable regression techniques which map multiple independent/predictor variables to a univariate dependent variable/outcome. The unfortunate cost of this convenience is loss of information in terms of statistical relationships present in the higherfidelity multidimensional representation. Our framework is aligned with the RDoC and proposes refocusing on disorder features, signs and symptoms. We have made a case for ‘quantising’ these highdimensional spaces by learning prototypes that represent datadriven locations in the space of patient signatures that provide meaningful structure. This converts highdimensional signature spaces into a tractable number of locations such that trajectories can be defined as events to be modelled using probability theory and statistics. While we discussed learning or estimating models of prototypes and trajectories using a statistical interpretation, our proposal can also be viewed within nonprobabilistic frameworks such as highdimensional function approximation or within the broader framework of manifold theory [112, 113].
Future work is needed to expand on how uncertainty in outcomes can be represented in our framework. In frequentist and Bayesian statistics, effect sizes on outcomes have confidence and credible intervals respectively. While there is no theoretical obstacle, evidence derived from trials usually hinges on dichotomised outcomes for ease of implementation, e.g. as treatment protocols. It is likely that consulting tables of statistically significant regression coefficients for a single outcome will not suffice, and benefit instead from the implementation of evidential reasoning algorithms such as Bayesian networks in Experiment 2.
Inference for discovery versus prediction
Whether used for predicting outcome or diagnosis, any statistical method including machine learning, multivariate statistics or the more familiar generalised linear models, is vulnerable to difficulties replicating results or model selection and validation on independent samples and this has hampered current research. For example, in the context of gene hunting, searching for neuroimaging markers of disease, or investigating the success or failure of an intervention, we seek an answer to one or more hypotheses while simultaneously, bolstering support for strong association or causality by controlling for confounds within the experimental design. The definition of validity in explanatory modelling [114] is therefore the discrepancy between the model’s prediction output with respect to the data it was fitted to. In this approach, hypotheses are tested (i.e. that the estimated regression coefficients on predictor variables are nonzero) with no explicit requirement that they generalise beyond the study data. Instead, replication and validation of findings proceeds by repeating the experiment (or trial), fitting similar models and retrospective comparison with the previous studies.
We contrast this to inference for prediction [114] and stratification where we are less concerned with discovery (explanation), and focus instead on how available information can be used to prospectively optimise treatment for the individual. In this case, the ‘gold standard’ for model validity is not its fit to collected data, but rather, how well the model predicts in the case of novel independent data sets. Our proposal thus comes with a health warning: inferences drawn from analyses such as those presented above using the CATIE trial data should not be used to retrospectively assert causality as inference for discovery. Instead, we should couch stratification and prediction in terms of probabilistic assignments of ‘weight’ given available information: e.g. “this patient has a 0.75 probability of a positive outcome on positive symptoms given treatment 1, and 0.60 with treatment 2 given they were close to prototype A, but not B”—see [115] for examples of similar reasoning with clinical examples. We propose that techniques suited to this task are likely to resemble efforts in the probabilistic reasoning community, such as graphical models that enable statistical structure in multivariate data to be discovered, modelled and then used for prediction [90, 96].
Methodological challenges: model validation and replication
We have argued that multivariate signatures represent the fundamental ‘unit’ for describing patients and that assessment or measure of outcome in the form of trajectories should respect this. To evaluate the efficacy of such a stratified psychiatry, any predictive model should be validated by its’ ability to predict on independent samples (cf. explanatory models) and this has been lacking in over three decades of published studies (see [116] for a review). Validation on truly independent samples is challenging because patientlevel data is either not suitably collected, or made available; of the existing studies similar to our proposed framework (reviewed in “Learning models of signatures” section), only two [64, 71] make use of independent samples and all rely on crossvalidation for model validation and selection to mitigate against overoptimistic results due to overfitting e.g. the biasvariance tradeoff [57] and inductive bias [117]. Model validation is beyond the scope of this paper, but is addressed in the literature on model selection more generally [118] and criteria for accepting any model (acquired by any statistical method) can be approached from a frequentist, information theoretic or Bayesian perspective [119]. The choice of criteria varies depending on the intended use of the algorithm and the definition of ‘associational’ versus ‘causal’ [120, 121].
As we have demonstrated in the systematic review, data collected using existing clinical trial methods rarely align with the proposals we describe, or those of others [3, 110]. To partially mitigate against this problem, the field is currently constrained to within sample crossvalidation [122, 123], whereby one part of a dataset is used to build a model, while the remainder is used for testing its predictive performance. Moreover, given small sample sizes of expensive biomarker data (e.g. neuroimaging and ‘omics’ data), crossvalidation is performed with recourse to its logical extreme (leaveoneout, LOOCV), but replication failure is unfortunately predictable given that these subsamples are not truly independent. Robust validation requires testing against novel large datasets with adequate resolution in patient level data and compatible measures, conserved across multiple time points; in other words, organised with stratification and prediction in mind from the outset. There are, to our knowledge, no such repositories for psychiatric research. A particularly useful proposal [124] is to build repositories of not only data, but trained models which can be applied to other data sets. This requires standardisation and interfaces that allow exchange of data and implementations of algorithms in a common format. In this way, stratified psychiatry, and our proposal which speaks to it, can progress from the exploratory phase to validation and prospective testing.
Conclusion
We have presented a concrete proposal in response to the growing calls and clear need for the realisation of a truly stratified psychiatry. Our approach meets this need by integrating principles and methods from the mature fields of multivariate statistics and probabilistic reasoning with an evolving nosology in psychiatry. We have drawn on schizophrenia as a particularly challenging and pertinent field of research and clinical practice where treatment resistant disease continues to plague the lives of patients and defy scientists to the extent that it is often considered a separate clinical entity. The thorny problem of treatment resistance may be rendered more tractable if it is addressed in the context of movement in multidimensional signature spaces. By preserving and exploiting heterogeneity and embracing the principles of stratified psychiatry, our work can begin to productively focus on defining trajectories for subgroups of individuals with quantifiably similar signatures.
References
 1.
World Health Organisation. The ICD10 classification of mental and behavioural disorders: clinical descriptions and diagnostic guidelines. 10th ed. Geneva: World Health Organisation; 1992.
 2.
American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders, 5th Edition (DSM5). Diagnostic Stat. Man. Ment. Disord. 4th Ed. TR. 2013.
 3.
Insel TR, Cuthbert BN, Whiteford HA, Collins FS, Varmus H, Insel T, et al. Brain disorders? Precisely. Science. 2015;348:499–500.
 4.
Cuthbert BN, Kozak MJ. Constructing constructs for psychopathology: the NIMH research domain criteria. J Abnorm Psychol. 2013;122:928–37.
 5.
Insel T, Cuthbert B, Garvey M, Heinssen R, Pine DS, Quinn K, et al. Research domain criteria (RDoC): toward a new classification framework for research on mental disorders. Am J Psychiatry. 2010;167:748–51.
 6.
Morris SE, Cuthbert BN. Research domain criteria: cognitive systems, neural circuits, and dimensions of behavior. Dialogues Clin Neurosci. 2012;14:29–37.
 7.
Sanislow CA, Pine DS, Quinn KJ, Kozak MJ, Garvey MA, Heinssen RK, et al. Developing constructs for psychopathology research: research domain criteria. J Abnorm Psychol. 2010;119:631–9.
 8.
Simmons JM, Quinn KJ. The NIMH research domain criteria (RDoC) Project: implications for genetics research. Mamm Genome. 2014;25:23–31.
 9.
Schumann G, Binder EB, Holte A, de Kloet ER, Oedegaard KJ, Robbins TW, et al. Stratified medicine for mental disorders. Eur Neuropsychopharmacol. 2014;24:5–50.
 10.
Chekroud AM, Zotti RJ, Shehzad Z, Gueorguieva R, Johnson MK, Trivedi MH, et al. Crosstrial prediction of treatment outcome in depression: a machine learning approach. Lancet Psychiatry. 2016;3:243–50.
 11.
Iniesta R, Malki K, Maier W, Rietschel M, Mors O, Hauser J, et al. Combining clinical variables to optimize prediction of antidepressant treatment outcomes. J Psychiatr Res. 2016;78:94–102.
 12.
Wallace ML, Frank E, Kraemer HC. A novel approach for developing and interpreting treatment moderator profiles in randomized clinical trials. JAMA Psychiatry. 2013;70:1241–7.
 13.
Ahn WY, Ramesh D, Moeller FG, Vassileva J. Utility of machinelearning approaches to identify behavioral markers for substance use disorders: impulsivity dimensions as predictors of current cocaine dependence. Front Psychiatry. 2016;7:1–11.
 14.
Ruderfer DM, Charney AW, Readhead B, Kidd BA, Kähler AK, Kenny PJ, et al. Polygenic overlap between schizophrenia risk and antipsychotic response: a genomic medicine approach. Lancet Psychiatry. 2016;3:350–7.
 15.
Meier MH, Caspi A, Reichenberg A, Keefe RSE, Fisher HL, Harrington H, et al. Neuropsychological decline in schizophrenia from the premorbid to the postonset period: evidence from a populationrepresentative longitudinal study. Am J Psychiatry. 2014;171:91–101.
 16.
Green MF, Kern RS, Heaton RK. Longitudinal studies of cognition and functional outcome in schizophrenia: implications for MATRICS. Schizophr Res. 2004;72:41–51.
 17.
Gøtzsche PC. Rational diagnosis and treatment: evidencebased clinical decision making. 4th ed. New York: Wiley; 2007.
 18.
Moriyama IM, Loy RM, RobbSmith AHTR. History of the statistical classification of diseases and causes of death. Hyattsville: National Center for Health Statistics; 2011.
 19.
Kawa S, Giordano J. A brief historicity of the Diagnostic and Statistical Manual of Mental Disorders: issues and implications for the future of psychiatric canon and practice. Philos Ethics Humanit Med. 2012;7:1.
 20.
Kinch MS, Patridge E. An analysis of FDAapproved drugs for psychiatric disorders. Drug Discov Today. 2015;20:292–5.
 21.
Berrios GE. Classifications in psychiatry: a conceptual history. Aust N Z J Psychiatry. 1999;33:145–60.
 22.
Amaral OB. Defining disease in the information age. PLoS Med. 2006;3:e317.
 23.
Burdick KE, Goldberg TE, Funke B, Bates JA, Lencz T, Kucherlapati R, et al. DTNBP1 genotype influences cognitive decline in schizophrenia. Schizophr Res. 2007;89:169–72.
 24.
Weickert TW, Goldberg TE, Mishara A, Apud JA, Kolachana BS, Egan MF, et al. CatecholOmethyltransferase val 108/158met genotype predicts working memory response to antipsychotic medications. Biol Psychiatry. 2004;56:677–82.
 25.
Tan HY, Callicott JH, Weinberger DR. Dysfunctional and compensatory prefrontal cortical systems, genes and the pathogenesis of schizophrenia. Cereb Cortex. 2007;17:i171–81.
 26.
Egan MF, Goldberg TE, Kolachana BS, Callicott JH, Mazzanti CM, Straub RE, et al. Effect of COMT Val108/158 Met genotype on frontal lobe function and risk for schizophrenia. Proc Natl Acad Sci USA. 2001;98:6917–22.
 27.
RebolloMesa I, Picchioni M, Shaikh M, Bramon E, Murray R, Toulopoulou T. COMT (Val(158/108)Met) genotype moderates the impact of antipsychotic medication on verbal IQ in twins with schizophrenia. Psychiatr Genet. 2011;21:98–105.
 28.
Potkin SG, Turner JA, Guffanti G, Lakatos A, Fallon JH, Nguyen DD, et al. A genomewide association study of schizophrenia using brain activation as a quantitative phenotype. Schizophr Bull. 2009;35:96–108.
 29.
Schizophrenia Working Group of the Psychiatric Genomics Consortium SWG of the PG, Ripke S, Neale BM, Corvin A, Walters JTR, Farh KH, Holmans PA, et al. Biological insights from 108 schizophreniaassociated genetic loci. Nature. 2014;511:421–7.
 30.
Treutlein J, Mühleisen TW, Frank J, Mattheisen M, Herms S, Ludwig KU, et al. Dissection of phenotype reveals possible association between schizophrenia and Glutamate Receptor Delta 1 (GRID1) gene promoter. Schizophr Res. 2009;111:123–30.
 31.
Wessman J, Paunio T, TuulioHenriksson A, Koivisto M, Partonen T, Suvisaari J, et al. Mixture model clustering of phenotype features reveals evidence for association of DTNBP1 to a specific subtype of schizophrenia. Biol Psychiatry. 2009;66:990–6.
 32.
Jablensky A. Subtyping schizophrenia: implications for genetic research. Mol Psychiatry. 2006;11(9):815–36. doi:10.1038/sj.mp.4001857.
 33.
Murray RM, Sham P, Van Os J, Zanelli J, Cannon M, McDonald C. A developmental model for similarities and dissimilarities between schizophrenia and bipolar disorder. Schizophr Res. 2004;71:405–16.
 34.
Walker J, Curtis V, Murray RM. Schizophrenia and bipolar disorder: similarities in pathogenic mechanisms but differences in neurodevelopment. Int Clin Psychopharmacol. 2002;17(Suppl 3):S11–9.
 35.
Demjaha A, MacCabe JH, Murray RM. How genes and environmental factors determine the different neurodevelopmental trajectories of schizophrenia and bipolar disorder. Schizophr Bull. 2012;38:209–14.
 36.
Barnow S, Arens EA, Sieswerda S, DinuBiringer R, Spitzer C, Lang S. Borderline personality disorder and psychosis: a review. Curr Psychiatry Rep. 2010;12:186–95.
 37.
Schroeder K, Fisher HL, Schäfer I. Psychotic symptoms in patients with borderline personality disorder: prevalence and clinical management. Curr Opin Psychiatry. 2013;26:113–9.
 38.
Glaser JP, Van Os J, Thewissen V, MyinGermeys I. Psychotic reactivity in borderline personality disorder. Acta Psychiatr Scand. 2010;121:125–34.
 39.
NishizonoMaher A, Ikuta N, Ogiso Y, Moriya N, Miyake Y, Minakawa K. Psychotic symptoms in depression and borderline personality disorder. J Affect Disord. 1993;28:279–85.
 40.
Annen S, Roser P, Brüne M. Nonverbal behavior during clinical interviews: similarities and dissimilarities among schizophrenia, mania, and depression. J Nerv Ment Dis. 2012;200:26–32.
 41.
Keshavan MS, Morris DW, Sweeney JA, Pearlson G, Thaker G, Seidman LJ, et al. A dimensional approach to the psychosis spectrum between bipolar disorder and schizophrenia: the SchizoBipolar Scale. Schizophr Res. 2011;133:250–4.
 42.
Jabben N, Arts B, Krabbendam L, Van Os J. Investigating the association between neurocognition and psychosis in bipolar disorder: further evidence for the overlap with schizophrenia. Bipolar Disord. 2009;11:166–77.
 43.
Green MF. Cognitive impairment and functional outcome in schizophrenia and bipolar disorder. J Clin Psychiatry. 2006;67:3–8.
 44.
Lichtenstein P, Yip BH, Björk C, Pawitan Y, Cannon TD, Sullivan PF, et al. Common genetic determinants of schizophrenia and bipolar disorder in Swedish families: a populationbased study. Lancet. 2009;373:234–9.
 45.
Craddock N, O’Donovan MC, Owen MJ. Psychosis genetics: modeling the relationship between schizophrenia, bipolar disorder and mixed (or “schizoaffective”) psychoses. Schizophr Bull. 2009;35:482–90.
 46.
Purcell SM, Wray NR, Stone JL, Visscher PM, O’Donovan MC, Sullivan PF, et al. Common polygenic variation contributes to risk of schizophrenia and bipolar disorder. Nature. 2009;460:748–52.
 47.
Reininghaus U, Böhnke JR, Hosang G, Farmer A, Burns T, McGuffin P, et al. Evaluation of the validity and utility of a transdiagnostic psychosis dimension encompassing schizophrenia and bipolar disorder. Br J Psychiatry. 2016;209:107–13.
 48.
Allardyce J, Suppes T, Van Os J. Dimensions and the psychosis phenotype. Int J Methods Psychiatr Res. 2007;16:S34–40.
 49.
Demjaha A, Morgan K, Morgan C, Landau S, Dean K, Reichenberg A, et al. Combining dimensional and categorical representation of psychosis: the way forward for DSMV and ICD11? Psychol Med. 2009;39:1943–55.
 50.
Goldberg D. The overlap between the common mental disorders—challenges for classification. Int Rev Psychiatry. 2012;24:549–55.
 51.
Djulbegovic B, Paul A, Uk F. From efficacy to effectiveness in the face of uncertainty. JAMA. 2011;305:2005–6.
 52.
Gong Q, Li L, Du M, PetterssonYeo W, Crossley N, Yang X, et al. Quantitative prediction of individual psychopathology in trauma survivors using restingstate FMRI. Neuropsychopharmacology. 2014;39:681–7.
 53.
Ziermans T, de Wit S, Schothorst P, Sprong M, van Engeland H, Kahn R, et al. Neurocognitive and clinical predictors of longterm outcome in adolescents at ultrahigh risk for psychosis: a 6year followup. PLoS ONE. 2014;9:e93994.
 54.
Young N. An introduction to Hilbert space. Cambridge: Cambridge University Press; 1988.
 55.
Semmes S. An introduction to analysis on metric spaces. N Am Math Soc. 2004;50:438–43.
 56.
Duda RO, Hart PE, Stork DG. Pattern classification. New York: Wiley; 2000.
 57.
Hastie T, Tibshirani R, Friedman J. The elements of statistical learning. 2nd ed. New York: Springer; 2009.
 58.
Kohonen T. Selforganizing maps. Berlin: Springer; 1995.
 59.
Fisher RA. The use of multiple measurements in taxonomic problems. Ann Eugen. 1936;7:179–88.
 60.
McCullagh P, Nelder JA. Generalized linear models. 2nd ed. Boca Raton: Chapman & Hall/CRC; 1989.
 61.
Cortes C, Vapnik V. Support vector networks. Mach Learn. 1995;20:273–97.
 62.
Rasmussen CE, Williams CKI. Gaussian processes for machine learning. Cambridge: The MIT Press; 2005.
 63.
Rodriguez A, Laio A. Clustering by fast search and find of density peaks. Science. 2014;344:1492–6.
 64.
Koutsouleris N, Kahn RS, Chekroud AM, Leucht S, Falkai P, Wobrock T, et al. Multisite prediction of 4week and 52week treatment outcomes in patients with firstepisode psychosis: a machine learning approach. Lancet Psychiatry. 2016;3:935–46.
 65.
GalatzerLevy IR, Karstoft KI, Statnikov A, Shalev AY. Quantitative forecasting of PTSD from early trauma responses: a machine learning application. J Psychiatr Res. 2014;59:68–76.
 66.
Orrù G, PetterssonYeo W, Marquand AF, Sartori G, Mechelli A. Using support vector machine to identify imaging biomarkers of neurological and psychiatric disease: a critical review. Neurosci Biobehav Rev. 2012;36:1140–52.
 67.
Koutsouleris N, Meisenzahl EM, Davatzikos C, Bottlender R, Frodl T, Scheuerecker J, et al. Use of neuroanatomical pattern classification to identify subjects in atrisk mental states of psychosis and predict disease transition. Arch Gen Psychiatry. 2009;66:700.
 68.
KhodayariRostamabad A, Hasey GM, MacCrimmon DJ, Reilly JP, de Bruin H. A pilot study to determine whether machine learning methodologies using pretreatment electroencephalography can predict the symptomatic response to clozapine therapy. Clin Neurophysiol. 2010;121:1998–2006.
 69.
KhodayariRostamabad A, Reilly JP, Hasey GM, de Bruin H, Maccrimmon DJ. A machine learning approach using EEG data to predict response to SSRI treatment for major depressive disorder. Clin Neurophysiol. 2013;124:1975–85.
 70.
de Wit S, Ziermans TB, Nieuwenhuis M, Schothorst PF, van Engeland H, Kahn RS, et al. Individual prediction of longterm outcome in adolescents at ultrahigh risk for psychosis: applying machine learning techniques to brain imaging data. Hum Brain Mapp 2016;38(2):704–14. doi:10.1002/hbm.23410.
 71.
Kessler RC, van Loo HM, Wardenaar KJ, Bossarte RM, Brenner LA, Cai T, et al. Testing a machinelearning algorithm to predict the persistence and severity of major depressive disorder from baseline selfreports. Mol Psychiatry. 2016;21:1366–71.
 72.
Berkowitz RL, Patel U, Ni Q, Parks JJ, Docherty JP. The impact of the clinical antipsychotic trials of intervention effectiveness (CATIE) on prescribing practices: an analysis of data from a large midwestern state. J Clin Psychiatry. 2012;73:498–503.
 73.
Turner EH, Knoepflmacher D, Shapley L, Dwan K, Altman D, Arnaiz J, et al. Publication bias in antipsychotic trials: an analysis of efficacy comparing the published literature to the US food and Drug Administration database. PLoS Med. 2012;9:e1001189.
 74.
Goldberg D. Should our major classifications of mental disorders be revised. Br J Psychiatry. 2010;196:255–6.
 75.
Schwarz E, Guest PC, Steiner J, Bogerts B, Bahn S. Identification of bloodbased molecular signatures for prediction of response and relapse in schizophrenia patients. Trans Psychiatry. 2012;2:e82.
 76.
Gaebel W, Riesbeck M. Are there clinically useful predictors and early warning signs for pending relapse? Schizophr Res. 2014;152:469–77.
 77.
Herz MI, Lamberti JS. Prodromal symptoms and relapse prevention in schizophrenia. Schizophr Bull. 1995;21:541–51.
 78.
Remington G, Foussias G, Agid O, Fervaha G, Takeuchi H, Hahn M. The neurobiology of relapse in schizophrenia. Schizophr Res. 2014;152:381–90.
 79.
Emsley R, Chiliza B, Asmal L, Harvey BH. The nature of relapse in schizophrenia. BMC Psychiatry. 2013;13:50.
 80.
Moore A, Derry S, Eccleston C, Kalso E. Expect analgesic failure; pursue analgesic success. Br Med J. 2013;346:2690.
 81.
Moore RA, Derry S, McQuay HJ, Straube S, Aldington D, Wiffen P, et al. Clinical effectiveness: an approach to clinical trial design more relevant to clinical practice, acknowledging the importance of individual differences. Pain. 2010;149:173–6.
 82.
McQuay HJ, Derry S, Moore RA, Poulain P, Legout V. Enriched enrolment with randomised withdrawal (EERW): time for a new look at clinical trial design in chronic pain. Pain. 2008;135:217–20.
 83.
Toth C, Mawani S, Brady S, Chan C, Liu C, Mehina E, et al. An enrichedenrolment, randomized withdrawal, flexibledose, doubleblind, placebocontrolled, parallel assignment efficacy study of nabilone as adjuvant in the treatment of diabetic peripheral neuropathic pain. Pain. 2012;153:2073–82.
 84.
Stroup TS, Lieberman JA, McEvoy JP, Davis SM, Swartz MS, Keefe RSE, et al. Results of phase 3 of the CATIE schizophrenia trial. Schizophr Res. 2009;107:1–12.
 85.
Stroup TS, Mcevoy JP, Swartz MS, Byerly MJ, Qlick ID, Canive JM, et al. The National Institute of Mental Health Clinical Antipsychotic Trials of Intervention Effectiveness (CATIE) project: schizophrenia trial design and protocol development. Schizophr Bull. 2003;29:15–31.
 86.
Chow SC. Adaptive clinical trial design. Annu Rev Med. 2014;65:405–15.
 87.
Chow SC, Chang M. Adaptive design methods in clinical trials—a review. Orphanet J Rare Dis. 2008;3:11.
 88.
Chang M, Chow SC, Pong A. Adaptive design in clinical research: issues, opportunities, and recommendations. J Biopharm Stat. 2006;16:299–309.
 89.
Papoulis A, Pillai SU. Probability, random variables, and stochastic processes. 4th ed. New York: McGrawHill; 2002.
 90.
Koller D, Friedman N. Probabilistic graphical models: principles and techniques. Cambridge: MIT Press; 2009.
 91.
Keefe RSE, Bilder RM, Harvey PD, Davis SM, Palmer BW, Gold JM, et al. Baseline neurocognitive deficits in the CATIE schizophrenia trial. Neuropsychopharmacology. 2006;31:2033–46.
 92.
Keefe RSE, Bilder RM, Davis SM, Harvey PD, Palmer BW, Gold JM, et al. Neurocognitive effects of antipsychotic medications in patients with chronic schizophrenia in the CATIE Trial. Arch Gen Psychiatry. 2007;64:633–47.
 93.
Hojsgaard S. Graphical independence networks with the gRain package for R. J Stat Softw. 2012;46:1–26.
 94.
Conklin HM, Curtis CE, Katsanis J, Iacono WG. Verbal working memory impairment in schizophrenia patients and their firstdegree relatives: evidence from the digit span task. Am J Psychiatry. 2000;157:275–7.
 95.
Lauritzen SL, Spiegelhalter DJ. Local computations with probabilities on graphical structures and their application to expert systems. J R Stat Soc Ser B. 1988;50:157–224.
 96.
Pearl J. Probabilistic reasoning in intelligent systems. Morgan Kauffman: San Mateo; 1988.
 97.
Roffman JL, Lamberti JS, Achtyes E, Macklin EA, Galendez GC, Raeke LH, et al. Randomized multicenter investigation of folate plus vitamin B12 supplementation in schizophrenia. JAMA Psychiatry. 2013;70:481–9.
 98.
Ayers KL, Cordell HJ. SNP selection in genomewide and candidate gene studies via penalized logistic regression. Genet Epidemiol. 2010;34:879–91.
 99.
Zou H, Hastie T. Regression shrinkage and selection via the elastic net, with applications to microarrays. J R Stat Soc Ser B. 2003;67:301–20.
 100.
Tibshirani R. Regression shrinkage and selection via the lasso: a retrospective. J R Stat Soc Ser B Stat Methodol. 2011;73:273–82.
 101.
Zou H, Hastie T. Regularization and variable selection via the elastic net. J R Stat Soc Ser B Stat Methodol. 2005;67:301–20.
 102.
Saeys Y, Inza I, Larrañaga P. A review of feature selection techniques in bioinformatics. Bioinformatics. 2007;23:2507–17.
 103.
Guyon I, Elisseeff A. An introduction to variable and feature selection. J Mach Learn Res. 2003;3:1157–82.
 104.
Bellman RE. Adaptive control processes: a guided tour. Princeton: Princeton University Press; 1961.
 105.
Hughes GF. On the mean accuracy of statistical pattern recognizers. IEEE Trans Inf Theory. 1968;14:55–63.
 106.
Wallwork RS, Fortgang R, Hashimoto R, Weinberger DR, Dickinson D. Searching for a consensus fivefactor model of the Positive and Negative Syndrome Scale for schizophrenia. Schizophr Res. 2012;137:246–50.
 107.
Daban C, Amado I, Baylé F, Gut A, Willard D, Bourdel MC, et al. Disorganization syndrome is correlated to working memory deficits in unmedicated schizophrenic patients with recent onset schizophrenia. Schizophr Res. 2003;61:323–4.
 108.
Liddle PF. The symptoms of chronic schizophrenia: a reexamination of the positivenegative dichotomy. Br J Psychiatry. 1987;151:145–51.
 109.
Nieuwenstein MR, Aleman A, de Haan EHF. Relationship between symptom dimensions and neurocognitive functioning in schizophrenia: a metaanalysis of WCST and CPT studies. J Psychiatr Res. 2001;35:119–25.
 110.
Schork NJ. Time for oneperson trials. Nature. 2015;520:609–11.
 111.
Borg I, Groenen PJF. Modern multidimensional scaling: theory and applications. 2nd ed. New York: Springer; 2005.
 112.
Tu LW. An introduction to manifolds. 2nd ed. New York: Springer; 2010.
 113.
De Silva V, Tenenbaum JB. Unsupervised learning of curved manifolds. In: Denison DD, Hansen MH, Holmes CC, Mallick B, Yu B, editors. Nonlinear estimation and classification. New York: Springer; 2003. p. 453–65.
 114.
Shmueli G. To explain or to predict? Stat Sci. 2010;25:289–310.
 115.
Goodman SN. Aligning statistical and scientific reasoning. Science. 2016;352:1180–1.
 116.
Marquand AF, Wolfers T, Mennes M, Buitelaar J, Beckmann CF. Beyond lumping and splitting: a review of computational approaches for stratifying psychiatric disorders. Biol Psychiatry Cogn Neurosci Neuroimaging. 2016;1:433–47.
 117.
Gordon DF, Des Jardins M. Evaluation and selection of biases in machine learning. Mach Learn J. 1995;20:1–17.
 118.
Burnham KP, Anderson DR. Model selection and multimodel inference: a practical informationtheoretic approach. 2nd ed. New York: Springer; 2002.
 119.
Vehtari A, Ojanen J. A survey of Bayesian predictive methods for model assessment, selection and comparison. Stat Surv. 2012;6:142–228.
 120.
Holland PW. Statistics and causal inference. J Am Stat Assoc. 1986;81:945.
 121.
Pearl J. Causal inference in statistics: an overview. Stat Surv. 2009;3:96–146.
 122.
Arlot S, Celisse A. A survey of crossvalidation procedures for model selection. Stat Surv. 2010;4:40–79.
 123.
Kohavi R. A Study of Crossvalidation and bootstrap for accuracy estimation and model selection. In: Proceedings of 14th international joint conference artificial intelligence. 1995. p. 1137–1143.
 124.
Young J, Kempton MJ, McGuire P. Using machine learning to predict outcomes in psychosis. Lancet Psychiatry. 2016;3:908–9.
 125.
R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3900051070. http://www.Rproject.org (2008).
 126.
Wickham H. ggplot2: elegant graphics for data analysis. New York: Springer; 2009.
 127.
Wilke CO. cowplot: streamlined plot theme and plot annotations for ‘ggplot2’. R package version 0.6.1. https://CRAN.Rproject.org/package=cowplot. (2016).
 128.
Wickham H. Reshaping data with the reshape package. J Stat Softw. 2007;21(12):1–20.
 129.
Højsgaard S. Graphical independence networks with the gRain package for R. J Stat Softw. 2012;46(10):1–26.
 130.
Venables WN, Ripley BD. Modern applied statistics with S. 4th ed. New York:Springer; 2002. ISBN 0387954570.
 131.
Csardi G, Nepusz T: The igraph software package for complex network research, InterJournal, Complex Systems 1695. 2006. http://igraph.org.
 132.
Kuhn M (coauthors: Wing J, Weston S, Williams A, Keefer C, Engelhardt A, Cooper T, Mayer Z, Kenkel B, The R Core Team, Benesty M, Lescarbeau R, Ziem A, Scrucca L, Tang Y, Candan C). caret: classification and regression training. R package version 6.064. https://CRAN.Rproject.org/package=caret. (2016).
 133.
Cain CK, McCue M, Bello I, Creedon T, Tang DI, Laska E, Goff DC. DCycloserine augmentation of cognitive remediation in schizophrenia. Schizophr Res. 2014;153:177–83. doi:10.1016/j.schres.2014.01.016.
 134.
Davis MC, Lee J, Horan WP, Clarke AD, McGee MR, Green MF, Marder SR. Effects of single dose intranasal oxytocin on social cognition in schizophrenia. Schizophr Res. 2013;147:393–7. doi:10.1016/j.schres.2013.04.023.
 135.
Kane JM, D'Souza DC, Patkar AA, Youakim JM, Tiller JM, Yang R, Keefe RS. Armodafinil as adjunctive therapy in adults with cognitive deficits associated with schizophrenia: a 4week, doubleblind, placebocontrolled study. J Clin Psychiatry. 2010;71:1475–81. doi:10.4088/JCP.09m05950gry.
 136.
Zhang XY, Liu L, Liu S, Hong X, Chen DC, Xiu MH, Kosten TR. Shortterm tropisetron treatment and cognitive and p50 auditory gating deficits in schizophrenia. Am J Psychiatry. 2012;169:974–81. doi:10.1176/appi.ajp.2012.11081289.
 137.
Egan MF, Zhao X, Gottwald R, HarperMozley L, Zhang Y, Snavely D, Lines C, Michelson D. Randomized crossover study of the histamine H3 inverse agonist MK0249 for the treatment of cognitive impairment in patients with schizophrenia. Schizophr Res. 2013;146:224–30. doi:10.1016/j.schres.2013.02.030.
 138.
Goff DC, Cather C, Gottlieb JD, Evins AE, Walsh J, Raeke L, Otto MW, Schoenfeld D, Green MF. Onceweekly dcycloserine effects on negative symptoms and cognition in schizophrenia: an exploratory study. Schizophr Res, 2008;106:320–7. doi:10.1016/j.schres.2008.08.012.
 139.
Buchanan RW, Keefe RS, Lieberman JA, Barch DM, Csernansky JG, et al. A randomized clinical trial of MK0777 for the treatment of cognitive impairments in people with schizophrenia. Biol Psychiatry. 2011;69(5):442–9. doi:10.1016/j.biopsych.2010.09.052.
 140.
Goff DC, Herz L, Posever T, Shih V, Tsai G, Henderson DC, Freudenreich O, Evins OE, Yovel I, Zhang H, Schoenfeld D. A sixmonth, placebocontrolled trial of dcycloserine coadministered with conventional antipsychotics in schizophrenia patients. Psychopharmacology. 2005;179(1):144–50. doi:10.1007/s0021300420322.
Authors’ contributions
DWJ and JP curated the trial data and performed the systematic review. DWJ conducted the analysis and wrote the first draft. DWJ and AA designed the data visualisation and with DK and SSS, wrote the final draft. SSS provided overall supervision of the project. All authors read and approved the final manuscript.
Acknowledgements
The authors gratefully acknowledge data used in the preparation of this article that resides in the NIHsupported NIMH Data Repositories. All analysis was conducted using the R statistics environment [125] using the following packages: ggplot2 [126], cowplot [127], reshape2 [128 http://www.jstatsoft.org/v21/i12/], gRain [129 http://www.jstatsoft.org/v46/i10/], MASS [130], igraph [131], and caret [132].
Competing interests
The authors declare they have no competing financial interests. Three of the authors (DWJ, DT, SSS) are practicing psychiatrists.
Availability of data and materials
The datasets generated during and/or analysed during the current study are available in the NIMH Data Archive repository, Study ID: N01 MH09000106, https://ndar.nih.gov/edit_collection.html?id=2081.
Consent for publication
No individually identifiable data is included.
Ethical approval
The study reanalyses data contained in the NIMH National Clinical Trials Database and was analysed according to an approved Data Use Certificate.
Funding
DWJ is funded by a National Institute of Health Research Integrated Academic Training Clinical Lectureship. S.S.S. is supported by a European Research Council Consolidator Award (Grant Number 311686) and the National Institute for Health Research (NIHR) Mental Health Biomedical Research Centre at South London and Maudsley NHS Foundation Trust and King’s College London. The funders had no role in study design, data collection, data analysis, data interpretation, or writing of the report.
Author information
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Stratified psychiatry
 Trials
 Methodology
 Multivariate