Home » date » 2010 » Dec » 29 »

paper (14)

*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: Wed, 29 Dec 2010 19:20:12 +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/29/t12936503917hg9xoxepvt6wbx.htm/, Retrieved Wed, 29 Dec 2010 20:19:55 +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/29/t12936503917hg9xoxepvt6wbx.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 «

9,3 14,2 17,3 23 16,3 18,4 14,2 9,1 5,9 7,2 6,8 8 14,3 14,6 17,5 17,2 17,2 14,1 10,4 6,8 4,1 6,5 6,1 6,3 9,3 16,4 16,1 18 17,6 14 10,5 6,9 2,8 0,7 3,6 6,7 12,5 14,4 16,5 18,7 19,4 15,8 11,3 9,7 2,9 0,1 2,5 6,7 10,3 11,2 17,4 20,5 17 14,2 10,6 6,1

Output produced by software:

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 seconds
R Server	'Herman Ole Andreas Wold' @ www.yougetit.org

Seasonal Decomposition by Loess - Parameters
Component	Window	Degree	Jump
Seasonal	561	0	57
Trend	19	1	2
Low-pass	13	1	2

Seasonal Decomposition by Loess - Time Series Components
t	Observed	Fitted	Seasonal	Trend	Remainder
1	9.3	7.27787551221527	0.0206500189642456	11.3014744688205	-2.02212448778473
2	14.2	13.8229943766453	3.05375084093265	11.5232547824221	-0.377005623354744
3	17.3	16.9881162911357	5.86684861284062	11.7450350960237	-0.311883708864313
4	23	25.6391023899078	8.42280575612615	11.938091853966	2.63910238990784
5	16.3	13.9900833196333	6.4787680684584	12.1311486119083	-2.30991668036671
6	18.4	20.1947915050885	4.30685145264669	12.2983570422648	1.79479150508853
7	14.2	15.4994976197142	0.434936907664548	12.4655654726212	1.29949761971422
8	9.1	8.84112869650511	-3.2451599649695	12.6040312684644	-0.258871303494889
9	5.9	6.29236914128705	-7.2348662055946	12.7424970643076	0.392369141287054
10	7.2	9.24471900437254	-7.52797701598902	12.6832580116165	2.04471900437254
11	6.8	7.3720690651996	-6.396088024125	12.6240189589254	0.572069065199594
12	8	7.81805337693193	-4.18052050067076	12.3624671237388	-0.18194662306807
13	14.3	16.4784346924835	0.0206500189642456	12.1009152885522	2.17843469248351
14	14.6	14.3114031059102	3.05375084093265	11.8348460531572	-0.288596894089807
15	17.5	17.5643745693973	5.86684861284062	11.5687768177621	0.0643745693973123
16	17.2	14.5745593233213	8.42280575612615	11.4026349205526	-2.62544067667874
17	17.2	16.6847389081985	6.4787680684584	11.2364930233431	-0.515261091801497
18	14.1	12.7594727321467	4.30685145264669	11.1336758152066	-1.34052726785329
19	10.4	9.33420448526537	0.434936907664548	11.0308586070701	-1.06579551473463
20	6.8	5.84079317473804	-3.2451599649695	11.0043667902315	-0.959206825261964
21	4.1	4.45699123220176	-7.2348662055946	10.9778749733928	0.356991232201757
22	6.5	9.51763148669913	-7.52797701598902	11.0103455292899	3.01763148669913
23	6.1	7.55327193893807	-6.396088024125	11.0428160851869	1.45327193893807
24	6.3	5.74605472505956	-4.18052050067076	11.0344657756112	-0.553945274940443
25	9.3	7.55323451500028	0.0206500189642456	11.0261154660355	-1.74676548499972
26	16.4	18.838333939098	3.05375084093265	10.9079152199694	2.43833393909797
27	16.1	15.5434364132561	5.86684861284062	10.7897149739033	-0.556563586743904
28	18	16.9786196581409	8.42280575612615	10.598574585733	-1.02138034185915
29	17.6	18.3137977339789	6.4787680684584	10.4074341975627	0.7137977339789
30	14	13.3642631472466	4.30685145264669	10.3288854001067	-0.635736852753443
31	10.5	10.3147264896847	0.434936907664548	10.2503366026508	-0.185273510315341
32	6.9	6.76913710495946	-3.2451599649695	10.27602286001	-0.130862895040536
33	2.8	2.53315708822532	-7.2348662055946	10.3017091173693	-0.266842911774685
34	0.7	-1.46205745039481	-7.52797701598902	10.3900344663838	-2.16205745039481
35	3.6	3.11772820872663	-6.396088024125	10.4783598153984	-0.482271791273366
36	6.7	6.96228908572193	-4.18052050067076	10.6182314149488	0.26228908572193
37	12.5	14.2212469665365	0.0206500189642456	10.7581030144993	1.72124696653647
38	14.4	14.8575247164862	3.05375084093265	10.8887244425812	0.457524716486169
39	16.5	16.1138055164963	5.86684861284062	11.0193458706631	-0.386194483503694
40	18.7	17.9590283877051	8.42280575612615	11.0181658561687	-0.740971612294887
41	19.4	21.3042460898672	6.4787680684584	11.0169858416744	1.90424608986722
42	15.8	16.3939873001345	4.30685145264669	10.8991612472188	0.593987300134518
43	11.3	11.3837264395723	0.434936907664548	10.7813366527632	0.0837264395722581
44	9.7	11.9919239800575	-3.2451599649695	10.653235984912	2.29192398005752
45	2.9	2.50973088853384	-7.2348662055946	10.5251353170608	-0.390269111466164
46	0.1	-2.72100995822562	-7.52797701598902	10.4489869742146	-2.82100995822562
47	2.5	1.02324939275649	-6.396088024125	10.3728386313685	-1.47675060724351
48	6.7	7.24593107457508	-4.18052050067076	10.3345894260957	0.545931074575082
49	10.3	10.2830097602129	0.0206500189642456	10.2963402208228	-0.0169902397870878
50	11.2	9.10086323633162	3.05375084093265	10.2453859227357	-2.09913676366838
51	17.4	18.7387197625108	5.86684861284062	10.1944316246486	1.33871976251077
52	20.5	22.4136462725701	8.42280575612615	10.1635479713038	1.91364627257008
53	17	17.3885676135827	6.4787680684584	10.1326643179589	0.388567613582682
54	14.2	13.9760947874313	4.30685145264669	10.117053759922	-0.223905212568688
55	10.6	10.6636198904504	0.434936907664548	10.1014432018851	0.0636198904503864
56	6.1	5.3556535490445	-3.2451599649695	10.089506415925	-0.744346450955499

