*Before running, ensure that the dataset has an ID variable representing the unique client ID, and a date variable.  
*All CCAPS items should be labeled CCAPS_01 through CCAPS_62. Syntax to convert variable names from the old CCAPS 70 is available upon request.

IF  (nvalid(CCAPS_01,CCAPS_03,CCAPS_06,CCAPS_07,CCAPS_08,CCAPS_11,CCAPS_12,
CCAPS_16,CCAPS_19,CCAPS_21,CCAPS_22,CCAPS_24,CCAPS_28,CCAPS_31,CCAPS_33,CCAPS_34,CCAPS_37,
CCAPS_38,CCAPS_39,CCAPS_42,CCAPS_45,CCAPS_48,CCAPS_50,CCAPS_54,CCAPS_55, 
CCAPS_58,CCAPS_61,CCAPS_62)=0) Is_CCAPS34=1. 
EXECUTE. 
RECODE Is_CCAPS34 (SYSMIS=0). 
EXECUTE. 
IF  (Is_CCAPS34=0) Is_CCAPS62=1. 
EXECUTE. 
RECODE Is_CCAPS62 (SYSMIS=0). 
EXECUTE. 
IF  (Is_CCAPS62=0) CCAPS_type=34*Is_CCAPS34+62*Is_CCAPS62. 
EXECUTE. 
COMPUTE CCAPS_type=34*Is_CCAPS34+62*Is_CCAPS62. 
EXECUTE. 


*COMPUTE NUMBER OF MISSING ITEMS PER SUBSCALE.AND TOTAL.
Compute Depression34_MISS = NMISS(CCAPS_09, CCAPS_10,
CCAPS_20, CCAPS_23, CCAPS_40, CCAPS_46).
EXECUTE .
Compute Anxiety34_MISS = NMISS(CCAPS_04, CCAPS_14,
CCAPS_17, CCAPS_18, CCAPS_27, CCAPS_30).
EXECUTE.
Compute Social_Anxiety34_MISS = NMISS(CCAPS_02,
CCAPS_35, CCAPS_41, CCAPS_44, CCAPS_47).
EXECUTE.
Compute Academics34_MISS = NMISS(CCAPS_15, CCAPS_51,
CCAPS_53, CCAPS_59).
EXECUTE.
Compute Eating34_MISS = NMISS(CCAPS_05, CCAPS_13,
CCAPS_25).
EXECUTE.
Compute Hostility34_MISS = NMISS(CCAPS_32, CCAPS_36,
CCAPS_43, CCAPS_52, CCAPS_57, CCAPS_60).
EXECUTE.
Compute Alcohol34_MISS = NMISS(CCAPS_26, CCAPS_29,
CCAPS_49, CCAPS_56).
EXECUTE.
Compute DI34_MISS = NMISS(CCAPS_04, CCAPS_09,
CCAPS_10, CCAPS_14, CCAPS_17, CCAPS_18, CCAPS_20,
CCAPS_23, CCAPS_27, CCAPS_30, CCAPS_36, CCAPS_40,
CCAPS_41, CCAPS_43, CCAPS_46, CCAPS_47, CCAPS_51,
CCAPS_52, CCAPS_53).
Execute.
Compute CCAPS34_Tmiss = Sum(Depression34_MISS, Anxiety34_MISS, Social_Anxiety34_MISS, Academics34_MISS, Eating34_MISS, Hostility34_MISS, Alcohol34_MISS).
EXECUTE.

*Compute VARIANCE of all items provided.
COMPUTE VARIANCE=VARIANCE(CCAPS_01 TO CCAPS_62). 
EXECUTE.

*Indicate whether the subscale or administration is invalid by (a) 33% of subscale missing or (b) more than 50% of total missing or (c) variance=0.
IF  (Depression34_MISS > 2) or (CCAPS34_Tmiss>17) or (variance=0) Invalid_depression = 1.
Execute.
If sysmis(invalid_depression) invalid_depression = 0.
Execute.
If (Anxiety34_MISS > 2) or (CCAPS34_Tmiss>17) or (variance=0) Invalid_anxiety = 1.
Execute.
If sysmis(Invalid_anxiety) Invalid_anxiety = 0.
Execute.
If (Social_Anxiety34_MISS > 1) or (CCAPS34_Tmiss>17) or (variance=0) Invalid_social_anxiety = 1.
Execute. 
If sysmis(Invalid_social_anxiety) Invalid_social_anxiety = 0.
Execute. 
If (Academics34_MISS>1) or (CCAPS34_Tmiss>17) or (variance=0) Invalid_academics = 1.
Execute.
If sysmis(invalid_academics) invalid_academics = 0.
Execute.
If (Eating34_MISS > 1) or (CCAPS34_Tmiss>17) or (variance=0) Invalid_eating = 1.
Execute.
If sysmis(invalid_eating) invalid_eating = 0.
Execute.
If (Hostility34_MISS > 2) or (CCAPS34_Tmiss>17) or (variance=0) Invalid_hostility = 1.
Execute.
If sysmis(invalid_hostility) invalid_hostility = 0.
Execute.
If (Alcohol34_MISS > 1) or (CCAPS34_Tmiss>17) or (variance=0) Invalid_alcohol = 1.
Execute.
If sysmis(invalid_alcohol) invalid_alcohol = 0.
Execute.
If (DI34_MISS > 6)  or (CCAPS34_Tmiss>17) or (variance=0) Invalid_DI = 1.
Execute.
If sysmis(Invalid_DI) invalid_DI = 0.
execute.

