* Encoding: UTF-8.
* Encoding: .

**************************This syntax will analysis the data to produce the results for the alcohol shopper lab study. 
**************************The syntax was created by KD (Bristol Univeristy) on 12/03/2020.
**************************The syntax was checked and edited by MP (Cambridge Univeristy) . Checked as of 10/06/2020.  

*The file "ShopperLabStudy_CodeForRecoding_SPSS(v24)_v1.2_20200303" must be run prior to running the below code. 

*Discussed clustering coding with RM (Bristol Univeristy) and MP (Cambridge Univeristy) May 2020. Decision to keep cluster but use GENLINMIXED command.
*Decision not to use transformed primary outcome. Instead use counts where possible (e.g.Totalnumberofalcoholicdrinksselected OF Totalnumberofdrinksselected).

*399 participants are included in the dataset. Missing data are labelled below as it varies by outcome variable. 

* ***********************************demographics*****************************************************

CTABLES
  /VLABELS VARIABLES=Gender Highestqualification Ethnicity BMIcategory Condition_1_3 DISPLAY=BOTH
  /TABLE Gender [C] + Highestqualification [C] + Ethnicity [C] + BMIcategory [C] BY Condition_1_3 
    [C][COUNT F40.0, COLPCT.COUNT PCT40.1]
  /CATEGORIES VARIABLES=Gender Highestqualification Ethnicity BMIcategory Condition_1_3 ORDER=A 
    KEY=VALUE EMPTY=INCLUDE
  /CRITERIA CILEVEL=95.

* Custom Tables.
CTABLES
  /VLABELS VARIABLES=Age BMI Condition_1_3 DISPLAY=BOTH
  /TABLE Age [S][MEAN, STDDEV] + BMI [S][MEAN, STDDEV] BY Condition_1_3 [C]
  /CATEGORIES VARIABLES=Condition_1_3 ORDER=A KEY=VALUE EMPTY=INCLUDE
  /CRITERIA CILEVEL=95.

* Custom Tables.
CTABLES
  /VLABELS VARIABLES=Weeklyconsumptioncategory Condition_1_3 DISPLAY=BOTH
  /TABLE Weeklyconsumptioncategory [C] BY Condition_1_3 [C][COUNT F40.0, COLPCT.COUNT PCT40.1]
  /CATEGORIES VARIABLES=Weeklyconsumptioncategory Condition_1_3 ORDER=A KEY=VALUE EMPTY=INCLUDE
  /CRITERIA CILEVEL=95.

* Custom Tables.
CTABLES
  /VLABELS VARIABLES=AUDITscore Weeklyalcoholconsumed Weeklyalcoholpurchased Condition_1_3 
    DISPLAY=BOTH
  /TABLE AUDITscore [S][MEAN, STDDEV] + Weeklyalcoholconsumed [S][MEAN, STDDEV] + 
    Weeklyalcoholpurchased [S][MEAN, STDDEV] BY Condition_1_3 [C]
  /CATEGORIES VARIABLES=Condition_1_3 ORDER=A KEY=VALUE EMPTY=INCLUDE
  /CRITERIA CILEVEL=95.

*******************************************descriptives of the outcomes******************************

CTABLES
  /VLABELS VARIABLES=Primary_outcome_proportion Mean_alcohol_units_1a 
    Numberofalcoholunitsselected_1b Totalspendonalcoholicdrinks_2a 
    Totalspendonalcoholicdrink_proportion_2b Condition_1_3 DISPLAY=BOTH
  /TABLE Primary_outcome_proportion [S][MEAN, STDDEV, VALIDN F40.0] + Mean_alcohol_units_1a 
    [S][MEAN, STDDEV, VALIDN F40.0] + Numberofalcoholunitsselected_1b [S][MEAN, STDDEV, VALIDN F40.0] + 
    Totalspendonalcoholicdrinks_2a [S][MEAN, STDDEV, VALIDN F40.0] + 
    Totalspendonalcoholicdrink_proportion_2b [S][MEAN, STDDEV, VALIDN F40.0] BY Condition_1_3 [C]
  /CATEGORIES VARIABLES=Condition_1_3 ORDER=A KEY=VALUE EMPTY=INCLUDE
  /CRITERIA CILEVEL=95.

