Leunketa esponentzial

Wikipedia(e)tik
Hona jo: nabigazioa, Bilatu
Leunketa esponentziala denbora serie bat leundu eta epe laburreko aurresanak egiteko erabiltzen da. Irudian (x ardatza: denbora; y ardatza: seriea) puntu gorriek eta lerro urdinak serie historikoa eta horren leunketa esponentziala adierazten dituzte: leunketa esponentziala aurreko balio guztiak haztatzen ditu balio leundua finkatzeko.

Leunketa esponentziala denbora serieak aztertu eta epe laburrerako aurresanak egiteko teknika estatistikoen multzo bat da, jasotako datu historikoetatik abiaturik errepikapenezko formulak erabiltzen dituena, leunketa-parametroak finkatu ondoren. Bere bertsio sinpleenean, batezbesteko higikorra besterik ez da baina datu historikoak zenbat atzerago izan orduan eta eta neurri txikiagoan haztatu egiten dituena. Esponentzial adjektiboak datu historikoei denboran atzera joan ahala eta funtzio esponentzialaren arabera gero eta garrantzi edo pisu txikiagoa ematen zaiela adierazten du. Bere abantaila nagusia sinpletasuna da: kalkuluaren aldetik teknika errazak dira, ezarri beharreko parametroak eta eskuratzen diren emaitzak aise interpretatzen dira eta datu kopuru txikiekin ere erabil daiteke. Leunketa esponentziala datu historikoetan bakarrik oinarritzen da eta aurrez ez du datuei buruz inongo eredu estatistikorik ezartzen, denbora serieei aplikatzen zaizkien eredu batzuetarako emaitzak hobezinak direla frogatu bada ere. Bhear bezala aplikatuz gero, emaitzak eta aurresanak fidagarriak direla egiaztatu da.[1]

Leunketa esponentzialeko teknikak denbora-seriearen izaerara egokitzen dira. Bertsio sinpleenean, leunketa esponentzial sinplea joera eta urtarokotasunik gabeko serieetarako erabiltzen da. Leunketa esponentzialeko teknika konplexuagoak asmatu dira, joera lineal eta ez-lineal kontuan hartu eta urtarokotasuna era batukorrean eta biderkakorrean gaineratzen dutenak, horretarako leunketa anizkoitzak burutuz joera eta urtarokotasun osagaietarako. Leunketa esponentziala epe laburraz haraindi egokiagoa izan dadin, joera moteldu egiten duten parametroak ere barneratu daitezke.

Leunketa esponentzial sinplea[aldatu | aldatu iturburu kodea]

x_t t aldiko datu historiko edo erreala eta s_t t aldiko leunketaren emaitza izanik, ikertzaileak finkatzen duen 0 < \alpha < 1 leunketa parametro baterako, honako hau da leunketa esponentzial sinplearen errepikapen-formula, aurreko leunketa-baliotik leunketa balioak eguneratu edo berritu egiten dituena:

s_t = \alpha x_{t} + (1-\alpha)s_{t-1}

Leunketa abiarazten duen s0 hasierako leunketa-balioa finkatzeko, datu historiko guztien batezbestekoa [1]edo serieko lehenengo balioen batezbestekoa[2] eta serie historikoko lehen balioa [3] (aurreragoko adibidean erabiltzen dena) proposatu dira.

Eskubiko adierazpeneko s_{t-1} leunketa-balioetarako formula aplikatuz:

s_{t-1} = \alpha x_{t-1} + (1-\alpha)s_{t-2}
s_{t-2} = \alpha x_{t-2} + (1-\alpha)s_{t-3}
\ldots\ldots\ldots\ldots

eta lehenengo formulan ordeztuz, leunketa-balioetarako formula hau eskuratzen da:

 s_t = \alpha x_t + \alpha (1- \alpha ) x_{t-1} + \alpha (1- \alpha )^2 x_ {t-2} + ... + \alpha (1- \alpha )^{t-1} x_1 + \alpha (1- \alpha )^{t} x_0
 = \alpha \sum_{i=0}^{t-1}{(1- \alpha )^i x_{t-i}}+ \alpha (1- \alpha )^{t} x_0