*computing 34 subscales for valid cases.
IF (Invalid_depression=0)  
Depression34 = MEAN(CCAPS_09, CCAPS_10,
CCAPS_20, CCAPS_23, CCAPS_40, CCAPS_46).
Execute.
IF (Invalid_anxiety=0)  
Anxiety34 = MEAN(CCAPS_04, CCAPS_14,
CCAPS_17, CCAPS_18, CCAPS_27, CCAPS_30).
Execute.
IF (Invalid_social_anxiety = 0)  
Social_Anxiety34 = MEAN(CCAPS_02,
CCAPS_35, CCAPS_41, CCAPS_44, CCAPS_47).
Execute.
IF (Invalid_academics =0)  
Academics34 = MEAN(CCAPS_15, CCAPS_51,
CCAPS_53, CCAPS_59).
Execute.
IF (Invalid_eating=0)  
Eating34 = MEAN(CCAPS_05, CCAPS_13,
CCAPS_25).
Execute.
IF (Invalid_hostility=0)  
Hostility34 = MEAN(CCAPS_32, CCAPS_36,
CCAPS_43, CCAPS_52, CCAPS_57, CCAPS_60).
Execute.
IF (Invalid_alcohol=0)  
Alcohol34 = MEAN(CCAPS_26, CCAPS_29,
CCAPS_49, CCAPS_56).
Execute.
IF (Invalid_DI=0) 
DI = Mean(CCAPS_04, CCAPS_09,
CCAPS_10, CCAPS_14, CCAPS_17, CCAPS_18, CCAPS_20,
CCAPS_23, CCAPS_27, CCAPS_30, CCAPS_36, CCAPS_40,
CCAPS_41, CCAPS_43, CCAPS_46, CCAPS_47, CCAPS_51,
CCAPS_52, CCAPS_53).
EXECUTE.

Sort cases by ID (A) Date (A).

*Create a new dataset with one line per client.  
DATASET DECLARE Center_Change.
SORT CASES BY ID.
AGGREGATE
  /OUTFILE='Center_Change'
  /PRESORTED
  /BREAK=ID
  /Depression34_first=FIRST(Depression34)
  /Anxiety34_first=FIRST(Anxiety34)
  /Social_Anxiety34_first=FIRST(Social_Anxiety34)
  /Academics34_first=FIRST(Academics34)
  /Eating34_first=FIRST(Eating34)
  /Hostility34_first=FIRST(Hostility34)
  /Alcohol34_first=FIRST(Alcohol34)
  /DI_first=FIRST(DI)
  /Depression34_last=LAST(Depression34)
  /Anxiety34_last=LAST(Anxiety34)
  /Social_Anxiety34_last=LAST(Social_Anxiety34)
  /Academics34_last=LAST(Academics34)
  /Eating34_last=LAST(Eating34)
  /Hostility34_last=LAST(Hostility34)
  /Alcohol34_last=LAST(Alcohol34)
  /DI_last=LAST(DI)
  /CCAPS_N=N.

************Open new dataset created by previous syntax and run the rest of the syntax on that datset***************

*Create a variable indicating points of change on each subscale from first CCAPS to last CCAPS.  
******Positive values indicate a reduction in symptoms, while negative values indicate an increase in symptoms*****.
COMPUTE Depression_Change=Depression34_first-Depression34_last.
COMPUTE Anxiety_Change=Anxiety34_first-Anxiety34_last.
COMPUTE SocialAnxiety_Change=Social_Anxiety34_first-Social_Anxiety34_last.
COMPUTE Academics_Change=Academics34_first-Academics34_last.
COMPUTE Eating_Change=Eating34_first-Eating34_last.
COMPUTE Hostility_Change=Hostility34_first-Hostility34_last.
COMPUTE Eating_Change=Eating34_first-Eating34_last.
COMPUTE Hostility_Change=Hostility34_first-Hostility34_last.
COMPUTE Alcohol_Change=Alcohol34_first-Alcohol34_last.
COMPUTE DI_Change=DI_first-DI_last.
Execute.

*Filter out people from the datset who only had 1 CCAPS, for whom change can't be calculated.
USE ALL. 
COMPUTE filter_$=(CCAPS_N>1). 
VARIABLE LABELS filter_$ 'CCAPS_N>1 (FILTER)'. 
VALUE LABELS filter_$ 0 'Not Selected' 1 'Selected'. 
FORMATS filter_$ (f1.0). 
FILTER BY filter_$. 
EXECUTE.

*Descriptives of change for each subscale.
DESCRIPTIVES VARIABLES=Depression_Change Anxiety_Change SocialAnxiety_Change Academics_Change 
    Eating_Change Hostility_Change Alcohol_Change DI_Change
  /STATISTICS=MEAN STDDEV VARIANCE MIN MAX.