CTABLES
  /VLABELS VARIABLES=Selectedalcohol_3 Condition_1_3 DISPLAY=BOTH
  /TABLE Selectedalcohol_3 [C] BY Condition_1_3 [C][COUNT F40.0, COLPCT.COUNT PCT40.1]
  /CATEGORIES VARIABLES=Selectedalcohol_3 Condition_1_3 ORDER=A KEY=VALUE EMPTY=INCLUDE
  /CRITERIA CILEVEL=95.


CTABLES
  /VLABELS VARIABLES=NEAMEANAVERAGE_IMAGEandTEXT_4a NEAMEANAVERAGE_TEXTONLY_4a DISPLAY=BOTH
  /TABLE BY NEAMEANAVERAGE_IMAGEandTEXT_4a [S][MEAN, STDDEV] + NEAMEANAVERAGE_TEXTONLY_4a [S][MEAN, 
    STDDEV]
  /CRITERIA CILEVEL=95.

CTABLES
  /VLABELS VARIABLES=Acceptability_IMAGEandTEXT_4b Acceptability_TEXTONLY_4b DISPLAY=BOTH
  /TABLE BY Acceptability_IMAGEandTEXT_4b [S][MEAN, STDDEV, VALIDN F40.0] + 
    Acceptability_TEXTONLY_4b [S][MEAN, STDDEV, VALIDN F40.0]
  /CRITERIA CILEVEL=95.

*****************************************Primary outcome - Proportion of drinks selected that are alcoholic (n = 392)************************************ 

*full factorial model first* 
*As per SAP model includes condition, weekday/weekend, morning/afternoon and the 2-way and 3-way interactions.
*counts used*
*with cluster as per SAP (statistical analysis plan)*
*not reported in manuscript because too many parameters for too little data*

GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Totalnumberofalcoholicdrinksselected TRIALS=FIELD(Totalnumberofdrinksselected)  
    OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=BINOMIAL LINK=LOGIT
  /FIXED  EFFECTS=Condition_1_3 Weekday_weekend Morning_afternoon Weekday_weekend*Morning_afternoon 
    Condition_1_3*Weekday_weekend Condition_1_3*Morning_afternoon 
    Condition_1_3*Weekday_weekend*Morning_afternoon USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001
  /EMMEANS TABLES=Condition_1_3 COMPARE=Condition_1_3 CONTRAST=PAIRWISE 
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD.ILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001.

*this  model includes three main effects of  condition, weekday/weekend, morning/afternoon only (with no interactions)*
*with cluster as per the SAP
*counts used*
*now a better fitted model (goodness of fit table reduced)* SE smaller*
*this is the best model (and so is reported in manuscript)*

GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Totalnumberofalcoholicdrinksselected TRIALS=FIELD(Totalnumberofdrinksselected)  
    OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=BINOMIAL LINK=LOGIT
  /FIXED  EFFECTS=Condition_1_3 Weekday_weekend Morning_afternoon USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /EMMEANS TABLES=Condition_1_3 COMPARE=Condition_1_3 CONTRAST=PAIRWISE 
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD.


*As a sensitivity analysis (not in the SAP) add in Weeklyalcoholconsumed as a covariate (as Weeklyalcoholconsumed was strongly related to outcome) (no interatcions).
*Model slightly improved (goodness of fit table + small decrease in SE)
*reported as a sensitivity analysis in manuscript*

GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Totalnumberofalcoholicdrinksselected TRIALS=FIELD(Totalnumberofdrinksselected)  
    OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=BINOMIAL LINK=LOGIT
  /FIXED  EFFECTS=Condition_1_3 Weekday_weekend Morning_afternoon Weeklyalcoholconsumed 
    USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001
  /EMMEANS TABLES=Condition_1_3 COMPARE=Condition_1_3 CONTRAST=PAIRWISE 
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD.

*As a sensitivity analysis (not in the SAP) add in age, gender, BMI and education as covariates. as well as weekly alcohol consumption. (no interatcions).
*Model not improved. Will be reported in manuscript as requested by reviewer. 

RECODE BMIcategory (1=6) (2=5) (3=4) (4=3) (5=2) (6=1) (ELSE=Copy) INTO BMIcategory_rev.
EXECUTE.

VALUE LABELS BMIcategory_rev
    1 "Morbidly obese"
    2 "Severely obese"
    3 "Obese"
    4 "Overweight"
    5 "Healthy weight"
    6 " Underweight".
