Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_variancereduction.wasp
Title produced by softwareVariance Reduction Matrix
Date of computationThu, 11 Dec 2014 09:15:50 +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/11/t1418289427zx31kcgeujtl3y0.htm/, Retrieved Thu, 16 May 2024 06:48:58 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=265644, Retrieved Thu, 16 May 2024 06:48:58 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact101
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Variance Reduction Matrix] [ARIMA n 1] [2014-12-11 09:15:50] [7de4f24d5c21ad7c83693f758b02221d] [Current]
Feedback Forum

Post a new message
Dataseries X:
12.9
12.2
12.8
7.4
6.7
12.6
14.8
13.3
11.1
8.2
11.4
6.4
10.6
12.0
6.3
11.3
11.9
9.3
9.6
10.0
6.4
13.8
10.8
13.8
11.7
10.9
16.1
13.4
9.9
11.5
8.3
11.7
9.0
9.7
10.8
10.3
10.4
12.7
9.3
11.8
5.9
11.4
13.0
10.8
12.3
11.3
11.8
7.9
12.7
12.3
11.6
6.7
10.9
12.1
13.3
10.1
5.7
14.3
8.0
13.3
9.3
12.5
7.6
15.9
9.2
9.1
11.1
13.0
14.5
12.2
12.3
11.4
8.8
14.6
12.6
13.0
12.6
13.2
9.9
7.7
10.5
13.4
10.9
4.3
10.3
11.8
11.2
11.4
8.6
13.2
12.6
5.6
9.9
8.8
7.7
9.0
7.3
11.4
13.6
7.9
10.7
10.3
8.3
9.6
14.2
8.5
13.5
4.9
6.4
9.6
11.6
11.1
4.35
12.7
18.1
17.85
16.6
12.6
17.1
19.1
16.1
13.35
18.4
14.7
10.6
12.6
16.2
13.6
18.9
14.1
14.5
16.15
14.75
14.8
12.45
12.65
17.35
8.6
18.4
16.1
11.6
17.75
15.25
17.65
16.35
17.65
13.6
14.35
14.75
18.25
9.9
16
18.25
16.85
14.6
13.85
18.95
15.6
14.85
11.75
18.45
15.9
17.1
16.1
19.9
10.95
18.45
15.1
15
11.35
15.95
18.1
14.6
15.4
15.4
17.6
13.35
19.1
15.35
7.6
13.4
13.9
19.1
15.25
12.9
16.1
17.35
13.15
12.15
12.6
10.35
15.4
9.6
18.2
13.6
14.85
14.75
14.1
14.9
16.25
19.25
13.6
13.6
15.65
12.75
14.6
9.85
12.65
19.2
16.6
11.2
15.25
11.9
13.2
16.35
12.4
15.85
18.15
11.15
15.65
17.75
7.65
12.35
15.6
19.3
15.2
17.1
15.6
18.4
19.05
18.55
19.1
13.1
12.85
9.5
4.5
11.85
13.6
11.7
12.4
13.35
11.4
14.9
19.9
11.2
14.6
17.6
14.05
16.1
13.35
11.85
11.95
14.75
15.15
13.2
16.85
7.85
7.7
12.6
7.85
10.95
12.35
9.95
14.9
16.65
13.4
13.95
15.7
16.85
10.95
15.35
12.2
15.1
17.75
15.2
14.6
16.65
8.1




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Sir Ronald Aylmer Fisher' @ fisher.wessa.net

\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 & 1 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ fisher.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=265644&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]1 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Ronald Aylmer Fisher' @ fisher.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=265644&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=265644&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 time1 seconds
R Server'Sir Ronald Aylmer Fisher' @ fisher.wessa.net







Variance Reduction Matrix
V(Y[t],d=0,D=0)11.5217Range15.6Trim Var.7.6577
V(Y[t],d=1,D=0)14.7492Range19.9Trim Var.9.53755
V(Y[t],d=2,D=0)43.4189Range33.55Trim Var.28.0486
V(Y[t],d=3,D=0)142.75Range62.65Trim Var.93.044
V(Y[t],d=0,D=1)17.353Range25.3Trim Var.10.3776
V(Y[t],d=1,D=1)31.162Range32.85Trim Var.20.36
V(Y[t],d=2,D=1)94.7259Range56.95Trim Var.59.7181
V(Y[t],d=3,D=1)318.797Range105.65Trim Var.201.334
V(Y[t],d=0,D=2)51.4767Range37.8Trim Var.32.3332
V(Y[t],d=1,D=2)94.639Range52.05Trim Var.61.6714
V(Y[t],d=2,D=2)291.002Range93.35Trim Var.184.164
V(Y[t],d=3,D=2)986.388Range169Trim Var.636.036

\begin{tabular}{lllllllll}
\hline
Variance Reduction Matrix \tabularnewline
V(Y[t],d=0,D=0) & 11.5217 & Range & 15.6 & Trim Var. & 7.6577 \tabularnewline
V(Y[t],d=1,D=0) & 14.7492 & Range & 19.9 & Trim Var. & 9.53755 \tabularnewline
V(Y[t],d=2,D=0) & 43.4189 & Range & 33.55 & Trim Var. & 28.0486 \tabularnewline
V(Y[t],d=3,D=0) & 142.75 & Range & 62.65 & Trim Var. & 93.044 \tabularnewline
V(Y[t],d=0,D=1) & 17.353 & Range & 25.3 & Trim Var. & 10.3776 \tabularnewline
V(Y[t],d=1,D=1) & 31.162 & Range & 32.85 & Trim Var. & 20.36 \tabularnewline
V(Y[t],d=2,D=1) & 94.7259 & Range & 56.95 & Trim Var. & 59.7181 \tabularnewline
V(Y[t],d=3,D=1) & 318.797 & Range & 105.65 & Trim Var. & 201.334 \tabularnewline
V(Y[t],d=0,D=2) & 51.4767 & Range & 37.8 & Trim Var. & 32.3332 \tabularnewline
V(Y[t],d=1,D=2) & 94.639 & Range & 52.05 & Trim Var. & 61.6714 \tabularnewline
V(Y[t],d=2,D=2) & 291.002 & Range & 93.35 & Trim Var. & 184.164 \tabularnewline
V(Y[t],d=3,D=2) & 986.388 & Range & 169 & Trim Var. & 636.036 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=265644&T=1

