close Warning: Can't synchronize with repository "(default)" (/var/svn/tolp does not appear to be a Subversion repository.). Look in the Trac log for more information.

Changes between Version 1 and Version 2 of tolKit


Ignore:
Timestamp:
Jul 8, 2015, 10:06:13 AM (10 years ago)
Author:
Pedro Gea
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • tolKit

    v1 v2  
    1010
    1111Incorpora las funciones {{{Estimate}}} y {{{CalcForecasting}}}.
     12
     13== Utilidades ==
     14
     15El paquete {{{tolKit}}} incorpora funcionalidades de TOL al entorno R de una manera amigable para el usuario de R,
     16evitando en la medida de lo posible la necesidad de conocer las particularidades de la programación en TOL.
     17
     18En su primera versión, este paquete facilita el acceso a las funciones de estimación y previsión de modelos
     19lineales con estructura ARIMA denominadas: {{{Estimate}}} y {{{CalcForecasting}}}.
     20
     21=== tolEstimate  ===
     22
     23La función {{{Estimate}}} de TOL se incorpora a R como {{{tolEstimate}}}.
     24
     25{{{tolEstimate}}} estima un modelo lineal cuyo ruido presenta estructura ARIMA. Es decir, un modelo ARIMA con regresores externos.
     26
     27Nótese que no se trata de un modelo ARIMAX habitual, ya que el bloque AR afecta tanto al output (variable observada) como a los inputs (regresores).
     28
     29Ejemplo:
     30{{{
     31#!ruby
     32# Se utiliza la serie USAccDeaths
     33ser <- as.Serie(USAccDeaths)
     34# Se descartan los últimos dos años de datos
     35# para compararlos con la previsión
     36ser.in = Ssub(ser, to=-25)
     37ser.out = Ssub(ser, from=-24)
     38plot(ser.in, ser.out, style="TOL")
     39}}}
     40
     41{{{
     42#!ruby
     43# Se estima un ARIMA(0,1,1)xSARIMA(12,0,1,1)
     44est <- tolEstimate(ser.in,
     45  dif=(1-B)*(1-B^12), MA=list(1-0.1*B, 1-0.1*B^12))
     46print(est)
     47plot(est, style="TOL")
     48acf(residuals(est))
     49}}}
     50
     51=== CalcForecasting ===
     52
     53La función {{{CalcForecasting }}} de TOL se incorpora a R como {{{tolCalcForecasting}}}.
     54
     55{{{tolCalcForecasting}}} calcula previsiones para un modelo previamente estimado con {{{tolEstimate}}}.
     56
     57Ejemplo:
     58{{{
     59#!ruby
     60fst <- tolCalcForecasting(est, 24)
     61print(fst)
     62plot(fst, actual=ser.out, style="TOL")
     63plot(fst, actual=ser.out, style="TOL",
     64  mar=c(6,4,2,2), legend.names=list())
     65}}}