[R-br] Duvida em estatística espacial - Kernel em duas dimensões levando em conta prevalências.
Augusto Ribas
ribas.aca em gmail.com
Quarta Fevereiro 25 18:21:45 BRT 2015
Ola.
Gostaria de saber se alguém poderia me sugerir algo para avaliar o risco,
ou aonde uma doença ocorre espacialmente, mais comumente.
Meus dados são constituídos de coordenadas geográficas de cidades, na forma
de latitude e longitude. Dai para cada cidade eu tenho o número de casos de
uma doença para um total de vaquinhas observadas. E calculo a prevalência
como os casos pelo total de vaquinhas observadas.
Então eu tenho esses dados, e gostaria de fazer um mapinha, mostrando
espacialmente se existe uma área que tem maior ocorrência, maior chance de
ocorrência. Apenas para explorar os dados a principio.
Segue um pequeno exemplo de como estão os dados.
#########################################################
library(MASS)
set.seed(1)
dados<-data.frame(lat=sample(1:10,20,replace=T),
lon=sample(1:10,20,replace=T),
casos=rpois(20,5),
total=rpois(20,60)
)
dados$prevalencia<-dados$casos/dados$total
#Eu gostaria de fazer uma figura mais ou menos como essa, mas levando em
consideração as prevalências.
contorno<-kde2d(dados$lat,dados$lon)
contour(contorno)
#############################################################
Alguém tem alguma sugestão? Eu não sei nem o que começar a ler exatamente
para
Agradeço desde ja a atenção.
--
Grato
Augusto C. A. Ribas
Site Pessoal: http://recologia.com.br/ <http://augustoribas.heliohost.org>
Github: https://github.com/Squiercg
Lattes: http://lattes.cnpq.br/7355685961127056
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150225/418ecfe9/attachment.html>
Mais detalhes sobre a lista de discussão R-br