Outliers in the gene expression data were detected using Grubbs test. Treatment by gender interactions for single genes were tested using two-way analysis of variance models (ANOVAs). Comparisons of mRNA expression levels within each sex to assess direction of change in males and females were performed by Welch two-sample t-tests. Pearson correlations were calculated to assess and visualize the relationship between expression levels of genes using multivariate analyses in the JMP 12 software (SAS Institute, Cary, NC). Open field locomotor activity experiments were analyzed using ANOVAs with repeated measures using a multivariate model followed by post-hoc pairwise t-test comparisons when appropriate to address the effect of THC at specific time points.
