* Encoding: UTF-8.

*********************************Tableware study (Size and shape of plates and size of wine glasses and bottles: impact on selection of food and alcohol) ******************************************

***This syntax will analyse the data for the Tableware study.
***The syntax was created by Katie De-loyde (Bristol University) on 05/08/2020
***The syntax has been checked and edited by RM (Bristol University) and MP (Cambridge University) (August 2020)

***Run  'Tableware_DataRecode_v1.0_2020716' code first. 

*************************************************************************************wide format database***********************************
    
***demographic table***
    
CTABLES
  /VLABELS VARIABLES=Sex Age White_Nonwhite Qualifications BMIcode BMI DISPLAY=BOTH
  /TABLE Sex [COUNT F40.0, COLPCT.COUNT PCT40.1] + Age [MEAN, STDDEV] + White_Nonwhite [COUNT 
    F40.0, COLPCT.COUNT PCT40.1] + Qualifications [COUNT F40.0, COLPCT.COUNT PCT40.1] + BMIcode [COUNT 
    F40.0, COLPCT.COUNT PCT40.1] + BMI [MEAN, STDDEV]
  /CATEGORIES VARIABLES=Sex White_Nonwhite Qualifications BMIcode ORDER=A KEY=VALUE EMPTY=INCLUDE 
    MISSING=EXCLUDE
  /CRITERIA CILEVEL=95.

****descriptive for all 6 plate groups***

CTABLES
  /VLABELS VARIABLES=Rice_C1_weight_difference Rice_C2_weight_difference Rice_C3_weight_difference 
    Rice_C4_weight_differece Rice_C5_weight_difference Rice_C6_weight_difference DISPLAY=BOTH
  /TABLE Rice_C1_weight_difference [MEAN, STDDEV, MINIMUM, MAXIMUM] + Rice_C2_weight_difference 
    [MEAN, STDDEV, MINIMUM, MAXIMUM] + Rice_C3_weight_difference [MEAN, STDDEV, MINIMUM, MAXIMUM] + 
    Rice_C4_weight_differece [MEAN, STDDEV, MINIMUM, MAXIMUM] + Rice_C5_weight_difference [MEAN, 
    STDDEV, MINIMUM, MAXIMUM] + Rice_C6_weight_difference [MEAN, STDDEV, MINIMUM, MAXIMUM]
  /CRITERIA CILEVEL=95.

****descriptive for all 6 wine groups***

CTABLES
  /VLABELS VARIABLES=Wine_C1_weight_difference_ML Wine_C2_weight_difference_ML 
    Wine_C3_weight_difference_ML Wine_C4_weight_difference_ML Wine_C5_weight_difference_ML 
    Wine_C6_weight_difference_ML DISPLAY=BOTH
  /TABLE Wine_C1_weight_difference_ML [MEAN, STDDEV] + Wine_C2_weight_difference_ML [MEAN, STDDEV] 
    + Wine_C3_weight_difference_ML [MEAN, STDDEV] + Wine_C4_weight_difference_ML [MEAN, STDDEV] + 
    Wine_C5_weight_difference_ML [MEAN, STDDEV] + Wine_C6_weight_difference_ML [MEAN, STDDEV]
  /CRITERIA CILEVEL=95.

***repeated measures plate***

GLM Rice_C1_weight_difference Rice_C2_weight_difference Rice_C3_weight_difference 
    Rice_C4_weight_differece Rice_C5_weight_difference Rice_C6_weight_difference
  /WSFACTOR=Plate 6 Polynomial 
  /METHOD=SSTYPE(3)
  /PLOT=PROFILE(Plate)
  /EMMEANS=TABLES(Plate) COMPARE ADJ(LSD)
  /PRINT=DESCRIPTIVE ETASQ PARAMETER 
  /CRITERIA=ALPHA(.05)
  /WSDESIGN=Plate.

***repeated measures wine***

