2 Exemplo de dados simulados

A função sim.citrus() implementa os modelos descritos.


  > args(sim.citrus)


  function (coords, p.death = 0.1, model = c("pdist", "edist",
      "hGRF", "clipGRF", "transGRF"), nsim = 1, ...)
  NULL

É necessário entrar com as coordenadas, a incidência, o modelo e os parâmetros adicionais do modelo escolhido. Na Figura 2 vizualiza-se alguns exemplos de dados simulados, comforme os comandos a seguir:


  > coo <- expand.grid(7 * (1:20), 3 * (1:40))
  > par(mfrow = c(5, 3), mar = c(2, 2, 2, 0.1), mgp = c(1, 0.3,
  +     0))
  > set.seed(123)
  > plot(sim.citrus(coo, p = 0.1, model = "edist", a1 = 0, a2 = 0.1),
  +     main = expression(paste("P. Exponencial, ", alpha, "2 = 0.1")),
  +     pch = 19)
  > set.seed(123)
  > plot(sim.citrus(coo, p = 0.1, model = "edist", a1 = 0, a2 = 0.5),
  +     main = expression(paste("P. Exponencial, ", alpha, "2 = 0.5")),
  +     pch = 19)
  > set.seed(123)
  > plot(sim.citrus(coo, p = 0.1, model = "edist", a1 = 0, a2 = 1.5),
  +     main = expression(paste("P. Exponencial, ", alpha, "2 = 1.5")),
  +     pch = 19)
  > set.seed(123)
  > plot(sim.citrus(coo, p = 0.1, model = "pdist", a1 = 0, a2 = 0.1),
  +     main = expression(paste("P. Potencia, ", alpha, "2 = 0.1")),
  +     pch = 19)
  > set.seed(123)
  > plot(sim.citrus(coo, p = 0.1, model = "pdist", a1 = 0, a2 = 0.5),
  +     main = expression(paste("P. Potencia, ", alpha, "2 = 0.5")),
  +     pch = 19)
  > set.seed(123)
  > plot(sim.citrus(coo, p = 0.1, model = "pdist", a1 = 0, a2 = 1.5),
  +     main = expression(paste("P. Potencia, ", alpha, "2 = 1.5")),
  +     pch = 19)
  > set.seed(123)
  > plot(sim.citrus(coo, p = 0.1, model = "hGRF", cov.pars = c(1,
  +     5)), main = expression(paste("GRF Hierarquico, ", phi,
  +     " =   5")), pch = 19)
  > set.seed(123)
  > plot(sim.citrus(coo, p = 0.1, model = "hGRF", cov.pars = c(1,
  +     10)), main = expression(paste("GRF Hierarquico, ", phi,
  +     " = 10")), pch = 19)
  > set.seed(123)
  > plot(sim.citrus(coo, p = 0.1, model = "hGRF", cov.pars = c(1,
  +     15)), main = expression(paste("GRF Hierarquico, ", phi,
  +     " = 15")), pch = 19)
  > set.seed(123)
  > plot(sim.citrus(coo, p = 0.1, model = "clipGRF", cov.pars = c(1,
  +     5)), main = expression(paste("GRF Cliped,", phi, " =   5")),
  +     pch = 19)
  > set.seed(123)
  > plot(sim.citrus(coo, p = 0.1, model = "clipGRF", cov.pars = c(1,
  +     10)), main = expression(paste("GRF Cliped, ", phi, " = 10")),
  +     pch = 19)
  > set.seed(123)
  > plot(sim.citrus(coo, p = 0.1, model = "clipGRF", cov.pars = c(1,
  +     15)), main = expression(paste("GRF Cliped, ", phi, " = 15")),
  +     pch = 19)
  > set.seed(123)
  > plot(sim.citrus(coo, p = 0.1, model = "transGRF", cov.pars = c(1,
  +     5)), main = expression(paste("GRF Tranformado, ", phi,
  +     " =   5")), pch = 19)
  > set.seed(123)
  > plot(sim.citrus(coo, p = 0.1, model = "transGRF", cov.pars = c(1,
  +     10)), main = expression(paste("GRF Tranformado, ", phi,
  +     " = 10")), pch = 19)
  > set.seed(123)
  > plot(sim.citrus(coo, p = 0.1, model = "transGRF", cov.pars = c(1,
  +     15)), main = expression(paste("GRF Tranformado, ", phi,
  +     " = 15")), pch = 19)


pict

Figura 2: Vizualização de dados simulados segundo três diferentes modelos (colunas) com variação dos parâmetros em dada modelo (linhas)