[R-es] como hacer grafico de un modelo setar

Victoria Yañez vyj20 en hotmail.com
Lun Sep 2 19:10:46 CEST 2013


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>


Más información sobre la lista de distribución R-help-es