Eu acredito que deva ter um comando que faça isso porque é uma espécie de média móvel cuja janela é a vizinhança. Aqui vai minha sugestão.

## Matriz que representa uma vizinhança em forma de grid. O valor na
## cédula é o número de eventos no setor.
x <- matrix(rpois(25, l=10), 5, 5)

## Vizinhança de nível. As coordenadas dos vizinhos do setor [2,2].
viz <- 1
i <- 2; j <- 2
outer((i-viz):(i+viz), (j-viz):(j+viz), paste)

## Matriz acrescida de borda com NA para evitar problema ao passar pelos
## setores de borda.
xB <- cbind(NA, rbind(NA, x, NA), NA)

## Aplicando a soma.
sum(xB[-viz:viz+i, -viz:viz+j], na.rm=TRUE)

## Aplicando à todos os índices.
S <- x
for(i in 1:nrow(x)){
    for(j in 1:ncol(x)){
        S[i,j] <- sum(xB[-viz:viz+i, -viz:viz+j], na.rm=TRUE)
    }
}

## Resultado.
S

À disposição.
Walmes.