Hi,
When we perform an Anova (for example), we can check normality by two ways: a numerical way via a S-W test and a graphical way via a QQ-plot of the residuals.
For homoscedasticity (who invented that term?), we only have a numerical test with a Levene's test. Would you consider adding a graphical way such as an S-L plot? See here for R details about that --> https://mgimond.github.io/ES218/Week07b.html
Happy holidays to everyone and a special thought to our Australian friends who are cooked in a continent-sized furnace
Keep up the good work!