Re: [R-br] grafico de função de autocorrelação
library(tseries) y = log(AirPassengers) acf(diff(y)) library(plotrix) draw.ellipse(0.5,0.2,0.5,0.2) (S,f,P) Allaman
Obrigada! 2013/2/22 Ivan Bezerra Allaman <ivanalaman@yahoo.com.br>:
library(tseries) y = log(AirPassengers) acf(diff(y))
library(plotrix) draw.ellipse(0.5,0.2,0.5,0.2)
(S,f,P) Allaman
_______________________________________________ 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.
-- Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743
Com pequenas alterações library(tseries) y = log(AirPassengers) radf=acf(diff(y)) li=-sd(diff(y))*1.65 ss=which(radf$acf<li | radf$acf>-li) library(plotrix) for (i in 2:length(ss)) draw.ellipse(radf$lag[ss[i]],ifelse(radf$acf[ss[i]]>0,-li,li),0.02,0.08,border="red") Sds 2013/2/22 Ivan Bezerra Allaman <ivanalaman@yahoo.com.br>:
library(tseries) y = log(AirPassengers) acf(diff(y))
library(plotrix) draw.ellipse(0.5,0.2,0.5,0.2)
(S,f,P) Allaman
_______________________________________________ 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.
--
Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743 _______________________________________________ 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.
Natália, Meu centavo na discussão, veja a função ellipse do pacote mixtools: http://www.inside-r.org/packages/cran/mixtools/docs/ellipse Acho que serve para alguma coisa! att, FH 2013/2/22 Robert Iquiapaza <rbali@ufmg.br>
Com pequenas alterações
library(tseries) y = log(AirPassengers) radf=acf(diff(y)) li=-sd(diff(y))*1.65 ss=which(radf$acf<li | radf$acf>-li) library(plotrix) for (i in 2:length(ss))
draw.ellipse(radf$lag[ss[i]],ifelse(radf$acf[ss[i]]>0,-li,li),0.02,0.08,border="red")
Sds
2013/2/22 Ivan Bezerra Allaman <ivanalaman@yahoo.com.br>:
library(tseries) y = log(AirPassengers) acf(diff(y))
library(plotrix) draw.ellipse(0.5,0.2,0.5,0.2)
(S,f,P) Allaman
_______________________________________________ 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.
--
Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743 _______________________________________________ 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.
_______________________________________________ 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.
participantes (4)
-
FHRB Toledo -
Ivan Bezerra Allaman -
Natalia Martins -
Robert Iquiapaza