psychology research

profileGhazi33
ContentServer.pdf

Timing Matters: Change Depends on the Stage of Treatment in Cognitive Behavioral Therapy for Panic Disorder With Agoraphobia

Andrew T. Gloster Technische Universität Dresden and University of Basel

Jens Klotsche Technische Universität Dresden and German Rheumatism

Research Centre Berlin

Alexander L. Gerlach University of Cologne

Alfons Hamm Ernst Moritz Arndt University of Greifswald

Andreas Ströhle Charité-Universitätsmedizin Berlin

Siegfried Gauggel University of Aachen

Tilo Kircher Philipps-University Marburg

Georg W. Alpers University of Mannheim

Jürgen Deckert University of Würzburg

Hans-Ulrich Wittchen Technische Universität Dresden

Objective: The mechanisms of action underlying treatment are inadequately understood. This study examined 5 variables implicated in the treatment of panic disorder with agoraphobia (PD/AG): cata- strophic agoraphobic cognitions, anxiety about bodily sensations, agoraphobic avoidance, anxiety sen- sitivity, and psychological flexibility. The relative importance of these process variables was examined across treatment phases: (a) psychoeducation/interoceptive exposure, (b) in situ exposure, and (c) generalization/follow-up. Method: Data came from a randomized controlled trial of cognitive behavioral therapy for PD/AG (n � 301). Outcomes were the Panic and Agoraphobia Scale (Bandelow, 1995) and functioning as measured in the Clinical Global Impression scale (Guy, 1976). The effect of process variables on subsequent change in outcome variables was calculated using bivariate latent difference score modeling. Results: Change in panic symptomatology was preceded by catastrophic appraisal and agoraphobic avoidance across all phases of treatment, by anxiety sensitivity during generalization/follow- up, and by psychological flexibility during exposure in situ. Change in functioning was preceded by agoraphobic avoidance and psychological flexibility across all phases of treatment, by fear of bodily symptoms during generalization/follow-up, and by anxiety sensitivity during exposure. Conclusions: The effects of process variables on outcomes differ across treatment phases and outcomes (i.e., symptomatology vs. functioning). Agoraphobic avoidance and psychological flexibility should be in- vestigated and therapeutically targeted in addition to cognitive variables.

Keywords: mechanism of action, panic disorder, agoraphobia, avoidance, psychological flexibility

This article was published Online First November 4, 2013. Andrew T. Gloster, Institute of Clinical Psychology and Psychotherapy,

Technische Universität Dresden, Dresden, Germany, and Department of Psy- chology, Division of Clinical Psychology and Epidemiology, University of Basel, Basel, Switzerland; Jens Klotsche, Institute of Clinical Psychology and Psychotherapy, Technische Universität Dresden, Dresden, Germany, and Ep- idemiology Unit, German Rheumatism Research Centre Berlin, Leibnitz In- stitute, Berlin, Germany; Alexander L. Gerlach, Department of Psychology, Clinical Psychology, and Psychotherapy, University of Cologne, Cologne, Germany; Alfons Hamm, Institute for Psychology, Physiological and Clinical Psychology/Psychotherapy, Ernst Moritz Arndt University of Greifswald, Greifswald, Germany; Andreas Ströhle, Department of Psychiatry, Charité- Universitätsmedizin Berlin, Berlin, Germany; Siegfried Gauggel, Institute for Medical Psychology and Medical Sociology, University of Aachen, Aachen, Germany; Tilo Kircher, Clinic for Psychiatry and Psychotherapy, Philipps-

University Marburg, Marburg, Germany; Georg W. Alpers, Clinical and Biological Psychology and Psychotherapy, University of Mannheim, Mann- heim, Germany; Jürgen Deckert, Department of Psychiatry, Psychosomatics, and Psychotherapy, University of Würzburg, Würzburg, Germany; Hans- Ulrich Wittchen, Institute of Clinical Psychology and Psychotherapy, Tech- nische Universität Dresden, Dresden, Germany.

This work is part of the German multicenter trial “Mechanisms of Action in CBT” (MAC). The MAC study is funded by the German Federal Ministry of Education and Research (BMBF; Project No. 01GV0615) as part of the BMBF Psychotherapy Research Funding Initiative. The randomized clinical trial (RCT) project was approved by the Ethics Committee of the Medical Faculty of the Technical University of Dresden, Germany (EK 164082006). The neuroimaging components were approved by the Ethics Committee of the Medical Faculty of the Rheinisch-Westfälische Hochschule University Aachen, Germany (EK 073/07). The experimental pharmacology study

T hi

s do

cu m

en t

is co

py ri

gh te

d by

th e

A m

er ic

an Ps

yc ho

lo gi

ca l

A ss

oc ia

tio n

or on

e of

its al

lie d

pu bl

is he

rs .

T hi

s ar

tic le

is in

te nd

ed so

le ly

fo r

th e

pe rs

on al

us e

of th

e in

di vi

du al

us er

an d

is no

t to

be di

ss em

in at

ed br

oa dl

y.

Journal of Consulting and Clinical Psychology © 2013 American Psychological Association 2014, Vol. 82, No. 1, 141–153 0022-006X/14/$12.00 DOI: 10.1037/a0034555

141

The efficacy of cognitive behavioral therapy (CBT) for panic disorder (PD) and agoraphobia (AG) is undeniable, yet the crucial task of elucidating the mechanisms of action lags behind. Far from an abstract theoretical concern, understanding the mechanisms of action of treatment would provide knowledge on how to generalize the principles utilized in treatment studies and offer hope for the sizable minority of patients that do not respond to current treat- ments (Hofmann & Smits, 2008) by amplifying those specific processes known to affect outcome. Toward this aim, conceptual and methodological clarity are crucial.

Much theoretical and empirical work has been devoted to ex- plaining the nature and treatment of PD and AG. Cognitive ac- counts suggest that catastrophic misinterpretation of bodily sensa- tions influences the etiology and maintenance of PD and AG (Clark, 1986). Cognitive therapy thus targets the content and frequency of associated thoughts through numerous methods. Anx- iety sensitivity, or the fear of anxiety and fear, has also been conceptualized as a risk factor associated with the subsequent onset of panic disorder (Ehlers, 1995) and with avoidance behav- iors (Zvolensky & Forsyth, 2002). Anxiety sensitivity is conse- quently considered an important therapeutic target for PD/AG (Smits, Powers, Cho, & Telch, 2004) achieved through various means including interoceptive exposure. However, the pernicious effects of both cognitive appraisals and anxiety sensitivity depend in part on how an individual attempts to regulate their negative affect (Kashdan, Zvolensky, & McLeish, 2008). This suggests that successful therapy must also target the way one interacts with these negative appraisals, beliefs, and emotions.

A common regulation strategy for these negative appraisals and emotions is avoidance. Indeed, agoraphobic avoidance, or the avoidance of feared situations, is a defining feature of agoraphobia (Chambless, Caputo, Jasin, Gracely, & Williams, 1985) even in the absence of PD (Wittchen, Gloster, Beesdo-Baum, Fava, & Craske, 2010; Wittchen et al., 2008). Although not always explicitly tar-

geted (Hofmann & Spiegel, 1999), reduction of agoraphobic avoidance is a common therapeutic target and is associated with successful outcome (Gloster, Wittchen, et al., 2011). Another regulatory strategy recently implicated in PD/AG is psychological flexibility. Psychological flexibility refers to the ability to mind- fully accept cognitions and emotions when doing so is useful for living a meaningful life (Bond et al., 2011). Similar to anxiety sensitivity, psychological flexibility is not exclusively relevant to PD/AG, yet it has been implicated in panic-related distress (Kar- ekla, Forsyth, & Kelly, 2004), baseline functioning in anxiety disorders (Gloster, Klotsche, Chaker, Hummel, & Hoyer, 2011), and treatment outcome (Forman, Herbert, Moitra, Yeomans, & Geller, 2007), and it is conceptually distinct from anxiety sensi- tivity in patients with PD/AG (Kämpfe et al., 2012).

Given that appraisal of anxiety symptoms, anxiety sensitivity, avoidance, and psychological flexibility are all associated with various aspects of PD/AG, it is important to understand to what degree some or all these constructs are active mechanisms for successful treatment outcome. Mediation analysis (Baron & Kenny, 1986) has emerged as one important analytical procedure for the critical testing of putative mechanisms of action in therapy (Kazdin, 2007), and a handful of formal mediation analyses have been conducted across variations of CBT for PD/AG (e.g., group vs. individual therapy). These studies provide positive evidence for the mediating or partially mediating role of cognitive content, cognitive appraisal, and self-efficacy (Casey, Newcombe, & Oei, 2005; Hofmann et al., 2007; Meulenbeek, Spinhoven, Smit, Van Balkom, & Cuijpers, 2010; Vögele et al., 2010) and anxiety sensitivity (Smits et al., 2004) in reducing the severity of panic disorder. Although an important step toward isolating active mech- anisms, these findings are limited by the fact that the assessment of target variables did not precede outcome assessments. That is, the process variables were tested concurrently with the outcome mea- sure (e.g., both measured pre–post). This lack of temporal order

was approved by the Ethics Committee of the state of Berlin, Germany (EudraCT: 2006-00-4860-29). The study was registered with the Interna- tional Standard Randomised Controlled Trial Number (ISRCTN): IS- RCTN80046034.

Principal investigators (PIs) with respective areas of responsibility in the MAC study are V. Arolt (Münster, Germany: overall MAC program coordi- nation), Hans-Ulrich Wittchen (Dresden, Germany: PI for the RCT and man- ual development), Alfons Hamm (Greifswald, Germany: PI for psychophysi- ology), Alexander L. Gerlach (Münster, Germany: PI for psychophysiology and panic subtypes), Andreas Ströhle (Berlin, Germany: PI for experimental pharmacology), Tilo Kircher (Marburg, Germany: PI for functional neuroim- aging), and Jürgen Deckert (Würzburg, Germany: PI for genetics). Additional site directors in the RCT component of the program are Georg W. Alpers (Würzburg, Germany), T. Fydrich, and L. Fehm (Berlin–Adlershof, Germany) as well as T. Lang (Bremen, Germany). All principle investigators take responsibility for the integrity of the respective study data and their compo- nents. All authors and coauthors had full access to all study data. Data analysis and manuscript preparation were completed by the authors and coauthors of this article, who take responsibility for its accuracy and content.

Acknowledgments and staff members by site are as follows: Greifswald, Germany (coordinating site for psychophysiology): Christiane Melzig, Jan Richter, Susan Richter, and Matthias von Rad; Berlin–Charité, Germany (coordinating center for experimental pharmacology): Harald Bruhn, Anja Siegmund, Meline Stoy, and Andre Wittmann; Berlin–Adlershof, Germany:

