From Karacayli et al. (2023) and from Ravoux et al. (2023) The sysetmatics matrix might need some work, since it is not diagonal