EXECUTE.

GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Totalnumberofalcoholicdrinksselected TRIALS=FIELD(Totalnumberofdrinksselected)  
    OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=BINOMIAL LINK=LOGIT
  /FIXED  EFFECTS=Condition_1_3 Weekday_weekend Morning_afternoon Weeklyalcoholconsumed Age Gender 
    BMIcategory_rev Highestqualification USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001
  /EMMEANS TABLES=Condition_1_3 COMPARE=Condition_1_3 CONTRAST=PAIRWISE 
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD.

**********************************************secondary outcomes********************************

***interactions fitted for all secondary ouctomes, but not reported in manuscript because too many parameters for too little data*

******Mean alcohol units purchased - secondary outcome 1a (n=392)***

*the model includes three main effects of  condition, weekday/weekend, morning/afternoon only (no interactions)*
*with cluster as per the SAP
*counts used*
*now a better fitted model (goodness of fit table reduced) compared to model with interactions* SE smaller*
*this is the best model (will be reported in manuscript)*

GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Mean_alcohol_units_1a TRIALS=NONE OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=NORMAL LINK=IDENTITY
  /FIXED  EFFECTS=Condition_1_3 Weekday_weekend Morning_afternoon USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001
  /EMMEANS TABLES=Condition_1_3 COMPARE=Condition_1_3 CONTRAST=PAIRWISE 
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD.

****************Total alcohol units purchased - secondary outcome 1b (n=396)***

* the model includes three main effects of  condition, weekday/weekend, morning/afternoon only (no interactions)*
*with cluster as per the SAP
*now a better fitted model (goodness of fit table reduced) compared to model with interactions* SE smaller*
*this is the best model (will be reported in manuscript)*

GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Numberofalcoholunitsselected_1b TRIALS=NONE OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=NORMAL LINK=IDENTITY
  /FIXED  EFFECTS=Condition_1_3 Weekday_weekend Morning_afternoon USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001
  /EMMEANS TABLES=Condition_1_3 COMPARE=Condition_1_3 CONTRAST=PAIRWISE 
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD.

***********************total spend on alcoholic drinks - secondary outcome 2a (n = 396)***************

* the model includes three main effects of  condition, weekday/weekend, morning/afternoon only*
*with cluster as per the SAP
*now a better fitted model (goodness of fit table reduced) compared to a model with interactions* SE smaller*
*this is the best model (will be reported in manuscript)*

GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Totalspendonalcoholicdrinks_2a TRIALS=NONE OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=NORMAL LINK=IDENTITY
  /FIXED  EFFECTS=Condition_1_3 Weekday_weekend Morning_afternoon USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001
  /EMMEANS TABLES=Condition_1_3 COMPARE=Condition_1_3 CONTRAST=PAIRWISE 
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD.

**********************Spend (£) on alcoholic drinks as a proportion of total spend  - secondary outcome 2b (n=392)*******************

* the model includes three main effects of  condition, weekday/weekend, morning/afternoon only*
*with cluster as per the SAP
*now a better fitted model (goodness of fit table reduced) compared to a model with interactions* SE smaller*
*this is the best model (will be reported in manuscript)*

GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Totalspendonalcoholicdrink_proportion_2b TRIALS=NONE OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=NORMAL LINK=IDENTITY
  /FIXED  EFFECTS=Condition_1_3 Weekday_weekend Morning_afternoon USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001
  /EMMEANS TABLES=Condition_1_3 CONTRAST=NONE 
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD.

*****************************Proportion of participants selecting at least one alcoholic beverage in the selection task - secondary outcome 3 (n = 396)*

CTABLES
  /VLABELS VARIABLES=Condition_1_3 Selectedalcohol_3 DISPLAY=BOTH
  /TABLE Condition_1_3 [C][COUNT F40.0, ROWPCT.COUNT PCT40.1] BY Selectedalcohol_3 [C]
  /CATEGORIES VARIABLES=Condition_1_3 Selectedalcohol_3 ORDER=A KEY=VALUE EMPTY=INCLUDE
  /CRITERIA CILEVEL=95.

GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Selectedalcohol_3 TRIALS=NONE OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=BINOMIAL LINK=LOGIT
  /FIXED  EFFECTS=Condition_1_3 Weekday_weekend Morning_afternoon USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001
  /EMMEANS TABLES=Condition_1_3 COMPARE=Condition_1_3 CONTRAST=PAIRWISE 
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD.