GLM Wine_C1_weight_difference_ML Wine_C2_weight_difference_ML Wine_C3_weight_difference_ML 
    Wine_C4_weight_difference_ML Wine_C5_weight_difference_ML Wine_C6_weight_difference_ML
  /WSFACTOR=wine 6 Polynomial 
  /METHOD=SSTYPE(3)
  /PRINT=DESCRIPTIVE ETASQ PARAMETER 
  /CRITERIA=ALPHA(.05)
  /WSDESIGN=wine.

***************************************************************************************************************long format database*******************************

***to check descriptive in long format***

CTABLES
  /VLABELS VARIABLES=Plate food_served_outcome DISPLAY=BOTH
  /TABLE Plate BY food_served_outcome [MEAN, STDDEV]
  /CATEGORIES VARIABLES=Plate ORDER=A KEY=VALUE EMPTY=INCLUDE
  /CRITERIA CILEVEL=95.

CTABLES
  /VLABELS VARIABLES=Wine wine_served_outcome DISPLAY=BOTH
  /TABLE Wine BY wine_served_outcome [MEAN, STDDEV]
  /CATEGORIES VARIABLES=Wine ORDER=A KEY=VALUE EMPTY=INCLUDE
  /CRITERIA CILEVEL=95.

***descriptive for primary analysis***

CTABLES
  /VLABELS VARIABLES=Plate_size Plate_shape food_served_outcome DISPLAY=BOTH
  /TABLE Plate_size + Plate_shape BY food_served_outcome [MEAN, STDDEV]
  /CATEGORIES VARIABLES=Plate_size Plate_shape ORDER=A KEY=VALUE EMPTY=INCLUDE
  /CRITERIA CILEVEL=95.

CTABLES
  /VLABELS VARIABLES=Wine_glass Wine_bottle wine_served_outcome DISPLAY=BOTH
  /TABLE Wine_glass + Wine_bottle BY wine_served_outcome [MEAN, STDDEV]
  /CATEGORIES VARIABLES=Wine_glass Wine_bottle ORDER=A KEY=VALUE EMPTY=INCLUDE
  /CRITERIA CILEVEL=95.

********************data analysis - primary outcome - plate study************************************

***lots of models were run: diagonal, Compound Symmetry: correlation metric, Heterogeneous compound symmetry, GENLIN, ect. But Unstructured showed the better diagnostics and so is reported in the manuscript (confirmed by RM and MP August 2020).

***plate mixed model (Unstructured) *** Interaction included***

MIXED food_served_outcome BY Plate_size Plate_shape
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Plate_size Plate_shape Plate_size*Plate_shape | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Plate_size*Plate_shape | SUBJECT(Participant_ID) COVTYPE(UN)
  /SAVE=RESID
  /EMMEANS=TABLES(Plate_size) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Plate_shape) COMPARE ADJ(LSD).

***plate mixed model (Unstructured) *** Interaction NOT included***REPORTED in paper as best model**

MIXED food_served_outcome BY Plate_size Plate_shape
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Plate_size Plate_shape | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Plate_size*Plate_shape | SUBJECT(Participant_ID) COVTYPE(UN)
  /SAVE=RESID
  /EMMEANS=TABLES(Plate_size) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Plate_shape) COMPARE ADJ(LSD).

***extra analysis where plate size is a covariate (not a fixed factor)

MIXED food_served_outcome BY Plate_shape WITH Plate_size
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Plate_shape Plate_size | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Plate_size*Plate_shape | SUBJECT(Participant_ID) COVTYPE(UN)
  /SAVE=RESID
  /EMMEANS=TABLES(Plate_shape) COMPARE ADJ(LSD).

***extra analysis where plate area is a covariate 

MIXED food_served_outcome BY Plate_shape WITH Plate_area
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Plate_shape Plate_area | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Plate_area*Plate_shape | SUBJECT(Participant_ID) COVTYPE(UN)
  /SAVE=RESID
  /EMMEANS=TABLES(Plate_shape) COMPARE ADJ(LSD).


