Ola pessoal.<br>Bem pra economizar tempo, antes de ficar falando minha historinha eu vou fazer a minha pergunta pra não tomar muito o tempo de vocês.<br>Como correlacionar series temporais com um atributo externo de cada
serie temporal e estimar o quanto esse atributo deixa a serie temporal
estável.<br>
Acho que essa é a pergunta que eu tenho que responder e ainda não sei como.<br><br><br>Estava tentando aprender sobre series temporais mas estou apanhando aqui com o R.<br>Basicamente eu gostaria de comparar varias series temporais com uma variável externa.<br>
A ideia é mais ou menos o seguinte, tenho dados da fenologia de varias especies de planta. <br>Por exemplo, todo mes eu vou la no campo e para varias arvores de cada especie eu olho se ha brotos e estimo qual a porcentagem da copa que tem esses brotos.<br>
Ai fico com uma matriz onde colunas são especies e linhas são o estado fenológico da planta.<br>Além disso, cada especie tem uma medida caracteristica dela, da especie, que é a densidade de madeira.<br>A principio, me parece que especies com densidade de madeira alta so dão brotos no período de seca, enquanto as especies que tem densidade de madeira baixa são bagunçadas, brotam na chuva, fora da chuva, estão nem ai pra paçoca.<br>
<br>Eu estava lendo o livro Analyzing Ecological Data (2007 by Zuur, Ieno & Smith), e pra esse tipo de coisa eles sugerem o uso de Min/max auto-correlation factor analysis (MAFA) e Dynamic factor analysis (DFA). Mas nos scripts do livro ele fala que ambos os métodos<br>
não eram implementados no R na época e ainda não achei no R eles.<br><br>Eu estou garimpando atrás de algo nesse sentido aqui <a href="http://cran.r-project.org/web/views/TimeSeries.html">http://cran.r-project.org/web/views/TimeSeries.html</a> <br>
Mas Time Series parece ser um campo muito amplo e eu to ficando perdido.<br><br>Alguem tem alguma sugestão de literatura para Time Series mais voltado pra ecologia, pra fazer relação de series temporais com variáveis externas, mas não tão matematico ja que eu não tenho muita facilidade com matematica?<br>
<br>Eu fiz um exemplo de dados de mais ou menos como estão as coisas.<br>O exemplo é uma matriz de series temporais chamada brota, um vetor chamado chuva que é a pluviosidade no mês e um vetor que são valores de densidade de madeira das especies:<br>
Eu estava a principio explorando graficamente tudo , se alguem tiver alguma sugestão.<br>Segue o CMR.<br><br><br>#Exemplo:<br>#Criando dados<br>#Densidade de madeira<br>set.seed(69)<br>dmad<-round(c(rnorm(5,mean=50),rnorm(5,mean=90)))<br>
names(dmad)<-c(paste("sp",1:10,sep=""))<br>dmad<br><br>#Chuva<br>chuva<-c(150,100,50,40,20,20,30,50,70,100,150,200,150,100,50,30,20,20,40,50,70,100,150,200)<br><br><br>#Padrão Anual<br>ydeb<-c(10,10,10,10,20,40,50,40,20,10,10,10)<br>
ydea<-rnorm(12,mean=mean(ydeb),sd=sd(ydeb))<br><br>#Amostras por 2 anos para especies com baixa e alta densidade de madeira<br>brota<-matrix(NA,nrow=24, ncol=10,dimnames=list(paste("mes",1:24,sep=""),paste("sp",1:10,sep="")))<br>
for (i in 1:5) {<br> brota[,i]<-round(c(ydeb+rnorm(12,mean=5,sd=5),ydeb+rnorm(12,mean=5,sd=5)),digits=2)<br> }<br>for (i in 6:10) {<br> brota[,i]<-round(c(ydea+rnorm(12,mean=5,sd=5),ydea+rnorm(12,mean=5,sd=5)),digits=2)<br>
}<br>#TRansformando em serie temporal<br>brota<-ts(brota)<br>#Observando as Series<br>plot(brota)<br><br>#Observando as series juntas, separando as densidade de madeira por cores<br>#Mais vermelho = mais densa, mais claro = menos densa<br>
cores<-heat.colors(10,alpha=1)<br>matplot(c(1:24),brota,type="l",lwd=2,lty=1,xlab="meses",ylab="brota",col=cores[order(dmad)])<br><br>#Observando a chuva junto no grafico anterior<br>graficobarras<-barplot(chuva/max(chuva)*100)<br>
matlines(graficobarras,brota,type="l",lwd=2,lty=1,xlab="meses",ylab="brota",col=cores[order(dmad)])<br><br>#Obsevando a partir de uma pca como são, mas acredito que isso esta errado fazer.<br>
brota.pca<-princomp(brota)<br>biplot(brota.pca)<br><br clear="all"><br>-- <br><div>Grato<br>Augusto C. A. Ribas</div>
<div> </div>
<div>Site Pessoal: <a href="http://augustoribas.heliohost.org" target="_blank">http://augustoribas.heliohost.org</a></div>
<div>Lattes: <a href="http://lattes.cnpq.br/7355685961127056" target="_blank">http://lattes.cnpq.br/7355685961127056</a><br></div><br>