***************************Post intervention ratings (within-subjects): NEA and acceptability. - secondary outcome 4a and 4b***

CTABLES
  /VLABELS VARIABLES=NEAMEANAVERAGE_TEXTONLY_4a NEAMEANAVERAGE_IMAGEandTEXT_4a 
    Acceptability_TEXTONLY_4b Acceptability_IMAGEandTEXT_4b DISPLAY=NAME
  /TABLE NEAMEANAVERAGE_TEXTONLY_4a [MEAN, STDDEV] + NEAMEANAVERAGE_IMAGEandTEXT_4a [MEAN, STDDEV] 
    + Acceptability_TEXTONLY_4b [MEAN, STDDEV] + Acceptability_IMAGEandTEXT_4b [MEAN, STDDEV]
  /CRITERIA CILEVEL=95.

GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Difference_NEA_4a TRIALS=NONE OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=NORMAL LINK=IDENTITY
  /FIXED  EFFECTS=Condition_1_3 USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD.

GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Difference_Acceptability_4b TRIALS=NONE OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=NORMAL LINK=IDENTITY
  /FIXED  EFFECTS=Condition_1_3 USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD.

***with added adjustment for order (Order_within_subjects)*** not used in manuscript as uncertain about order variable***

GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Difference_NEA_4a TRIALS=NONE OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=NORMAL LINK=IDENTITY
  /FIXED  EFFECTS=Condition_1_3 Order_within_subjects USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD
  /SAVE PEARSON_RESIDUALS(PearsonResidual).

GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Difference_Acceptability_4b TRIALS=NONE OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=NORMAL LINK=IDENTITY
  /FIXED  EFFECTS=Condition_1_3 Order_within_subjects USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD
  /SAVE PEARSON_RESIDUALS(PearsonResidual).

**********************************Sensitivity analysis  of the primary outcome**********************************

***Sensitivity analysis 1 (Participants from the two HWL groups that answered ‘no’ to the manipulation check (did you notice any labels on the products?) will be excluded for a sensitivity analysis of the primary outcome)***

IF  ((Postshoppingtaskquestion_1  = 1) and (Condition_1_3 = 1)) Sensitivity_analysis1= 0.
EXECUTE.
IF  ((Postshoppingtaskquestion_1  = 1) and (Condition_1_3 = 2)) Sensitivity_analysis1= 0.
EXECUTE.
IF  ((Postshoppingtaskquestion_1  = 1) and (Condition_1_3 = 3)) Sensitivity_analysis1= 1.
EXECUTE.
IF (Postshoppingtaskquestion_1  = 2) Sensitivity_analysis1= 1.
EXECUTE.

IF (ParticipantID = 1003 or ParticipantID = 1002 or ParticipantID = 1001) Sensitivity_analysis1 = -999.

VALUE LABELS Sensitivity_analysis1
0 "not included"
1 "included".
EXECUTE.

MISSING VALUES Sensitivity_analysis1 (-999).

***HOWEVER NC wanted this repeated using a different variable***
***Sensitivity analysis 1 (Participants from the two HWL groups that answered ‘no’ to the manipulation check (did you notice any labels on the products?) will be excluded for a sensitivity analysis of the primary outcome)***

If (Sensitivity_analysis1 = 3) and ((Condition_1_3 = 2)) Sensitivity_analysis1b= 0.
If (Sensitivity_analysis1 = 3) and ((Condition_1_3 = 1)) Sensitivity_analysis1b= 0.
If (Sensitivity_analysis1 = 3) and ((Condition_1_3 = 3)) Sensitivity_analysis1b= 1.


VALUE LABELS Sensitivity_analysis1b
0 "not included"
1 "included".
EXECUTE.

TEMPORARY.
SELECT IF Sensitivity_analysis1b = 1. 
GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Totalnumberofalcoholicdrinksselected TRIALS=FIELD(Totalnumberofdrinksselected)  
    OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=BINOMIAL LINK=LOGIT
  /FIXED  EFFECTS=Condition_1_3 Weekday_weekend Morning_afternoon USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001
  /EMMEANS TABLES=Condition_1_3 COMPARE=Condition_1_3 CONTRAST=PAIRWISE 
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD.

