Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_decompose.wasp
Title produced by softwareClassical Decomposition
Date of computationSun, 14 Dec 2014 20:14:31 +0000
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2014/Dec/14/t14185880996g2507wdarsvo33.htm/, Retrieved Thu, 16 May 2024 23:19:02 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=267840, Retrieved Thu, 16 May 2024 23:19:02 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact106
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Classical Decomposition] [KUL paper CD vrouw] [2014-12-14 20:06:08] [bb1b6762b7e5624d262776d3f7139d34]
-    D  [Classical Decomposition] [KUL paper cd 3 te...] [2014-12-14 20:11:05] [bb1b6762b7e5624d262776d3f7139d34]
- R  D      [Classical Decomposition] [KUL paper CD vrouw 2] [2014-12-14 20:14:31] [8568a324fefbb8dbb43f697bfa8d1be6] [Current]
- R           [Classical Decomposition] [KUL paper CD tes...] [2014-12-14 20:23:08] [bb1b6762b7e5624d262776d3f7139d34]
- R             [Classical Decomposition] [KUL paper CD tes...] [2014-12-14 20:25:37] [bb1b6762b7e5624d262776d3f7139d34]
- RM            [Spectral Analysis] [KUL paper Spectra...] [2014-12-14 20:30:41] [bb1b6762b7e5624d262776d3f7139d34]
- R               [Spectral Analysis] [KUL paper spectra...] [2014-12-14 20:36:27] [bb1b6762b7e5624d262776d3f7139d34]
- RM              [Classical Decomposition] [KUL paper CD test ] [2014-12-14 20:40:24] [bb1b6762b7e5624d262776d3f7139d34]
Feedback Forum

Post a new message
Dataseries X:
7.5
NA
6.5
NA
NA
NA
8.5
NA
NA
NA
NA
NA
NA
5
2.5
5
NA
3.5
NA
4
NA
NA
4.5
NA
NA
NA
NA
7
NA
5.5
2.5
5.5
NA
NA
NA
NA
4.5
NA
NA
5
NA
NA
NA
NA
NA
4.5
NA
NA
7.5
NA
NA
NA
NA
NA
NA
NA
0
NA
3.5
NA
NA
6
1.5
NA
3.5
NA
4
NA
NA
6
5
5.5
3.5
NA
6.5
6.5
NA
7
3.5
NA
4
7.5
4.5
NA
3.5
NA
NA
4.5
2.5
7.5
NA
NA
NA
3
NA
3.5
NA
NA
NA
NA
NA
4.5
NA
NA
NA
2.5
7
0
1
3.5
5.5
NA
NA
NA
NA
NA
8.5
NA
NA
10
NA
8.5
9
NA
NA
NA
NA
NA
NA
NA
NA
7.5
NA
NA
NA
NA
NA
NA
9
NA
NA
NA
NA
NA
8
9
NA
7
5.5
NA
2
NA
NA
8.5
NA
NA
NA
9
7.5
6
10.5
NA
8
NA
10.5
NA
9.5
NA
7.5
5
NA
10
NA
NA
NA
NA
NA
10
NA
3
6
7
NA
7
NA
8
10
5.5
6
NA
NA
NA
NA
9.5
8
NA
5.5
7
9
8
NA
NA
6
8
NA
9
NA
NA
9.5
NA
NA
NA
5
7
8
NA
NA
NA
NA
8
8.5
3.5
NA
NA
10.5
8.5
8
NA
NA
9.5
9
10
NA
NA
NA
NA
6.5
NA
NA
NA
6
4
NA
10.5
NA
NA
8.5
NA
7
NA
NA
5
NA
8.5
NA
9.5
NA
1.5
6
NA
NA
7.5
NA
NA
9
NA
8.5
7
NA
NA
9.5
NA
8
9.5
NA
8
9




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time0 seconds
R Server'Sir Ronald Aylmer Fisher' @ fisher.wessa.net
R Engine error message
Error in decompose(x, type = par1) : 
  time series has no or less than 2 periods
Execution halted

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input & view raw input (R code)  \tabularnewline
Raw Output & view raw output of R engine  \tabularnewline
Computing time & 0 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ fisher.wessa.net \tabularnewline
R Engine error message & 
Error in decompose(x, type = par1) : 
  time series has no or less than 2 periods
Execution halted
\tabularnewline \hline \end{tabular} %Source: https://freestatistics.org/blog/index.php?pk=267840&T=0

[TABLE]
[ROW][C]Summary of computational transaction[/C][/ROW]
[ROW][C]Raw Input[/C][C]view raw input (R code) [/C][/ROW]
[ROW][C]Raw Output[/C][C]view raw output of R engine [/C][/ROW]
[ROW][C]Computing time[/C][C]0 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Ronald Aylmer Fisher' @ fisher.wessa.net[/C][/ROW]
[ROW][C]R Engine error message[/C][C]
Error in decompose(x, type = par1) : 
  time series has no or less than 2 periods
Execution halted
[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=267840&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=267840&T=0

As an alternative you can also use a QR Code:  

The GUIDs for individual cells are displayed in the table below:

Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time0 seconds
R Server'Sir Ronald Aylmer Fisher' @ fisher.wessa.net
R Engine error message
Error in decompose(x, type = par1) : 
  time series has no or less than 2 periods
Execution halted



Parameters (Session):
par1 = additive ; par2 = 1 ;
Parameters (R input):
par1 = additive ; par2 = 1 ;
R code (references can be found in the software module):
par2 <- '1'
par1 <- 'additive'
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,signif(m$trend[i]+m$seasonal[i],6)) else a<-table.element(a,signif(m$trend[i]*m$seasonal[i],6))
a<-table.element(a,signif(m$trend[i],6))
a<-table.element(a,signif(m$seasonal[i],6))
a<-table.element(a,signif(m$random[i],6))
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')