* Encoding: UTF-8.
*  ======================================================================
  
*  *Greaves, L. M., Barlow, F. K., *Lee, C. H. J., *Matika, C. M., 
*  *Wang, Q., *Lindsay, C-J., *Case, C. J. B., *Sengupta, N. K., *Huang, Y., 
*  *Cowie, L. J., *Stronge, S., *Storey, M., *De Souza, L., *Manuela, S., 
*  *Hammond, M. D., *Milojev, P., *Townrow, C. S., *Muriwai, E., *Satherley, N., 
*  *Fraser, G., West-Newman, T., Houkamau, C., Bulbulia, J., Osborne, D., 
*  Wilson, M. S., & Sibley, C. G. (2017). The diversity and prevalence of 
*  sexual orientation self-labels in a New Zealand national sample. 
*  Archives of Sexual Behavior, 46, 1325-1336.

*  ======================================================================

* Using the NZAVS Time 5 SPSS data file.
* Applying the sample weighting procedure implemented in Sibley (2015) to weight sample by gender, ethnicity, and region.

WEIGHT BY w.GendEthnicRegion.T05.

* creating the sexual orientation variable to analyse with the 6 categories - heterosexual, lesbian/gay, bisexual, bicurious, 
    pansexual/open, asexual.

RECODE SexualOrientationL2.T05 (10=10) (20=20) (30=30) (35=35) (40=40) (50=50) (ELSE=SYSMIS) INTO 
    SexualOrientationL2NoMissing.T05.
EXECUTE.

* analyses including pairwise comparisons for sexual orientation at each level of gender, and gender at each level of sexual orientation.

UNIANOVA Age.T05 BY Gender.T05 SexualOrientationL2NoMissing.T05
  /METHOD=SSTYPE(3)
  /INTERCEPT=INCLUDE
  /POSTHOC=SexualOrientationL2NoMissing.T05(BONFERRONI) 
  /EMMEANS=TABLES(OVERALL) 
  /EMMEANS=TABLES(Gender.T05) 
  /EMMEANS=TABLES(SexualOrientationL2NoMissing.T05) 
  /EMMEANS=TABLES(Gender.T05*SexualOrientationL2NoMissing.T05) COMPARE (Gender.T05) ADJ(LSD)
 /EMMEANS=TABLES(Gender.T05*SexualOrientationL2NoMissing.T05) COMPARE (SexualOrientationL2NoMissing.T05) ADJ(LSD)
 /PRINT=ETASQ DESCRIPTIVE
  /CRITERIA=ALPHA(.05)
  /DESIGN=Gender.T05 SexualOrientationL2NoMissing.T05 Gender.T05*SexualOrientationL2NoMissing.T05.