Irene Schulz; Münster, Germany (overall MAC program coordination, genet- ics, and functional neuroimaging): Andreas Behnken, Katharina Domschke, Adrianna Ewert, Carsten Konrad, Bettina Pfleiderer, and Peter Zwanzger; Münster, Germany (coordinating site for psychophysiology and subtyping): Judith Eidecker, Swantje Koller, Fred Rist, and Anna Vossbeck-Elsebusch; Marburg/Aachen, Germany (coordinating center for functional neuroimaging): Barbara Drüke, Sonja Eskens, Thomas Forkmann, Siegfried Gauggel, Susan Gruber, Andreas Jansen, Thilo Kellermann, Isabelle Reinhardt, and Nina Vercamer-Fabri; Dresden, Germany (coordinating site for data collection, analysis, and the RCT): Franziska Einsle, Christine Fröhlich, Andrew T. Gloster, Christina Hauke, Simone Heinze, Michael Höfler, Ulrike Lueken, Peter Neudeck, Stephanie Preiß, and Dorte Westphal; Würzburg Psychiatry Department, Würzburg, Germany (coordinating center for genetics): Andreas Reif; Würzburg Psychology Department, Würzburg, Germany: Julia Dürner, Hedwig Eisenbarth, Antje B. M. Gerdes, Harald Krebs, Paul Pauli, Silvia Schad, and Nina Steinhäuser; Bremen, Germany: Veronika Bamann, Sylvia Helbig-Lang, Anne Kordt, Pia Ley, Franz Petermann, and Eva-Maria Schröder. Additional support was provided by the coordinating center for clinical studies in Dresden, Germany (KKS Dresden): Xina Grählert and Marko Käppler. Correspondence concerning this article should be addressed to Andrew T. Gloster, University of Basel, Division of Clinical Psychology and Epide- miology, Missionsstrasse 62A, CH-4055 Basel, Switzerland. E-mail: [email protected]

T hi

s do

cu m

en t

is co

py ri

gh te

d by

th e

A m

er ic

an Ps

yc ho

lo gi

ca l

A ss

oc ia

tio n

or on

e of

its al

lie d

pu bl

is he

rs .

T hi

s ar

tic le

is in

te nd

ed so

le ly

fo r

th e

pe rs

on al

us e

of th

e in

di vi

du al

us er

an d

is no

t to

be di

ss em

in at

ed br

oa dl

y.

142 GLOSTER ET AL.

hinders interpretation because it is unclear if the outcome variables also influence the process variables and because such designs leave open the possibility that process variables exert their effects at different points during the therapy.

Longitudinal temporal testing of putative mechanisms within the course of treatment for PD/AG has been examined only in a few studies. In one such study, 12 patients diagnosed with PD/AG completed daily diaries for 30 weeks (Bouchard et al., 2007). CBT was administered in groups and emphasized either cognitive or exposure interventions. Daily dairies were used to assess beliefs about the consequences of panic, self-efficacy to control panic attacks in the face of bodily sensations and catastrophic thoughts, and anticipatory anxiety about having a panic attack that day. During the course of therapy, all 12 patients recorded changes in their beliefs and level of self-efficacy prior to recording changes in anticipatory anxiety, irrespective of condition. Despite the small sample size, this study demonstrated with temporal sensitivity that changes in cognitive variables preceded change in other aspects of symptomatology for all patients, though the magnitude of change differed across patients.

To our knowledge, only a few further studies temporally exam- ined whether salient process variables preceded subsequent change in panic-related outcomes. Using cognitive therapy and guided mastery—both administered in a group format—Hoffart (1995) examined the relevance of self-efficacy, catastrophic beliefs, and perceived control of thoughts on subsequent fear in a behavioral avoidance test (BAT). Results from the 46 patients indicated that change in self-efficacy was the strongest and most consistent predictor of subsequent change in fear. A second study (Teach- man, Marker, & Clerkin, 2010) examined whether catastrophic misinterpretations subsequently affected various facets of panic symptomatology. Panic control treatment was administered in a group format to 43 patients. Using bivariate difference score modeling analysis, results indicated that change in catastrophic misinterpretation predicted subsequent change in panic symptom- atology. The reverse pattern (i.e., symptomatology predicting sub- sequent change in catastrophic misinterpretation) was not consis- tently found except for distress/apprehension. A final study examined cognitive process variables during the first phase of treatment (i.e., 4 weeks) in 41 patients diagnosed with PD/AG (Meuret, Rosenfield, Seidel, Bhaskara, & Hofmann, 2010) across two distinct treatment conditions. Process variables were opera- tionalized as a composite score of questionnaires that measure anxiety sensitivity/fear about the consequences of panic (i.e., Anx- iety Sensitivity Index [ASI; Peterson & Reiss, 1993] and Body Sensation Questionnaire [BSQ; Chambless, Caputo, Bright, & Gallagher, 1984], respectively) in addition to perceived control. During the phase of treatment examined, patients received either cognitive treatment or capnometry-assisted respiratory training. So designed, the specificity of the cognitive process variables could be tested across relatively pure intervention conditions. Indeed, results suggested that cognitions were bidirectionally associated with changes in panic severity only in the cognitive training condition, whereas perceived control was bidirectionally associ- ated with panic symptom change in both conditions. This excellent study included information only from the first half of treatment (4 weeks), however, thus limiting information about how mecha- nisms unfold over the full course of treatment or generalize fol- lowing treatment. Taken together, these studies provide strong

support for the role of cognitively oriented process variables defined as catastrophic misinterpretations and self-efficacy in the prediction of subsequent change in symptomatology.

To our knowledge, no other PD/AG relevant process variables than those discussed above (i.e., feared consequences/anxiety sen- sitivity and self-efficacy/perceived control to cope with panic) have been tested longitudinally. Examination and direct compari- son of other variables implicated in the treatment of PD/AG such as avoidance behavior and psychological flexibility is a crucial step in the process of understanding the mechanisms of treatment (Kazdin, 2007).

In addition to expanding the scope of process variables under investigation, outcome variables also need to be expanded. To date, all process studies examined the effect on panic-related variables, but only one examined how putative mechanisms affect other outcomes such as functioning (Smits et al., 2004). It remains an open question whether the mechanisms of action involved in symptom reduction are identical in importance and sequence to those involved in other treatment targets. Social, occupational, and psychological functioning are certainly related to symptomatology, yet it is a broader measuring stick. Indeed, the impetus for patients to seek therapy may be primarily related to functioning and in our quest to better understand mechanisms care should be taken not to reduce patients to their symptomatology.

The purpose of the present study was to investigate the degree to which five process variables affect treatment outcome across the active and follow-up phases of a standardized CBT for PD/AG. Toward this end, the process variables were examined across phases of therapy for two outcomes: severity of PD/AG symptom- atology and overall functioning. The process variables were ex- amined longitudinally using bivariate latent difference score mod- eling to determine the relative effects of the process variables at different points in the therapeutic process. So doing, the relative importance of the process variables were examined for their rela- tionship to different components of the therapy. We hypothesized that the variables would differentially predict subsequent symp- toms as a function of treatment phase (i.e., psychoeducation, functional analysis, interoceptive exposure [pre-treatment to inter- mediate assessment following the fourth session]; exposure in situ, anticipatory anxiety and specified interoceptive exposure [in- termediate assessment to post-treatment following the 12th session]; and generalization period with two booster sessions that reviewed progress, helped set goals, and addressed diffi- culties [post-treatment to 6-month follow-up period]) and out- come variable (panic and agoraphobia symptoms vs. function- ing). Specifically, we predicted that (a) cognitive appraisal would predict subsequent change in panic symptoms, but not functioning, only during Sessions 1– 4 because these sessions addressed psychoeducation and engaged the patient intellectu- ally; (b) anxiety elicited by bodily symptoms would predict subsequent change in panic symptom, but not functioning, only during Sessions 1– 4 because these sessions introduced intero- ceptive exposure; (c) anxiety sensitivity would predict subse- quent change panic symptoms during Sessions 1–12 because of interoceptive exposure during the first four sessions and expo- sure in situ during second half of treatment and functioning during treatment because improvement in function is likely related to a new relationship with the anxiety; (d) avoidance behavior would predict subsequent change in both panic symp-

T hi

s do

cu m

en t

is co

py ri

gh te

d by

th e

A m

er ic

an Ps

yc ho

lo gi

ca l

A ss

oc ia

tio n

or on

e of

its al

lie d

pu bl

is he

rs .

T hi

s ar

tic le

is in

te nd

ed so

le ly

fo r

th e

pe rs

on al

us e

of th

e in

di vi

du al

us er

an d

is no

t to

be di

ss em

in at

ed br

oa dl

y.

143TIMING MATTERS

toms and functioning across Sessions 1–12 and the follow-up period because avoidance was directly and intensively targeted in the therapy as a maintaining factor; and (e) and psychological flexibility would predict subsequent change in panic symptoms during Sessions 5–12 because it is believed to facilitate expo- sure and functioning during all phases because it is closely tied to functioning.

Method

Design

Data were collected within the Mechanisms of Action for CBT (MAC) study. The MAC study was a multicenter, randomized controlled trial for patients with PD/AG. The methods and main outcomes of the study were published elsewhere (Gloster et al., 2009; Gloster, Wittchen, et al., 2011). The MAC study was ap- proved by the internal review board of all relevant institutions. The current study included all patients (n � 369), but the longitudinal analyses were limited only to those patients who received treat- ment (n � 301). Thus, the n � 68 waitlist patients were excluded from this set of analyses.

Participants

All patients met Diagnostic and Statistical Manual of Mental Disorders (4th ed., text rev.; DSM–IV–TR; American Psychiatric Association, 2000) diagnostic criteria for PD with AG, scored �18 on the Hamilton Anxiety Scale (HAM-A; Shear et al., 2001), and scored �4 on the Clinical Global Impression scale (CGI; Guy, 1976). Other current comorbid diagnoses, including unipolar de- pression and other anxiety disorders, were allowed unless they were of primary clinical concern. Over 90% of the sample had at least one comorbid condition, with nearly half the sample diag- nosed with two or more mental disorders in addition to PD/AG. The most frequent comorbid conditions were specific phobia (n � 214; 71.1%), harmful use/abuse of alcohol (n � 135; 44.9%), social phobia (n � 126; 42.4%), and major depression (n � 118; 39.2%). As such, this sample can be considered both relatively severe and representative of patients seen in clinical practice. All patients were free from psychopharmacological medication. Ex- tensive details about inclusion and exclusion criteria have been previously published (Gloster et al., 2009; Gloster, Wittchen, et al., 2011).

The 301 patients in this study had a mean age of 35.5 years (SD � 10.7). A majority of patients were women (n � 228; 75.8%), and 131 (43.5%) had at least some higher education. Nearly one third were married (n � 98; 32.7%), half were single (n � 165; 55.0%), and the rest were divorced or widowed. Consistent with the demo- graphic characteristics of the population from which these data were sampled, all participants were of Caucasian origin.

Treatment

Patients received a 12-session manualized treatment protocol (Lang, Helbig-Lang, Westphal, Gloster, & Wittchen, 2011), im- plemented over 6 weeks, and followed by two booster sessions. Sixty-three certified therapists, all of whom were either advanced graduate students or post-docs, administered treatment. All thera-

pists went through a thorough training and certification procedure. Treatment integrity, training, randomization, and further design issues are published elsewhere (Gloster, Wittchen, et al., 2011).

