Population inference: How to adjust for covariates

In case you need to adjust for covariates in the population model:

from fmristats import *

model = PopulationModel(sample, 'status + smoking + sex + cohort')
result = model.fit()