[R-br] grafico de função de autocorrelação

FHRB Toledo fernandohtoledo em gmail.com
Sexta Fevereiro 22 13:59:38 BRT 2013


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 em 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 em 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 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.
>>
>>
>>
>> --
>>
>> 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 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.
>>
>
>
> _______________________________________________
> 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/20130222/3023c354/attachment.html>


Mais detalhes sobre a lista de discussão R-br