***Sensitivity analysis 2 (Participants identified by the onsite research team 
*(i.e. those that supplied nonsensical responses, did not adequately comply with the protocol, or did not engage with the task or follow instructions), will be excluded for a sensitivity analysis of the primary outcome)***

IF (ParticipantID NE 165) OR (ParticipantID NE 264) OR (ParticipantID NE 268) OR (ParticipantID NE 281) OR (ParticipantID NE 357) OR (ParticipantID NE 360) Sensitivity_analysis2 = 1.
IF (ParticipantID = 165) Sensitivity_analysis2 = 0.
IF (ParticipantID = 264) Sensitivity_analysis2 = 0.
IF (ParticipantID = 268) Sensitivity_analysis2 = 0.
IF (ParticipantID = 281) Sensitivity_analysis2 = 0.
IF (ParticipantID = 357) Sensitivity_analysis2 = 0.
IF (ParticipantID = 360) Sensitivity_analysis2 = 0.


VALUE LABELS Sensitivity_analysis2
0 "not included"
1 "included".
EXECUTE.

SELECT IF Sensitivity_analysis2 = 1. 
GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Totalnumberofalcoholicdrinksselected TRIALS=FIELD(Totalnumberofdrinksselected)  
    OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=BINOMIAL LINK=LOGIT
  /FIXED  EFFECTS=Condition_1_3 Weekday_weekend Morning_afternoon USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001
  /EMMEANS TABLES=Condition_1_3 COMPARE=Condition_1_3 CONTRAST=PAIRWISE 
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD.

*********************************************************************************************excluding outliers***********************

******Mean alcohol units purchased - secondary outcome 1a*** < 9.80. 14 excluded. 

TEMPORARY.
SELECT IF Mean_alcohol_units_1a < 9.80.
GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Mean_alcohol_units_1a TRIALS=NONE OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=NORMAL LINK=IDENTITY
  /FIXED  EFFECTS=Condition_1_3 Weekday_weekend Morning_afternoon USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001
  /EMMEANS TABLES=Condition_1_3 CONTRAST=NONE 
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD.

****************Total alcohol units purchased - secondary outcome 1b*** 7 excluded.

TEMPORARY.
SELECT IF Numberofalcoholunitsselected_1b < 36.
GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Numberofalcoholunitsselected_1b TRIALS=NONE OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=NORMAL LINK=IDENTITY
  /FIXED  EFFECTS=Condition_1_3 Weekday_weekend Morning_afternoon USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001
  /EMMEANS TABLES=Condition_1_3 COMPARE=Condition_1_3 CONTRAST=PAIRWISE 
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD.

***********************total spend on alcoholic drinks - secondary outcome 2a*************** 22 excluded

TEMPORARY.
SELECT IF Totalspendonalcoholicdrinks_2a < 25.
GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Totalspendonalcoholicdrinks_2a TRIALS=NONE OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=NORMAL LINK=IDENTITY
  /FIXED  EFFECTS=Condition_1_3 Weekday_weekend Morning_afternoon USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001
  /EMMEANS TABLES=Condition_1_3 COMPARE=Condition_1_3 CONTRAST=PAIRWISE 
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD.

**********************Spend (£) on alcoholic drinks as a proportion of total spend  - secondary outcome 2b******************* exclude 47 0s

TEMPORARY.
SELECT IF Totalspendonalcoholicdrink_proportion_2b > 0.
GENLINMIXED
  /DATA_STRUCTURE SUBJECTS=ParticipantID
  /FIELDS TARGET=Totalspendonalcoholicdrink_proportion_2b TRIALS=NONE OFFSET=NONE
  /TARGET_OPTIONS DISTRIBUTION=NORMAL LINK=IDENTITY
  /FIXED  EFFECTS=Condition_1_3 Weekday_weekend Morning_afternoon USE_INTERCEPT=TRUE
  /RANDOM EFFECTS=Cluster USE_INTERCEPT=FALSE SUBJECTS=ParticipantID 
    COVARIANCE_TYPE=VARIANCE_COMPONENTS 
  /BUILD_OPTIONS TARGET_CATEGORY_ORDER=ASCENDING INPUTS_CATEGORY_ORDER=ASCENDING MAX_ITERATIONS=100 
    CONFIDENCE_LEVEL=95 DF_METHOD=RESIDUAL COVB=ROBUST PCONVERGE=0.000001(ABSOLUTE) SCORING=0 
    SINGULAR=0.000000000001
  /EMMEANS TABLES=Condition_1_3 CONTRAST=NONE 
  /EMMEANS_OPTIONS SCALE=ORIGINAL PADJUST=LSD.

