[R-br] plot de uma timeseries univariada com seu modelo SARIMA fitted?

[carlos pombo sonderblohm]

Bom dia pessoal,
estou a tentar fazer um plot de um modelo SARIMA de uma serie temporal>

library(stats)
library(astsa)
library(graphics)

SL <-c(-0.289220360456518, -0.412287982738895, -0.30186163626253,
            -0.325463729174269, -0.950943918812091, -0.621873409193565,
-0.454370839103369,
            -0.201717177298859, -0.0125484487133019, 0.37501627306213,
0.860628415115297,
            0.573001978891618, 0.949171439258297, 0.884107931496505,
0.321528581746018,
            0.197039928762375, -0.0511602427005644, -0.14797533453817,
0.0997648324025207,
            -0.356603216398346, -0.486741440877119, -0.796433339558246,
-0.797717898751611,
            -1.03991357760201, -0.272391281348932, -0.280840105493716,
-0.220566038983021,
            -0.092982026032348, -0.737941509253829, -0.25720443208744,
0.281594104947966,
            0.295827677265899, 0.381173023222877, 0.410685495501984,
0.0935538548738117,
            0.320870181634448, -0.176916211372866, -0.434468255568435,
-0.177142353391016,
            -0.479723718529784, -1.18506069524595, -0.845078218877219,
-0.466484365108015,
            -0.378310531861291, -0.568334035779147, -0.756612458676533,
-0.612289956791277,
            -0.595067441206468, 0.227790700535617, -0.328474103547927,
0.0382893789935963,
            -0.107472883393398, 0.104836312989005, 0.206908642082212,
0.566811114425553,
            0.226484501566576, 0.268795384938629, 0.247658319827575,
0.728293948754543,
            0.704287103883026, 0.665057803685019, 1.01752383083208,
1.34104381438303,
            1.10724539596799, 0.850746327229352, 0.746239914170796,
0.665531282137068,
            0.469573210434806, -0.0626525410819481, -0.294702780105021,
-0.160850370184094,
            0.511318161080905)

SL <- ts(SL, start = c (1995, 1), freq = 12) #transformação numa ts

ar1 <-sarima(SL, 1,0,0,0,0,0,12) #o modelo seleccionado

plot(SL, type ='o')) # o plot da serie original
lines(fitted(ar1),col=2) # estou a tentar adicionar o modelo fitted a meus
dados

Mas não consigo fazer o grafico pretendido, quero fazer que o grafico tenha
os pontos o linhas da serie original e no mesmo grafico, o modelo ajustado
para ver como bate com os dados originais,
alguma dica?

Meu R es:
R version 2.13.1 (2011-07-08)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: x86_64-pc-mingw32/x64 (64-bit)
SO, windows 7

Obrigado

Carlos


-- 
Carlos A. Pombo Sonderblohm
PhD Student on Marine Science (Fisheries)
Faculdade de Ciências e Tecnología
Universidade do Algarve,
Campus de Gambelas
8005-139 Faro
Portugal
Tef. 289 800 905 ext. 7605
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121018/6aef4ad3/attachment.html>