Trial summary
The clinical trial was conducted in five primary care facilities in the city of Santiago – Chile and compared two interventions in people with risky alcohol use (i.e., Alcohol Use Disorders Identification Test (AUDIT) total score between 8 and 15)[7]. One group received an informative pamphlet (n = 168), and the other received the same pamphlet plus a HT-delivered BI based on the contents from the pamphlet (n = 174). The BI was based on motivational interview and is similar in contents and duration to a standard intervention for risky alcohol use.
Our outcomes were the AUDIT risk category (primary outcome), and the AUDIT-total and AUDIT-C scores (i.e., sum of items one to three), measured 6 months after enrollment. For practical reasons, the only outcome measurement was a one-year retrospective AUDIT. This minimal outcome collection allowed protocol implementation without changing much the way the program operates in the real world.
Previous analysis
All analysis were conducted on participants that completed follow-up (n = 294) using Chi-squared and T-tests accordingly. Then, mixed-effects regressions were used to adjust for demographic variables as fixed-effects (age, sex, and educational level) and the health center as random effects. See the main paper for a detailed description of the trial[4].
Current analysis
We used Bayesian inference to estimate the same mixed-effects regression models specified in the original analysis. The alcohol risk status (i.e., primary outcome), was modeled with a Bernoulli distribution with fixed effects for all the covariates, and random intercepts for health centers. The model for the primary outcome is presented in Eq. 1.
Risky drinking ~ Bernoulli ( q )
\(\text{log}\left(\frac{q}{1-q}\right)\) = \(\beta\)1 + \(\beta\)2 Group + \(\beta\)3 Sex + \(\beta\)4 Education(complete) + \(\beta\)5 Education(superior) + \(\beta\)6 Age + C Center
\(\beta\) [1−6] ~ Normal ( 0, 1 )
C ~ Normal ( 0, sigma_c )
sigma_c ~ Normal ( 0, 1 )
AUDIT-total and AUDIT-C scores (i.e., secondary outcomes), were modeled with normal distributions with fixed effects for all covariates and random intercepts for health centers. The model for the secondary outcomes is presented in Eq. 2.
AUDIT-score ~ Normal ( mu, sigma )
mu = \(\beta\)1 + \(\beta\)2 Group + \(\beta\)3 Sex + \(\beta\)4 Education(complete) + \(\beta\)5 Education(superior) + \(\beta\)6 Age + \(\beta\)7 AUDIT-score(baseline) + C Center
\(\beta\) [1−7] ~ Normal ( 0, 1 )
sigma = Normal ( 0, 1 )
C ~ Normal ( 0, sigma_c )
sigma_c ~ Normal ( 0, 1 )
Hamiltonian Markov chain Monte Carlo was used for Bayesian inference. For each model, 50,000 iterations were run with 25,000 warm-up iterations in four chains. The analyses were run using R with R-Stan version 2.21.0 on a Mac mini M1.