1 Modelos

O modelo beta-binomial inclui um parâmetro de agregação, porém este parâmetro não apresenta dependência espacial e não pode ser utilizado para simulação. O modelo geoestatístico modela a dependência espacial por um processo estocástico subjacente e pode ser utilizado para simular dados com dependência espacial. Também há modelos para dados de incidência de doenças de plantas baseados na pressão infectiva exercida por uma planta doente nas plantas sadias (?) e (?).

Para simular dados de incidência de plantas, foi adaptado a função grf do pacote geoR (?) e foram implementados dois modelos de pressão infectiva baseada em distância. O resultado da função grf é um valor contínuo e para obter a incidência é feita uma transformação logito. Os modelos de pressão infectiva implementados tem dois parâmetros: α1, que controla a incidência, e α2 que controla a agregação. Em um modelo a pressão infectiva é proporcional a uma potência da distância α1 + d-2*α2 e no outro, é proporcional a um exponencial da distância α1 + exp (-d * α2).

Na Figura 1, pode-se vizualizar a pressão infectiva em função da distância para alguns parâmetros de ambos os modelos, e a correlação segundo uma função de correlação exponencial, a mais utilizada no modelo geoestatístico.


  > par(mfrow = c(1, 1), mar = c(3, 3, 3, 0.3), mgp = c(2, 1,
  +     0))
  > plot(function(d) 0 + exp(-0.5 * d), 0, 20, lty = 3, ylab = "",
  +     xlab = "Distância")
  > plot(function(d) 0 + exp(-1 * d), 0, 20, ad = T, col = 2,
  +     lty = 3)
  > plot(function(d) 0 + exp(-2 * d), 0, 20, ad = T, col = 4,
  +     lty = 3)
  > plot(function(d) 0 + d^(-2 * 0.5), 0, 20, ad = T, col = 1,
  +     lty = 2)
  > plot(function(d) 0 + d^(-2 * 1), 0, 20, ad = T, col = 2,
  +     lty = 2)
  > plot(function(d) 0 + d^(-2 * 2), 0, 20, ad = T, col = 4,
  +     lty = 2)
  > plot(function(d) exp(-d/10), 0, 20, ad = T, col = 1)
  > plot(function(d) exp(-d/5), 0, 20, ad = T, col = 2)
  > plot(function(d) exp(-d/3), 0, 20, ad = T, col = 3)
  > legend(10, 1, c(expression(paste("corr exponencial, ", phi,
  +     "  = 10")), expression(paste("corr exponencial, ", phi,
  +     "  =   5")), expression(paste("corr exponencial, ", phi,
  +     "  =   3")), expression(paste("pres potencia, ", alpha,
  +     "2  = 0.5")), expression(paste("pres potencia, ", alpha,
  +     "2  = 1")), expression(paste("pres potencia, ", alpha,
  +     "2  = 2")), expression(paste("pres exponencial, ", alpha,
  +     "2  = 0.5")), expression(paste("pres exponencial, ",
  +     alpha, "2  = 1")), expression(paste("pres exponencial, ",
  +     alpha, "2  = 2")), ), lty = rep(1:3, each = 3), col = 1:3)


pict

Figura 1: Valores de correlação e pressão infectiva para diferentes modelos e parâmetros de agregação.


Nota-se que o modelo de pressão infectiva potência, a pressão infectiva decai mais rapidamente em função da distância, que no modelo com pressão infectiva exponencial. Nota-se também que a função de correlação do modelo geoestatístico é equivalente à pressão infectiva exponecial.