The treatment was highly efficacious (Gloster, Wittchen, et al., 2011) and consisted of three phases: (1) psychoeducation, individ- ualized behavioral analysis, rationale for exposure, interoceptive exposure exercises (Sessions 1–4); (2) standardized in-situ expo- sure exercises, anticipatory anxiety, individualized in situ exposure exercises (Sessions 5–12); and (3) the generalization period through the 6-month follow-up assessment. The study had two active treatment groups that varied only with respect to the imple- mentation of a single component (in situ exposure with [T�] vs. without [T�] the therapist present), but not content. No relaxation exercises, breathing retraining, or explicit logical empiricism and disconfirmation of thoughts were undertaken in either group.

Assessment

Measures were assessed pre-treatment, at the intermediate point in treatment (between the fourth and fifth sessions), post-treatment (after the 12th session), and at the 6-month follow-up, which occurred 6 months after the post-treatment assessment. The one exception is the Anxiety Sensitivity Index, which was not mea- sured at the intermediate assessment.

PD/AG symptomatology and clinical functioning. Panic and Agoraphobia Scale (PAS; Bandelow, 1995). The

PAS is a patient self-report, 13-item questionnaire that measures the severity of panic attacks, avoidance, anticipatory anxiety, dis- ability, and worries about health. All items are scored from 0 to 4. Scores on the PAS have good reliability and are sensitive to change (Bandelow, 1995; Gloster, Wittchen, et al., 2011). The internal consistency of the PAS in this sample was � � .86.

CGI–Severity Subscale–Functioning Item (Guy, 1976). The CGI is a clinician-rated scale that measures the overall severity of a disorder, with scores that range between 1 (no disorder) and 7 (among the most severely ill patients). The scale normally queries for information across the facets of panic symptoms, anxiety, anticipatory anxiety, avoidance, and overall functional level before making the global rating. Scores on the CGI are sensitive to change in panic treatment (Barlow, Gorman, Shear, & Woods, 2000; Gloster, Wittchen, et al., 2011). For this study, we only used the one item measuring overall functioning in order to maximize conceptual distinctness from the PAS.

Process variables. Agoraphobic Cognitions Questionnaire (ACQ; Chambless et

al., 1984). The ACQ is a 14-item self-report questionnaire that measures the frequency of catastrophic beliefs about the possible consequences of experienced anxiety and panic. Each item is rated on a 5-point scale ranging from 1 (never) to 5 (always). The ACQ has sound psychometrics and is a standard assessment in PD/AG research (Zgourides, Warren, & Englert, 1989). The internal con- sistency of the ACQ in this sample was � � .74.

Bodily Sensations Questionnaire (BSQ; Chambless et al., 1984). The BSQ is a 17 item self-report questionnaire that mea- sures the degree of anxiety elicited by body sensations. Each item is rated on a 5-point scale ranging from 1 (not at all) to 5 (extremely). The BSQ has sound psychometrics and is a standard assessment in PD/AG research (Zgourides et al., 1989). The inter- nal consistency of the BSQ in this sample was � � .87.

T hi

s do

cu m

en t

is co

py ri

gh te

d by

th e

A m

er ic

an Ps

yc ho

lo gi

ca l

A ss

oc ia

tio n

or on

e of

its al

lie d

pu bl

is he

rs .

T hi

s ar

tic le

is in

te nd

ed so

le ly

fo r

th e

pe rs

on al

us e

of th

e in

di vi

du al

us er

an d

is no

t to

be di

ss em

in at

ed br

oa dl

y.

144 GLOSTER ET AL.

Mobility Inventory (MI; Chambless et al., 1985). The MI is a self-report questionnaire that measures the degree to which 27 situations are avoided. Items are scored from 1 (never avoid the situation) to 5 (always avoid the situation), with the mean of all items as the total score. Scores of the MI are highly reliable and sensitive to change (Chambless et al., 1985; Gloster, Wittchen, et al., 2011). For this study, only the ratings for the “alone” subscale are utilized. The internal consistency of the MI in this sample was � � .93.

Anxiety Sensitivity Inventory (ASI; Peterson & Reiss, 1993). The ASI is a 16-item self-report questionnaire that measures beliefs about potential harmful consequences of anxiety related symptoms. Each item is rated on a 5-point scale ranging from 0 (very little) to 4 (very much). The ASI has demonstrated sound psychometrics and is associated with various indices of PD/AG and other anxiety disorders (Rodriguez, Bruce, Pagano, Spencer, & Keller, 2004). The internal consistency of the ASI in this sample was � � .86.

Acceptance and Action Questionnaire–II (AAQ-II; Bond et al., 2011). The AAQ-II is a 7-item self-report questionnaire that measures psychological flexibility. Each item is rated on a 7-point scale ranging from 1 (never true) to 7 (always true). The AAQ-II has demonstrated sound psychometrics and is associated with various indices of PD/AG and other anxiety disorders, with good discriminant validity (Bond et al., 2011; Gloster, Klotsche, et al., 2011). The internal consistency of the AAQ-II in this sample was � � .94.

Statistical Analyses

Lower level mediation analyses were conducted to investigate the association of the five process variables ACQ, BSQ, MI, ASI, and AAQ-II and the change in treatment outcomes PAS and CGI over time (Kenny, Korchmaros, & Bolger, 2003) as a preliminary step in data analyses. The associations were estimated by multi- level linear mixed models with time as predictor variable, the five process variables ACQ, BSQ, MI, ASI, and AAQ-II as time varying covariates, and PAS and CGI as outcomes (Kenny et al., 2003; Singer & Willet, 2003).

Latent difference score (LDS) models provide a tool where change and individual differences in change are represented in the model (Selig & Preacher, 2009). LDS incorporates features of latent growth curve modeling and cross-lagged regression models. We only shortly describe our analytic strategy, a detailed presen- tation of the theory can be found in McArdle and Nesselroade (1994) or Hawley, Ho, Zuroff, and Blatt (2006). We evaluated different univariate LDS models for the change of PAS, CGI, ACQ, BSQ, ASI, AAG-II, and MI over time for investigating the nature of change in a first step. The latent change in a repeatedly observed score Y in an individual n at time t can be expressed by

�y�t�n � y�t�n � y�t � 1�n � �ysyn � �yy�t � 1�n, (1) where the observed score Y(t)n can be decomposed into a true score y(t)n and a measurement error en with a mean of zero and a positive variance. The latent change in Y is the sum of two components in Equation 1, an additive (�n syn) and a proportional (�y y(t � 1)n) change component. The coefficient syn corresponds to an intercept in the equation, which may vary across individuals and is constant over time. The � coefficient is a factor loading and

fixed to one for model identification purposes. The coefficient �y represents the proportional effect of the previous latent variable on the change rate. We compared univariate LDS models for time- invariant and time-varying proportional coefficients �y as well as a no change score LDS model (�y(t)n � 0, �y � �y � 0 in Equation 1) for each considered score. The univariate LDS models where combined to establish bivariate LDS in a second step. Bivariate LDS models provide an appealing feature for investigat- ing whether one score is the leading indicator of change in the other variable. A coupling parameter � is included into the equa- tions of two univariate LDS models representing the effect of one score on the rate of change in the other. The bivariate LDS model with another score z(t) at time t can be written by

�y�t�n � y�t�n � y�t � 1�n � �ysyn � �yy�t � 1�n � �zz�t � 1�n �z�t�n � z�t�n � z�t � 1�n � �zszn � �zz�t � 1�n � �yy�t � 1�n.

(2)

The relationship between the two dual change LDS models is given by the components �z z(t � 1)n and �y y(t � 1)n besides the additive and proportional change components. The subsequent latent change in one variable is predicted by the other variable occurring earlier in time in case of coupling between the two univariate LDS models. We investigated different patterns of cou- pling between two univariate LDS models by restricting the path coefficients in the models. The analyses included models with (a) no coupling (�z � 0 and �y � 0) between the two series, (b) unidirectional coupling exists in which one variables predicts later change in the other and vice versa (�z � 0 and �y � 0 or �z � 0 and �y � 0), and (c) bidirectional coupling exists between the two scores (�z � 0 and �y � 0). We also compared models with time-invariant and time-varying coupling coefficients �z and �y. Whenever the final model indicated that more than one � coeffi- cient (one per phase of treatment) per process variable was signif- icant, the coefficients were tested for significant differences. The third step of our analyses concerns the hypotheses whether treat- ment condition (T� vs. T�) predicts the subsequent rate of change in the studied variables over the treatment process. Treatment condition is added by the term ( TX) in Equation 2. All path coefficients are reported as unstandardized coefficients. The pa- rameters of the LDS models were estimated in Mplus Version 6.1 (Muthén & Muthén, 2011). We used the full the full-information maximum likelihood estimator due to missing data in some cases. This approach ensures the use of all available data for parameter estimation. Thus, also patients with incomplete sessions were incorporated into analyses.

Results

Baseline Values

The mean and standard deviations for the outcome variables and process variables at baseline and post-treatment are displayed in Table 1. The correlations between variables at baseline and post- treatment are likewise displayed in Table 1.

Lower Level Mediation Models

The association of the five process variables ACQ, BSQ, ASI, MI, AAQ-II, and the treatment outcomes PAS and CGI were

T hi

s do

cu m

en t

is co

py ri

gh te

d by

th e

A m

er ic

an Ps

yc ho

lo gi

ca l

A ss

oc ia

tio n

or on

e of

its al

lie d

pu bl

is he

rs .

T hi

s ar

tic le

is in

te nd

ed so

le ly

fo r

th e

pe rs

on al

us e

of th

e in

di vi

du al

us er

an d

is no

t to

be di

ss em

in at

ed br

oa dl

y.

145TIMING MATTERS

investigated by lower level mediation models. This preliminary step was conducted in view of the existing literature. ACQ, BSQ, MI, ASI, and AAQ-II all partially mediated the treatment out- comes of PAS and CGI as indicated by a significant mediated effect in the mediator analyses (available upon request). However, lower level mediation models are inadequate to show sequencing across time.

The waitlist (WL) group reported only negligible pre-treatment to post-treatment changes and was significantly worse than both treatment groups at post-treatment (see Gloster, Wittchen, et al., 2011). Nevertheless, the WL group was tested here using lower level models. As expected, the WL group did not demonstrate any meditational effects. As no meaningful change was observed in this group, predicting change was not possible, and this group was excluded from further longitudinal analyses below.

Univariate Latent Difference Score Models

The change in PAS, CGI, ACQ, BSQ, MI, ASI, and AAQ-II was investigated by univariate LDS models including the no change model and the two dual change models with both time- varying and time-invariant proportional effects �(t). The no change LDS models consistently resulted in a poor model fit (standardized root-mean-square residual [SRMR] ranges from .26 for ASI to .42 for AAQ-II). The univariate LDS models including time-varying proportional effects �(t) substantially improved model fit compared to time-invariant proportional effects in all analyzed models. The model fit of the univariate LDS models can be considered to be acceptable for modeling the change in the seven variables over time by a dual change model with time- varying proportional effects (comparative fit index [CFI] ranges from .86 for CGI to .98 for ASI; Tucker–Lewis index [TLI] ranges from .84 for MI to .98 for AAQ-II; root-mean-square error of

approximation [RMSEA] ranges from .08 for AAQ-II and MI to .14 for CGI; SRMR ranges from .04 for ASI to .10 for MI). Unstandardized parameter estimates for the proportional effects �(t) were statistically significant (ps ranging from .001 to .043) except for the BSQ and MI. The latent BSQ (�1 � �.12, p � .35) and MI (�1 � �.15, p � .10) at baseline assessment did not significantly predict the subsequent rate of change. Detailed infor- mation is reported in Table 2 about model fit and parameter estimates for the additive and proportional change components.

