Recovery from depression often demonstrates a nonlinear pattern of treatment response, where the largest reduction in symptoms is observed early followed by smaller improvements. This study investigated whether this exponential pattern could model the antidepressant response to repetitive transcranial magnetic stimulation (TMS). Symptom ratings from 97 patients treated with TMS for depression were collected at baseline and after every five sessions. A nonlinear mixed-effects model was constructed using an exponential decay function. This model was also applied to group-level data from several published clinical trials of TMS for treatment-resistant depression. These nonlinear models were compared to corresponding linear models. In our clinical sample, response to TMS was well modeled with the exponential decay function, yielding significant estimates for all parameters and demonstrating superior fit compared to a linear model. Similarly, when applied to multiple studies comparing TMS modalities as well as to previously identified treatment response trajectories, the exponential decay models yielded consistently better fits compared to linear models. These results demonstrate that the antidepressant response to TMS follows a nonlinear pattern of improvement that is well modeled with an exponential decay function. This modeling offers a simple and useful framework to inform clinical decisions and future studies.
Article
Modeling the Antidepressant Treatment Response to Transcranial Magnetic Stimulation using an Exponential Decay Function
https://doi.org/10.21203/rs.3.rs-1874941/v1
This work is licensed under a CC BY 4.0 License
Journal Publication
published 02 May, 2023
You are reading this latest preprint version
Major depressive disorder is the leading cause of disability worldwide [1] and represents a common risk factor for poor clinical outcomes. Unfortunately, many individuals with major depression do not respond to conventional pharmacological interventions [2]. Repetitive transcranial magnetic stimulation (TMS) offers a noninvasive nonpharmacological treatment option that is effective for treatment-resistant depression [3]. This therapeutic approach uses repetitive pulsed magnetic fields to induce depolarization in brain regions that modulate neural networks and alter pathological patterns of brain activity seen in depression [4]. TMS lacks the systematic side effects associated with pharmacological antidepressant treatments and is generally well tolerated and safe. Furthermore, there is emerging literature supporting the use of TMS for a wide range of neuropsychiatric disorders for which pharmacotherapy offers limited efficacy (reviewed in [5,6]).
Despite this therapeutic promise, TMS remains an underutilized psychiatric treatment. Challenges limiting more widespread adoption of TMS include the time commitment for patients and clinicians and variable responses to treatment. A standard course of TMS requires daily in-clinic treatments for thirty sessions, usually followed by a taper series [3]. Standard treatment delivery limits the number of patients that can be treated on a single device to 60–80 per year, making TMS a limited resource that can only be administered to a small percentage of individuals with treatment-resistant depression. While recent studies have suggested that faster responses might be obtained with additional daily sessions, this literature remains mixed, and the time commitments are still substantial (e.g., [7–9]). Furthermore, response to TMS is variable, as some patients respond rapidly and robustly, whereas others demonstrate slower responses or no response at all [10,11].
Several studies have investigated the potential of utilizing neuroimaging techniques to develop biomarkers of TMS response that could identify patients more likely to respond [12–14]. However, these assessments remain exploratory research tools that are challenging to implement and have not yet demonstrated clinical utility [15]. There have also been efforts to use early symptomatic improvement to predict a later response to TMS [16,17], but findings have been mixed when examined in different settings [18]. To this end, a generalized model of TMS treatment response could inform the testing and development of measurement-based care strategies to optimize outcomes and identify patients more likely to respond to additional TMS treatments, thus increasing the percentage of successful courses of treatment.
Pharmacological antidepressant treatments have been shown to exhibit a nonlinear pattern of treatment response, with relatively large improvements in symptoms in the first weeks followed by smaller but continued improvements approaching a plateau [19–21]. This decrease in depressive symptoms has been modeled using an exponential decay function across multiple medications [22–24]. The exponential decay model suggests a defined relationship between early clinical response and eventual treatment outcomes. However, it is unclear if this model of treatment response can be extrapolated to other modalities, such as TMS for treatment-resistant depression.
The objectives of this study were to investigate whether an exponential decay function could model the treatment response to TMS using a clinical cohort sample and a broad range of recent TMS studies for treatment-resistant depression from around the world, and to test whether this model could predict clinical outcomes based on early clinical response. This line of inquiry aims to develop a framework for understanding the pattern of recovery from depression with TMS to inform clinical decisions and future protocols.
Overview
We report on data from our retrospective clinical sample of 97 patients undergoing TMS for treatment resistant depression, as well as novel modeling analyses of data from published clinical trials of TMS protocols. Primary outcome measures were longitudinal ratings of depressive symptoms acquired during TMS protocols for treatment resistant depression. The primary analyses utilize nonlinear mixed-effects models (NLME) to apply the exponential decay function to individual-level and group-level data of TMS treatment response and compare these nonlinear models against linear models. Additionally, we assess the predictive utility of rearranging the exponential decay function to predict treatment outcomes based on early response.
Participants and Data Acquisition
De-identified symptom rating data were collected from a naturalistic sample of 97 patients (mean age 53.76, SD = 13.20; 14% female) treated with clinical TMS at the VA Providence Healthcare System (Providence, RI, USA). Stimulation protocols followed clinical practice, with most patients receiving 10 Hz TMS delivered to the left dorsolateral prefrontal cortex. TMS was administered every weekday for up to 30 sessions, followed by six taper sessions. Symptom ratings were assessed using the Patient Health Questionnaire 9 (PHQ-9) [25] at baseline and after every five TMS treatment sessions. The PHQ-9 is a standard clinical rating scale for depression that provides adequate information about symptom severity and minimizes rater burden, with established psychometrics that facilitate its comparison with other depression rating scales. Scores range from 0 to 27, with higher scores indicative of more severe depression.
Group-level symptom rating data were also collected from published clinical trials investigating various TMS protocols. These studies included 1) a randomized noninferiority study comparing Theta Burst Stimulation (TBS) (n = 193) to standard 10 Hz stimulation (n = 192) [26], 2) a randomized controlled trial comparing standard and high-dose protocols of left-sided 10 Hz stimulation (standard n = 59, high-dose n = 91) and right-sided 1 Hz stimulation (standard n = 57, high-dose n = 93) [27], 3) a randomized controlled trial comparing supra- (120% motor threshold (MT)) and sub-threshold (80% MT) accelerated bilateral TBS (suprathreshold n = 87, subthreshold n = 93) to standard left-sided 10 Hz stimulation (n = 72) [9], 4) and an unblinded trial of accelerated, high-dose TBS using functional MRI guided targeting (n = 20) [7]. These studies contribute heterogenous data incorporating a variety of TMS paradigms, supporting the generalizability and relevance of our analysis. Additional group-level symptom rating data were collected from a secondary analysis of the THREE-D Study [26] in which unique group-level treatment response trajectories were identified in subjects receiving either left-sided 10 Hz TMS or left-sided TBS (n = 388) [10]. Primary outcome measures in these studies were the Hamilton Depression Rating Scale (HAM-D) [28] (17 item and 6 item), and the Quick Inventory of Depressive Symptoms Clinician-Rated (QIDS-C) [29]. Longitudinal group-level mean symptom rating data were obtained from the included tables [9,27] or extracted from figures using http://www.graphreader.com/ [7,10,26].
Statistical Analyses
The TMS treatment response was modeled using the exponential decay function in Eq. 1:
in which symptom ratings (D) at time t are described using the magnitude of total response (A), decaying at a time constant (B), and approaching a minimum value after treatment completion (C) as shown in Figure 1. In this model, parameters A and C are in the units of the symptom rating scale (i.e., score) and parameter B is time in weeks.
NLME models were constructed for our naturalistic sample receiving clinical TMS using Eq. 1. Model parameters were treated as random effects at the patient level and as fixed effects at the group level. A model treating all parameters (A, B, and C) as random effects was compared to the simpler model fitting only parameters A and C as random effects. The simpler model provided an equivalent fit as assessed by the likelihood ratio test (LRT) and lower Akaike information criterion (AIC) and Bayesian information criterion (BIC) values, and this simpler model was chosen. This NLME model was compared to a corresponding linear mixed-effect (LME) model in which slope and intercept were treated as random effects at the patient level using the AIC, BIC, and LRT. The predictive utility of this model was assessed using leave-one-out cross-validation (LOOCV), estimating the time constant (B) on a subset of the sample and using this estimate to calculate predicted treatment outcomes (C) for the left-out individual based on symptom rating scores at baseline and after five or ten TMS sessions.
Similar NLME models were constructed using the exponential decay function in Eq. 1 for three published group-level studies comparing various TMS protocols [9,26,27] with parameters A and C treated as random effects at the treatment group level. As above, these simpler models provided equivalent fits when compared to models using all parameters as random effects. These NLME models were subsequently compared to corresponding LME models in which the slope and intercept were treated as random effects at the group level using AIC, BIC, and LRT.
The exponential decay model was also fit at the group level using nonlinear least squares for each treatment group in the above studies [9,26,27] as well as the unblinded trial of accelerated, high-dose TBS [7]. Corresponding linear models were compared to nonlinear exponential decay models using AIC and BIC.
Finally, a NLME model was constructed using the exponential decay function for the unique group-level symptom trajectories previously identified in 388 subjects receiving either left-sided 10 Hz TMS or left-sided TBS [10]. Here, a model incorporating all parameters (A, B, C) as random effects yielded significant LRT and lower AIC and BIC values compared with a simpler NLME using only parameters A and C as random effects. This NLME model was then compared to a corresponding LME model in which slope and intercept were considered random effects at the group level using AIC, BIC, and LRT.
All mixed effect models were fit using maximum likelihood as implemented in the Nonlinear Mixed-Effects Models library in R [30].
Ethics declaration
This study was performed in accordance with the Declaration of Helsinki. The procedures for gathering and analyzing de-identified clinical data were approved as quality improvement by the VA Palo Alto/Stanford Institutional Review Board [31]. Additional group-level analyses were based on publicly available data from previously published studies [7,9,10,26,27].
Modeling Treatment Response to TMS in a Clinical Sample
The NLME model demonstrated that the TMS treatment response was well modeled with the exponential decay function, yielding significant estimates for model parameters A, B, and C (all p < 0.0001). Each individual contributed 2 to 9 (median = 6) longitudinal measurements of depressive symptoms for a total of 562 observations. The exponential model captured a wide range of individual treatment response trajectories, as illustrated in Fig. 2. The magnitude of total response (A) was estimated to be a clinically meaningful 5.8-point drop in PHQ-9 scores with a time constant (B) of 1.2 weeks or 6 TMS treatments. When compared to a corresponding LME, the exponential decay model displayed lower AIC and BIC values and a significant likelihood ratio (LRT = 63.2, p < 0.0001), suggesting that the NLME model based on the exponential decay function is a better fit.
Predicting Treatment Outcomes to TMS based on Early Response
In subsets containing complete longitudinal data from this naturalistic sample, LOOCV yielded consistent estimates of time constant B (mean 1.21 weeks, sd 0.023). By rearranging Eq. 1 using this time constant estimate (B) and symptom scores at baseline and after one week (n = 90) or two weeks (n = 82) of treatment, we calculated predicted values of C for the left-out patients, as shown in Eq. 2, constrained to scores of 0 to 27, inclusive.
Predicted C values yielded significant correlations with PHQ-9 scores at the end of treatment, accounting for 38–59% of the variance in final scores, using PHQ-9 scores at baseline and after one week (adjusted R2 = 0.383, p < 0.0001) or two weeks (adjusted R2 = 0.585, p < 0.0001) of daily TMS sessions.
Modeling Group-Level Treatment Response Patterns to TMS
The NLME model was extended to group-level data from three published studies comparing multiple TMS modalities for treatment-resistant depression [9,26,27], including left-sided high-frequency (10 Hz) repetitive TMS, standard and accelerated TBS, and right-sided, low-frequency (1 Hz) stimulation. These NLME models using the exponential decay function yielded significant estimates for all model parameters (all p < 0.005).
In a randomized noninferiority study comparing TBS to 10 Hz stimulation, the “THREE-D” study [26], the NLME model estimated the magnitude of total response (A) to be a drop of 12.2 points on the HAM-D with a time constant of 2.3 weeks. When compared to the corresponding LME model, the NLME model yielded a significant likelihood ratio (LRT = 18.1, p < 0.0001) and lower AIC (49 vs 65) and BIC values (53 vs 69), demonstrating the NLME model is a better fit.
Similarly, when this exponential decay model was applied to a randomized controlled trial comparing left-sided 10 Hz stimulation and right-sided 1 Hz stimulation [27], estimates of the treatment response (A) demonstrated an average 13.6 point drop on the HAM-D with a time constant B of 1.2 weeks. This NLME model demonstrated a significant likelihood ratio (LRT = 22.7, p < 0·0001) and lower AIC (81 vs. 102) and BIC (88 vs. 108) values compared to the corresponding LME model, indicating the nonlinear model is a better fit.
Finally, when applied to a randomized controlled trial comparing accelerated bilateral TBS above (120%) and below (80%) MT to standard left-sided 10 Hz stimulation [9], the NLME model yielded estimates of a total 6.0 point drop on the QIDS-C [29] with a time constant of 0.7 weeks. When compared to the corresponding LME model, the NLME model yielded a significant likelihood ratio (LRT = 44.1, p < 0.0001) and lower AIC (30 vs 72) and BIC (35 vs 76) values, suggesting the exponential decay model was a superior fit.
The exponential decay model was also fit at the individual group level using nonlinear least squares for each group in the above studies [9,26,27], an unblinded accelerated, high dose TBS study [7] and our naturalistic clinical sample. Although these individual group fits did not always yield significant estimates for all model parameters based on limited degrees of freedom, when the AIC and BIC values from these nonlinear fits were compared to corresponding linear models of the same group-level data, exponential decay models consistently yielded lower AIC and BIC values for all TMS samples studied (Fig. 3).
Modeling Unique Group-Level Response Trajectories to TMS
When the exponential decay model was applied to the unique group-level symptom trajectories identified in the secondary analysis of the THREE-D study [10], the NLME model was able to describe each group trajectory in terms of parameters A, B, and C, (Eq. 1, Fig. 4). The three groups of responders demonstrated similar estimates of the magnitude of treatment response (parameter A), with reductions of 13.0, 13.9 and 14.5 points on the HAM-D for “Lower baseline symptoms, linear response,” “Higher baseline symptoms, linear response,” and “Rapid response” groups, respectively. However, these groups differed more markedly in their group estimates of the time constant B, with “Rapid response” at 1.1 weeks, followed by “Lower baseline symptoms, linear response” at 2.4 weeks, and “Higher baseline symptoms, linear response” at 3.8 weeks. These differences correspond to the speed at which these groups approach their minimum depressive symptom rating scores during TMS treatment. In contrast, the “Nonresponse” group was associated with a minimal estimate of total treatment response (parameter A) of a 2.7-point reduction on the HAM-D and a long estimate for the time constant B of 7.5 weeks. This NLME model was associated with a significant likelihood ratio (LRT = 30.7, p < 0.0001) and lower AIC (121 vs. 143) and BIC (134 vs. 151) values compared to a corresponding LME model, indicating that these trajectories were better modeled with the nonlinear exponential decay function.
These results provide substantial evidence that the antidepressant response to TMS demonstrates a nonlinear pattern of improvement that follows an exponential decay function. This pattern of early large reductions in depressive symptom ratings followed by smaller improvements is similar to the trajectory seen with pharmacological depression treatments [19–23]. This pattern is further replicated at the group level in multiple clinical TMS trials utilizing various stimulation protocols and clinical rating scales representing over 1000 individuals [7,9,26,27]. This study also demonstrated that this model can provide a parsimonious method to describe a variety of treatment response trajectories to TMS, including those previously identified using group-based trajectory modeling [10]. Finally, this study provides evidence that the exponential decay model can be applied to define and predict individual treatment response outcomes based on early clinical response.
An advantage of this nonlinear modeling approach is that it provides fitted parameters that have real-world meaning. In this model, A is the total magnitude of response as measured by a point reduction in clinical rating scales, B is the time constant of the response representing the speed of recovery as a function of time, and C is the minimum expected symptom rating at the end of treatments. This simple quantification of the magnitude, speed, and outcome of TMS treatment response provides a comprehensive estimate of the treatment course that can be readily translated into common metrics of clinical outcomes for depression treatment; categorical treatment response (50% reduction in symptoms) would occur when A is greater than or equal to C, and remission would occur when C is below a rating scale threshold (e.g., PHQ9 < 5). Additionally, this model incorporates the relationship between the speed of improvement and length of TMS treatment, providing a formula that relates early treatment response to eventual treatment outcomes. At time B, one would expect the magnitude of treatment response to reflect ~ 63% (1–1/e) of the total response. At time 3xB, this would reflect 95% (1–3/e) of the total expected response. While most of the groups included in this study only incorporated time points during TMS treatment, the study by Chen et al. [9] also included longitudinal measurements weeks after treatment, implying estimates of C may extend to post-treatment outcomes.
Using our clinical TMS sample, we demonstrated that one could use a population average value for the time constant B to predict individual treatment outcomes based on the response at weeks one or two. This work is consistent with previous studies demonstrating that early TMS treatment response can predict later treatment outcomes [16,17]. However, it is unlikely a single time constant will describe all individuals’ treatment response patterns adequately. Additional studies using well-sampled longitudinal measurements integrating this variability in speed of recovery are needed to further elucidate the utility of this model at the individual level. At the group-level, we found evidence that the time constant of recovery may differ based on the population, treatment protocol, and the depression rating scale used. Additionally, applying this model to the group trajectories identified by Kaster et al. [10] suggests that variability in the speed of recovery reflected by this time constant is a key difference defining the unique group and, by extrapolation, individual treatment response trajectories. For an individual undergoing TMS, one could estimate this time constant using nonlinear least squares to fit the exponential decay model to a target outcome (C) to determine whether a patient could benefit from additional TMS treatments given slower early response in concordance with current clinical recommendations [31].
The limitations of this work are its use of available datasets and retrospective analysis; future research will be needed to test whether these findings can be used prospectively to determine optimal treatment courses. Furthermore, while TMS, given its resource-intensive nature, is an appropriate context for developing this model, novel pharmacotherapies and psychotherapies could also benefit from application of the model. Whether and how these models differ from placebo response will be an important area of future inquiry. These limitations aside, these data clearly demonstrate the nonlinear trajectory of symptom improvement with TMS for treatment resistant depression across many different protocols and illustrate methods to predict treatment outcomes based on early response. This modeling offers a simple and useful framework to describe the antidepressant response to TMS that could inform clinical decisions and future studies.
Funding:
This work was supported by the National Institute of Mental Health (R25 MH101076, R01 MH120126), National Institute of General Medical Science (U54GM115677, P20 GM130452), Department of Veterans Affairs (IK2CX002115), and VA RR&D Center for Neurorestoration and Neurotechnology.
Acknowledgments: We thank Michelle R. Madore, William D. Rooney, Linda L. Carpenter, Benjamin D. Greenberg, Robert M. Swift, and the Career Development Institute for Psychiatry. The opinions herein represent those of the authors and not the U.S. Department of Veterans Affairs.
Author contributions:
Conceptualization: YAB, AZ, LHP, NSP
Formal Analysis: YAB
Methodology: YAB, AZ
Supervision: NSP
Validation: AZ
Visualization: YAB
Writing – original draft: YAB, NSP
Writing – review & editing: YAB, AZ, MCB, LMW, LHP, NSP
Data and materials availability: Data from the clinical sample and analysis code are available from the corresponding author on reasonable request. Group-level data are available in the main text and in the text of previously published studies [7,9,10,26,27].
Competing interests:
Authors declare that they have no competing interests.
- Friedrich, M.J. Depression is the leading cause of disability around the world. JAMA317(15), 1517; 10.1001/jama.2017.3826 (2017).
- Rush, A.J., et al. Acute and longer-term outcomes in depressed outpatients requiring one or several treatment steps: a STAR*D report. Am J Psychiatry163(11), 1905-17; 10.1176/ajp.2006.163.11.1905 (2006).
- O'Reardon, J.P., et al. Efficacy and safety of Transcranial Magnetic Stimulation in the acute treatment of major depression: a multisite randomized controlled trial. Biol Psychiatry62(11), 1208-16;10.1016/j.biopsych.2007.01.018 (2007).
- Liston, C., et al. Default mode network mechanisms of Transcranial Magnetic Stimulation in depression. Biol Psychiatry76(7), 517-26; 10.1016/j.biopsych.2014.01.023 (2014).
- Cosmo, C., Zandvakili, A., Petrosino, N.J., Berlow, Y.A. & Philip, N.S. Repetitive Transcranial Magnetic Stimulation for treatment-resistant depression: Recent critical advances in patient care. Curr Treat Options Psychiatry8(2), 47-63; 10.1007/s40501-021-00238-y (2021).
- George, M.S. Whither TMS: A one-trick pony or the beginning of a neuroscientific revolution? Am J Psychiatry176(11), 904-10; 10.1176/appi.ajp.2019.19090957 (2019).
- Cole, E.J., et al. Stanford accelerated intelligent neuromodulation therapy for treatment-resistant depression. Am J Psychiatry177(8), 716-26.; 10.1176/appi.ajp.2019.19070720 (2020).
- Blumberger, D.M., et al. A randomized sham controlled comparison of once vs twice-daily intermittent theta burst stimulation in depression: A Canadian rTMS treatment and biomarker network in depression (CARTBIND) study. Brain Stimul14(6), 1447-55; 10.1016/j.brs.2021.09.003 (2021).
- Chen, L., et al. Accelerated theta burst stimulation for the treatment of depression: A randomised controlled trial. Brain Stimul14(5), 1095-105;10.1016/j.brs.2021.07.018 (2021).
- Kaster, T.S., et al. Trajectories of response to dorsolateral prefrontal rTMS in major depression: A THREE-D Study. Am J Psychiatry176(5), 367-75; 10.1176/appi.ajp.2018.18091096 (2019).
- Yip, A.G., George, M.S., Tendler, A., Roth, Y., Zangen, A. & Carpenter, L.L. 61% of unmedicated treatment resistant depression patients who did not respond to acute TMS treatment responded after four weeks of twice weekly deep TMS in the Brainsway pivotal trial. Brain Stimul10(4), 847-9; 10.1016/j.brs.2017.02.013 (2017).
- Drysdale, et al. Resting-state connectivity biomarkers define neurophysiological subtypes of depression. Nat Med23(1), 28-38; 10.1038/nm.4246 (2017).
- Zandvakili, A., et al. Use of machine learning in predicting clinical response to transcranial magnetic stimulation in comorbid posttraumatic stress disorder and major depression: A resting state electroencephalography study. J Affect Disord252, 47-54; 10.1016/j.jad.2019.03.077 (2019).
- Williams, et al. Identifying response and predictive biomarkers for transcranial magnetic stimulation outcomes: Protocol and rationale for a mechanistic study of functional neuroimaging and behavioral biomarkers in veterans with Pharmacoresistant depression. BMC Psychiatry21(1), 35; 10.1186/s12888-020-03030-z (2021).
- Fitzgerald, P.B. Targeting repetitive transcranial magnetic stimulation in depression: Do we really know what we are stimulating and how best to do it? Brain Stimul14(3), 730-6; 10.1016/j.brs.2021.04.018 (2021).
- Feffer, K., et al. Early symptom improvement at 10 sessions as a predictor of rTMS treatment outcome in major depression. Brain Stimul11(1), 181-9; 10.1016/j.brs.2017.10.010 (2018).
- Mondino, M., et al. Predicting treatment response to 1Hz rTMS using early self-rated clinical changes in major depression. Brain Stimul13(6), 1603-5; 10.1016/j.brs.2020.10.004 (2020).
- Beck, Q.M., Tirrell, E., Fukuda, A.M., Kokdere, F. & Carpenter, L.L.. Can early treatment response serve as a predictor of antidepressant outcome of repetitive transcranial magnetic stimulation? Brain Stimul13(2), 420-1; 10.1016/j.brs.2019.12.002 (2020).
- Taylor, M.J., Freemantle, N., Geddes, J.R. & Bhagwagar, Z. Early onset of selective serotonin reuptake inhibitor antidepressant action: Systematic review and meta-analysis. Arch Gen Psychiatry63(11), 1217-23; 10.1001/archpsyc.63.11.1217 (2006).
- Lam, R.W. Onset, time course and trajectories of improvement with antidepressants. Eur Neuropsychopharmacol22 Suppl 3, S492-8; 10.1016/j.euroneuro.2012.07.005 (2012).
- Posternak, M.A. & Zimmerman, M. Is there a delay in the antidepressant effect? A meta-analysis. J Clin Psychiatry66(2), 148-58; 10.4088/jcp.v66n0201 (2005).
- Priest, R.G., et al. Recovery from depressive illness does fit an exponential model. J Clin Psychopharmacol16(6), 420-4; 10.1097/00004714-199612000-00002 (1996).
- Boyer, W. Antidepressant response follows an exponential decay curve: Replication and extension. J Clin Psychopharmacol17(5), 426-7;10.1097/00004714-199710000-00015 (1997).
- Friede, T., Kieser, M. & Miller, F. Modeling the recovery from depressive illness by an exponential model with mixed effects. Methods Inf Med39(1), 12-5; 10.1055/s-0038-1634261 (2000).
- Kroenke, K., Spitzer, R.L. & Williams, J.B. The PHQ-9: Validity of a brief depression severity measure. J Gen Intern Med16(9), 606-13;10.1046/j.1525-1497.2001.016009606.x (2001).
- Blumberger, D.M., et al. Effectiveness of theta burst versus high-frequency repetitive transcranial magnetic stimulation in patients with depression (THREE-D): A randomised non-inferiority trial. Lancet391(10131), 1683-92;10.1016/S0140-6736(18)30295-2 (2018).
- Fitzgerald, P.B., et al. A pragmatic randomized controlled trial exploring the relationship between pulse number and response to repetitive transcranial magnetic stimulation treatment in depression. Brain Stimul13(1), 145-52; 10.1016/j.brs.2019.09.001 (2020).
- Hamilton, M. A rating scale for depression. J Neurol Neurosurg Psychiatry23, 56-62; 10.1136/jnnp.23.1.56 (1960).
- Rush, A.J., et al. The 16-item quick inventory of depressive symptomatology (QIDS), clinician rating (QIDS-C), and self-report (QIDS-SR): A psychometric evaluation in patients with chronic major depression. Biol Psychiatry54(5), 573-83; 10.1016/s0006-3223(02)01866-8 (2003).
- Pinheiro, J., Bates, D., DebRoy, S., Sarkar, D. & R Core Team. NLME: Linear and nonlinear mixed effects models. https://CRAN.R-project.org/package=nlme (2019).
- Madore, M. R. et al. Prefrontal Transcranial Magnetic Stimulation for Depression in US Military Veterans - A Naturalistic Cohort Study in the Veterans Health Administration. J Affect Disord, doi:10.1016/j.jad.2021.10.025 (2021).
- Perera, T., et al. The clinical TMS society consensus review and treatment recommendations for TMS therapy for major depressive disorder. Brain Stimul 9(3), 336-46; 10.1016/j.brs.2016.03.010 (2016).
No competing interests reported.