Como exibir o número de grupos em tabelas de regressão fixest usando modelsummary?

Estou usando o fixestpacote para regressões de efeitos fixos e modelsummarypara criar tabelas. Gostaria que as tabelas mostrassem o número de grupos (por exemplo, número de valores únicos para as variáveis ​​de efeito fixo). Aqui está um exemplo mínimo:

library(fixest)
library(modelsummary)

# Example fixed effects regression
reg1 <- feols(Sepal.Length ~
                Sepal.Width + Petal.Length + Petal.Width |
                Species,
              data = iris
)

# Create summary table
msummary(reg1)

Isso produz uma tabela legal, mas eu gostaria de adicionar o número de grupos no meu efeito fixo (neste exemplo, seriam 3 espécies). Eu tentei usar gof_mappara incluir essas informações, mas não parece extrair o número de grupos por padrão:

msummary(reg1, 
         gof_map = c("nobs", "r.squared", "adj.r.squared", "n"))

Existe uma maneira de fazer modelsummaryexibir o número de grupos para cada efeito fixo no meu modelo? Eu preferiria uma solução que funcionasse com a modelsummaryfuncionalidade padrão, mas estou aberto a funções personalizadas, se necessário.

Observação : na minha pesquisa atual, estou trabalhando com dados em painel, onde o número de grupos é importante para interpretar os resultados.