Bivariate Latent Difference Score Models

The parameter estimates and the model fit indices are reported in Table 3 and Table 4 for the final bivariate LDS models. The final models were selected based on considering a combination of the Bayesian information criterion (BIC), CFI, TLI, RMSEA, and SRMR.

Panic and agoraphobia symptoms. Five bivariate LDS mod- els were conducted for evaluating the coupling between the uni- variate series of PAS and ACQ, BSQ, ASI, MI, and AAQ-II as reported in Table 3. We compared four models for each variable combination: (a) no coupling exists, (b) unidirectional coupling from PAS to the process variable, (c) unidirectional coupling from the process variable to PAS, and (d) bidirectional coupling be- tween PAS and process variable. Given our results, the bivariate LDS models including unidirectional coupling from ACQ to PAS (SRMR � .05), BSQ to PAS (SRMR � .06), and ASI to PAS (SRMR � .03) resulted in best model fit. Latent ACQ significantly predicts later change in PAS for baseline to intermediate (�1,ACQ_PAS), intermediate to post (�2,ACQ_PAS), and post to follow-up (�3,ACQ_PAS) assessments. The coupling coefficient �2,ACQ_PAS significantly differs from �3,ACQ_PAS, �

2(1) � 12.78, p .001, indicating that the strongest association existed for latent

Table 1 Distribution Outcome and Process Variables at Baseline and Post-Assessment and Correlations Between Variables

Variables M (SD) PAS CGI ACQ BSQ MI ASI AAQ-II

Correlations between measures at baseline Outcome variables

PAS 27.8 (9.8) 1 CGI 4.54 (0.90) .44 1

Process variables ACQ 2.18 (0.57) .37 .15 1 BSQ 48.0 (12.4) .35 .12 .59 1 MI 2.98 (0.81) .55 .40 .17 .26 1 ASI 31.4 (11.5) .40 .12 .53 .53 .22 1 AAQ-II 45.3 (10.2) .19 .16 .36 .28 .14 .50 1

Correlations between measures at post Outcome variables

PAS 14.4 (9.3) 1 CGI 3.04 (1.25) .64 1

Process variables ACQ 1.63 (0.46) .58 .32 1 BSQ 34.5 (11.9) .57 .36 .63 1 MI 1.96 (0.87) .56 .49 .43 .36 1 ASI 16.6 (10.8) .65 .41 .63 .71 .37 1 AAQ-II 52.3 (10.0) .48 .32 .48 .43 .28 .57 1

Note. All correlations are significant at the 5% level. PAS � Panic and Agoraphobia Scale; CGI � Clinical Global Impression scale; ACQ � Agoraphobic Cognitions Questionnaire; BSQ � Bodily Sensations Questionnaire; MI � Mobility Inventory; ASI � Anxiety Sensitivity Inventory; AAQ-II � Acceptance and Action Questionnaire–II.

T hi

s do

cu m

en t

is co

py ri

gh te

d by

th e

A m

er ic

an Ps

yc ho

lo gi

ca l

A ss

oc ia

tio n

or on

e of

its al

lie d

pu bl

is he

rs .

T hi

s ar

tic le

is in

te nd

ed so

le ly

fo r

th e

pe rs

on al

us e

of th

e in

di vi

du al

us er

an d

is no

t to

be di

ss em

in at

ed br

oa dl

y.

146 GLOSTER ET AL.

ACQ predicts change in PAS at the interval post to follow-up assessment. Latent BSQ did not predict later change in PAS over time. Latent ASI at post assessment predicted later change in PAS (�3,ASI_PAS � .28, p � .02). The longitudinal association of PAS and MI was best modeled by a bivariate LDS model with time- invariant coupling coefficients (SRMR � .06). Latent MI pre- dicted later change in PAS (�1,MI_PAS � �2,MI_PAS � �3,MI_PAS � 11.0, p .001), suggesting the effect of MI was not different across treatment phases. Notably, latent PAS also predicted later change in MI (�1,PAS_MI � �2,PAS_MI � �3,PAS_MI � .13, p .001). The bivariate LDS model with bidirectional coupling and time- varying coupling coefficients for PAS and AAQ-II achieved best model fit (SRMR � .04). Latent AAQ-II predicted later change in PAS in the interval intermediate assessment and post assessment (�3,AAQ-II_PAS � .13, p � .03).

Clinical functioning. Bivariate LDS models were applied for investigating the associations of CGI and the five process variables over time as shown in Table 4. The final models included bidi- rectional coupling coefficients. The alternative bivariate LDS models for CGI and ACQ including no coupling, unidirectional coupling, and bidirectional coupling resulted in an acceptable (RMSEA ranges from .12 to .13, SRMR ranges from .17 to .18, CFI is .9, and TLI ranges from .81 to .85). The LDS model for CGI and ACQ with the closest fit included coupling coefficients that were constraint to be equal over time. Neither latent ACQ nor CGI predicted later latent change in the other variable. Latent BSQ at post assessment significantly predicted later change in CGI (�3,BSQ_CGI � .20, p � .03). In contrast, change in CGI did not predict later change in BSQ. The best fitting bivariate LDS model for ASI and CGI (SRMR � .04) included bidirectional time- varying coupling coefficients for ASI predicting later change in CGI and time-invariant coupling for CGI predicting later change in ASI. Latent ASI at post assessment predicted later change in CGI (�2,ASI_CGI � .06, p .001). Notably, change in CGI significantly predicted later change in ASI over time (�1,CGI_ASI � �2,CGI_ASI � �3,CGI_ASI � 5.79, p .001). The longitudinal association of CGI and MI was best modeled by a bivariate LDS model with time- invariant coupling coefficients (SRMR � .18), suggesting the effect of MI is not significantly different across treatment phases. It is notable that even the best fitting model did not result in a consistent good model fit. The indices CFI (CFI � .90) and TLI (TLI � .85) suggested an acceptable model fit, whereas RMSEA (RMSEA � .14) and SRMR (SRMR � .18) suggested a poor model fit. Latent MI predicted later change in CGI (�1,MI_CGI � �2,MI_CGI � �3,MI_CGI � 2.58, p .001) and latent CGI also predicted later change in MI (�1,CGI_MI � �2,CGI_MI � �3,CGI_MI � .96, p .001). The bivariate LDS model for AAQ-II and CGI including time-varying coupling coefficients resulted in an acceptable model fit (SRMR � .09). Latent AAQ-II predicted later change in CGI over time (�1,AAQ-II_CGI � .40, p � .008; �2,AAQ-II_CGI � .43, p � .01; �3,AAQ-II_CGI � .47, p .014) and vice versa (�1,CGI_AAQ-II � 9.4, p � .004; �2,CGI_AAQ-II � 8.5, p � .011; �3,CGI_AAQ-II � 6.1, p � .028). Although the three coupling coefficients for latent AAQ-II predicting later change in CGI differed, these differences were not significantly different throughout treatment.

Treatment condition. We tested whether the treatment con- dition in our study (T� vs. T�) predicted the rate of change in outcome and process variables. We added treatment condition to the best fitting bivariate LDS model as presented in Table 3. ForT

ab le

2 U

ni va

ri at

e L

D S

M od

el s

fo r

P A

S, C

G I,

A C

Q ,

B SQ

, A

SI ,

M I,

an d

A A

Q -I

I

Pa ra

m et

er s

an d

fi t

in di

ce s

PA S

C G

I A

C Q

B SQ

A SI

M I

A A

Q -I

I

A dd

iti ve

co ef

fi ci

en ts

E [s

] (S

E );

p va

lu e

27 .8

(0 .6

); p

.0

01 5.

2 (0

.1 );

p

.0 01

2. 2

(0 .0

3) ;p

.0

01 48

.0 (0

.7 );

p

.0 01

31 .3

(0 .7

); p

.0

01 3.

0 (0

.0 5)

;p

.0 01

45 .3

(0 .6

); p

.0

01 �

2 (s

) 67

.8 5

0. 32

0. 25

10 2.

10 12

2. 48

0. 65

77 .7

9 Pr

op or

tio na

l co

ef fi

ci en

ts �

1 (S

E );

p va

lu e

� .2

2 (.

09 );

p �

.0 21

� .4

6 (.

15 );

p �

.0 02

� .3

8 (.

11 );

p �

.0 01

� .1

2 (.

13 );

p �

.3 48

— †

� .1

5 (.

09 );

p �

.1 02

� .2

5 (.

12 );

p �

.0 37

� 2

(S E

); p

va lu

e �

.4 3

(. 11

); p

.0

01 �

.3 3

(. 16

); p

� .0

43 �

.5 6

(. 12

); p

.0

01 �

.3 4

(. 13

); p

� .0

09 �

.8 2

(. 06

); p

.0

01 a

� .3

7 (.

10 );

p

.0 01

� .3

9 (.

16 );

p �

.0 15

� 3

(S E

); p

va lu

e �

.5 0

(. 17

); p

� .0

04 �

.5 0

(. 21

); p

� .0

17 �

.5 7

(. 15

); p

.0

01 �

.2 0

(. 15

); p

� .0

24 �

.8 2

(. 11

); p

.0

01 �

.3 9

(. 13

); p

� .0

03 �

.2 4

(. 11

); p

� .0

29 G

oo dn

es s

of fi

t pa

ra m

et er

s N

o. 9

9 9

9 8

9 9

B IC

7, 45

9 3,

06 6

1, 10

5 7,

82 7

5, 67

0 1,

88 0

7, 37

1 �

2 (d

f) ;

p va

lu e

21 .8

(5 );

p �

.0 01

70 .8

(5 );

p �

.0 02

20 .3

(5 );

p �

.0 01

22 .4

(5 );

p

.0 01

18 5.

3 (3

); p

.0

01 86

.0 (5

); p

.0

01 14

.6 (5

); p

� .0

1 C

FI .9

5 .8

6 .9

7 .9

6 .9

8 .8

7 .9

6 T

L I

.9 5

.8 6

.9 6

.9 5

.9 4

.8 4

.9 5

R M

SE A

.1 1

.1 4

.1 0

.1 1

.1 1

.0 8

.0 8

SR M

R .0

6 .0

9 .0

5 .0

8 .0

4 .1

0 .0

9

N ot

e. L

D S

� la

te nt

di ff

er en

ce sc

or e;

PA S

� Pa

ni c

an d

A go

ra ph

ob ia

Sc al

e; C

G I

� C

lin ic

al G

lo ba

l Im

pr es

si on

sc al

e; A

C Q

� A

go ra