***************************************************************************************************************data analysis - primary outcome - wine study***************************************************************

***lots of models were run: diagonal, Compound Symmetry: correlation metric, Heterogeneous compound symmetry, GENLIN, ect. But Unstructured showed the better diagnostics and so is reported in the manuscript (confirmed by RM and MP August 2020).

***wine mixed model 3a (Unstructured) *** Interaction included***

MIXED wine_served_outcome BY Wine_bottle Wine_glass
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Wine_bottle Wine_glass Wine_bottle*Wine_glass | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Wine_bottle*Wine_glass | SUBJECT(Participant_ID) COVTYPE(UN)
  /SAVE=RESID
  /EMMEANS=TABLES(Wine_bottle) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Wine_glass) COMPARE ADJ(LSD).

***wine mixed model 3a (Unstructured) *** Interaction NOT included*** REPORTED in paper as best model**

MIXED wine_served_outcome BY Wine_bottle Wine_glass
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Wine_bottle Wine_glass | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Wine_bottle*Wine_glass | SUBJECT(Participant_ID) COVTYPE(UN)
  /SAVE=RESID
  /EMMEANS=TABLES(Wine_bottle) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Wine_glass) COMPARE ADJ(LSD).

***extra analysis where glass size is a covariate (not a fixed factor)

MIXED wine_served_outcome BY Wine_bottle WITH Wine_glass
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Wine_bottle Wine_glass | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Wine_bottle*Wine_glass | SUBJECT(Participant_ID) COVTYPE(UN)
  /EMMEANS=TABLES(Wine_bottle) COMPARE ADJ(LSD).

***extra analysis where glass volume is a covariate

MIXED wine_served_outcome BY Wine_bottle WITH Glass_volume
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Wine_bottle Glass_volume | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Wine_bottle*Glass_volume | SUBJECT(Participant_ID) COVTYPE(UN)
  /EMMEANS=TABLES(Wine_bottle) COMPARE ADJ(LSD).

*******************************************************primary outcome models for both studies with adjustments for gender and education***

MIXED food_served_outcome BY Plate_size Plate_shape Sex Education
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Plate_size Plate_shape Sex Education | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Plate_size*Plate_shape | SUBJECT(Participant_ID) COVTYPE(UN)
  /EMMEANS=TABLES(Plate_size) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Plate_shape) COMPARE ADJ(LSD).

MIXED wine_served_outcome BY Wine_bottle Wine_glass Sex Education
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Wine_bottle Wine_glass Sex Education | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Wine_bottle*Wine_glass | SUBJECT(Participant_ID) COVTYPE(UN)
  /EMMEANS=TABLES(Wine_bottle) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Wine_glass) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Sex) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Education) COMPARE ADJ(LSD).


**************************************************trying with interactions (for plate only)***

MIXED food_served_outcome BY Plate_size Plate_shape Sex Education
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Plate_size Plate_shape Sex Education Plate_size*Sex Plate_shape*Sex Plate_shape*Education 
    Plate_size*Education | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Plate_size*Plate_shape | SUBJECT(Participant_ID) COVTYPE(UN)
  /EMMEANS=TABLES(Plate_size) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Plate_shape) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Sex) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Education) COMPARE ADJ(LSD).

***********************************************primary outcome models for both studies with adjustments for study order***

MIXED food_served_outcome BY Plate_size Plate_shape Sequence_rice_or_wine_first
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Plate_size Plate_shape Sequence_rice_or_wine_first | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Plate_size*Plate_shape | SUBJECT(Participant_ID) COVTYPE(UN)
  /SAVE=RESID
  /EMMEANS=TABLES(Plate_size) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Plate_shape) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Sequence_rice_or_wine_first) COMPARE ADJ(LSD).

