R help es - Sep 2013 - como hacer grafico de un modelo setar

Hola,
Estoy intentando realizar un ajuste a un modelo setar con R y me surgen 
problemas a la hora de representar el modelo setar sin la constante. El 
código es el siguiente:

# Estimacion TAR(2;2,2) con delta 1:
mod.setar <- setar(x, m = 2, mL = 2, mH = 2, thDelay = 1)
mod.setar
summary(mod.setar)

mod.setarc<-setar(x, m=2, mL=2, mH=2, thDelay=1, include=c("none"))
##eliminamos la constante del modelo
mod.setarc
summary(mod.setarc)

#Ajustamos el modelo para los primeros 104 valores de la serie
set.seed(10)
mod.test <- list()
x.train <- window(x, end = 1924)
x.test <- window(x, start = 1925)
mod.test[["linear"]] <- linear(x.train, m = 2)
mod.test[["setar"]] <- setar(x.train, m = 2, thDelay = 1)
mod.test[["setarc"]] <- setar(x.train, m = 2, thDelay = 1, 
include=c("none"))
mod.test[["nnet"]] <- nnetTs(x.train, m = 2, size = 3, control = 
list(maxit = 1e+05))
mod.test[["aar"]] <- aar(x.train, m = 2)

#ahora la predicción:
frc.test <- lapply(mod.test, predict, n.ahead = 10)

#Dibujamos en una misma gráfica todas las predicciones que hemos 
calculado para los modelos estimados:
plot(x.test, ylim = range(x))
for (i in 1:length(frc.test))
     lines(frc.test[[i]], lty = i+1, col = i + 1)
     legend(1925, 2.4, lty=1:(length(frc.test)+1), 
col=1:(length(frc.test)+1), legend=c("observed", names(frc.test)))

y el resultado que obtengo es:

en el que como pueden ver noaparece la representacion del modelo setarc, 
que es el setar sin la constante.
¿Alguien ha trabajado con esto y me puede ayudar?

Muchas gracias de antemano por la atención.
Reciban un cordial saludo.


------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20130902/18d5f6fb/attachment-0001.html>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: djbfedab.jpg
Type: image/jpeg
Size: 86214 bytes
Desc: no disponible
URL:
<https://stat.ethz.ch/pipermail/r-help-es/attachments/20130902/18d5f6fb/attachment-0001.jpg>

Victoria,

es probable que las predicciones del modelo setarc estén fuera del  
área del gráfico (definido por ylim).
Un saludo. Olivier
--  
____________________________________

Olivier G. Nuñez
Email: onunez@unex.es
http://matematicas.unex.es/~onunez
Tel : +34 663 03 69 09
Departamento de Matemáticas
Universidad de Extremadura

____________________________________





El 02/09/2013, a las 19:10, Victoria Yañez escribió:
> Hola,
> Estoy intentando realizar un ajuste a un modelo setar con R y me  
> surgen problemas a la hora de representar el modelo setar sin la  
> constante. El código es el siguiente:
>
> # Estimacion TAR(2;2,2) con delta 1:
> mod.setar <- setar(x, m = 2, mL = 2, mH = 2, thDelay = 1)
> mod.setar
> summary(mod.setar)
>
> mod.setarc<-setar(x, m=2, mL=2, mH=2, thDelay=1,
include=c("none"))
> ##eliminamos la constante del modelo
> mod.setarc
> summary(mod.setarc)
>
> #Ajustamos el modelo para los primeros 104 valores de la serie
> set.seed(10)
> mod.test <- list()
> x.train <- window(x, end = 1924)
> x.test <- window(x, start = 1925)
> mod.test[["linear"]] <- linear(x.train, m = 2)
> mod.test[["setar"]] <- setar(x.train, m = 2, thDelay = 1)
> mod.test[["setarc"]] <- setar(x.train, m = 2, thDelay = 1,
include=c
> ("none"))
> mod.test[["nnet"]] <- nnetTs(x.train, m = 2, size = 3, control
=
> list(maxit = 1e+05))
> mod.test[["aar"]] <- aar(x.train, m = 2)
>
> #ahora la predicción:
> frc.test <- lapply(mod.test, predict, n.ahead = 10)
>
> #Dibujamos en una misma gráfica todas las predicciones que hemos  
> calculado para los modelos estimados:
> plot(x.test, ylim = range(x))
> for (i in 1:length(frc.test))
>     lines(frc.test[[i]], lty = i+1, col = i + 1)
>     legend(1925, 2.4, lty=1:(length(frc.test)+1), col=1:(length 
> (frc.test)+1), legend=c("observed", names(frc.test)))
>
> y el resultado que obtengo es:
>
> en el que como pueden ver noaparece la representacion del modelo  
> setarc, que es el setar sin la constante.
> ¿Alguien ha trabajado con esto y me puede ayudar?
>
> Muchas gracias de antemano por la atención.
> Reciban un cordial saludo.
>
>
> <djbfedab.jpg>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es

	[[alternative HTML version deleted]]

Quitar la constante del modelo, no parece una buena estrategia para  
estos datos.
En cualquier caso, si quieres visualizar las predicciones del modelo  
setarc con las de los demás modelos,
prueba sustituir la linea

plot(x.test, ylim = range(x))

por

plot(x.test, ylim = range( union( x , frc.test[["setarc"]] )))

Un saludo. Olivier
--  
____________________________________

Olivier G. Nuñez
Email: onunez en iberstat.es
Tel : +34 663 03 69 09
Web: http://www.iberstat.es

____________________________________




El 02/09/2013, a las 21:03, Victoria Yañez escribió:
> y sabrías decirme  cómo puedo solucionar esto, por favor?