How do I change colors in a visualization of a two-factor ANOVA interaction?

head(ToothGrowth)
library(HH)
attach(ToothGrowth)
AOV<-aov(len ~ suppdose)
summary(AOV)
interaction_result <-interaction2wt(len ~ supp
dose,col=c("red","blue","black"))


I tried col() but could only change the color of the box diagram, I wanted to change the color of the upper left and lower right interaction lines