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 computationSat, 06 Dec 2014 11:51:44 +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/06/t1417866732jdhp6y0et25glmt.htm/, Retrieved Thu, 16 May 2024 22:20:27 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=263592, Retrieved Thu, 16 May 2024 22:20:27 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact107
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Variance Reduction Matrix] [variance reductio...] [2014-12-06 11:51:44] [e4bec374a19c70fe4499af2adad38eb7] [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
6,3
11,3
11,9
9,3
9,6
10
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
9,7
10,8
10,3
10,4
12,7
9,3
11,8
5,9
11,4
13
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
13,3
9,3
12,5
7,6
15,9
9,2
9,1
11,1
13
14,5
12,2
12,3
11,4
8,8
14,6
12,6
13
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
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'Herman Ole Andreas Wold' @ wold.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 & 'Herman Ole Andreas Wold' @ wold.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=263592&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]'Herman Ole Andreas Wold' @ wold.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=263592&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=263592&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'Herman Ole Andreas Wold' @ wold.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)14.7492Range19.9Trim Var.9.53755
V(Y[t],d=1,D=1)43.4189Range33.55Trim Var.28.0486
V(Y[t],d=2,D=1)142.75Range62.65Trim Var.93.044
V(Y[t],d=3,D=1)494.164Range117.05Trim Var.318.455
V(Y[t],d=0,D=2)43.4189Range33.55Trim Var.28.0486
V(Y[t],d=1,D=2)142.75Range62.65Trim Var.93.044
V(Y[t],d=2,D=2)494.164Range117.05Trim Var.318.455
V(Y[t],d=3,D=2)1763.79Range215.8Trim Var.1133.58

\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) & 14.7492 & Range & 19.9 & Trim Var. & 9.53755 \tabularnewline
V(Y[t],d=1,D=1) & 43.4189 & Range & 33.55 & Trim Var. & 28.0486 \tabularnewline
V(Y[t],d=2,D=1) & 142.75 & Range & 62.65 & Trim Var. & 93.044 \tabularnewline
V(Y[t],d=3,D=1) & 494.164 & Range & 117.05 & Trim Var. & 318.455 \tabularnewline
V(Y[t],d=0,D=2) & 43.4189 & Range & 33.55 & Trim Var. & 28.0486 \tabularnewline
V(Y[t],d=1,D=2) & 142.75 & Range & 62.65 & Trim Var. & 93.044 \tabularnewline
V(Y[t],d=2,D=2) & 494.164 & Range & 117.05 & Trim Var. & 318.455 \tabularnewline
V(Y[t],d=3,D=2) & 1763.79 & Range & 215.8 & Trim Var. & 1133.58 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=263592&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]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=1,D=1)[/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=2,D=1)[/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=3,D=1)[/C][C]494.164[/C][C]Range[/C][C]117.05[/C][C]Trim Var.[/C][C]318.455[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=2)[/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=1,D=2)[/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=2,D=2)[/C][C]494.164[/C][C]Range[/C][C]117.05[/C][C]Trim Var.[/C][C]318.455[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=2)[/C][C]1763.79[/C][C]Range[/C][C]215.8[/C][C]Trim Var.[/C][C]1133.58[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=263592&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=263592&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)14.7492Range19.9Trim Var.9.53755
V(Y[t],d=1,D=1)43.4189Range33.55Trim Var.28.0486
V(Y[t],d=2,D=1)142.75Range62.65Trim Var.93.044
V(Y[t],d=3,D=1)494.164Range117.05Trim Var.318.455
V(Y[t],d=0,D=2)43.4189Range33.55Trim Var.28.0486
V(Y[t],d=1,D=2)142.75Range62.65Trim Var.93.044
V(Y[t],d=2,D=2)494.164Range117.05Trim Var.318.455
V(Y[t],d=3,D=2)1763.79Range215.8Trim Var.1133.58



Parameters (Session):
par1 = 1 ;
Parameters (R input):
par1 = 1 ;
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()