Bom dia,
Tenho um banco de dados:

    regiao semana notificados confirmados sexo local
1      1      1          82                              45    M     A
2      1      2         105                             32    F     A
3      1      3         134                             71    M     A
4      1      4          80                              44    M     A
5      1      5         102                             32    M     A
6      1      6         108                             46    F     A

Formado por 50 regiões, a variável confirmados apresenta missings (NA).
Preciso substituir os missings pela média de casos confirmados por região.
Calculei a média por região por esse comando, sem considerar os missings:

media <- bancos%>%
  group_by(regiao)%>%
  summarise(mean(na.exclude(confirmados)))
as.data.frame(media)

Com esse resultado:

 regiao mean(na.exclude(confirmados))
1       1                      50.15152
2       2                      51.52604
3       3                      69.90769
4       4                      45.18557
5       5                      61.65816
6       6                      57.79688
7       7                      47.30526
8       8                      69.71795
9       9                      68.40704
10     10                      62.25521
11     11                      58.00000

Só que não estou conseguindo trocar os missings por essas médias e preciso de ajuda nisso.
Desde já agradeço,
Suzana