Leunketa esponentzialerako haztapenak serie historikoan: leunketa esponentzialean, serie historikoko balioak α, α(1-α)2, α(1-α)3, ... koefiziente beherakorren arabera (0<α<1) haztatzen dira leunketa balioa batezbesteko batez kalkulatzean; α handia denean (gorriz), azken balioari haztapen handia ematen zaio eta beraz, hori kontuan hartzen da batezbestekoa kalkulatzean; α txikia denean (urdinez), berriz, azken balioak haztapen handienekoa izan arren, ez da besteak baino askoz ere nuerri handiagoan haztatzen.

Era horretan, leunketa esponentziala batezbesteko aritmetiko haztatu moduan agertzen da, non datu historikoei ematen zaizkien haztapenek segida geometriko bati jarraitzen dioten. Segida geometriko bat funtzio esponentzialaren bertsio diskretua delako ematen zaio leunketa esponentzial izena. Ohartzekoa da haztapen-koefiziente horiek gero eta txikiagoak direla, datuetan zenbat eta atzerago joan: datu berriak zaharrak baino neurri handiagoan haztatu edo kontuan hartzen dira.

Aurresanak[aldatu | aldatu iturburu kodea]

Leunketa esponentziala aurresanerako erabiltzen denean, t+1 aldirako \hat{x}_{t+1} aurreikuspen moduan t aldiko leunketaren emaitza ezartzen da:

\hat{x}_{t+1}=s_t

Eskuarki, etorkizuneko aurresanetarako ere kalkulaturiko azken leunketa balioa hartzen da:

\hat{x}_{t+k}=s_t\ ,\ t=1,2,\ldots.

Aurreko adierazpen honetatik leunketa esponentziala honela ere adieraz daiteke, aurresana egitean sortzen den errorea erabiliz:

\hat{x}_{t+1}=s_t =\alpha \cdot x_{t} + (1-\alpha) \cdot s_{t-1}=
=\alpha \cdot x_{t} + (1-\alpha) \cdot \hat{x}_{t}=
=\hat{x}_{t} + \alpha \cdot (x_t-\hat{x}_{t})=
=\hat{x}_{t} + \alpha \cdot e_t

Horrela, leunketa esponentzialaren arabera hurrengo aldirako aurresana aurreko aurresanean bertan sortutako errorearen zati bat gehituz kalkulatzen da.[4]

α leunketa parametroaren interpretazioa eta finkapena[aldatu | aldatu iturburu kodea]

α leunketa parametroa [0,1] tarteko balioak hartzen dituen balio bat da, subjektiboki aukeratu edo optimotasun-irizpide bati jarraiki finkatzen dena. Zenbat eta txikiagoa izan, leunketa orduan eta nabarmenagoa da, hau da, joerarik gabeko denbora-serie batek erakuts ditzakeen gorabehera puntualak nabariago leundu eta joera orduan eta era garbiagoan azalaraziko du. Leunketa parametroa handia denean, berriz, leunketaren balioak azken behaketa edo datuen mendean geratzen dira neurri handiago batean eta horrenbestez nabarmenago islatuko dira gorabehera puntualak leunketaren seriean. Adibidez, α>0.9 balioak erabiltzen direnean, leunketa balioek denbora seriea atzerapen txiki batez ia errepikatu besterik ez dute egiten.

Aurresanei dagokienean, α balioa handia denean, erabiltzailea bereziki azken balioetan oinarritzen da bere aurresanak egiteko; α balio txikietarako, aitzitik, serieko balio guztiak modu orekatuago batez hartzen ditu kontuan, leunketa esponentzialean zenbat eta atzerago joan, haztapenak orduan eta txikiagoak direla kontuan harturik betiere. Adibidez, salmenten aurresanak egiteko erabiltzen dituen saltzaile batek, α handia finkatuko du, bereziki azken egunetako salmentei begiratzen badie; atzerago joan eta aspaldiko salmentak ere hartzen baditu kontuan, α txikia finkatu beharko du.

