Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_variancereduction.wasp
Title produced by softwareVariance Reduction Matrix
Date of computationWed, 17 Dec 2008 12:52:13 -0700
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2008/Dec/17/t122954356948ywunyx4etor5n.htm/, Retrieved Sun, 19 May 2024 05:34:47 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=34529, Retrieved Sun, 19 May 2024 05:34:47 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsVan Dooren Leen
Estimated Impact175
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Variance Reduction Matrix] [paper] [2008-12-17 19:52:13] [d175f84d503eb4f2a43145d5e67795b5] [Current]
Feedback Forum

Post a new message
Dataseries X:
98,6
98,0
106,8
96,6
100,1
107,7
91,5
97,8
107,4
117,5
105,6
97,4
99,5
98,0
104,3
100,6
101,1
103,9
96,9
95,5
108,4
117,0
103,8
100,8
110,6
104,0
112,6
107,3
98,9
109,8
104,9
102,2
123,9
124,9
112,7
121,9
100,6
104,3
120,4
107,5
102,9
125,6
107,5
108,8
128,4
121,1
119,5
128,7
108,7
105,5
119,8
111,3
110,6
120,1
97,5
107,7
127,3
117,2
119,8
116,2
111,0
112,4
130,6
109,1
118,8
123,9
101,6
112,8
128,0
129,6
125,8
119,5
115,7
113,6
129,7
112,0
116,8
127,0
112,1
114,2
121,1
131,6
125,0
120,4
117,7
117,5
120,6
127,5
112,3
124,5
115,2
104,7
130,5




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time0 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135

\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 & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=34529&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]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=34529&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=34529&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'Gwilym Jenkins' @ 72.249.127.135







Variance Reduction Matrix
V(Y[t],d=0,D=0)105.336633941094Range40.1Trim Var.80.8696620628857
V(Y[t],d=1,D=0)122.697022216914Range48.4Trim Var.79.2745377898223
V(Y[t],d=2,D=0)332.875567765568Range77.1Trim Var.229.509604938272
V(Y[t],d=3,D=0)989.184830212235Range139Trim Var.666.891897151899
V(Y[t],d=0,D=1)122.697022216914Range48.4Trim Var.79.2745377898223
V(Y[t],d=1,D=1)332.875567765568Range77.1Trim Var.229.509604938272
V(Y[t],d=2,D=1)989.184830212235Range139Trim Var.666.891897151899
V(Y[t],d=3,D=1)3100.67493615935Range235Trim Var.2061.74388510224
V(Y[t],d=0,D=2)332.875567765568Range77.1Trim Var.229.509604938272
V(Y[t],d=1,D=2)989.184830212235Range139Trim Var.666.891897151899
V(Y[t],d=2,D=2)3100.67493615935Range235Trim Var.2061.74388510224
V(Y[t],d=3,D=2)10055.1667803030Range449.5Trim Var.6506.68660839161

\begin{tabular}{lllllllll}
\hline
Variance Reduction Matrix \tabularnewline
V(Y[t],d=0,D=0) & 105.336633941094 & Range & 40.1 & Trim Var. & 80.8696620628857 \tabularnewline
V(Y[t],d=1,D=0) & 122.697022216914 & Range & 48.4 & Trim Var. & 79.2745377898223 \tabularnewline
V(Y[t],d=2,D=0) & 332.875567765568 & Range & 77.1 & Trim Var. & 229.509604938272 \tabularnewline
V(Y[t],d=3,D=0) & 989.184830212235 & Range & 139 & Trim Var. & 666.891897151899 \tabularnewline
V(Y[t],d=0,D=1) & 122.697022216914 & Range & 48.4 & Trim Var. & 79.2745377898223 \tabularnewline
V(Y[t],d=1,D=1) & 332.875567765568 & Range & 77.1 & Trim Var. & 229.509604938272 \tabularnewline
V(Y[t],d=2,D=1) & 989.184830212235 & Range & 139 & Trim Var. & 666.891897151899 \tabularnewline
V(Y[t],d=3,D=1) & 3100.67493615935 & Range & 235 & Trim Var. & 2061.74388510224 \tabularnewline
V(Y[t],d=0,D=2) & 332.875567765568 & Range & 77.1 & Trim Var. & 229.509604938272 \tabularnewline
V(Y[t],d=1,D=2) & 989.184830212235 & Range & 139 & Trim Var. & 666.891897151899 \tabularnewline
V(Y[t],d=2,D=2) & 3100.67493615935 & Range & 235 & Trim Var. & 2061.74388510224 \tabularnewline
V(Y[t],d=3,D=2) & 10055.1667803030 & Range & 449.5 & Trim Var. & 6506.68660839161 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=34529&T=1

