Categorical covariate with two way anova in R

Question

This is a random sample of my data set:

structure(list(DTI_ID = structure(c(31L, 241L, 84L, 298L, 185L, 
269L, 198L, 24L, 286L, 177L, 228L, 158L, 57L, 293L, 2