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]]