Horrela, α balioa finkatzeko irizpide nagusiak seriearen egonkortasunari erreparatzen dio: seriea egonkorra bada eta aparteko aldaketarik izaten ez badu, α txikia finkatu behar da. Berriz, seriearen parametroak maiz aldatzen badira, egokiagoa izango da α balio handia finkatzea, aldaketei azkar erantzuteko aurresanetan. Balioa noiz den txikia edo handia zehaztearren, α parametroarako gehienetan erabiltzen diren balioak 0.01-0.30 bitartekoak izaten direla aipatu da, 0.1 balioaren inguruko parametroa ohikoa izanik. [5]

Irizpide orokor horiek zehaztearren, serieko ex post edo atzerako aurresanetarako doitasun-neurriak erabili dira leunketa-parametroak finkatzeko, hala nola batezbesteko errore koadratikoa,

\text{BEK}=\frac{\sum_te_t^2}{n}=\frac{\sum_t(x_t-\hat{x}_t)^2}{n} \,,

eta batezbesteko errore absolutua

\text{BEA}=\frac{\sum_t|e_t|}{n}=\frac{\sum_t|x_t-\hat{x}_t|}{n} \,.

Hala ere, nabarmendu behar da neurri hauek minimotuz eskuratzen diren α parametro balioak ex post aurresanen doitasuna soilik neurtzen dutela eta ondorioz ez dutela ziurtatzen ex ante edo ondorengo aurresanen doitasuna.[6][1]

Beste alde batetik, froga daiteke leunketa esponentziala α=1-d parametroa S=\sum_{j=0}^{\infty}d^{j+1}(x_{t-j}-s_t)^2 adierazpena minimotzen duen balioa dela.[7]

Adibidea[aldatu | aldatu iturburu kodea]

Honako taula honetan serie historiko bateko datuetarako bi leunketa esponentzial, leunketa parametro banarekin, zehazten dira, dagozkien aurresanekin batera. Eskubian, leunketaren abiarazpena eta lehenengo leunketa-balioen eta aurresanen kalkulua azaltzen da:

Denbora

(t)

Serie historikoa

(x_t)

\alpha=0.3 \alpha=0.8
Leunketa balioak

(s_t)

Aurresanak

(\hat{x}_t)

Leunketa balioak

(s_t)

Aurresanak

(\hat{x}_t)

0 14 14 - 14 -
1 12 13.40 14 12.40 14
2 13 13.28 13.40 12.88 12.40
3 15 13.80 13.28 14.58 12.88
4 14 13.86 13.80 14.12 14.58
5 22 16.30 13.86 20.42 14.12
6 13 15.31 16.30 14.48 20.42
7 14 14.92 15.31 14.10 14.48
8 12 14.04 14.92 12.42 14.10
9 (etorkizuna) - - 14.04 - 12.42
10 (etorkizuna) - - 14.04 - 12.42
\cdots - - 14.04 - 12.42


\text{Abiarazpena:}\ s_0=x_0=14

\text{Parametroa:}\ \alpha=0.3


s_1=\alpha x_1 + (1-\alpha)s_0=0.3 \times 12 + 0.7 \times 14=13.40;

\hat{x}_1=s_0=14


s_2=\alpha x_2 + (1-\alpha)s_1=0.3 \times 13 + 0.7 \times 13.40=13.28;

\hat{x}_2=s_1=13.40


s_3=\alpha x_3 + (1-\alpha)s_2=0.3 \times 15 + 0.7 \times 13.28=13.80;

\hat{x}_3=s_2=13.28


\cdots\cdots\cdots\cdots\cdots\cdots\cdots\cdots

