Geoestatística_dúvida...
Você pode tentar usar um modelo espaço-temporal de variograma, como o de Gneiting, disponível no R. Nesse caso não é necessário que você tenha os dados para toda a série de tempo. Uma boa referência é: Gneiting, T., Genton, M. G. and Guttorp, P. (2007). Geostatistical space-time models, stationarity, separability and full symmetry.<http://www.statsnetbase.com/books/6408/C5939_ch04.pdf>In Finkenstadt, B., Held, L. and Isham, V. (eds.), Statistical Methods for Spatio-Temporal Systems,<http://books.google.de/books?id=j-rtW1nVuCoC>Chapman & Hall/CRC, Boca Raton, pp. 151-175. []´s Em 11 de novembro de 2012 21:56, andrebvs <andrebvs@bol.com.br> escreveu:
Olá pessoal!
Estou iniciando os estudos em Geoestatística e não tenho muito conhecimento na área, então, sei que alguns colegas aqui trabalham na área, e se possível, gostaria de saber o seguinte:
Tenho alguns dados de totais anuais de precipitação pluviométrica de um determinado Estado, os quais estão organizados com as respectivas coordenadas em graus (Latitude, Longitude e Altitude). Cada coordenada (com o total pluviométrico), é referente à cada um dos 46 municipios do Estado. A pergunta é: o que dá pra fazer de análise com esses dados dentro da geoestatística, tendo em vista que alguns municipios vão do ano de 1911 à 1990, outros vão de 1920 a 1985, ou seja, há um certo desbalanceamento dos anos de um municipio para outro.
desde já agradeço!
*Att.* *André*
_______________________________________________ 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.
Pessoal, Estou tentando reproduzir um exemplo de modelagem espacial no INLA com dados de geoestatística, do artigo de Blangiardo et al (2012). O código em R está disponível e o exemplo 3 (SwissRainfall.R). Ao rodar o modelo aparece a seguinte mensagem:
mod = inla(formula, + data=inla.stack.data(stack, spde=spde),family="gaussian", + control.predictor=list(A=inla.stack.A(stack), compute=TRUE), + control.compute=list(cpo=TRUE, dic=TRUE), + keep=FALSE, verbose=TRUE) Erro em .local(x, i, j, ..., value) : not-yet-implemented 'Matrix[<-' method Calls: inla -> inla.predictor.section -> [<- -> [<- -> .local
Alguém tem idéia do quê possa estar acontecendo? Att, Taynãna Pesquisadora em Saúde Pública - ENSP/FIOCRUZ
Taynana Isto me parece um erro interno da função sugiro que voce coloque um exemplo reproduzível e mande para a equipe do INLA On Wed, 14 Nov 2012, Taynãna César Simões wrote:
Pessoal,
Estou tentando reproduzir um exemplo de modelagem espacial no INLA com dados de geoestatística, do artigo
de Blangiardo et al (2012). O código em R está disponível e o exemplo 3 (SwissRainfall.R).
Ao rodar o modelo aparece a seguinte mensagem:
mod = inla(formula, + data=inla.stack.data(stack, spde=spde),family="gaussian", + control.predictor=list(A=inla.stack.A(stack), compute=TRUE), + control.compute=list(cpo=TRUE, dic=TRUE), + keep=FALSE, verbose=TRUE) Erro em .local(x, i, j, ..., value) : not-yet-implemented 'Matrix[<-' method Calls: inla -> inla.predictor.section -> [<- -> [<- -> .local
Alguém tem idéia do quê possa estar acontecendo?
Att, Taynãna Pesquisadora em Saúde Pública - ENSP/FIOCRUZ
Oi Taynãna, qual a versão do seu INLA. inla.version() Eu consegui executar sem problemas na versão de 18 de agosto Elias T. Krainski
________________________________ De: Taynãna César Simões <tay_estatistica@yahoo.com.br> Para: "r-br@listas.c3sl.ufpr.br" <r-br@listas.c3sl.ufpr.br> Enviadas: Quarta-feira, 14 de Novembro de 2012 10:12 Assunto: [R-br] INLA Geoestatística
Pessoal,
Estou tentando reproduzir um exemplo de modelagem espacial no INLA com dados de geoestatística, do artigo
de Blangiardo et al (2012). O código em R está disponível e o exemplo 3 (SwissRainfall.R).
Ao rodar o modelo aparece a seguinte mensagem:
mod = inla(formula, + data=inla.stack.data(stack, spde=spde),family="gaussian", + control.predictor=list(A=inla.stack.A(stack), compute=TRUE), + control.compute=list(cpo=TRUE, dic=TRUE), + keep=FALSE, verbose=TRUE) Erro em .local(x, i, j, ..., value) : not-yet-implemented 'Matrix[<-' method Calls: inla -> inla.predictor.section -> [<- -> [<- -> .local
Alguém tem idéia do quê possa estar acontecendo?
Att, Taynãna Pesquisadora em Saúde Pública - ENSP/FIOCRUZ _______________________________________________ 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.
Olá André, o que eu faço geralmente é primeiro fazer um controle de qualidade de dados, eliminando ou criando um variável indicando se o dado é suspeito ou confiável. Segundo é determinar a disponibilidade de dados por estação e o período de dados de cada estação. Uma forma efetiva de fazer isso é através de um grafico do tipo levelplot (pacote lattice) on de x é o tempo e y a estação e z é uma variável (0,1) que indica se há dados ou não naquela estação. Assim você tem uma visão geral da disponibilidade de dados das estações e o maior período de dados simultâneos entre as estações. Um exemplo é o gráfico fornecido pela função "dwi" (Days with Information) do pacote hydroTSM. Se seu objetivo é obter séries médias espaciais (por estado, bacia hidrográfica, etc) é importante considerar a variação do número de estações disponíveis com o tempo. [ ]'s Jonatan Tatsch 2012/11/11 andrebvs <andrebvs@bol.com.br>
Olá pessoal!
Estou iniciando os estudos em Geoestatística e não tenho muito conhecimento na área, então, sei que alguns colegas aqui trabalham na área, e se possível, gostaria de saber o seguinte:
Tenho alguns dados de totais anuais de precipitação pluviométrica de um determinado Estado, os quais estão organizados com as respectivas coordenadas em graus (Latitude, Longitude e Altitude). Cada coordenada (com o total pluviométrico), é referente à cada um dos 46 municipios do Estado. A pergunta é: o que dá pra fazer de análise com esses dados dentro da geoestatística, tendo em vista que alguns municipios vão do ano de 1911 à 1990, outros vão de 1920 a 1985, ou seja, há um certo desbalanceamento dos anos de um municipio para outro.
desde já agradeço!
*Att.* *André*
_______________________________________________ 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.
-- #------------------------------------------------------------------------------# # Jônatan Dupont Tatsch # # Climate and Biosphere Laboratory # # Department of Atmospheric Sciences # # Institute of Astronomy, Geophysics and Atmospheric Sciences # # University of São Paulo # # Rua do Matão, 1226 # # Cid. Universitária, São Paulo, SP, Brazil, CEP: 05508-090 # # Phone:+55 11 3091-4772, Fax:+55 11 3091-4714 # # http://jonatandupont.weebly.com/index.html # #------------------------------------------------------------------------------#
participantes (6)
-
andrebvs -
Elias T. Krainski -
Jônatan -
Luis Paulo Braga -
Paulo Justiniano -
Taynãna César Simões