[R-br] Gerar mapa de distribuição por UF

[Gilbert Queiroz]

Então, Elias,
Eu montei a matriz da seguinte forma:
- Uma coluna com os municípios;
- Em seguidam uma coluna com o Nr de casos/Município
- Depois, umas 6 colunas com variáveis ambientais.

E fui repetindo desta forma para cada ano (empilhando) de 2002 a 2011...

Minha ideia é montar um modelo que tenha Nr de casos/mun = f( variáveis ambientais)

Depois, quero jogar isto num shapefile usando o TerraView para mostra a evolução dos casos naqueles anos, por município.

É possível???

Abs.





________________________________
 De: Elias Teixeira Krainski <eliaskrainski em yahoo.com.br>
Para: r-br em listas.c3sl.ufpr.br 
Enviadas: Domingo, 3 de Março de 2013 17:35
Assunto: Re: [R-br] Gerar mapa de distribuição por UF
 

> Como eu entro com os dados no R para obter uma modelagem de dados
> espaço-temporais????
> Tenho um arquivo onde nas colunas estão as variáveis que quero
> trabalhar e nas linhas os municípios. Só que eu repito estes
> municípios para os anos de 2002 a 2011. É assim mesmo????

Que tipo de modelo vc tem em mente? 

Eu tenho usado INLA, não sei como trabalham outros pacotes. 

No INLA trabalhamos com indices. No seu caso, é só 
'empilhar' os dados. Assim vc pode ter um vetor de 
dados (em vez de matriz), um índice para área e 
outro para tempo.

Na classe dos modelos dinâmicos, veja exemplos em 
modelos dinâmicos espaço temporais em 
http://www.sciencedirect.com/science/article/pii/S0167947311003999
e maiores detalhes da implementação em 
http://www.r-inla.org/examples/case-studies/ruiz-cardenas-kraonski-and-rue-2010

O INLA esta em constante desenvolvimento e atualmente 
é possível simplificar um nível nos códigos nesses links. 

Por exemplo, usando a extratégia de definir matriz 
de precisão pelo produto de Kronecker, por exemplo, 
entre um efeito espacial e um temporal. 

Na classe de modelos de interação espaço temporais 
do Held vc pode usar a estratégia de Kronecker. 
Para interação do tipo 4, você pode fazer simplesmente 

f(i, model='besag', graph=grafo, group=t, 
    control.group=list(model='rw1'))

em que 'i' indica as áreas e 't' o tempo.




_______________________________________________
R-br mailing list
R-br em 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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130304/092a99be/attachment.html>