Rodrigo, aumente a variância em cov.pars = c(56, 4) Att. Elias T. Krainski
________________________________ De: Rodrigo Sant'Ana <rodrigo.gringo@gmail.com> Para: r-br@listas.c3sl.ufpr.br Enviadas: Quarta-feira, 1 de Fevereiro de 2012 17:26 Assunto: [R-br] Gaussian Random Fields - grf::geoR
Pessoal,
estou com uma pequena dúvida na simulação de um campo Gaussiano aleatório, com a função grf do pacote geoR. Acredito que vocês poderão me ajudar e/ou até me corrigir.
A idéia central é simular um campo onde os valores dos dados representaram densidades de uma determinada espécie, para posterior determinação da densidade total da área em questão. Para tanto, estou replicando a seguir o código executado por mim.
### set.seed(3)
dados <- grf(400, grid = "reg", xlims = c(0, 20), ylims = c(0, 20), cov.pars = c(56, 4), cov.model = "gaussian", nugget = 1, mean = 57.8)
summary(dados1) $n [1] 400
$coords.summary x y min 0 0 max 20 20
$distances.summary min max 1.052632 28.284271
$data.summary Min. 1st Qu. Median Mean 3rd Qu. Max. 44.06 55.45 59.33 59.36 63.71 71.96
$borders.summary Var1 Var2 min 0 0 max 20 20
$others [1] "cov.model" "nugget" "cov.pars" "kappa" "lambda" [6] "aniso.pars" "method" ".Random.seed" "messages" "call"
attr(,"class") [1] "summary.geodata"
=======
Dúvida: Ao solicitar o summary do objeto "dados", observa-se que o amplitude dos valores ($data.summary) gerados pela função grf vai de 44.06 a 71.96. Entretanto, gostaria de simular algo mais próximo da realidade das espécies de animais alvo desta simulação, onde esta amplitude deveria que ser maior, com valores minimos mais próximos a zero e os valores máximos poderiam ser os mesmos (próximos á 80).
Seria possível aumentar a amplitude dos valores gerados para os dados ou estou devaneando demais?
P.S.: A idéia não é induzir a geração de valores aleatórios, o intuito é aumentar a amplitude dos dados ok! Desculpem-me se ficou um pouco tendenciosa a explanação a cima, mas a ideia era meramente explicar minha dúvida.
Obrigado pela atenção e ajuda de todos.
Abs _______________ Rodrigo Sant'Ana Oceanógrafo
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.