Is it possible to get an overall measure of variance explained by different levels of a categorical predictor, without using a reference level? Let's take a model like this:
mod <- lm(Petal.Width ~ Species, data = iris)

Is there any way to know how much variance is explained by each level of Species overall? By overall I mean NOT just relative to a reference level.
I know that emmeans can compare coefficients of each level of a categorical predictor to the overall mean. Is something like this possible for amount of variance explained?

The difficulty here is that a level of a categrical variable is not a variable --- it is an outcome of the variable.  So this is a bit like asking if it is possible to get a measure of variance explained for a continuous variable for the outcome $X=3$.
While it is an odd question, I suppose that the answer is yes, it is possible (though it would no longer be the same model).  What you would do is to remove the initial categorical/continuous variable and replace it with a single indicator variable for the outcome of interest.  You can then use ANOVA to measure the variance explained by this indicator variable.  This would still compare the outcome to a reference level, but the reference level would be all cases where that outcome does not occur.  Of course, this is not the same model as you are using when you include the full categorical variable (which I'm assuming has ore than two outcomes), but that would be the appropriate comparison.