Serie historikoa (beltzez) eta horren gainean egindako leunketa esponentziala bi, parametro desberdinekin. Berdez agertzen den leunketa α parametro txiki batekin dago, azken balioei haztapen gutxiago ematen diena; horregatik, muturrekotzat jo daitekeen 22 balioa gertatzean, ondorengo leunketa balioek ez dute horren eragina nabarmenki jasotzen. Gorriz agertzen den leunketa α parametro handi batekin burutu da eta beraz, azkenik gertatzen diren balioak nabarmenago jasotzen ditu, 22 balioaren kasuan ikus daitekeen bezala, baina joera ez du hain garbiro azalarazten. Beste alde batetik, leunketa esponentzialak serie historikoarekiko duen atzerapena ere ikus daiteke.

Leunketa moldatzaileak[aldatu | aldatu iturburu kodea]

Oinarrizko formuletan α parametroa konstantea bada ere, serie historikoaren bilakaerara eta aldiro sortutako aurresanen errorera egokitu edo moldatu asmoz parametroaren balioa aldi batetik bestera aldatzen duten leunketa modatzaileak ere proposatu dira. α parametro moldatzaileak ezartzeko irizpideak arrazoizkoak izan arren, leunketa moldatzaileak ez dira beti parametro finkoko leunketak baino hobeak. Areago, kasu batzuetan aurresanetan desegonkortasunak sortzen dituztela egiaztatu da.[8]

Gehien erabiltzen den leunketa moldatzaileak parametro finkoko leunketaren formulak erabiltzen ditu, baina α parametro finkoaren ordez aldiro aldatzen den αt parametroa ezarriz:[9]

\alpha_t=\Bigg\vert\frac{A_t}{M_t}\Bigg\vert
A_t=\phi e_t+(1-\phi)A_{t-1}
M_t=\phi |e_t|+(1-\phi)M_{t-1}

\phi parametroaren ohiko balioa 0.2 da baina erabiltzaileak beste balio bat finka dezake.

Moldaketak aurresanen erorrearen bilakaera eta horretan sor daitezkeen alborapena edo errore sistematikoa hautematen duten seinaleen jarraipen batean oinarritzen dira gehienetan: arestiko prozeduran A_t balioak dira aurresanetan burutzen den errore sistematikoa adierazten dutenak.

Leunketa esponentzialerako metodo estandarrak[aldatu | aldatu iturburu kodea]

Leunketa esponentzial sinplea joerarik gabeko serieetarako erabiltzen da soilik. Horretaz gainera, joera eta urtarokotasuna kontuan hartzen duten leunketa metodoak ere garatu dira. Metodo guztiak bateratzen dituen nomenklatura hau garatu da, leunketak barnehartzen dituen denbora-osagaien eta horien izaeraren arabera:[10]

Joera osagaia Urtarokotasuna
Ez ((Ingelesez) none) Batukorra ((Ingelesez) additive) Biderkakorra ((Ingelesez) multiplicative)
Ez ((Ingelesez) none) NN NA NM
Batukorra ((Ingelesez) additive) AN AA AM
Biderkakorra ((Ingelesez) multiplicative) MN MA MM
Batukor moteldua ((Ingelesez) damped additive) DAN DAA DAM
Biderkakor moteldua ((Ingelesez) damped multiplicative) DMN DMA DMM

Joera batukorra denean, beste denbora osagaiekin duen elkarregintza modu batukorrean gertatzen dela (hau da, balio bat gehituz) esan nahi du. Biderkakorra denean, berriz, elkarregintza hau modu biderkakorrean gertatzen dela esan nahi du (balio bat bidertuz). Joera motelduko leunketa esponentziala joera iraunkorra ez denean eta aurresanak epe laburraz haraindi egiten direnean da egokia. Zehatzago, oraingo unetik aldendu ahala, gaineratzen den joera osagaia aldatzen da, \phi parametro baten arabera: 0<\phi<1 ezartzen bada, joera osagaia moteltzen doa; \phi>1 denean, berriz, joeraren hazkundea esponentziala da. Urtarokotasunaren batukortasuna eta biderkakortasuna era berean interpretatu behar da. Jarraian, metodo bakoitzeko formulak zehazten dira:[11]


