<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Hola,<br>
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:<br>
<br>
# Estimacion TAR(2;2,2) con delta 1:<br>
mod.setar <- setar(x, m = 2, mL = 2, mH = 2, thDelay = 1) <br>
mod.setar<br>
summary(mod.setar)<br>
<br>
mod.setarc<-setar(x, m=2, mL=2, mH=2, thDelay=1,
include=c("none")) ##eliminamos la constante del modelo<br>
mod.setarc<br>
summary(mod.setarc)<br>
<br>
#Ajustamos el modelo para los primeros 104 valores de la serie<br>
set.seed(10)<br>
mod.test <- list()<br>
x.train <- window(x, end = 1924)<br>
x.test <- window(x, start = 1925)<br>
mod.test[["linear"]] <- linear(x.train, m = 2)<br>
mod.test[["setar"]] <- setar(x.train, m = 2, thDelay = 1)<br>
mod.test[["setarc"]] <- setar(x.train, m = 2, thDelay = 1,
include=c("none"))<br>
mod.test[["nnet"]] <- nnetTs(x.train, m = 2, size = 3, control =
list(maxit = 1e+05))<br>
mod.test[["aar"]] <- aar(x.train, m = 2)<br>
<br>
#ahora la predicción:<br>
frc.test <- lapply(mod.test, predict, n.ahead = 10)<br>
<br>
#Dibujamos en una misma gráfica todas las predicciones que hemos
calculado para los modelos estimados:<br>
plot(x.test, ylim = range(x))<br>
for (i in 1:length(frc.test)) <br>
lines(frc.test[[i]], lty = i+1, col = i + 1)<br>
legend(1925, 2.4, lty=1:(length(frc.test)+1),
col=1:(length(frc.test)+1), legend=c("observed", names(frc.test)))<br>
<br>
y el resultado que obtengo es:<br>
<img src="cid:part1.05040906.09090308@hotmail.com" alt=""><br>
en el que como pueden ver noaparece la representacion del modelo
setarc, que es el setar sin la constante.<br>
¿Alguien ha trabajado con esto y me puede ayudar?<br>
<br>
Muchas gracias de antemano por la atención.<br>
Reciban un cordial saludo.<br>
<br>
<br>
</body>
</html>