Making use of car::Anova seems fine. To handle power of continuous variables (keeping track of terms), parts of this code need to be rewritten. The problem is the use of attr(terms(mod), “factors”). In ffmanova this matrix has be reformulated by the use of ffmanova::fixModelMatrix (see the documentation). (Unfortunately ModelMatrix has another meaning that model.matrix because the code was originally written in Matlab.)