Charts produced by software:

http://www.freestatistics.org/blog/date/2010/Dec/29/t12936503917hg9xoxepvt6wbx/1dn131293650411.png (open in new window)

http://www.freestatistics.org/blog/date/2010/Dec/29/t12936503917hg9xoxepvt6wbx/1dn131293650411.ps (open in new window)

http://www.freestatistics.org/blog/date/2010/Dec/29/t12936503917hg9xoxepvt6wbx/2dn131293650411.png (open in new window)

http://www.freestatistics.org/blog/date/2010/Dec/29/t12936503917hg9xoxepvt6wbx/2dn131293650411.ps (open in new window)

http://www.freestatistics.org/blog/date/2010/Dec/29/t12936503917hg9xoxepvt6wbx/35w0o1293650411.png (open in new window)

http://www.freestatistics.org/blog/date/2010/Dec/29/t12936503917hg9xoxepvt6wbx/35w0o1293650411.ps (open in new window)

http://www.freestatistics.org/blog/date/2010/Dec/29/t12936503917hg9xoxepvt6wbx/45w0o1293650411.png (open in new window)

http://www.freestatistics.org/blog/date/2010/Dec/29/t12936503917hg9xoxepvt6wbx/45w0o1293650411.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 <- as.numeric(par1) #seasonal period

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')

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.

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