R - problema con group_by

group_by me agrupa los datos de una tabla solo dándome un resultado al hacer la mediana, aunque yo lo quiero ordenado por provincias:

library(tidyverse)
library(pxR)
library(knitr)

library(dplyr)

censo2021 <- pxR::read.px("https://www.ine.es/jaxi/files/tpx/es/px/55243.px") |> as_tibble()
censo2021=censo2021|>select(-1)
censo2021=censo2021|>dplyr::rename(nacimiento=3,municipio=2,pob=value)
censo2021=censo2021|>separate(municipio,c("cod","municipio"),sep=" ",extra="merge")
censo2021=censo2021|>mutate(prov=as.numeric(cod)%/%1000)

censo2021|>filter(Sexo=="Ambos sexos",nacimiento%in%c("TOTAL","España"),!is.na(prov))|>
pivot_wider(names_from = nacimiento, values_from = pob)|>mutate(extranj=TOTAL-España)|>
group_by(prov)|>summarise(mediana=median(extranj))|>arrange(-mediana)