2 Exemplo de dados simulados
A função sim.citrus() implementa os modelos descritos.
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)