[TABLE]
[ROW][C]Variance Reduction Matrix[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=0)[/C][C]105.336633941094[/C][C]Range[/C][C]40.1[/C][C]Trim Var.[/C][C]80.8696620628857[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=0)[/C][C]122.697022216914[/C][C]Range[/C][C]48.4[/C][C]Trim Var.[/C][C]79.2745377898223[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=0)[/C][C]332.875567765568[/C][C]Range[/C][C]77.1[/C][C]Trim Var.[/C][C]229.509604938272[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=0)[/C][C]989.184830212235[/C][C]Range[/C][C]139[/C][C]Trim Var.[/C][C]666.891897151899[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=1)[/C][C]122.697022216914[/C][C]Range[/C][C]48.4[/C][C]Trim Var.[/C][C]79.2745377898223[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=1)[/C][C]332.875567765568[/C][C]Range[/C][C]77.1[/C][C]Trim Var.[/C][C]229.509604938272[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=1)[/C][C]989.184830212235[/C][C]Range[/C][C]139[/C][C]Trim Var.[/C][C]666.891897151899[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=1)[/C][C]3100.67493615935[/C][C]Range[/C][C]235[/C][C]Trim Var.[/C][C]2061.74388510224[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=2)[/C][C]332.875567765568[/C][C]Range[/C][C]77.1[/C][C]Trim Var.[/C][C]229.509604938272[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=2)[/C][C]989.184830212235[/C][C]Range[/C][C]139[/C][C]Trim Var.[/C][C]666.891897151899[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=2)[/C][C]3100.67493615935[/C][C]Range[/C][C]235[/C][C]Trim Var.[/C][C]2061.74388510224[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=2)[/C][C]10055.1667803030[/C][C]Range[/C][C]449.5[/C][C]Trim Var.[/C][C]6506.68660839161[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=34529&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=34529&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)105.336633941094Range40.1Trim Var.80.8696620628857
V(Y[t],d=1,D=0)122.697022216914Range48.4Trim Var.79.2745377898223
V(Y[t],d=2,D=0)332.875567765568Range77.1Trim Var.229.509604938272
V(Y[t],d=3,D=0)989.184830212235Range139Trim Var.666.891897151899
V(Y[t],d=0,D=1)122.697022216914Range48.4Trim Var.79.2745377898223
V(Y[t],d=1,D=1)332.875567765568Range77.1Trim Var.229.509604938272
V(Y[t],d=2,D=1)989.184830212235Range139Trim Var.666.891897151899
V(Y[t],d=3,D=1)3100.67493615935Range235Trim Var.2061.74388510224
V(Y[t],d=0,D=2)332.875567765568Range77.1Trim Var.229.509604938272
V(Y[t],d=1,D=2)989.184830212235Range139Trim Var.666.891897151899
V(Y[t],d=2,D=2)3100.67493615935Range235Trim Var.2061.74388510224
V(Y[t],d=3,D=2)10055.1667803030Range449.5Trim Var.6506.68660839161



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(x,lag=1,differences=smalld)
if (bigd > 0) myx <- diff(myx,lag=par1,differences=bigd)
a<-table.element(a,var(myx))
a<-table.element(a,'Range',header=TRUE)
a<-table.element(a,max(myx)-min(myx))
a<-table.element(a,'Trim Var.',header=TRUE)
smyx <- sort(myx)
sn <- length(smyx)
a<-table.element(a,var(smyx[smyx>quantile(smyx,0.05) & smyxa<-table.row.end(a)
}
}
a<-table.end(a)
table.save(a,file='mytable.tab')