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.

Opened 15 years ago

Closed 15 years ago

#940 closed defect (fixed)

AIA returns series with ratio applied

Reported by: cperales Owned by: Víctor de Buen Remiro
Priority: highest Milestone:
Component: Kernel Version: 2.0.1
Severity: blocker Keywords: AIA PulseOut DifEq
Cc: atorre@…, jconde@…

Description (last modified by Víctor de Buen Remiro)

AIA returns series with ratio applied instead pulses when you ask it to find pulses.
The InputDef's description is correct (i.e. @InputDef(0.8000483238833983, RationExpand(y2009m11d14,C,PulseOut))) but the compiled serie inside is the difEq(Ratio,Pulse).

¿Is this correct?

i.e. code:

Set AIA(
  SubSer(Gaussian(0,0.5),y2007,y2010),
  Ration (1-0.8056076456602815*B^7)/(1-0.773530260494357*B+0.02932819395975552*B^2-B^7+0.773530260494357*B^8-0.02932819395975552*B^9),
  SetOfRation(PulseOut),
  4,
  20
);

Thank you.

Change History (3)

comment:1 Changed 15 years ago by Víctor de Buen Remiro

Description: modified (diff)
Status: newaccepted

comment:2 Changed 15 years ago by Víctor de Buen Remiro

Efectivamente había un desliz y se usaba el cociente de polinomios correspondiente al modelo en lugar del correspondiente a la forma del outlier encontrado. El error había pasado desapercibido porque lo que se suele hacer es usar la descripción porque es más manejable como texto que como serie a la hora de manipular modelos.

Independientemente del error que ya está resuelto, debo avisar de que el test propuesto no sirve como chequeo de la función AIA desde el punto de vista estadístico, puesto que hay que pasarle una serie de residuos compatible con la estructura de intervención y el modelo indicados. Por ejemplo, este otro test crea una serie de residuos correspondiente a un output con un pulso de 5 sigmas

Real sigma = 0.5;
Polyn ar = 1-0.8056076456602815*B^7;
Polyn ma = 1-0.773530260494357*B+0.02932819395975552*B^2-B^7+0.773530260494357*B^8-0.02932819395975552*B^9;

Serie res = 
    SubSer(Gaussian(0,sigma),y2010,y2011)+
  5*sigma*SubSer(RationExpand(y2010m07,C,ar/ma),y2010,y2011);

Set AIA(
  res,
  ar/ma,
  SetOfRation(PulseOut),
  4,
  20
);

comment:3 Changed 15 years ago by Víctor de Buen Remiro

Resolution: fixed
Status: acceptedclosed

(In [2405]) Fixes #940

Note: See TracTickets for help on using tickets.