<p dir="ltr">Veja o pacote spatstat para mapear padrões. Supondo q a doenca pode acontecer em toda area de estudo.</p>
<div class="gmail_quote">Em 25/02/2015 18:21, "Augusto Ribas" <<a href="mailto:ribas.aca@gmail.com">ribas.aca@gmail.com</a>> escreveu:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>Ola.<br><br></div>Gostaria de saber se alguém poderia me sugerir algo para avaliar o risco, ou aonde uma doença ocorre espacialmente, mais comumente.<br><br></div>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.<br><br></div>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.<br><br></div>Segue um pequeno exemplo de como estão os dados.<br><div><br><br>#########################################################<br><div><div>library(MASS)<br>set.seed(1)<br>dados<-data.frame(lat=sample(1:10,20,replace=T),<br>                  lon=sample(1:10,20,replace=T),<br>                  casos=rpois(20,5),<br>                  total=rpois(20,60)<br>                  )<br>dados$prevalencia<-dados$casos/dados$total<br><br></div><div>#Eu gostaria de fazer uma figura mais ou menos como essa, mas levando em consideração as prevalências.<br></div><div>contorno<-kde2d(dados$lat,dados$lon)<br>contour(contorno)<br>#############################################################<br><br></div><div>Alguém tem alguma sugestão? Eu não sei nem o que começar a ler exatamente para <br></div><div>Agradeço desde ja a atenção.<br></div><div><br><div><div>-- <br><div><div dir="ltr"><div>Grato<br>Augusto C. A. Ribas</div>
<div> </div>
<div>Site Pessoal: <a href="http://recologia.com.br/" target="_blank">http://recologia.com.br/</a><a href="http://augustoribas.heliohost.org" target="_blank"></a></div><div>Github: <a href="https://github.com/Squiercg" target="_blank">https://github.com/Squiercg</a></div>
<div>Lattes: <a href="http://lattes.cnpq.br/7355685961127056" target="_blank">http://lattes.cnpq.br/7355685961127056</a><br></div></div></div>
</div></div></div></div></div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div>