GAMMs: Implementing a 3-way interaction with categorical variables (i.e. 2 x 2 x 2) using binary difference smooths

Question

I am attempting to construct a model which conceptually includes an interaction of three 2-level categorical factors, following the guidance of Wieling’s 2018 GAMMs tutorial