ph ob

ic C

og ni

tio ns

Q ue

st io

nn ai

re ;

B SQ

� B

od ily

Se ns

at io

ns Q

ue st

io nn

ai re

; A

SI �

A nx

ie ty

Se ns

iti vi

ty In

ve nt

or y;

M I

� M

ob ili

ty In

ve nt

or y;

A A

Q -I

I �

A cc

ep ta

nc e

an d

A ct

io n

Q ue

st io

nn ai

re –I

I; �

1 ,

� 2 ,

an d

� 3

� di

st in

ct tim

e- va

ry in

g pr

op or

tio na

l ch

an ge

co ef

fi ci

en ts

;N o.

� nu

m be

r of

m od

el pa

ra m

et er

s; B

IC �

B ay

es ia

n in

fo rm

at io

n cr

ite ri

on ;C

FI �

co m

pa ra

tiv e

fi ti

nd ex

;T L

I �

T uc

ke r–

L ew

is in

de x;

R M

SE A

� ro

ot -m

ea n-

sq ua

re er

ro r

of ap

pr ox

im at

io n;

SR M

R �

st an

da rd

iz ed

ro ot

-m ea

n- sq

ua re

re si

du al

; —

† �

th e

pa ra

m et

er is

no t

es tim

at ed

du e

to m

is si

ng A

SI at

in te

rm ed

ia te

as se

ss m

en t.

a Pr

op or

tio na

l ch

an ge

co ef

fi ci

en t

� 2

re fe

rs to

th e

in te

rv al

ba se

lin e

to po

st -a

ss es

sm en

t.

T hi

s do

cu m

en t

is co

py ri

gh te

d by

th e

A m

er ic

an Ps

yc ho

lo gi

ca l

A ss

oc ia

tio n

or on

e of

its al

lie d

pu bl

is he

rs .

T hi

s ar

tic le

is in

te nd

ed so

le ly

fo r

th e

pe rs

on al

us e

of th

e in

di vi

du al

us er

an d

is no

t to

be di

ss em

in at

ed br

oa dl

y.

147TIMING MATTERS

example, treatment condition was established as a predictor for the rate of change in ACQ in the bivariate LDS model for ACQ and PAS with unidirectional coupling (latent ACQ predicts later change in PAS). The goodness-of-fit parameters indicate a good model fit, �2(23) � 41.2, p � .01; CFI � .98; TLI � .97; RMSEA � .05; SRMR � .05. However, treatment condition did not significantly predict change in ACQ throughout treatment ( 1 � �0.02, p � .88; 2 � 0.21, p � .79; 3 � 1.60, p � .07). A similar pattern was found for the other bivariate LDS models with treatment condition as an additional explanatory variable.

Discussion

This study examined the mechanism of action in CBT for panic disorder with agoraphobia using temporally sensitive bivariate latent difference score modeling in a large sample of 301 patients. Importantly, the current study found evidence for the temporal specificity of process-outcome effects over the course of therapy and differences between outcome measures. Change in a compre- hensive measure of panic and agoraphobic symptoms (PAS) was consistently predicted across all phases of treatment by previous values of catastrophic appraisal (ACQ) and agoraphobic avoidance (MI). In contrast to the unidirectional relationship from ACQ to

subsequent scores on the PAS, the relation between MI and PAS was bidirectional (i.e., scores on the PAS also predicted subsequent change on the MI). During the second phase of treatment (i.e., exposure in situ), scores on psychological flex- ibility (AAQ-II) predicted subsequent change in the PAS at post-treatment. This suggests that psychological flexibility is particularly relevant during the phase of treatment that patients are asked to face their fears. Further changes on the PAS during the 6-month follow-up period were unidirectionally associated with scores on the ASI at post-treatment. Fear of bodily symp- toms (BSQ) did not predict subsequent change in the PAS during any stage of therapy.

Change in global functioning (CGI) presented a somewhat dif- ferent picture. Scores in both agoraphobic avoidance (MI) and psychological flexibility (AAQ-II) predicted subsequent change in functioning across all phases of treatment. Likewise, scores on the CGI predicted subsequent changes in the MI and AAQ-II during these phases (bidirectional relations). This suggests that both avoidance and psychological flexibility are strongly related to functioning across the therapy and follow-up periods and are complexly intertwined with functioning. In addition, scores on anxiety sensitivity (ASI) at baseline predicted subsequent change

Table 3 Bivariate LDS Models for Outcome PAS and Process Variables ACQ, BSQ, ASI, MI, and AAQ-II

Parameters and fit indices ACQ BSQ ASI MI AAQ-II

Additive coefficients E[�o] (SE); p value 27.81 (0.56); p .001 27.81 (0.57); p .001 27.78 (0.57); p .001 27.88 (0.55); p .001 27.77 (0.57); p .001 �2(�o) 67.10 69.76 54.44 55.47 65.99 E[�p] (SE); p value 2.18 (0.03); p .001 48.00 (0.71); p .001 31.35 (0.67); p .001 2.98 (0.05); p .001 45.34 (0.59); p .001 �2(�p) 0.24 103.00 98.22 0.59 73.96

Proportional coefficients �1,PAS (SE); p value �.44 (.11); p .001 �.85 (.11); p .001 —

† �1.24 (.21); p .001 .02 (.22); p � .942 �2,PAS (SE); p value �.80 (.15); p .001 �1.23 (.12); p .001 �1.01 (.09); p .001 �1.54 (.23); p .001 �.12 (.25); p � .644 �3,PAS (SE); p value �1.20 (.22); p .001 �1.80 (.23); p .001 �1.40 (.17); p .001 �1.67 (.26); p .001 �.32 (.27); p � .227

�1,pv (SE); p value �.38 (.11); p � .001 �.18 (.10); p � .069 — † �.79 (.24); p � .001 .03 (.22); p � .891

�2,pv (SE); p value �.55 (.12); p .001 �.41 (.11); p .001 �.82 (.06); p .001 a �.43 (.23); p � .067 .04 (.22); p � .865

�3,pv (SE); p value �.56 (.15); p .001 �.28 (.14); p � .044 �.28 (.11); p .001 .20 (.25); p � .413 �.11 (.23); p � .637 Cross-lag coefficients

�1,PAS_pv (SE); p value — †† —†† —† .13 (.02); p .001 �.27 (.31); p � .386

�2,PAS_pv (SE); p value — †† —†† —†† .13 (.02); p .001 �.04 (.33); p � .896

�3,PAS_pv (SE); p value — †† —†† —†† .13 (.02); p .001 .08 (.36); p � .836

�1,pv_PAS (SE); p value 15.20 (4.68); p � .001 .46 (.38); p � .665 — † 11.02 (2.31); p .001 .10 (.22); p � .657

�2,pv_PAS (SE); p value 17.18 (5.43); p � .002 .53 (.34); p � .706 .17 (.09); p � .092 a 11.02 (2.31); p .001 .13 (.27); p � .031

�3,pv_PAS (SE); p value 22.94 (7.06); p � .001 .35 (.36); p � .685 .28 (.14); p � .022 11.02 (2.31); p .001 .12 (.22); p � .583 Goodness of fit parameters

No. 26 26 24 25 29 BIC 8,339 15,066 10,919 9,082 14,753 �2 (df); p value 38.1 (18); p .001 35.8 (18); p .001 7.04 (4); p .134 79.5 (19); p .001 3.3 (15); p .001 CFI .98 .98 1.00 .95 .98 TLI .97 .97 .98 .93 .97 RMSEA .06 .06 .05 .10 .06 SRMR .05 .06 .03 .06 .04

Note. LDS � latent difference score; PAS � Panic and Agoraphobia Scale; ACQ � Agoraphobic Cognitions Questionnaire; BSQ � Bodily Sensations Questionnaire; ASI � Anxiety Sensitivity Inventory; MI � Mobility Inventory; AAQ-II � Acceptance and Action Questionnaire–II; �1, �2, and �3 � distinct time-varying proportional change coefficients; �1,PAS_pv, �2,PAS_pv, �3,PAS_pv � distinct coupling coefficients for latent PAS predicting later change in process variable; �1,pv_PAS, �2,pv_PAS, �3,pv_PAS � distinct coupling coefficients for latent process variable predicting later change in PAS; No. � number of model parameters; BIC � Bayesian information criterion; CFI � comparative fit index; TLI � Tucker–Lewis index; RMSEA � root-mean-square error of approximation; SRMR � standardized root-mean-square residual; —† � the parameter is not estimated due to missing ASI at intermediate assessment; —†† � the parameter was not estimated because during the process of model building better model fit was obtained by excluding the parameter. a Proportional change coefficient �2 refers to the interval baseline to post-assessment.

T hi

s do

cu m

en t

is co

py ri

gh te

d by

th e

A m

er ic

an Ps

yc ho

lo gi

ca l

A ss

oc ia

tio n

or on

e of

its al

lie d

pu bl

is he

rs .

T hi

s ar

tic le

is in

te nd

ed so

le ly

fo r

th e

pe rs

on al

us e

of th

e in

di vi

du al

us er

an d

is no

t to

be di

ss em

in at

ed br

oa dl

y.

148 GLOSTER ET AL.

in functioning from pre-treatment to post-treatment, as did func- tioning predict subsequent change in anxiety sensitivity (bi- directional relation). These bi-directional relationships may be similar to those observed by Teachman et al. (2010) with the variables distress/apprehension. The only unidirectional relation with functioning was observed for scores on the fear of bodily symptoms (BSQ) at post-treatment for subsequent change in func- tioning between post-treatment and follow-up. Catastrophic ap- praisal (ACQ) did not predict subsequent change in functioning at any point during the study.

Taken together, our hypotheses were partially supported. In partial contrast to our hypothesis, cognitive appraisal predicted subsequent change in panic and agoraphobia symptoms across all time points and not just during Sessions 1–4. Consistent with our hypothesis, cognitive appraisal did not predict subsequent change in functioning. Contrary to our hypothesis, anxiety elicited by bodily symptoms did not predict subsequent change in panic symptoms at any time point and did predict subsequent functioning at the follow up assessment. Contrary to our hypothesis, anxiety sensitivity predicted subsequent panic and agoraphobia symptoms only during the follow-up period. However, consistent with our hypothesis, anxiety sensitivity was related to functioning during the treatment phase. Consistent with our hypotheses, avoidance

behavior was related to subsequent change in panic symptoms and functioning across all phases. Also consistent with our hypothesis, psychological flexibility predicted subsequent change in panic symptoms during Sessions 5–12 and functioning across all time points.

This research builds on previous studies, all of which used panic and/or agoraphobic symptoms as an outcome variable. Although some of these studies largely lacked prospective temporal designs that measured process variables and outcome measures longitudi- nally (for exceptions, see Meuret et al., 2010; Teachman et al., 2010), results from these studies help piece together the puzzle of the processes relevant for effective treatment. Indeed, our results are consistent with the reliable finding that measures of one’s appraisal of symptoms (e.g., ACQ, BSQ, and ASI) mediated or partially mediated outcome (Casey et al., 2005; Hofmann et al., 2007; Meulenbeek et al., 2010; Meuret et al., 2010; Smits et al., 2004; Vögele et al., 2010). This finding was also found in studies that used different analytical frameworks such as time series analysis (Bouchard et al., 2007) structural equation modeling (Hoffart, Sexton, Hedley, & Martinsen, 2008), and bivariate latent difference score modeling (Teachman et al., 2010).