[TABLE]
[ROW][C]Variance Reduction Matrix[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=0)[/C][C]11.5217[/C][C]Range[/C][C]15.6[/C][C]Trim Var.[/C][C]7.6577[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=0)[/C][C]14.7492[/C][C]Range[/C][C]19.9[/C][C]Trim Var.[/C][C]9.53755[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=0)[/C][C]43.4189[/C][C]Range[/C][C]33.55[/C][C]Trim Var.[/C][C]28.0486[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=0)[/C][C]142.75[/C][C]Range[/C][C]62.65[/C][C]Trim Var.[/C][C]93.044[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=1)[/C][C]17.353[/C][C]Range[/C][C]25.3[/C][C]Trim Var.[/C][C]10.3776[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=1)[/C][C]31.162[/C][C]Range[/C][C]32.85[/C][C]Trim Var.[/C][C]20.36[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=1)[/C][C]94.7259[/C][C]Range[/C][C]56.95[/C][C]Trim Var.[/C][C]59.7181[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=1)[/C][C]318.797[/C][C]Range[/C][C]105.65[/C][C]Trim Var.[/C][C]201.334[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=2)[/C][C]51.4767[/C][C]Range[/C][C]37.8[/C][C]Trim Var.[/C][C]32.3332[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=2)[/C][C]94.639[/C][C]Range[/C][C]52.05[/C][C]Trim Var.[/C][C]61.6714[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=2)[/C][C]291.002[/C][C]Range[/C][C]93.35[/C][C]Trim Var.[/C][C]184.164[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=2)[/C][C]986.388[/C][C]Range[/C][C]169[/C][C]Trim Var.[/C][C]636.036[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=265644&T=1

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

As an alternative you can also use a QR Code:  

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

Variance Reduction Matrix
V(Y[t],d=0,D=0)11.5217Range15.6Trim Var.7.6577
V(Y[t],d=1,D=0)14.7492Range19.9Trim Var.9.53755
V(Y[t],d=2,D=0)43.4189Range33.55Trim Var.28.0486
V(Y[t],d=3,D=0)142.75Range62.65Trim Var.93.044
V(Y[t],d=0,D=1)17.353Range25.3Trim Var.10.3776
V(Y[t],d=1,D=1)31.162Range32.85Trim Var.20.36
V(Y[t],d=2,D=1)94.7259Range56.95Trim Var.59.7181
V(Y[t],d=3,D=1)318.797Range105.65Trim Var.201.334
V(Y[t],d=0,D=2)51.4767Range37.8Trim Var.32.3332
V(Y[t],d=1,D=2)94.639Range52.05Trim Var.61.6714
V(Y[t],d=2,D=2)291.002Range93.35Trim Var.184.164
V(Y[t],d=3,D=2)986.388Range169Trim Var.636.036



Parameters (Session):
par1 = 12 ;
Parameters (R input):
par1 = 12 ;
R code (references can be found in the software module):
par1 <- as.numeric(par1)
n <- length(x)
sx <- sort(x)
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Variance Reduction Matrix',6,TRUE)
a<-table.row.end(a)
for (bigd in 0:2) {
for (smalld in 0:3) {
mylabel <- 'V(Y[t],d='
mylabel <- paste(mylabel,as.character(smalld),sep='')
mylabel <- paste(mylabel,',D=',sep='')
mylabel <- paste(mylabel,as.character(bigd),sep='')
mylabel <- paste(mylabel,')',sep='')
a<-table.row.start(a)
a<-table.element(a,mylabel,header=TRUE)
myx <- x
if (smalld > 0) myx <- diff(myx,lag=1,differences=smalld)
if (bigd > 0) myx <- diff(myx,lag=par1,differences=bigd)
a<-table.element(a,signif(var(myx), digits=6))
a<-table.element(a,'Range',header=TRUE)
a<-table.element(a,signif(max(myx)-min(myx), digits=6))
a<-table.element(a,'Trim Var.',header=TRUE)
smyx <- sort(myx)
sn <- length(smyx)
a<-table.element(a,signif(var(smyx[smyx>quantile(smyx,0.05) & smyxa<-table.row.end(a)
}
}
a<-table.end(a)
table.save(a,file='mytable.tab')
bitmap(file='pic0.png')
op <- par(mfrow=c(2,2))
plot(x,type='l',xlab='time',ylab='value',main='d=0, D=0')
plot(diff(x,lag=1,differences=1),type='l',xlab='time',ylab='value',main='d=1, D=0')
plot(diff(x,lag=par1,differences=1),type='l',xlab='time',ylab='value',main='d=0, D=1')
plot(diff(diff(x,lag=1,differences=1),lag=par1,differences=1),type='l',xlab='time',ylab='value',main='d=1, D=1')
par(op)
dev.off()