Statistical analysis
psychotic-like reactions
According to our power calculation, at 80% power and Bonferroni adjusted alpha <0.008, a sample size of The effect of THC was determined by comparing outcome scores from the baseline visit with those following administration with THC alone (0:1) using paired t-tests. For the primary analysis, we used linear mixed models to assess the effect of varying the CBD:THC ratio on delayed recall on the HVLT-R. The four CBD:THC ratios (0:1, 1:1, 2:1, 3:1) were included as a fixed effect, with participant as a random effect to account for the dependency between repeated measures. All 6 contrasts were of interest (0:1 vs 1:1, 0:1 vs 2:1, 0:1 vs 3:1, 1:1 vs 2:1, 1:1 vs 3:1, 2:1 vs 3:1) and alpha was set according to the results of our power calculation at For pharmacokinetics, VAS scores and physiological outcomes, both peak effects (0 min for pharmacokinetic and physiological outcomes) and area under the curve (AUC) were investigated. For the AUC analyses, values were baseline corrected before using the spline method using the bayestestR package (version 0.7.5.1) [The relationships between both inhalation time and coughing with peak plasma THC and CBD, in addition with their respective AUCs, were assessed using Pearson’s correlation coefficients.We additionally categorised clinically significant psychotic-like reactions as increases in PANSS scores from baseline of ≥3 points, as in previous studies due to floor effects [Multiple imputation chain equations (MICE) were used to impute missing values in pharmacokinetic, cognitive, pleasurable, and physiological outcomes using the mice package (version 3.13.0) [All analyses were conducted using R version 3.5.3. lme4 (version 1.1-26) [