The longitudinal analysis used in the current study advance our understanding of when and to some degree how the consistent

Table 4 Bivariate LDS Models for Outcome CGI and Process Variables ACQ, BSQ, ASI, MI, and AAQ-II

Parameters and fit indices ACQ BSQ ASI MI AAQ-II

Additive coefficients E[�o] (SE); p value 5.17 (0.05); p .001 5.18 (0.05); p .001 5.17 (0.05); p .001 5.18 (0.05); p .001 5.17 (0.05); p .001 �2(�o) 0.32 0.02 1.15 0.26 0.31 E[�p] (SE); p value 2.18 (0.03); p .001 47.98 (0.70); p .001 31.37 (0.66); p .001 2.99 (0.05); p .001 45.33 (0.59); p .001 �2(�p) 0.25 94.69 55.3 0.56 74.05

Proportional coefficients �1,CGI (SE); p value �.50 (.23); p � .025 �2.84 (.94); p � .003 —

† �1.47 (.37); p .001 �1.98 (.65); p � .002 �2,CGI (SE); p value �.37 (.23); p � .109 �3.39 (1.14); p � .003 �1.36 (.16); p � .001 �1.63 (.37); p .001 �2.14 (.68); p � .002 �3,CGI (SE); p value �.53 (.26); p � .037 �3.34 (1.16); p � .004 �1.52 (.17); p � .001 �1.52 (.40); p .001 �1.44 (.54); p � .007

�1,pv (SE); p value �.43 (.15); p � .004 1.43 (0.93); p � .123 — † �1.03 (.35); p � .003 �.48 (.34); p � .155

�2,pv (SE); p value �.59 (.15); p .001 1.53 (1.04); p � .141 .00 (.13); p � .982 a �.80 (.36); p � .025 �.50 (.25); p � .047

�3,pv (SE); p value �.56 (.17); p � .001 1.73 (1.16); p � .136 �.12 (.10); p � .252 �.70 (.37); p � .058 �.89 (.40); p � .027 Cross-lag coefficients

�1,CGI_pv (SE); p value .07 (.07); p � .308 �1.58 (6.17); p � .086 — † .96 (.24); p .001 9.36 (3.27); p � .004

�2,CGI_pv (SE); p value .07 (.07); p � .308 14.01 (7.50); p � .062 5.79 (1.63); p .001 .96 (.24); p .001 8.52 (3.35); p � .011 �3,CGI_pv (SE); p value .07 (.07); p � .308 12.86 (7.63); p � .092 5.79 (1.63); p .001 .96 (.24); p .001 6.13 (2.79); p � .028

�1,pv_CGI (SE); p value �.23 (.61); p � .710 .13 (.08); p � .109 — † 2.58 (0.49); p .001 .40 (.15); p � .008

�2,pv_CGI (SE); p value �.23 (.61); p � .710 .11 (.08); p � .191 .06 (.02); p � .001 a 2.58 (0.49); p .001 .43 (.17); p � .010

�3,pv_CGI (SE); p value �.23 (.61); p � .710 .20 (.09); p � .031 .05 (.05); p � .684 2.58 (0.49); p .001 .47 (.19); p � .014 Goodness of fit parameters

No. 25 29 23 26 26 BIC 4,125 10,817 7,967 4,761 10,385 �2 (df); p value 97.6 (19); p .001 45.3 (15); p .001 13.2 (4); p .011 117.6 (18); p .001 6.7 (18); p .001 CFI .90 .96 .98 .90 .96 TLI .85 .92 .91 .85 .93 RMSEA .12 .08 .09 .14 .08 SRMR .18 .05 .04 .18 .09

Note. LDS � latent difference score; CGI � Clinical Global Impression scale; ACQ � Agoraphobic Cognitions Questionnaire; BSQ � Bodily Sensations Questionnaire; ASI � Anxiety Sensitivity Inventory; MI � Mobility Inventory; AAQ-II � Acceptance and Action Questionnaire–II; �1, �2, and �3 � distinct time-varying proportional change coefficients; �1,CGI_pv, �2,CGI_pv, �3,CGI _pv � distinct coupling coefficients for latent Panic and Agoraphobia Scale (PAS) predicting later change in process variable; �1,pv_ CGI, �2,pv_ CGI, �3,pv_ CGI � distinct coupling coefficients for latent process variable predicting later change in PAS; No. � number of model parameters; BIC � Bayesian information criterion; CFI � comparative fit index; TLI � Tucker–Lewis index; RMSEA � root-mean-square error of approximation; SRMR � standardized root-mean-square residual; —† � the parameter is not estimated due to missing ASI at intermediate assessment. a Proportional change coefficient �2 refers to the interval baseline to post-assessment.

T hi

s do

cu m

en t

is co

py ri

gh te

d by

th e

A m

er ic

an Ps

yc ho

lo gi

ca l

A ss

oc ia

tio n

or on

e of

its al

lie d

pu bl

is he

rs .

T hi

s ar

tic le

is in

te nd

ed so

le ly

fo r

th e

pe rs

on al

us e

of th

e in

di vi

du al

us er

an d

is no

t to

be di

ss em

in at

ed br

oa dl

y.

149TIMING MATTERS

finding that one’s appraisals of panic and agoraphobic symptoms mediates outcome. Consistent with previous studies, results of the current study suggest that panic symptomatology is affected by one’s catastrophic beliefs (ACQ) during all phases of treatment. These analyses also suggest that this effect is strongest during the generalization phase. Similarly, agoraphobic avoidance is associ- ated with subsequent change in PAS across treatment. Interest- ingly, psychological flexibility seems to affect panic symptomatol- ogy during the in situ exposure phase of treatment. This would suggest that exposure in situ requires a patient to engage with the feared stimuli in a flexible manner and take steps to reduce avoidance behavior (see Gloster, Hummel, Lyudmirskaya, Hauke, & Sonntag, 2012). During the follow-up period, agoraphobic avoidance, catastrophic cognitions, and fear of fear are the salient process variables. In sum, whereas cognitive variables do affect panic and agoraphobic-related outcome, the present results suggest that not all cognitive variables predict outcome, and which cogni- tive variables are the most salient predictors depends on the phase of treatment. However, given their exploratory nature these find- ings clearly require replication before firm conclusions can be drawn.

We also found evidence for the process of change in two variables not previously tested longitudinally. First, the degree of self-reported situational agoraphobic avoidance (MI) was most consistently associated with the reduction in panic and agorapho- bic symptoms and functioning. The bidirectional relation suggests a complex relation between these variables, likely due in part to a partial overlap of the constructs. It is important to note that two previous mediation studies included agoraphobic avoidance in their analyses but treated it as a dependent variable (Meulenbeek et al., 2010; Vögele et al., 2010). We treated agoraphobic avoidance as a potential process variable, however, because the treatment conceptualized avoidance and safety behaviors as a maintaining variable and directly targeted them (both those publically observ- able and those only observable to the patient). As such, the mobility inventory captured one class of this behavior. The final process variable associated with changes in the outcome was psychological flexibility (AAQ-II). This variable is not specific to panic and agoraphobia and the AAQ-II does not contain any words specifically referring to panic or agoraphobia. Instead, it is a broader construct that measures the degree to which one can mindfully accept thoughts and emotions while engaging in one’s life when it is important to do so. As such, it is theoretically consistent that psychological flexibility was associated with change in panic and agoraphobia symptoms only during the phase of treatment that concentrated on exposure in situ but not the phase that concentrated on psychoeducation (Gloster et al., 2012): drop- ping subtle avoidance behaviors and mindfully accepting associ- ated thoughts and emotions promotes change. It is likewise theo- retically consistent that psychological flexibility was consistently related to subsequent change in functioning: promotion of psycho- logical flexibility increases one’s ability to engage with that which is important to the patient.

This study also expanded the examination of process variables on the outcome of global functioning. In addition to adding infor- mation about how the putative process variables affect a broader target, testing the process variables against the CGI also served as a test of specificity for the process variables. The variables asso- ciated with the change in global functioning were agoraphobic

avoidance and psychological flexibility across all treatment phases, anxiety sensitivity during the active phase of treatment, and fear of bodily symptoms during the generalization phase of treatment. Interestingly, the cognitive appraisal process variables that have consistently been found to be associated with change in panic and agoraphobic symptoms were no longer significantly related to global functioning in the longitudinal models and only anxiety sensitivity was related to change in functioning during the active treatment phase. This, therefore, partially supports and is simultaneously in partial contrast to Smits et al. (2004), who found that anxiety sensitivity statistically mediated functioning in a cross-sectional analysis. Differences may have resulted from the timing of measurements (concurrent measurement at pre and post vs. longitudinal), measurement format (clinician judgment in the present study vs. questionnaire in Smits et al., 2004), differences in the treatment, or a combination of these factors. Once again, there is a critical need for replications before the processes that lead to change in global functioning can be established. Results clearly point to crucial importance of testing across various definitions of outcome and especially of expanding beyond purely symptom- based definitions. If replicated, these results suggest that different processes are involved in the change of symptomatology and functioning across the various treatment phases.

Treatment group (T� vs. T�) did not contribute to the expla- nation of relation between processes and outcome and were not included in the final models. This suggests that despite the slight advantage seen by the T� group in outcome (Gloster, Wittchen, et al., 2011), both treatment variants seem to work through similar processes. This is not surprising as both treatment variants had identical content and differed only with respect to the therapist’s presence during exposure in situ. It remains a possibility that the presence of the therapist may have facilitated the dropping of safety behaviors or offered more intense guidance, but the sum total of such effects—if they do indeed exist—are not strong enough to be detected by these analyses.

By linking process-outcome effects with specific phases and elements of treatment, we are in a stronger position to tie together results from outcome trials with current theories about the mech- anisms that underlie treatment. For example, inhibitory learning that promotes tolerance of anxiety and develops competing non- threat expectancies and that can be generalized across contexts is believed to be a crucial mechanism in exposure therapy (Arch & Craske, 2008). This study, then, shows with temporal fidelity that some of the therapeutic techniques and processes are involved at different time points during the therapy and may point to a spec- ification of what is involved in the processes of inhibitory learning. That is, during the intensive exposure in situ phase of treatment, cognitive attribution, agoraphobic avoidance, and psychological flexibility are associated with changes in PD/AG severity, whereas changes in global functioning are associated with agoraphobic avoidance, psychological flexibility, and anxiety sensitivity. Al- though clearly in need of replication in other variations of CBT for PD/AG, this type of analysis aids in the understanding of treatment processes at specific level.

This study needs to be interpreted in the light of several limi- tations. First, although consistent with previous studies, the pro- cess variables examined in this study were assessed using ques- tionnaires are limited by the retrospective recall bias inherent in questionnaires. Future studies using additional methodologies

T hi

s do

cu m

en t

is co

py ri

