Caro João, boa noite!
Sem os dados pra rodar o código fica difícil, mas arrisco um palpite:
Não está completo. Tente acrescentar ou alterar no seu código...
colors <-brewer.pal(9, "YlOrRd") ### a paleta com nove classes que você criou
colors; display.brewer.pal(9, "YlOrRd") ### visualizando a paleta
colors[10]="#D3D3D3" ### adiciona uma décima cor (cinza) na paleta (para os NA's)
### não passe findInterval() diretamente, aloje em uma nova variável
tmp <- findInterval(wld.aux$ValueMi, brks, all.inside=T) ### vai distribuir valores de 1:9 de acordo com seus 'brks' calculados
### se você listar 'tmp' deverá ver a atribuição de cores (1:9) e NA's
head(tmp); table(tmp)
### troca os NA's por 10 (nova cor da paleta)
tmp[which(is.na(tmp))] = 10
### muda o comando original usando 'tmp' no lugar de findInterval()
plot(wld.aux, col=colors[tmp], axes=F, bg=NA, mai=c(0,2,2,20))