COMPUTE position_feature = word_pos * 10 + sequence_length. COMPUTE form_feature=stress*10+syllables. EXECUTE. Variable Level position_feature (Nominal). Variable Level form_feature (Nominal). DESCRIPTIVES VARIABLES=pitchv_PC1 pitchv_PC2 duration /STATISTICS=MEAN STDDEV VARIANCE MIN MAX KURTOSIS SKEWNESS. MEANS TABLES=pitchv_PC1 pitchv_PC2 duration BY Sex /CELLS MEAN COUNT STDDEV /STATISTICS ANOVA LINEARITY. # The GLM GLM ampv_PC1 ampv_PC2 pitchv_PC1 pitchv_PC2 duration BY position_feature form_feature Sex Dialect WITH sentence_pos /METHOD=SSTYPE(3) /INTERCEPT=INCLUDE /EMMEANS=TABLES(OVERALL) WITH(sentence_pos=MEAN) /EMMEANS=TABLES(Sex) WITH(sentence_pos=MEAN) /EMMEANS=TABLES(Dialect) WITH(sentence_pos=MEAN) /PRINT=DESCRIPTIVE ETASQ OPOWER /CRITERIA=ALPHA(.05) /DESIGN=position_feature form_feature Sex Dialect sentence_pos. # The same GLM but using Participant instead of Sex and Dialect GLM ampv_PC1 ampv_PC2 pitchv_PC1 pitchv_PC2 duration BY position_feature form_feature Participant WITH sentence_pos /METHOD=SSTYPE(3) /INTERCEPT=INCLUDE /EMMEANS=TABLES(OVERALL) WITH(sentence_pos=MEAN) /EMMEANS=TABLES(Participant) WITH(sentence_pos=MEAN) /PRINT=DESCRIPTIVE ETASQ OPOWER /CRITERIA=ALPHA(.05) /DESIGN=position_feature form_feature sentence_pos participant. # Now without position_feature GLM ampv_PC1 ampv_PC2 pitchv_PC1 pitchv_PC2 duration BY form_feature word_pos sequence_length Sex Dialect WITH sentence_pos /METHOD=SSTYPE(3) /INTERCEPT=INCLUDE /PLOT=PROFILE(word_pos*sequence_length) /EMMEANS=TABLES(OVERALL) WITH(sentence_pos=MEAN) /EMMEANS=TABLES(word_pos*sequence_length) WITH(sentence_pos=MEAN) /PRINT=DESCRIPTIVE ETASQ OPOWER /CRITERIA=ALPHA(.05) /DESIGN=form_feature word_pos sequence_length Sex Dialect sentence_pos. # Now without form_feature GLM ampv_PC1 ampv_PC2 pitchv_PC1 pitchv_PC2 duration BY stress syllables position_feature Sex Dialect WITH sentence_pos /METHOD=SSTYPE(3) /INTERCEPT=INCLUDE /PLOT=PROFILE(stress*syllables) /EMMEANS=TABLES(OVERALL) WITH(sentence_pos=MEAN) /EMMEANS=TABLES(stress) WITH(sentence_pos=MEAN) /EMMEANS=TABLES(syllables) WITH(sentence_pos=MEAN) /PRINT=DESCRIPTIVE ETASQ OPOWER /CRITERIA=ALPHA(.05) /DESIGN=position_feature stress syllables Sex Dialect sentence_pos. # Graph PC1 intonation GGRAPH /GRAPHDATASET NAME="graphdataset" VARIABLES=word_pos MEAN(pitchv_PC1)[name="MEAN_pitchv_PC1"] sequence_length MISSING=LISTWISE REPORTMISSING=NO /GRAPHSPEC SOURCE=INLINE. BEGIN GPL SOURCE: s=userSource(id("graphdataset")) DATA: word_pos=col(source(s), name("word_pos"), unit.category()) DATA: MEAN_pitchv_PC1=col(source(s), name("MEAN_pitchv_PC1")) DATA: sequence_length=col(source(s), name("sequence_length"), unit.category()) GUIDE: axis(dim(1), label("word_pos")) GUIDE: axis(dim(2), label("Mean pitchv_PC1")) GUIDE: legend(aesthetic(aesthetic.color.interior), label("sequence_length")) SCALE: linear(dim(2), include(0)) ELEMENT: line(position(word_pos*MEAN_pitchv_PC1), color.interior(sequence_length), missing.wings()) END GPL. # Graph PC2 intonation GGRAPH /GRAPHDATASET NAME="graphdataset" VARIABLES=word_pos MEAN(pitchv_PC2)[name="MEAN_pitchv_PC2"] sequence_length MISSING=LISTWISE REPORTMISSING=NO /GRAPHSPEC SOURCE=INLINE. BEGIN GPL SOURCE: s=userSource(id("graphdataset")) DATA: word_pos=col(source(s), name("word_pos"), unit.category()) DATA: MEAN_pitchv_PC2=col(source(s), name("MEAN_pitchv_PC2")) DATA: sequence_length=col(source(s), name("sequence_length"), unit.category()) GUIDE: axis(dim(1), label("word_pos")) GUIDE: axis(dim(2), label("Mean pitchv_PC2")) GUIDE: legend(aesthetic(aesthetic.color.interior), label("sequence_length")) SCALE: linear(dim(2), include(0)) ELEMENT: line(position(word_pos*MEAN_pitchv_PC2), color.interior(sequence_length), missing.wings()) END GPL. # Checking for correlation between pitch and amplitude PARTIAL CORR /VARIABLES=ampv_PC1 ampv_PC2 pitchv_PC1 pitchv_PC2 BY position_feature form_feature sentence_pos /SIGNIFICANCE=TWOTAIL /STATISTICS=DESCRIPTIVES /MISSING=LISTWISE. # GLM for the influence of the control word IF (word_name <> 'Venlo') Control=0. IF (word_name = 'Venlo') Control=1. EXECUTE. COMPUTE filter_$=(stress = 1 AND syllables = 2). VARIABLE LABELS filter_$ 'stress = 1 AND syllables = 2 (FILTER)'. VALUE LABELS filter_$ 0 'Not Selected' 1 'Selected'. FORMATS filter_$ (f1.0). FILTER BY filter_$. EXECUTE. GLM ampv_PC1 ampv_PC2 pitchv_PC1 pitchv_PC2 duration BY position_feature Sex Dialect Control WITH sentence_pos /METHOD=SSTYPE(3) /INTERCEPT=INCLUDE /EMMEANS=TABLES(OVERALL) WITH(sentence_pos=MEAN) /EMMEANS=TABLES(Control) WITH(sentence_pos=MEAN) /PRINT=DESCRIPTIVE ETASQ OPOWER /CRITERIA=ALPHA(.05) /DESIGN=position_feature Sex Dialect sentence_pos Control. # T-test of component scores of starting words against sequence length USE ALL. COMPUTE filter_$=(word_pos = 1). VARIABLE LABELS filter_$ 'word_pos = 1 (FILTER)'. VALUE LABELS filter_$ 0 'Not Selected' 1 'Selected'. FORMATS filter_$ (f1.0). FILTER BY filter_$. EXECUTE. T-TEST GROUPS=sequence_length(2.5) /MISSING=ANALYSIS /VARIABLES=pitchv_PC1 pitchv_PC2 duration /CRITERIA=CI(.95).