gh te

d by

th e

A m

er ic

an Ps

yc ho

lo gi

ca l

A ss

oc ia

tio n

or on

e of

its al

lie d

pu bl

is he

rs .

T hi

s ar

tic le

is in

te nd

ed so

le ly

fo r

th e

pe rs

on al

us e

of th

e in

di vi

du al

us er

an d

is no

t to

be di

ss em

in at

ed br

oa dl

y.

150 GLOSTER ET AL.

(e.g., ecological momentary assessment, physiological variables; e.g., Domschke et al., 2011; Kircher et al., 2013; Richter et al., 2012), with different sources of method variance, are clearly needed. Second, although we broke down the effects of time across our treatment, the effects of time and the treatment components that occurred during that period of the treatment cannot be par- celed apart. Third, the ASI was not administered during the inter- mediate assessment. Further, the original ASI was utilized in this study. Subsequent versions of the ASI have expanded the measure and emphasized its multidimensional aspects (Taylor & Cox, 1998; Taylor et al., 2007). Although all versions of the ASI target the overarching concept of anxiety sensitivity, results from this study do not inform about dimensions of anxiety sensitivity as accentuated in more recent versions of the ASI. Likewise, these results cannot speak to the taxonic structure of the ASI. Fourth, although agoraphobic avoidance was revealed to be of core rele- vance in these analyses, other subtle aspects of avoidance—such as cognitive avoidance, utilization of safety signals, and so forth— were not specifically assessed, and therefore the relevance of these and other unassessed factors could not be modeled. Fifth, the examined process variables as well as the outcome variables are not without overlap. Whereas this is not unique to this study, construct overlap is extremely difficult to avoid in psychological research. To test the robustness of these results, we modeled several variations (i.e., with and without inclusion of the avoidance subscale on the PAS) and did not find any noticeable effect on the pattern of results. Sixth, it should be noted that not even sophis- ticated statistical analyses such as bivariate latent difference score modeling can establish the theoretical concepts, processes, and theories under investigation. Instead, statistical analysis is one approach to examining the process-outcome relations (Kenny et al., 2003; McArdle & Nesselroade, 1994). Finally, although the study from which these data are derived was partly designed to facilitate these types of analyses (i.e., assessment strategy) and had significantly more power than previous studies, the study did not randomize across the theoretical concepts under consideration. As such, the results should be considered post hoc in nature, and appropriate caution should be used in their interpretation.

Using bivariate latent difference score modeling, this study contributed to the understanding of processes underlying treatment in several ways. First, we replicated the importance of attribution variables consistently implicated in the process-outcome relation- ships in previous studies. Second, we expanded the list of process variables to agoraphobic avoidance and psychological flexibility. Third, we found clear evidence for the differentiation of medita- tional effects across outcomes (symptomatology vs. functioning). Finally, and most importantly, we found evidence that putative process variables are associated with changes in outcomes differ- ently at different stages in the treatment. Increasing the time resolution under investigation allows for a better understanding of how processes unfold over time by overcoming a limitation of cross-sectional data—namely, that they leave open the possibility that multiple constructs are relevant, but that they exert their effect at different points of time during therapy. The results in this study and similar studies have the potential to augment the effects of our current treatment and help therapists better deliver the treatments. The results point to specific processes at work and the timing of these processes. If substantiated within and across disorders, re- sults like these may help the sizeable minority of patients who do

not respond (Hofmann & Smits, 2008) and/or potentially improve the long-term prospects of patients, which is currently unclear (Durham, Higgins, Chambers, Swan, & Dow, 2012).

References

American Psychiatric Association. (2000). Diagnostic and statistical man- ual of mental disorders (4th ed., text rev.). Washington, DC: Author.

Arch, J. J., & Craske, M. G. (2008). Acceptance and commitment therapy and cognitive behavioral therapy for anxiety disorders: Different treat- ments, similar mechanisms? Clinical Psychology: Science and Practice, 15, 263–279. doi:10.1111/j.1468-2850.2008.00137.x

Bandelow, B. (1995). Assessing the efficacy of treatments for panic dis- order and agoraphobia: II. The Panic and Agoraphobia Scale. Interna- tional Clinical Psychopharmacology, 10, 73– 82. doi:10.1097/ 00004850-199506000-00003

Barlow, D. H., Gorman, J. M., Shear, M. K., & Woods, S. W. (2000). Cognitive-behavioral therapy, imipramine, or their combination for panic disorder: A randomized controlled trial. JAMA: Journal of the American Medical Association, 283, 2529–2536. doi:10.1001/jama.283 .19.2529

Baron, R. M., & Kenny, D. A. (1986). The moderator–mediator variable distinction in social psychological research: Conceptual, strategic, and statistical considerations. Journal of Personality and Social Psychology, 51, 1173–1182. doi:10.1037/0022-3514.51.6.1173

Bond, F. W., Hayes, S. C., Baer, R. A., Carpenter, K. M., Guenole, N., Orcutt, H. K., . . . Zettle, R. D. (2011). Preliminary psychometric properties of the Acceptance and Action Questionnaire–II: A revised measure of psychological inflexibility and experiential avoidance. Be- havior Therapy, 42, 676–688.

Bouchard, S., Gauthier, J., Nouwen, A., Ivers, H., Vallières, A., Simard, S., & Fournier, T. (2007). Temporal relationship between dysfunctional beliefs, self-efficacy and panic apprehension in the treatment of panic disorder with agoraphobia. Journal of Behavior Therapy and Experi- mental Psychiatry, 38, 275–292. doi:10.1016/j.jbtep.2006.08.002

Casey, L. M., Newcombe, P. A., & Oei, T. P. S. (2005). Cognitive mediation of panic severity: The role of catastrophic misinterpretation of bodily sensations and panic self-efficacy. Cognitive Therapy and Re- search, 29, 187–200. doi:10.1007/s10608-005-3164-3

Chambless, D. L., Caputo, G. C., Bright, P., & Gallagher, R. (1984). Assessment of fear of fear in agoraphobics: The Body Sensations Ques- tionnaire and the Agoraphobic Cognitions Questionnaire. Journal of Consulting and Clinical Psychology, 52, 1090–1097. doi:10.1037/0022- 006X.52.6.1090

Chambless, D. L., Caputo, G. C., Jasin, S. E., Gracely, E. J., & Williams, C. (1985). The Mobility Inventory for Agoraphobia. Behaviour Re- search and Therapy, 23, 35–44. doi:10.1016/0005-7967(85)90140-8

Clark, D. M. (1986). A cognitive approach to panic. Behaviour Research and Therapy, 24, 461–470. doi:10.1016/0005-7967(86)90011-2

Domschke, K., Reif, A., Weber, H., Richter, J., Hohoff, C., Ohrmann, P., . . . Deckert, J. (2011). Neuropeptide S receptor gene—Converging evidence for a role in panic disorder. Molecular Psychiatry, 16, 938– 948. doi:10.1038/mp.2010.81

Durham, R. C., Higgins, C., Chambers, J. A., Swan, J. S., & Dow, M. G. (2012). Long-term outcome of eight clinical trials of CBT for anxiety disorders: Symptom profile of sustained recovery and treatment-resistant groups. Journal of Affective Disorders, 136, 875–881. doi:10.1016/j.jad .2011.09.017

Ehlers, A. (1995). A 1-year prospective study of panic attacks: Clinical course and factors associated with maintenance. Journal of Abnormal Psychology, 104, 164–172. doi:10.1037/0021-843X.104.1.164

Forman, E. M., Herbert, J. D., Moitra, E., Yeomans, P. D., & Geller, P. A. (2007). A randomized controlled effectiveness trial of acceptance and

T hi

s do

cu m

en t

is co

py ri

gh te

d by

th e

A m

er ic

an Ps

yc ho

lo gi

ca l

A ss

oc ia

tio n

or on

e of

its al

lie d

pu bl

is he

rs .

T hi

s ar

tic le

is in

te nd

ed so

le ly

fo r

th e

pe rs

on al

us e

of th

e in

di vi

du al

us er

an d

is no

t to

be di

ss em

in at

ed br

oa dl

y.

151TIMING MATTERS

commitment therapy and cognitive therapy for anxiety and depression. Behavior Modification, 31, 772–799. doi:10.1177/0145445507302202

Gloster, A. T., Hummel, K. V., Lyudmirskaya, I., Hauke, T., & Sonntag, R. F. (2012). Aspects of exposure therapy in acceptance and commit- ment therapy. In P. Neudeck, & H.-U. Wittchen (Eds.), Exposure ther- apy: Rethinking the model—Refining the method (pp. 127–152). doi: 10.1007/978-1-4614-3342-2_8

Gloster, A. T., Klotsche, J., Chaker, S., Hummel, K. V., & Hoyer, J. (2011). Assessing psychological flexibility: What does it add above and beyond existing constructs? Psychological Assessment, 23, 970–982. doi:10.1037/a0024135

Gloster, A. T., Wittchen, H.-U., Einsle, F., Hofler, M., Lang, T., Helbig- Lang, S., . . . Arolt, V. (2009). Mechanism of action in CBT (MAC): Methods of a multi-center randomized controlled trial in 369 patients with panic disorder and agoraphobia. European Archives of Psychiatry and Clinical Neuroscience, 259, 155–166. doi:10.1007/s00406-009- 0065-6

Gloster, A. T., Wittchen, H.-U., Einsle, F., Lang, T., Helbig-Lang, S., Fydrich, T., . . . Arolt, V. (2011). Psychological treatment for panic disorder with agoraphobia: A randomized controlled trial to examine the role of therapist-guided exposure in situ in CBT. Journal of Consulting and Clinical Psychology, 79, 406–420. doi:10.1037/a0023584

Guy, W. (1976). ECDEU assessment manual for psychopharmacology. Rockville, MD: U.S. Department of Health, Education, and Welfare.

Hawley, L. L., Ho, M. H., Zuroff, D. C., & Blatt, S. J. (2006). The relationship of perfectionism, depression, and therapeutic alliance during treatment for depression: Latent difference score analysis. Journal of Consulting and Clinical Psychology, 74, 930–942.

Hoffart, A. (1995). Cognitive mediators of situational fear in agoraphobia. Journal of Behavior Therapy and Experimental Psychiatry, 26, 313– 320. doi:10.1016/0005-7916(95)00043-7

Hoffart, A., Sexton, H., Hedley, L., & Martinsen, E. (2008). Mechanisms of change in cognitive therapy for panic disorder with agoraphobia. Journal of Behavior Therapy and Experimental Psychiatry, 39, 262– 275. doi:10.1016/j.jbtep.2007.07.006

Hofmann, S. G., Meuret, A. E., Rosenfield, D., Suvak, M. K., Barlow, D. H., Gorman, J. M., . . . Woods, S. W. (2007). Preliminary evidence for cognitive mediation during cognitive-behavioral therapy of panic disorder. Journal of Consulting and Clinical Psychology, 75, 374–379. doi:10.1037/0022-006X.75.3.374

