2 Número médio de vizinhos doentes

Na análise do número de vizinhos próximos, pode-se usar a função neigh.test. Nesta função deve-se entrar com os dados em formato geodata, o número de simulações, o(s) código(s) que identifica(m) plantas doentes e o raio de vizinhança. Na Figura 2 vizualiza-se os resultados.


  > args(neigh.test)


  function (obj, NMC = 99, death = 1, ray = 30, evaluation = 1,
      edge.cor = TRUE)
  NULL


  > neit <- neigh.test(v303.geo, NMC = 199, death = 1:3, evaluation = 1:3)


  test evaluation: 1 2 3


  > summary(neit)


  Results for 199 Monte Carlo simulations!
  Observed: 3.61 4.12 10.1
  Randoms:
            Av1  Av2  Av3
  Min.    0.983 1.34 5.73
  1st Qu. 1.770 2.11 6.47
  Median  2.060 2.39 6.77
  Mean    2.100 2.43 6.79
  3rd Qu. 2.410 2.70 7.07
  Max.    3.610 4.14 8.48
  P-value: 0.01 0.01 0.005


  > par(mfrow = c(2, 3), mar = c(2, 2, 2, 0.1), mgp = c(1.2,
  +     0.5, 0))
  > hist(neit, main = "Número de Vizinhos", evaluation = 1:3)
  > plot(neit, main = "Número de Vizinhos", evaluation = 1:3)


pict

Figura 2: Vizualização dos resultados do teste de Monte Carlo para o número médio de vizinhos doentes.