<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">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.<br>
<br><span style="font-family:courier new,monospace">## Matriz que representa uma vizinhança em forma de grid. O valor na<br>## cédula é o número de eventos no setor.<br>x <- matrix(rpois(25, l=10), 5, 5)<br><br>## Vizinhança de nível. As coordenadas dos vizinhos do setor [2,2].<br>
viz <- 1<br>i <- 2; j <- 2<br>outer((i-viz):(i+viz), (j-viz):(j+viz), paste)<br><br>## Matriz acrescida de borda com NA para evitar problema ao passar pelos<br>## setores de borda.<br>xB <- cbind(NA, rbind(NA, x, NA), NA)<br>
<br>## Aplicando a soma.<br>sum(xB[-viz:viz+i, -viz:viz+j], na.rm=TRUE)<br><br>## Aplicando à todos os índices.<br>S <- x<br>for(i in 1:nrow(x)){<br>    for(j in 1:ncol(x)){<br>        S[i,j] <- sum(xB[-viz:viz+i, -viz:viz+j], na.rm=TRUE)<br>
    }<br>}<br><br>## Resultado.<br>S<br></span><br>À disposição.<br>Walmes.<br></div></div>