Hofmann, S. G., & Smits, J. A. J. (2008). Cognitive-behavioral therapy for adult anxiety disorders: A meta-analysis of randomized placebo- controlled trials. Journal of Clinical Psychiatry, 69, 621–632. doi: 10.4088/JCP.v69n0415

Hofmann, S. G., & Spiegel, D. A. (1999). Panic control treatment and its applications. Journal of Psychotherapy Practice & Research, 8, 3–11.

Kämpfe, C. K., Gloster, A. T., Wittchen, H.-U., Helbig-Lang, S., Lang, T., Gerlach, A. L., . . . Deckert, J. (2012). Experiential avoidance and anxiety sensitivity in patients with panic disorder and agoraphobia: Do both constructs measure the same? International Journal of Clinical and Health Psychology, 12, 5–22.

Karekla, M., Forsyth, J. P., & Kelly, M. M. (2004). Emotional avoidance and panicogenic responding to a biological challenge procedure. Behav- ior Therapy, 35, 725–746. doi:10.1016/S0005-7894(04)80017-0

Kashdan, T. B., Zvolensky, M. J., & McLeish, A. C. (2008). Anxiety sensitivity and affect regulatory strategies: Individual and interactive risk factors for anxiety-related symptoms. Journal of Anxiety Disorders, 22, 429–440. doi:10.1016/j.janxdis.2007.03.011

Kazdin, A. E. (2007). Mediators and mechanisms of change in psycho- therapy research. Annual Review of Clinical Psychology, 3, 1–27. doi: 10.1146/annurev.clinpsy.3.022806.091432

Kenny, D. A., Korchmaros, J. D., & Bolger, N. (2003). Lower level mediation in multilevel models. Psychological Methods, 8, 115–128. doi:10.1037/1082-989X.8.2.115

Kircher, T., Arolt, V., Jansen, A., Pyka, M., Reinhardt, I., Kellermann, T., . . . Straube, B. (2013). Effect of cognitive behavioural therapy on neural correlates of fear conditioning in panic disorder. Biological Psychiatry, 73, 93–101. doi:10.1016/j.biopsych.2012.07.026

Lang, T., Helbig-Lang, S., Westphal, D., Gloster, A. T., & Wittchen, H. U. (2011). Expositionsbasierte Therapie der Panikstörung mit Agorapho- bie: Ein Behandlungsmanual [Exposure-based therapy for panic disor- der with agoraphobia: A treatment manual]. Göttingen, Germany: Hogrefe.

McArdle, J. J., & Nesselroade, J. R. (1994). Using multivariate data to structure developmental change. In S. H. Cohen & H. W. Reese (Eds.), Life-span developmental psychology: Methodological contributions (pp. 223–267). Hillsdale, NJ: Erlbaum.

Meulenbeek, P., Spinhoven, P., Smit, F., Van Balkom, A., & Cuijpers, P. (2010). Cognitive mediation of panic reduction during an early inter- vention for panic. Acta Psychiatrica Scandinavica, 122, 20–29. doi: 10.1111/j.1600-0447.2009.01530.x

Meuret, A. E., Rosenfield, D., Seidel, A., Bhaskara, L., & Hofmann, S. G. (2010). Respiratory and cognitive mediators of treatment change in panic disorder: Evidence for intervention specificity. Journal of Con- sulting and Clinical Psychology, 78, 691–704. doi:10.1037/a0019552

Muthén, B., & Muthén, L. (2011). Mplus Version 6.1 software user’s guide, 1998–2008. Los Angeles, CA: Muthén & Muthén.

Peterson, R., & Reiss, S. (1993). Anxiety Sensitivity Index—Revised test manual. Columbus, OH: IDS Publishing.

Richter, J., Hamm, A. O., Pané-Farré, C. A., Gerlach, A. L., Gloster, A. T., Wittchen, H. U., . . . Arolt, V. (2012). Dynamics of defensive reactivity in patients with panic disorder and agoraphobia: Implications for the etiology of panic disorder. Biological Psychiatry, 72, 512–520. doi: 10.1016/j.biopsych.2012.03.035

Rodriguez, B. F., Bruce, S. E., Pagano, M. E., Spencer, M. A., & Keller, M. P. (2004). Factor structure and stability of the Anxiety Sensitivity Index in a longitudinal study of anxiety disorder patients. Behaviour Research and Therapy, 42, 79–91. doi:10.1016/S0005-7967(03)00074-3

Selig, J. P., & Preacher, K. J. (2009). Mediation models for longitudinal data in developmental research. Research in Human Development, 6, 144–164.

Shear, M. K., Vander Bilt, J., Rucci, P., Endicott, J., Lydiard, B., Otto, M. W., . . . Frank, D. M. (2001). Reliability and validity of a structured interview guide for the Hamilton Anxiety Rating Scale (SIGH-A). Depression and Anxiety, 13, 166–178.

Singer, J. D., & Willett, J. B. (2003). Applied longitudinal data analysis. Oxford, England: Oxford University Press. doi:10.1093/acprof:oso/ 9780195152968.001.0001

Smits, J. A. J., Powers, M. B., Cho, Y. R., & Telch, M. J. (2004). Mechanism of change in cognitive-behavioral treatment of panic disor- der: Evidence for the fear of fear mediational hypothesis. Journal of Consulting and Clinical Psychology, 72, 646–652. doi:10.1037/0022- 006X.72.4.646

Taylor, S., & Cox, B. J. (1998). An expanded Anxiety Sensitivity Index: Evidence for a hierarchic structure in a clinical sample. Journal of Anxiety Disorders, 12, 463–483. doi:10.1016/S0887-6185(98)00028-0

Taylor, S., Zvolensky, M. J., Cox, B. J., Deacon, B., Heimberg, R. G., Ledley, D. R., . . . Cardenas, S. J. (2007). Robust dimensions of anxiety sensitivity: Development and initial validation of the Anxiety Sensitivity Index—3. Psychological Assessment, 19, 176–188. doi:10.1037/1040- 3590.19.2.176

Teachman, B. A., Marker, C. D., & Clerkin, E. M. (2010). Catastrophic misinterpretations as a predictor of symptom change during treatment for panic disorder. Journal of Consulting and Clinical Psychology, 78, 964–973. doi:10.1037/a0021067

Vögele, C., Ehlers, A., Meyer, A. H., Frank, M., Hahlweg, K., & Margraf, J. (2010). Cognitive mediation of clinical improvement after intensive

T hi

s do

cu m

en t

is co

py ri

gh te

d by

th e

A m

er ic

an Ps

yc ho

lo gi

ca l

A ss

oc ia

tio n

or on

e of

its al

lie d

pu bl

is he

rs .

T hi

s ar

tic le

is in

te nd

ed so

le ly

fo r

th e

pe rs

on al

us e

of th

e in

di vi

du al

us er

an d

is no

t to

be di

ss em

in at

ed br

oa dl

y.

152 GLOSTER ET AL.

exposure therapy of agoraphobia and social phobia. Depression and Anxiety, 27, 294–301. doi:10.1002/da.20651

Wittchen, H.-U., Gloster, A. T., Beesdo-Baum, K., Fava, G. A., & Craske, M. G. (2010). Agoraphobia: A review of the diagnostic classificatory position and criteria. Depression and Anxiety, 27, 113–133. doi:10.1002/ da.20646

Wittchen, H.-U., Nocon, A., Beesdo, K., Pine, D. S., Hoefler, M., Lieb, R., & Gloster, A. T. (2008). Agoraphobia and panic. Psychotherapy and Psychosomatics, 77, 147–157. doi:10.1159/000116608

Zgourides, G. D., Warren, R., & Englert, M. E. (1989). Further evidence of construct validity for the Agoraphobic Cognitions Questionnaire and the

Body Sensations Questionnaire. Psychological Reports, 64, 590. doi: 10.2466/pr0.1989.64.2.590

Zvolensky, M. J., & Forsyth, J. P. (2002). Anxiety sensitivity dimensions in the prediction of body vigilance and emotional avoidance. Cognitive Therapy and Research, 26, 449–460. doi:10.1023/A:1016223716132

Received November 4, 2011 Revision received June 12, 2013

Accepted July 29, 2013 �

Members of Underrepresented Groups: Reviewers for Journal Manuscripts Wanted

If you are interested in reviewing manuscripts for APA journals, the APA Publications and Communications Board would like to invite your participation. Manuscript reviewers are vital to the publications process. As a reviewer, you would gain valuable experience in publishing. The P&C Board is particularly interested in encouraging members of underrepresented groups to participate more in this process.

If you are interested in reviewing manuscripts, please write APA Journals at [email protected]. Please note the following important points:

• To be selected as a reviewer, you must have published articles in peer-reviewed journals. The experience of publishing provides a reviewer with the basis for preparing a thorough, objective review.

• To be selected, it is critical to be a regular reader of the five to six empirical journals that are most central to the area or journal for which you would like to review. Current knowledge of recently published research provides a reviewer with the knowledge base to evaluate a new submission within the context of existing research.

• To select the appropriate reviewers for each manuscript, the editor needs detailed information. Please include with your letter your vita. In the letter, please identify which APA journal(s) you are interested in, and describe your area of expertise. Be as specific as possible. For example, “social psychology” is not sufficient—you would need to specify “social cognition” or “attitude change” as well.

• Reviewing a manuscript takes time (1–4 hours per manuscript reviewed). If you are selected to review a manuscript, be prepared to invest the necessary time to evaluate the manuscript thoroughly.

APA now has an online video course that provides guidance in reviewing manuscripts. To learn more about the course and to access the video, visit http://www.apa.org/pubs/authors/review- manuscript-ce-video.aspx.

T hi

s do

cu m

en t

is co

py ri

gh te

d by

th e

A m

er ic

an Ps

yc ho

lo gi

ca l

A ss

oc ia

tio n

or on

e of

its al

lie d

pu bl

is he

rs .

T hi

s ar

tic le

is in

te nd

ed so

le ly

fo r

th e

pe rs

on al

us e

of th

e in

di vi

du al

us er

an d

is no

t to

be di

ss em

in at

ed br

oa dl

y.

153TIMING MATTERS

  • Timing Matters: Change Depends on the Stage of Treatment in Cognitive Behavioral Therapy for Pan ...
    • Method
      • Design
      • Participants
      • Treatment
      • Assessment
        • PD/AG symptomatology and clinical functioning
          • Panic and Agoraphobia Scale (PAS; Bandelow, 1995)
          • CGI–Severity Subscale–Functioning Item (Guy, 1976)
        • Process variables
          • Agoraphobic Cognitions Questionnaire (ACQ; Chambless et al., 1984)
          • Bodily Sensations Questionnaire (BSQ; Chambless et al., 1984)
          • Mobility Inventory (MI; Chambless et al., 1985)
          • Anxiety Sensitivity Inventory (ASI; Peterson & Reiss, 1993)
          • Acceptance and Action Questionnaire–II (AAQ-II; Bond et al., 2011)
      • Statistical Analyses
    • Results
      • Baseline Values
      • Lower Level Mediation Models
      • Univariate Latent Difference Score Models
      • Bivariate Latent Difference Score Models
        • Panic and agoraphobia symptoms
        • Clinical functioning
        • Treatment condition
    • Discussion
    • References