MIXED wine_served_outcome BY Wine_bottle Wine_glass Sequence_rice_or_wine_first
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Wine_bottle Wine_glass Sequence_rice_or_wine_first | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Wine_bottle*Wine_glass | SUBJECT(Participant_ID) COVTYPE(UN)
  /SAVE=RESID
  /EMMEANS=TABLES(Wine_bottle) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Wine_glass) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Sequence_rice_or_wine_first) COMPARE ADJ(LSD).

*************************************primary outcome model with adjustments for plate sequence***

MIXED food_served_outcome BY Plate_size Plate_shape time_rice
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Plate_size Plate_shape time_rice | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Plate_size*Plate_shape | SUBJECT(Participant_ID) COVTYPE(UN)
  /SAVE=RESID
  /EMMEANS=TABLES(Plate_size) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Plate_shape) COMPARE ADJ(LSD)
 /EMMEANS=TABLES(time_rice) COMPARE ADJ(LSD).

*******************************************primary outcome model with adjustments for glass sequence***

MIXED wine_served_outcome BY Wine_bottle Wine_glass time_wine
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Wine_bottle Wine_glass time_wine | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Wine_bottle*Wine_glass | SUBJECT(Participant_ID) COVTYPE(UN)
  /SAVE=RESID
  /EMMEANS=TABLES(Wine_bottle) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Wine_glass) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(time_wine) COMPARE ADJ(LSD).


********************************************Sensitivity analysis 1 (using primary model) *********************************

TEMPORARY.
SELECT IF Studyaim_guessed = 2.
MIXED food_served_outcome BY Plate_size Plate_shape
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Plate_size Plate_shape | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Plate_size*Plate_shape | SUBJECT(Participant_ID) COVTYPE(UN)
  /SAVE=RESID
  /EMMEANS=TABLES(Plate_size) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Plate_shape) COMPARE ADJ(LSD).

TEMPORARY.
SELECT IF Studyaim_guessed = 2.
MIXED wine_served_outcome BY Wine_bottle Wine_glass
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Wine_bottle Wine_glass | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Wine_bottle*Wine_glass | SUBJECT(Participant_ID) COVTYPE(UN)
  /SAVE=RESID
  /EMMEANS=TABLES(Wine_bottle) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Wine_glass) COMPARE ADJ(LSD).

********************************************Sensitivity analysis 2 (using primary model)*********************************

TEMPORARY.
SELECT IF Sensitivty_analysis2 = 1.
MIXED food_served_outcome BY Plate_size Plate_shape
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Plate_size Plate_shape | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Plate_size*Plate_shape | SUBJECT(Participant_ID) COVTYPE(UN)
  /SAVE=RESID
  /EMMEANS=TABLES(Plate_size) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Plate_shape) COMPARE ADJ(LSD).

TEMPORARY.
SELECT IF Sensitivty_analysis2 = 1.
MIXED wine_served_outcome BY Wine_bottle Wine_glass
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Wine_bottle Wine_glass | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Wine_bottle*Wine_glass | SUBJECT(Participant_ID) COVTYPE(UN)
  /SAVE=RESID
  /EMMEANS=TABLES(Wine_bottle) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Wine_glass) COMPARE ADJ(LSD).

***************************************analysis excluding 1 outlier (ID 98) - plate only**************************

TEMPORARY.
SELECT IF Participant_ID NE 98.
MIXED food_served_outcome BY Plate_size Plate_shape
  /CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, 
    ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001, ABSOLUTE)
  /FIXED=Plate_size Plate_shape | SSTYPE(3)
  /METHOD=REML
  /PRINT=SOLUTION
  /REPEATED=Plate_size*Plate_shape | SUBJECT(Participant_ID) COVTYPE(UN)
  /SAVE=RESID
  /EMMEANS=TABLES(Plate_size) COMPARE ADJ(LSD)
  /EMMEANS=TABLES(Plate_shape) COMPARE ADJ(LSD).