Home » date » 2010 » Dec » 10 »

Blog2

*The author of this computation has been verified*
R Software Module: /rwasp_decomposeloess.wasp (opens new window with default values)
Title produced by software: Decomposition by Loess
Date of computation: Fri, 10 Dec 2010 19:49:25 +0000
 
Cite this page as follows:
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2010/Dec/10/t1292010632x22awpyur4ygr03.htm/, Retrieved Fri, 10 Dec 2010 20:50:38 +0100
 
BibTeX entries for LaTeX users:
@Manual{KEY,
    author = {{YOUR NAME}},
    publisher = {Office for Research Development and Education},
    title = {Statistical Computations at FreeStatistics.org, URL http://www.freestatistics.org/blog/date/2010/Dec/10/t1292010632x22awpyur4ygr03.htm/},
    year = {2010},
}
@Manual{R,
    title = {R: A Language and Environment for Statistical Computing},
    author = {{R Development Core Team}},
    organization = {R Foundation for Statistical Computing},
    address = {Vienna, Austria},
    year = {2010},
    note = {{ISBN} 3-900051-07-0},
    url = {http://www.R-project.org},
}
 
Original text written by user:
 
IsPrivate?
No (this computation is public)
 
User-defined keywords:
 
Dataseries X:
» Textbox « » Textfile « » CSV «
42.33600 42.14710 40.25640 39.18980 39.13170 38.15070 38.27070 39.13350 40.12190 41.28450 42.57490 43.90190 43.18350 43.61880 44.76240 45.19720 44.38810 43.55520 43.56780 44.21350 45.14510 45.80790 42.32820 37.89990 34.79640 35.21440 36.37270 36.25020 36.82610 36.77230 36.90420 37.04940 36.82590 36.13570 36.03000 35.79270 35.91740 35.40080 35.17230 34.92110 35.02920 34.77390 34.89990 34.90540 34.56800 34.40600 34.45780 34.73160 34.26020 33.88490 34.05490 34.27550 34.13930 34.15870 34.53860 33.79870 33.49730 33.68020 34.32840 34.15380 33.91840 34.32620 34.77500 35.01190 34.55130 34.69510 35.47300 35.97940 36.47890 36.39100 36.67040 37.41620 37.11850 36.30010 35.70020 35.58590 35.67770 35.24080 34.80160 34.43890 34.98810 36.06800 36.35660 36.12540 34.87100 35.21990 34.33900 33.82340 34.51990 35.53000 35.79660 33.84840 33.98710 34.11610 33.82350 32.45400 31.87740 31.11500 31.04360 30.88680 31.30010 30.05070 28.6 etc...
 
Output produced by software:


Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time3 seconds
R Server'RServer@AstonUniversity' @ vre.aston.ac.uk


Seasonal Decomposition by Loess - Parameters
ComponentWindowDegreeJump
Seasonal14310144
Trend911
Low-pass511


Seasonal Decomposition by Loess - Time Series Components
tObservedFittedSeasonalTrendRemainder
142.33642.31278857737760.072249927957877742.2869614946645-0.0232114226224027
242.147142.76570601350870.073585951422975441.45490803506830.618606013508703
340.256439.9118835280161-0.054865476521413640.6557819485053-0.344516471983887
439.189838.5614333062827-0.063264012312790639.8814307060301-0.62836669371729
539.131739.1847998651201-0.027706440574196439.10630657545410.0530998651200747
638.150737.3925589496880.072249927957877738.8365911223542-0.758141050312034
738.270737.49064635823510.073585951422975438.9771676903419-0.78005364176488
839.133538.8487892463857-0.054865476521413639.4730762301357-0.28471075361432
940.121939.9833997859176-0.063264012312790640.3236642263952-0.13850021408242
1041.284541.2643159571442-0.027706440574196441.33239048343-0.0201840428557958
1142.574942.86475621995980.072249927957877742.21279385208230.289856219959844
1243.901944.79289973614660.073585951422975442.93731431243040.890999736146611
1343.183542.8544166404819-0.054865476521413643.5674488360395-0.329083359518073
1443.618843.2751236168803-0.063264012312790644.0257403954325-0.343676383119686
1544.762445.3039467436662-0.027706440574196444.2485596969080.541546743666174
1645.197246.00149358836030.072249927957877744.32065648368180.804293588360295
1744.388144.42576167420170.073585951422975444.27685237437530.0376616742017077
1843.555242.9854332363591-0.054865476521413644.1798322401623-0.569766763640871
1943.567842.9893023258767-0.063264012312790644.2095616864361-0.57849767412327
2044.213544.1177176062186-0.027706440574196444.3369888343556-0.0957823937813771
2145.145146.13486613592130.072249927957877744.08308393612080.989766135921343
2245.807948.51579816707180.073585951422975443.02641588150522.70789816707184
2342.328243.4765371349983-0.054865476521413641.23472834152311.14833713499833
2437.899936.6935604958602-0.063264012312790639.1695035164526-1.20633950413979
2534.796432.2621080268624-0.027706440574196437.3583984137118-2.53429197313757
2635.214434.09099370317330.072249927957877736.2655563688689-1.12340629682674
2736.372736.64725629117280.073585951422975436.02455775740430.274556291172757
2836.250236.2675036683949-0.054865476521413636.28776180812650.0173036683948737
2936.826137.1276439388981-0.063264012312790636.58782007341470.301543938898121
3036.772336.8109857705189-0.027706440574196436.76132067005530.0386857705188817
3136.904236.91416888389470.072249927957877736.82198118814750.00996888389466477
3237.049437.27818743968310.073585951422975436.74702660889390.228787439683103
3336.825937.1244697223615-0.054865476521413636.58219575415990.298569722361464
3436.135735.9697209739234-0.063264012312790636.3649430383894-0.16597902607662
3536.0335.9690994123162-0.027706440574196436.118607028258-0.0609005876837756
3635.792735.63882504152040.072249927957877735.8743250305217-0.153874958479555
3735.917436.10221053098780.073585951422975435.65900351758920.184810530987846
3835.400835.396897211605-0.054865476521413635.4595682649164-0.00390278839498137
3935.172335.1404659554672-0.063264012312790635.2673980568456-0.031834044532836
4034.921134.7854787405446-0.027706440574196435.0844277000296-0.135621259455355
4135.029235.02245563015240.072249927957877734.9636944418897-0.00674436984756
4234.773934.57602422385210.073585951422975434.898189824725-0.197875776147946
4334.899935.0271504330823-0.054865476521413634.82751504343910.127250433082338
4434.905435.1404698579972-0.063264012312790634.73359415431560.23506985799721
4534.56834.5115022936343-0.027706440574196434.6522041469399-0.0564977063657253
4634.40634.15899130885870.072249927957877734.5807587631835-0.247008691141346
4734.457834.35889142814550.073585951422975434.4831226204315-0.0989085718544658
4834.731635.1429515494802-0.054865476521413634.37511392704120.411351549480194
4934.260234.2945940692097-0.063264012312790634.28906994310310.0343940692096751
5033.884933.5863519914617-0.027706440574196434.2111544491125-0.298548008538297
5134.054933.90084058538040.072249927957877734.1367094866617-0.154059414619567
5234.275534.34089105685310.073585951422975434.1365229917240.065391056853052
5334.139334.1400049324688-0.054865476521413634.19346054405260.000704932468764241
5434.158734.2113402701552-0.063264012312790634.16932374215760.05264027015523
5534.538635.0511875036286-0.027706440574196434.05371893694560.512587503628559
5633.798733.56801972098730.072249927957877733.9571303510548-0.230680279012667
5733.497333.00090081129570.073585951422975433.9201132372813-0.496399188704309
5833.680233.5065477283203-0.054865476521413633.9087177482011-0.173652271679714
5934.328434.7591433183714-0.063264012312790633.96092069394140.430743318371391
6034.153834.2383942838317-0.027706440574196434.09691215674250.0845942838316631
6133.918433.49801495635260.072249927957877734.2665351156896-0.420385043647443
6234.326234.16300434755490.073585951422975434.4158097010221-0.163195652445104
6334.77535.0574737167933-0.054865476521413634.54739175972810.282473716793326
6435.011935.3865488546183-0.063264012312790634.70051515769450.374648854618265
6534.551334.2382070234051-0.027706440574196434.8920994171691-0.313092976594938
6634.695134.18618290010190.072249927957877735.1317671719402-0.508917099898127
6735.47335.42658846311430.073585951422975435.4458255854627-0.0464115368856852
6835.979436.1911660062963-0.054865476521413635.82249947022510.211766006296351
6936.478936.8126205640526-0.063264012312790636.20844344826020.333720564052634
7036.39136.2616970624147-0.027706440574196436.5480093781595-0.129302937585273
7136.670436.51233656467540.072249927957877736.7562135073667-0.158063435324557
7237.416237.98011547736450.073585951422975436.77869857121250.563915477364532
7337.118537.6431864670836-0.054865476521413636.64867900943780.524686467083626
7436.300136.2673145338049-0.063264012312790636.3961494785079-0.0327854661950582
7535.700235.373126188466-0.027706440574196436.0549802521082-0.327073811534042
7635.585935.38969682338650.072249927957877735.7098532486556-0.196203176613473
7735.677735.86455826020090.073585951422975435.41725578837610.186858260200879
7835.240835.3521398748504-0.054865476521413635.1843256016710.11133987485043
7934.801634.6065076445512-0.063264012312790635.0599563677616-0.19509235544885
8034.438933.7952013015963-0.027706440574196435.1103051389779-0.643698698403718
8134.988134.58080643275980.072249927957877735.3231436392823-0.407293567240202
8236.06836.48973962043750.073585951422975435.57267442813950.421739620437499
8336.356637.0570098239384-0.054865476521413635.7110556525830.700409823938429
8436.125436.6569067945178-0.063264012312790635.6571572177950.53150679451776
8534.87134.4239558776893-0.027706440574196435.3457505628849-0.447044122310665
8635.219935.45659413809480.072249927957877734.91095593394730.236694138094791
8734.33933.96552341618490.073585951422975434.6388906323921-0.373476583815126
8833.823433.0527870500059-0.054865476521413634.6488784265155-0.77061294999411
8934.519934.3555243902719-0.063264012312790634.7475396220409-0.164375609728097
9035.5336.3139216555183-0.027706440574196434.77378478505590.783921655518341
9135.796636.78178271684180.072249927957877734.73916735520030.985182716841827
9233.848433.03899895872340.073585951422975434.5842150898536-0.809401041276573
9333.987133.8041931140649-0.054865476521413634.2248723624565-0.182906885935118
9434.116134.5649209557291-0.063264012312790633.73054305658370.448820955729104
9533.823534.4418201397887-0.027706440574196433.23288630078550.6183201397887
9632.45432.16974704510980.072249927957877732.6660030269324-0.284252954890228
9731.877431.6317280625770.073585951422975432.049485986-0.245671937423008
9831.11530.7449730575058-0.054865476521413631.5398924190157-0.370026942494246
9931.043630.9459286121199-0.063264012312790631.2045354001929-0.0976713878800624
10030.886830.9351145839727-0.027706440574196430.86619185660150.0483145839726511
10131.300132.14917947811470.072249927957877730.37877059392740.849079478114724
10230.050730.30752370273020.073585951422975429.72029034584680.256823702730244
10328.679928.4779838225137-0.054865476521413628.9366816540077-0.201916177486268
10427.643827.2133708959647-0.063264012312790628.1374931163481-0.430429104035273
10527.239426.9825893197495-0.027706440574196427.5239171208247-0.256810680250538
10626.854926.48897370778010.072249927957877727.148576364262-0.365926292219875
10727.015827.02398449843480.073585951422975426.93402955014220.00818449843480806
10826.918827.0576924065435-0.054865476521413626.83477306997790.13889240654348
10926.496626.2725316979561-0.063264012312790626.7839323143567-0.224068302043928
11026.78526.9471536728714-0.027706440574196426.65055276770280.162153672871359
11126.839827.23164896514120.072249927957877726.3757011069010.391848965141168
11226.447826.77821017696550.073585951422975426.04380387161150.330410176965511
11325.172824.6581546791545-0.054865476521413625.7423107973669-0.5146453208455
11424.878424.2804200217292-0.063264012312790625.5396439905836-0.597979978270796
11525.401525.3043899297022-0.027706440574196425.5263165108719-0.0971100702977488
11625.771625.73401417455920.072249927957877725.7369358974829-0.0375858254407859
11726.118126.17728337709690.073585951422975425.98533067148010.0591833770968968
11826.396926.8680791537682-0.054865476521413625.98058632275320.471179153768194
11926.657127.7124214768121-0.063264012312790625.66504253550061.05532147681215
12025.183925.1998805207479-0.027706440574196425.19562591982630.0159805207478563
12123.839422.86131245386010.072249927957877724.745237618182-0.97808754613989
12223.861923.12993559261880.073585951422975424.5202784559582-0.73196440738122
12324.258123.8581642091803-0.054865476521413624.7129012673411-0.399935790819676
12425.109825.0932978689205-0.063264012312790625.1895661433923-0.0165021310794842
12526.161726.6770127246505-0.027706440574196425.67409371592370.515312724650471
12626.808727.39225524004950.072249927957877726.15289483199260.58355524004951
12725.657724.63197210722130.073585951422975426.6098419413557-1.02572789277869
12827.098227.2145868083865-0.054865476521413627.03667866813490.116386808386533
12927.466527.5537947927977-0.063264012312790627.44246921951510.087294792797703
13028.28928.8777309933091-0.027706440574196427.72797544726510.588730993309063
13128.693329.49853945906270.072249927957877727.81581061297940.805239459062744
13227.240126.65984380864020.073585951422975427.7467702399368-0.580256191359808
13327.279827.1287268991303-0.054865476521413627.4857385773911-0.151073100869713
13427.640428.2573338148066-0.063264012312790627.08673019750620.616933814806568
13526.83127.0236296045273-0.027706440574196426.66607683604690.19262960452734
13626.246726.15063708158030.072249927957877726.2705129904619-0.0960629184197401
13725.221224.34204705314480.073585951422975426.0267669954322-0.879152946855214
13825.365324.7663152059851-0.054865476521413626.0191502705363-0.598984794014893
13926.259226.4587692471786-0.063264012312790626.12289476513420.199569247178619
14027.227928.2642003650266-0.027706440574196426.21930607554761.03630036502661
14126.331526.23103161141060.072249927957877726.3597184606316-0.100468388589441
14225.898125.24352670437960.073585951422975426.4790873441975-0.65457329562043
14326.689826.8662888535133-0.054865476521413626.56817662300820.176488853513256
 
Charts produced by software:
http://www.freestatistics.org/blog/date/2010/Dec/10/t1292010632x22awpyur4ygr03/1zq301292010560.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/10/t1292010632x22awpyur4ygr03/1zq301292010560.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/10/t1292010632x22awpyur4ygr03/2rz2l1292010560.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/10/t1292010632x22awpyur4ygr03/2rz2l1292010560.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/10/t1292010632x22awpyur4ygr03/3rz2l1292010560.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/10/t1292010632x22awpyur4ygr03/3rz2l1292010560.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/10/t1292010632x22awpyur4ygr03/4k91o1292010560.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/10/t1292010632x22awpyur4ygr03/4k91o1292010560.ps (open in new window)


 
Parameters (Session):
par1 = 12 ; par2 = periodic ; par3 = 0 ; par5 = 1 ; par7 = 1 ; par8 = FALSE ;
 
Parameters (R input):
par1 = 12 ; par2 = periodic ; par3 = 0 ; par5 = 1 ; par7 = 1 ; par8 = FALSE ;
 
R code (references can be found in the software module):
par1 <- 5
if (par2 != 'periodic') par2 <- as.numeric(par2) #s.window
par3 <- as.numeric(par3) #s.degree
if (par4 == '') par4 <- NULL else par4 <- as.numeric(par4)#t.window
par5 <- as.numeric(par5)#t.degree
if (par6 != '') par6 <- as.numeric(par6)#l.window
par7 <- as.numeric(par7)#l.degree
if (par8 == 'FALSE') par8 <- FALSE else par9 <- TRUE #robust
nx <- length(x)
x <- ts(x,frequency=par1)
if (par6 != '') {
m <- stl(x,s.window=par2, s.degree=par3, t.window=par4, t.degre=par5, l.window=par6, l.degree=par7, robust=par8)
} else {
m <- stl(x,s.window=par2, s.degree=par3, t.window=par4, t.degre=par5, l.degree=par7, robust=par8)
}
m$time.series
m$win
m$deg
m$jump
m$inner
m$outer
bitmap(file='test1.png')
plot(m,main=main)
dev.off()
mylagmax <- nx/2
bitmap(file='test2.png')
op <- par(mfrow = c(2,2))
acf(as.numeric(x),lag.max = mylagmax,main='Observed')
acf(as.numeric(m$time.series[,'trend']),na.action=na.pass,lag.max = mylagmax,main='Trend')
acf(as.numeric(m$time.series[,'seasonal']),na.action=na.pass,lag.max = mylagmax,main='Seasonal')
acf(as.numeric(m$time.series[,'remainder']),na.action=na.pass,lag.max = mylagmax,main='Remainder')
par(op)
dev.off()
bitmap(file='test3.png')
op <- par(mfrow = c(2,2))
spectrum(as.numeric(x),main='Observed')
spectrum(as.numeric(m$time.series[!is.na(m$time.series[,'trend']),'trend']),main='Trend')
spectrum(as.numeric(m$time.series[!is.na(m$time.series[,'seasonal']),'seasonal']),main='Seasonal')
spectrum(as.numeric(m$time.series[!is.na(m$time.series[,'remainder']),'remainder']),main='Remainder')
par(op)
dev.off()
bitmap(file='test4.png')
op <- par(mfrow = c(2,2))
cpgram(as.numeric(x),main='Observed')
cpgram(as.numeric(m$time.series[!is.na(m$time.series[,'trend']),'trend']),main='Trend')
cpgram(as.numeric(m$time.series[!is.na(m$time.series[,'seasonal']),'seasonal']),main='Seasonal')
cpgram(as.numeric(m$time.series[!is.na(m$time.series[,'remainder']),'remainder']),main='Remainder')
par(op)
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Seasonal Decomposition by Loess - Parameters',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Component',header=TRUE)
a<-table.element(a,'Window',header=TRUE)
a<-table.element(a,'Degree',header=TRUE)
a<-table.element(a,'Jump',header=TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Seasonal',header=TRUE)
a<-table.element(a,m$win['s'])
a<-table.element(a,m$deg['s'])
a<-table.element(a,m$jump['s'])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Trend',header=TRUE)
a<-table.element(a,m$win['t'])
a<-table.element(a,m$deg['t'])
a<-table.element(a,m$jump['t'])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Low-pass',header=TRUE)
a<-table.element(a,m$win['l'])
a<-table.element(a,m$deg['l'])
a<-table.element(a,m$jump['l'])
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Seasonal Decomposition by Loess - Time Series Components',6,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'t',header=TRUE)
a<-table.element(a,'Observed',header=TRUE)
a<-table.element(a,'Fitted',header=TRUE)
a<-table.element(a,'Seasonal',header=TRUE)
a<-table.element(a,'Trend',header=TRUE)
a<-table.element(a,'Remainder',header=TRUE)
a<-table.row.end(a)
for (i in 1:nx) {
a<-table.row.start(a)
a<-table.element(a,i,header=TRUE)
a<-table.element(a,x[i])
a<-table.element(a,x[i]+m$time.series[i,'remainder'])
a<-table.element(a,m$time.series[i,'seasonal'])
a<-table.element(a,m$time.series[i,'trend'])
a<-table.element(a,m$time.series[i,'remainder'])
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable1.tab')
 





Copyright

Creative Commons License

This work is licensed under a Creative Commons Attribution-Noncommercial-Share Alike 3.0 License.

Software written by Ed van Stee & Patrick Wessa


Disclaimer

Information provided on this web site is provided "AS IS" without warranty of any kind, either express or implied, including, without limitation, warranties of merchantability, fitness for a particular purpose, and noninfringement. We use reasonable efforts to include accurate and timely information and periodically update the information, and software without notice. However, we make no warranties or representations as to the accuracy or completeness of such information (or software), and we assume no liability or responsibility for errors or omissions in the content of this web site, or any software bugs in online applications. Your use of this web site is AT YOUR OWN RISK. Under no circumstances and under no legal theory shall we be liable to you or any other person for any direct, indirect, special, incidental, exemplary, or consequential damages arising from your access to, or use of, this web site.


Privacy Policy

We may request personal information to be submitted to our servers in order to be able to:

  • personalize online software applications according to your needs
  • enforce strict security rules with respect to the data that you upload (e.g. statistical data)
  • manage user sessions of online applications
  • alert you about important changes or upgrades in resources or applications

We NEVER allow other companies to directly offer registered users information about their products and services. Banner references and hyperlinks of third parties NEVER contain any personal data of the visitor.

We do NOT sell, nor transmit by any means, personal information, nor statistical data series uploaded by you to third parties.

We carefully protect your data from loss, misuse, alteration, and destruction. However, at any time, and under any circumstance you are solely responsible for managing your passwords, and keeping them secret.

We store a unique ANONYMOUS USER ID in the form of a small 'Cookie' on your computer. This allows us to track your progress when using this website which is necessary to create state-dependent features. The cookie is used for NO OTHER PURPOSE. At any time you may opt to disallow cookies from this website - this will not affect other features of this website.

We examine cookies that are used by third-parties (banner and online ads) very closely: abuse from third-parties automatically results in termination of the advertising contract without refund. We have very good reason to believe that the cookies that are produced by third parties (banner ads) do NOT cause any privacy or security risk.

FreeStatistics.org is safe. There is no need to download any software to use the applications and services contained in this website. Hence, your system's security is not compromised by their use, and your personal data - other than data you submit in the account application form, and the user-agent information that is transmitted by your browser - is never transmitted to our servers.

As a general rule, we do not log on-line behavior of individuals (other than normal logging of webserver 'hits'). However, in cases of abuse, hacking, unauthorized access, Denial of Service attacks, illegal copying, hotlinking, non-compliance with international webstandards (such as robots.txt), or any other harmful behavior, our system engineers are empowered to log, track, identify, publish, and ban misbehaving individuals - even if this leads to ban entire blocks of IP addresses, or disclosing user's identity.


FreeStatistics.org is powered by