*********************************************************************************************Bootstrapping*********************

PRESERVE.
SET RNG=MT MTINDEX=20200612.
SHOW RNG.
BOOTSTRAP
  /SAMPLING METHOD=SIMPLE
  /VARIABLES TARGET=Mean_alcohol_units_1a INPUT=Condition_1_3 Weekday_weekend Morning_afternoon   
  /CRITERIA CILEVEL=95 CITYPE=PERCENTILE  NSAMPLES=1000
  /MISSING USERMISSING=EXCLUDE.
MIXED Mean_alcohol_units_1a BY Condition_1_3 Weekday_weekend Morning_afternoon
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Condition_1_3 Weekday_weekend Morning_afternoon | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /EMMEANS=TABLES(Condition_1_3) COMPARE ADJ(LSD).
RESTORE.

PRESERVE.
SET RNG=MT MTINDEX=20200612.
SHOW RNG.
BOOTSTRAP
  /SAMPLING METHOD=SIMPLE
  /VARIABLES TARGET=Numberofalcoholunitsselected_1b INPUT=Condition_1_3 Weekday_weekend 
    Morning_afternoon   
  /CRITERIA CILEVEL=95 CITYPE=PERCENTILE  NSAMPLES=1000
  /MISSING USERMISSING=EXCLUDE.
MIXED Numberofalcoholunitsselected_1b BY Condition_1_3 Weekday_weekend Morning_afternoon
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Condition_1_3 Weekday_weekend Morning_afternoon | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /EMMEANS=TABLES(Condition_1_3) COMPARE ADJ(LSD).
RESTORE.

PRESERVE.
SET RNG=MT MTINDEX=20200612.
SHOW RNG.
BOOTSTRAP
  /SAMPLING METHOD=SIMPLE
  /VARIABLES TARGET=Totalspendonalcoholicdrinks_2a INPUT=Condition_1_3 Weekday_weekend 
    Morning_afternoon   
  /CRITERIA CILEVEL=95 CITYPE=PERCENTILE  NSAMPLES=1000
  /MISSING USERMISSING=EXCLUDE.
MIXED Totalspendonalcoholicdrinks_2a BY Condition_1_3 Weekday_weekend Morning_afternoon
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Condition_1_3 Weekday_weekend Morning_afternoon | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /EMMEANS=TABLES(Condition_1_3) COMPARE ADJ(LSD).
RESTORE.

PRESERVE.
SET RNG=MT MTINDEX=20200612.
SHOW RNG.
BOOTSTRAP
  /SAMPLING METHOD=SIMPLE
  /VARIABLES TARGET=Totalspendonalcoholicdrink_proportion_2b INPUT=Condition_1_3 Weekday_weekend 
    Morning_afternoon   
  /CRITERIA CILEVEL=95 CITYPE=PERCENTILE  NSAMPLES=1000
  /MISSING USERMISSING=EXCLUDE.
MIXED Totalspendonalcoholicdrink_proportion_2b BY Condition_1_3 Weekday_weekend Morning_afternoon
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Condition_1_3 Weekday_weekend Morning_afternoon | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /EMMEANS=TABLES(Condition_1_3) COMPARE ADJ(LSD).
RESTORE.

PRESERVE.
SET RNG=MT MTINDEX=20200612.
SHOW RNG.
BOOTSTRAP
  /SAMPLING METHOD=SIMPLE
  /VARIABLES TARGET=Selectedalcohol_3 INPUT=Condition_1_3 Weekday_weekend Morning_afternoon  
  /CRITERIA CILEVEL=95 CITYPE=PERCENTILE  NSAMPLES=1000
  /MISSING USERMISSING=EXCLUDE.
LOGISTIC REGRESSION VARIABLES Selectedalcohol_3
  /METHOD=ENTER Condition_1_3 Weekday_weekend Morning_afternoon 
  /CONTRAST (Condition_1_3)=Indicator
  /PRINT=CI(95)
  /CRITERIA=PIN(0.05) POUT(0.10) ITERATE(20) CUT(0.5).
RESTORE.