Home » date » 2010 » Dec » 16 »

Classical decomposition

*The author of this computation has been verified*
R Software Module: /rwasp_decompose.wasp (opens new window with default values)
Title produced by software: Classical Decomposition
Date of computation: Thu, 16 Dec 2010 15:47:36 +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/16/t1292514372t8o1qzv9aud9lod.htm/, Retrieved Thu, 16 Dec 2010 16:46:13 +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/16/t1292514372t8o1qzv9aud9lod.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 «
15561600 14917500 14805920 16958000 17605000 17131200 18474600 17286700 18574400 18056000 19701600 19061700 19681900 34521200 19922700 20177900 19759900 23076700 22532000 22029400 22587000 23256600 22680300 21916400 19640200 18813100 18730000 18154700 17848800 18077500 17133100 16602600 15878900 15789100 15422000 14661400 15879200 14339300 13169600 14528900 13375800 12309900 11933900 10061900 12609600 11156500 12187200 11284300 10177000 10970720 10820680 11492390 14573750 13992820 14727070 15685360 16736210 17950180 17002730 17415160 17929810 17865790 19202360 19085000 18188880 18466410 18520400 20025500 20636100 20672000 22589100 21864800 22750100 22548746 21325495 21556563 21415269 20401054 19062253 19085706 19279967 18552045 17800733 17142490 17593173 17633859 17336613 17008347 17951965 14520929 16941217 15436824 14744261 14248004 11540953 12881661 15185757 13554339 13575106 12238400 13303614 14151478 1417 etc...
 
Output produced by software:


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


Classical Decomposition by Moving Averages
tObservationsFitTrendSeasonalRandom
115561600NANA0.900288132417252NA
214917500NANA0.863765434895422NA
314805920NANA0.907841152495603NA
416958000NANA0.965843848011875NA
517605000NANA1.02986775797667NA
617131200NANA1.10379595379534NA
718474600NANA1.04416963049344NA
817286700NANA1.06729684004503NA
918574400NANA1.09645404680177NA
1018056000NANA1.08985937772229NA
1119701600NANA1.16842145088156NA
1219061700NANA1.15977935690514NA
1319681900NANA1.10282328990127NA
1434521200NANA1.1166612212153NA
1519922700NANA1.1184734102911NA
1620177900NANA1.2070943792417NA
1719759900NANA1.2428079398952NA
1823076700NANA1.24325565288817NA
1922532000NANA1.35440701915374NA
2022029400NANA1.30730071355346NA
2122587000NANA1.35733431410796NA
2223256600NANA1.3435497453056NA
2322680300NANA1.26775501927876NA
2421916400NANA1.27697403913775NA
2519640200NANA1.26398949541412NA
2618813100NANA1.2001969497807NA
271873000019088645.916518416977511.251.124348889270710.981211558007471
281815470018502329.953049516959121.53846151.090995775405470.981211558007472
291784880018190572.516537916949472.21153851.073223536963850.981211558007472
301807750018423651.711471516936477.11538461.087808969123570.981211558007471
311713310017461168.144810616915886.44230771.032234887858970.981211558007472
321660260016920510.021013816915407.51.000301649310770.981211558007472
331587890016182952.463630816909129.71153850.957054132276710.981211558007472
341578910016091432.95464516896212.30769230.952369244751920.981211558007472
351542200015717303.647866916891249.51923080.930499761428110.981211558007472
361466140014942139.521646716883222.59615380.8850288762437270.981211558007472
371587920016183258.208092916876593.26923080.958917356715470.981211558007472
381433930014613871.884182216872016.92307690.8661603381984450.981211558007472
391316960013421774.2264912168578850.7961718938343240.981211558007472
401452890014807102.384223416689156.05769230.8872289487279720.981211558007472
411337580013631922.586768116521299.42307690.8251120107252030.981211558007472
421230990012545612.512960516506350.38461540.760047631404550.981211558007472
431193390012162412.7871404165133100.7365217989089030.981211558007472
441006190010254567.343695516498612.88461540.6215411813957810.981211558007471
451260960012851051.230588916476039.80769230.7799842304695730.981211558007472
461115650011370126.971043116475006.15384620.6901440196663430.981211558007472
471218720012420563.027965516474991.73076920.7539040523321430.981211558007472
481128430011500374.111893716469753.71153850.6982723793760620.981211558007471
491017700010371871.302317616449920.45192310.6305119427556280.981211558007472
501097072011180789.617152516433433.50961540.6803684458643730.981211558007472
511082068011027876.62017916447041.50961540.6705082256727930.981211558007472
521149239011712448.662281716479378.26923080.7107336497124090.981211558007471
531457375014852811.181305916497841.79807690.9002881324172520.981211558007472
541399282014260757.413431816509988.51923080.8637654348954220.981211558007472
551472707015009066.984398316532701.72115380.9078411524956030.981211558007472
561568536015985706.519654116551025.8750.9658438480118750.981211558007472
571673621017056678.412946916562008.35576921.029867757976670.981211558007471
581795018018293893.761760316573619.15384621.103795953795340.981211558007471
591700273017328301.793068116595293.79807691.044169630493440.981211558007471
601741516017748629.088068116629515.25961541.067296840045030.981211558007471
611792981018273133.712784416665663.06730771.096454046801770.981211558007472
621786579018207887.844574216706639.60576921.089859377722290.981211558007472
631920236019570050.76356216749136.83653851.168421450881560.981211558007472
641908500019450443.529992216770813.70192311.159779356905140.981211558007472
651818888018537164.438763716808825.68269231.102823289901270.981211558007472
661846641018820008.640643616853821.26923081.11666122121530.981211558007471
671852040018875032.452337816875709.58653851.11847341029110.981211558007471
682002550020408952.41864616907503.48076921.20709437924170.981211558007472
692063610021031244.313820916922360.75961541.24280793989520.981211558007472
702067200021067831.734451116945695.50961541.243255652888170.981211558007472
712258910023021640.762030216997579.33653851.354407019153740.981211558007472
722186480022283471.71572317045406.22115381.307300713553460.981211558007472
732275010023185723.623352117081807.61538461.357334314107960.981211558007472
742254874622980514.055286217104326.90384621.34354974530560.981211558007472
752132549521733839.991963917143564.53846151.267755019278760.981211558007471
762155656321969332.529851717204212.34615381.276974039137750.981211558007472
772141526921825333.003096317267020.8751.263989495414120.981211558007472
782040105420791697.604365817323571.44230771.20019694978070.981211558007471
7919062253NANA1.12434888927071NA
8019085706NANA1.09099577540547NA
8119279967NANA1.07322353696385NA
8218552045NANA1.08780896912357NA
8317800733NANA1.03223488785897NA
8417142490NANA1.00030164931077NA
8517593173NANA0.95705413227671NA
8617633859NANA0.95236924475192NA
8717336613NANA0.93049976142811NA
8817008347NANA0.885028876243727NA
8917951965NANA0.95891735671547NA
9014520929NANA0.866160338198445NA
9116941217NANA0.796171893834324NA
9215436824NANA0.887228948727972NA
9314744261NANA0.825112010725203NA
9414248004NANA0.76004763140455NA
9511540953NANA0.736521798908903NA
9612881661NANA0.621541181395781NA
9715185757NANA0.779984230469573NA
9813554339NANA0.690144019666343NA
9913575106NANA0.753904052332143NA
10012238400NANA0.698272379376062NA
10113303614NANA0.630511942755628NA
10214151478NANA0.680368445864373NA
10314172009NANA0.670508225672793NA
10414022320NANA0.710733649712409NA
 
