[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