Rodrigo, aumente a variância em cov.pars = c(56, 4) Att. Elias T. Krainski
>________________________________ > De: Rodrigo Sant'Ana <[email protected]> >Para: [email protected] >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 >[email protected] >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. > > _______________________________________________ R-br mailing list [email protected] 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.
