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/

Github: https://github.com/Squiercg

Lattes: http://lattes.cnpq.br/7355685961127056

_______________________________________________
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.