Eredua Errepikapen-formula Errepikapen-formula Aurresanak
NN S_t=\alpha X_t+(1-\alpha)S_{t-1} S_t=S_{t-1}+\alpha e_t \hat{X}_{t+k}=S_t
AN[12] S_t=\alpha X_t+(1-\alpha)(S_{t-1}+T_{t-1})

T_t=\gamma (S_{t}-S_{t-1})+(1-\gamma)T_{t-1}

S_t=S_{t-1}+T_{t-1}+\alpha e_t

T_t=T_{t-1}+\alpha \gamma e_t

\hat{X}_{t+k}=S_t+kT_t
MN[12] S_t=\alpha X_t+(1-\alpha)(S_{t-1}R_{t-1})

R_t=\gamma (S_{t}/S_{t-1})+(1-\gamma)R_{t-1}

S_t=S_{t-1}R_{t-1}+\alpha e_t

R_t=(R_{t-1}+\alpha \gamma e_t)/S_{t-1}

\hat{X}_{t+k}=S_tR_t^k
DAN S_t=\alpha X_t+(1-\alpha)(S_{t-1}+\phi T_{t-1})

T_t=\gamma (S_{t}-S_{t-1})+(1-\gamma)\phi T_{t-1}

S_t=S_{t-1}+\phi T_{t-1}+\alpha e_t

T_t=\phi T_{t-1}+\alpha \gamma e_t

\hat{X}_{t+k}=S_t+\sum_{i=1}^k\phi^iT_t
DMN S_t=\alpha X_t+(1-\alpha)(S_{t-1}R_{t-1}^{\phi})

R_t=\gamma (S_{t}/S_{t-1})+(1-\gamma)R_{t-1}^{\phi}

S_t=S_{t-1}R_{t-1}^{\phi}+\alpha e_t

R_t=(R_{t-1}^{\phi}+\alpha \gamma e_t)/S_{t-1}

\hat{X}_{t+k}=S_tR_t^{\sum_{i=1}^k\phi^i}
NA S_t=\alpha(X_t-I_{t-p})+(1-\alpha)S_{t-1}

I_t=\delta (X_{t}-S_{t})+(1-\delta)I_{t-p}

S_t=S_{t-1}+\alpha e_t

I_t=I_{t-p}+\delta(1-\alpha)e_t

\hat{X}_{t+k}=S_t+I_{t-p+k}
AA S_t=\alpha(X_t-I_{t-p})+(1-\alpha)(S_{t-1}+T_{t-1})

T_t=\gamma (S_{t}-S_{t-1})+(1-\gamma)T_{t-1}

I_t=\delta (X_{t}-S_{t})+(1-\delta)I_{t-p}

S_t=S_{t-1}+T_{t-1}+\alpha e_t

T_t=T_{t-1}+\alpha \gamma e_t

I_t=I_{t-p}+\delta(1-\alpha)e_t

\hat{X}_{t+k}=S_t+kT_t+I_{t-p+k}
AM S_t=\alpha\Bigg(\frac{X_t}{I_{t-p}}\Bigg)+(1-\alpha)(S_{t-1}+T_{t-1})

T_t=\gamma (S_{t}-S_{t-1})+(1-\gamma)T_{t-1}

I_t=\delta \Bigg(\frac{X_{t}}{S_{t}}\Bigg)+(1-\delta)I_{t-p}

S_t=S_{t-1}+T_{t-1}+\frac{\alpha e_t}{I_{t-p}}

T_t=T_{t-1}+\frac{\alpha \gamma e_t}{I_{t-p}}

I_t=I_{t-p}+\frac{\delta(1-\alpha)e_t}{S_t}