Charts produced by software:
http://www.freestatistics.org/blog/date/2010/Dec/16/t1292514372t8o1qzv9aud9lod/111fc1292514452.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/16/t1292514372t8o1qzv9aud9lod/111fc1292514452.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/16/t1292514372t8o1qzv9aud9lod/2usfx1292514452.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/16/t1292514372t8o1qzv9aud9lod/2usfx1292514452.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/16/t1292514372t8o1qzv9aud9lod/3usfx1292514452.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/16/t1292514372t8o1qzv9aud9lod/3usfx1292514452.ps (open in new window)


http://www.freestatistics.org/blog/date/2010/Dec/16/t1292514372t8o1qzv9aud9lod/4mkei1292514452.png (open in new window)
http://www.freestatistics.org/blog/date/2010/Dec/16/t1292514372t8o1qzv9aud9lod/4mkei1292514452.ps (open in new window)


 
Parameters (Session):
par1 = multiplicative ; par2 = 1 ;
 
Parameters (R input):
par1 = multiplicative ; par2 = 4 ;
 
R code (references can be found in the software module):
par2 = 52
par2 <- as.numeric(par2)
x <- ts(x,freq=par2)
m <- decompose(x,type=par1)
m$figure
bitmap(file='test1.png')
plot(m)
dev.off()
mylagmax <- length(x)/2
bitmap(file='test2.png')
op <- par(mfrow = c(2,2))
acf(as.numeric(x),lag.max = mylagmax,main='Observed')
acf(as.numeric(m$trend),na.action=na.pass,lag.max = mylagmax,main='Trend')
acf(as.numeric(m$seasonal),na.action=na.pass,lag.max = mylagmax,main='Seasonal')
acf(as.numeric(m$random),na.action=na.pass,lag.max = mylagmax,main='Random')
par(op)
dev.off()
bitmap(file='test3.png')
op <- par(mfrow = c(2,2))
spectrum(as.numeric(x),main='Observed')
spectrum(as.numeric(m$trend[!is.na(m$trend)]),main='Trend')
spectrum(as.numeric(m$seasonal[!is.na(m$seasonal)]),main='Seasonal')
spectrum(as.numeric(m$random[!is.na(m$random)]),main='Random')
par(op)
dev.off()
bitmap(file='test4.png')
op <- par(mfrow = c(2,2))
cpgram(as.numeric(x),main='Observed')
cpgram(as.numeric(m$trend[!is.na(m$trend)]),main='Trend')
cpgram(as.numeric(m$seasonal[!is.na(m$seasonal)]),main='Seasonal')
cpgram(as.numeric(m$random[!is.na(m$random)]),main='Random')
par(op)
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Classical Decomposition by Moving Averages',6,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'t',header=TRUE)
a<-table.element(a,'Observations',header=TRUE)
a<-table.element(a,'Fit',header=TRUE)
a<-table.element(a,'Trend',header=TRUE)
a<-table.element(a,'Seasonal',header=TRUE)
a<-table.element(a,'Random',header=TRUE)
a<-table.row.end(a)
for (i in 1:length(m$trend)) {
a<-table.row.start(a)
a<-table.element(a,i,header=TRUE)
a<-table.element(a,x[i])
if (par1 == 'additive') a<-table.element(a,m$trend[i]+m$seasonal[i]) else a<-table.element(a,m$trend[i]*m$seasonal[i])
a<-table.element(a,m$trend[i])
a<-table.element(a,m$seasonal[i])
a<-table.element(a,m$random[i])
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.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