Boa noite R's, estou buscando, homericamente, gerar uma distribuição de densidades multivariada em uma área irregular, porém estou dando murro em ponta de faca. Acredito que devo estar errando em algo simples, e fácil. O intuito é simular a distribuição de densidade de um peixe em uma área conhecida, para que posteriormente, estas densidades sejam alvos de simulações amostrais. Abaixo segue o CMR de minha tentativa! ###---------------------------------- Simulação Densidade Normal Multivariada # Carregando os dados em anexo! load("R-BR.RData") ## A área de interesse deveria ser esta plotada a seguir plot(bor, type = "l") # Separação dos vetores Lat e Long do objeto dados long <- dados$long lat <- dados$lat ## Utilizando um kernel de Densidade da Normal Multivariada # Função de amostragem estratificada! system.time( d <- kde2d(long, lat, n = 100) ) ## Função criada para gerar as figuras em perspectiva (persp) e contorno ## (contour) de mesma densidade (biomassa). kde2dplot <- function(d, # Densidades 2D computada pela # função kde2d ncol=50, # Número de cores zlim=c(0,max(z)), # Limites da coordenada z nlevels=20, # vide nlevels em contour theta=30, # vide theta em persp phi=60) # vide phi em persp { z <- d$z # Recebe as densidades (biomassa) nrz <- nrow(z) # Recebe o número de linhas de z ncz <- ncol(z) # Recebe o número de colunas de z ## O objeto couleurs determina o esquema de cores a ser utilizado, sendo ## as opções: rainbow, heat.colors, terrain.colors, topo.colors, ## cm.colors. couleurs <- tail(terrain.colors(trunc(1.4 * ncol)), ncol) fcol <- couleurs[trunc(z/zlim[2]*(ncol-1)) + 1] dim(fcol) <- c(nrz, ncz) fcol <- fcol[-nrz, -ncz] par(mar = c(5, 5, 2, 2)) image(d, col = couleurs, xlab = "Longitude", ylab = "Latitude") contour(d, add = T, nlevels = nlevels, xlab = "Longitude", ylab = "Latitude") box() } ## Salva a figura de saída da simulação gerada kde2dplot(d) polygon(bor) ##----------------------------------------------------------------------------------------------------- Se alguém tiver uma alternativa, realmente salvaria uma alma! Brincadeiras a parte, queria muito agradece-los pela atenção. ______________ Rodrigo Sant'Ana Oceanógrafo