\hat{X}_{t+k}=(S_t+kT_t)I_{t-p+k}
NM S_t=\alpha(X_t/I_{t-p})+(1-\alpha)S_{t-1}

I_t=\delta (X_{t}/S_{t})+(1-\delta)I_{t-p}

S_t=S_{t-1}+\alpha e_t/I_{t-p}

I_t=I_{t-p}+\delta(1-\alpha)e_t/S_t

\hat{X}_{t+k}=S_tI_{t-p+k}
MA S_t=\alpha(X_t-I_{t-p})+(1-\alpha)S_{t-1}R_{t-1}

I_t=\delta (X_{t}-S_{t})+(1-\delta)I_{t-p}

R_t=\gamma (S_{t}/S_{t-1})+(1-\gamma)R_{t-1}

S_t=S_{t-1}R_{t-1}+\alpha e_t

I_t=I_{t-p}+\delta(1-\alpha)e_t

R_t=R_{t-1}+\alpha \gamma e_t / S_{t-1}

\hat{X}_{t+k}=S_tR_t^m+I_{t-p+k}

Erreferentziak[aldatu | aldatu iturburu kodea]

  1. a b c (Ingelesez)   Gardner, Everette S. (1985), «Exponential Smoothing: The State of the Art», Journal of Forecasting, http://www.bauer.uh.edu/gardner/Exp-Sm-1985.pdf .
  2. (Ingelesez)   Silver, Edward A.; Pyke, David F.; Peterson, Rein (1998), Inventory Management and Production Planning and Scheduling, 91. orrialdea .
  3. (Ingelesez)   Gupta, Surendra, Adaptive Response Rate Exponential Smoothing, 24. orrialdea, http://es.scribd.com/doc/52194214/20/Adaptive-Response-Rate-Exponential-Smoothing .
  4. (Ingelesez)   Lee, Cheng F.; Lee, John C.; Lee, Alice C. (2000), Statistics for Business and Financial Economics, 871. orrialdea, http://books.google.es/books?id=hnT5hL5YdCYC&printsec=frontcover&hl=es&source=gbs_ge_summary_r&cad=0#v=onepage&q&f=false .
  5. (Ingelesez)   Silver, Edward A.; Pyke, David F.; Peterson, Rein (1998), op. cit., 106-107. orrialdeak .
  6. (Gaztelaniaz)   Castro Zuluaga, Carlos A.; Botero Escobar, Sara C. (2012), Metodología para la selección del parámetro alpha en el modelo de Suavización Exponencial: Un enfoque empírico, http://www.laccei.org/LACCEI2012-Panama/StudentPapers/SP130.pdf .
  7. (Ingelesez)   Silver, Edward A.; Pyke, David F.; Peterson, Rein (1998), op. cit., 134. orrialdea .
  8. (Ingelesez)   Silver, Edward A.; Pyke, David F.; Peterson, Rein (1998), op. cit., 120-121. orrialdeak .
  9. (Ingelesez)   Taylor, James W. (2004), «Smooth Transition Exponential Smoothing», Journal of Forecasting, http://users.ox.ac.uk/~mast0315/STES.pdf .
  10. (Ingelesez)   Hyndman, Rob J.; Koehler, Anne B.; Snyder, Ralph D.; Grose, Simone (2002), «A state space framework for automatic forecasting using exponential smoothing methods», International Journal of Forecasting, http://www.forecasters.org/pdfs/ijf/HKSG.pdf .
  11. Denbora serieari dagokion aldagaia eta horren osagaiak maiuskulaz azaltzen dira, aurreko formuletan ez bezala.
  12. a b Charles C. Holtek proposatutako eredua da ondorengoa. Robert Goodell Brownen bertsioan parametro bakarra erabiltzen da bi formuletan.

Kanpo loturak[aldatu | aldatu iturburu kodea]

Commonsen badira fitxategi gehiago, gai hau dutenak: Leunketa esponentziala