Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_variancereduction.wasp
Title produced by softwareVariance Reduction Matrix
Date of computationFri, 18 Dec 2009 06:02:36 -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/2009/Dec/18/t1261141393rxqyn82zq3xffht.htm/, Retrieved Sat, 27 Apr 2024 07:15:07 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=69300, Retrieved Sat, 27 Apr 2024 07:15:07 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact118
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Variance Reduction Matrix] [] [2009-12-18 13:02:36] [02cb93c9d037d32bf77dfc632a3a9fbe] [Current]
Feedback Forum

Post a new message
Dataseries X:
100.00
110.08
123.28
106.30
114.60
110.22
117.82
138.84
131.13
153.26
148.76
110.75
132.00
138.55
126.50
119.47
152.49
152.46
134.08
162.33
130.46
142.83
148.39
122.88
125.86
133.31
140.97
110.30
123.04
125.99
112.24
136.10
111.86
109.63
135.75
114.39
121.79
101.33
147.01
113.01
101.89
117.85
128.68
117.27
121.55
109.18
127.88
104.92
100.93
116.79
107.50
109.18
128.31
83.80
93.16
103.99
106.33
97.55
111.71
105.43
103.62




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=69300&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=69300&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=69300&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)284.75978704918Range78.53Trim Var.176.891435486212
V(Y[t],d=1,D=0)343.730037175141Range90.19Trim Var.221.134402271139
V(Y[t],d=2,D=0)1044.20540280538Range145.82Trim Var.687.728341291727
V(Y[t],d=3,D=0)3455.97152065941Range263.33Trim Var.2387.50122986425
V(Y[t],d=0,D=1)394.449066326531Range81.75Trim Var.251.620169324474
V(Y[t],d=1,D=1)477.167742154255Range98.49Trim Var.255.650681997677
V(Y[t],d=2,D=1)1370.80627816836Range181.94Trim Var.656.50952804878
V(Y[t],d=3,D=1)4331.92222164251Range320.94Trim Var.2116.91743942308
V(Y[t],d=0,D=2)876.50835Range116.28Trim Var.568.597440530304
V(Y[t],d=1,D=2)1163.05140658730Range157.22Trim Var.612.234044758065
V(Y[t],d=2,D=2)3373.50416336135Range289.22Trim Var.1636.65642494624
V(Y[t],d=3,D=2)10795.5768178253Range526.48Trim Var.5378.66290678161

\begin{tabular}{lllllllll}
\hline
Variance Reduction Matrix \tabularnewline
V(Y[t],d=0,D=0) & 284.75978704918 & Range & 78.53 & Trim Var. & 176.891435486212 \tabularnewline
V(Y[t],d=1,D=0) & 343.730037175141 & Range & 90.19 & Trim Var. & 221.134402271139 \tabularnewline
V(Y[t],d=2,D=0) & 1044.20540280538 & Range & 145.82 & Trim Var. & 687.728341291727 \tabularnewline
V(Y[t],d=3,D=0) & 3455.97152065941 & Range & 263.33 & Trim Var. & 2387.50122986425 \tabularnewline
V(Y[t],d=0,D=1) & 394.449066326531 & Range & 81.75 & Trim Var. & 251.620169324474 \tabularnewline
V(Y[t],d=1,D=1) & 477.167742154255 & Range & 98.49 & Trim Var. & 255.650681997677 \tabularnewline
V(Y[t],d=2,D=1) & 1370.80627816836 & Range & 181.94 & Trim Var. & 656.50952804878 \tabularnewline
V(Y[t],d=3,D=1) & 4331.92222164251 & Range & 320.94 & Trim Var. & 2116.91743942308 \tabularnewline
V(Y[t],d=0,D=2) & 876.50835 & Range & 116.28 & Trim Var. & 568.597440530304 \tabularnewline
V(Y[t],d=1,D=2) & 1163.05140658730 & Range & 157.22 & Trim Var. & 612.234044758065 \tabularnewline
V(Y[t],d=2,D=2) & 3373.50416336135 & Range & 289.22 & Trim Var. & 1636.65642494624 \tabularnewline
V(Y[t],d=3,D=2) & 10795.5768178253 & Range & 526.48 & Trim Var. & 5378.66290678161 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=69300&T=1

[TABLE]
[ROW][C]Variance Reduction Matrix[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=0)[/C][C]284.75978704918[/C][C]Range[/C][C]78.53[/C][C]Trim Var.[/C][C]176.891435486212[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=0)[/C][C]343.730037175141[/C][C]Range[/C][C]90.19[/C][C]Trim Var.[/C][C]221.134402271139[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=0)[/C][C]1044.20540280538[/C][C]Range[/C][C]145.82[/C][C]Trim Var.[/C][C]687.728341291727[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=0)[/C][C]3455.97152065941[/C][C]Range[/C][C]263.33[/C][C]Trim Var.[/C][C]2387.50122986425[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=1)[/C][C]394.449066326531[/C][C]Range[/C][C]81.75[/C][C]Trim Var.[/C][C]251.620169324474[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=1)[/C][C]477.167742154255[/C][C]Range[/C][C]98.49[/C][C]Trim Var.[/C][C]255.650681997677[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=1)[/C][C]1370.80627816836[/C][C]Range[/C][C]181.94[/C][C]Trim Var.[/C][C]656.50952804878[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=1)[/C][C]4331.92222164251[/C][C]Range[/C][C]320.94[/C][C]Trim Var.[/C][C]2116.91743942308[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=2)[/C][C]876.50835[/C][C]Range[/C][C]116.28[/C][C]Trim Var.[/C][C]568.597440530304[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=2)[/C][C]1163.05140658730[/C][C]Range[/C][C]157.22[/C][C]Trim Var.[/C][C]612.234044758065[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=2)[/C][C]3373.50416336135[/C][C]Range[/C][C]289.22[/C][C]Trim Var.[/C][C]1636.65642494624[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=2)[/C][C]10795.5768178253[/C][C]Range[/C][C]526.48[/C][C]Trim Var.[/C][C]5378.66290678161[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=69300&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=69300&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)284.75978704918Range78.53Trim Var.176.891435486212
V(Y[t],d=1,D=0)343.730037175141Range90.19Trim Var.221.134402271139
V(Y[t],d=2,D=0)1044.20540280538Range145.82Trim Var.687.728341291727
V(Y[t],d=3,D=0)3455.97152065941Range263.33Trim Var.2387.50122986425
V(Y[t],d=0,D=1)394.449066326531Range81.75Trim Var.251.620169324474
V(Y[t],d=1,D=1)477.167742154255Range98.49Trim Var.255.650681997677
V(Y[t],d=2,D=1)1370.80627816836Range181.94Trim Var.656.50952804878
V(Y[t],d=3,D=1)4331.92222164251Range320.94Trim Var.2116.91743942308
V(Y[t],d=0,D=2)876.50835Range116.28Trim Var.568.597440530304
V(Y[t],d=1,D=2)1163.05140658730Range157.22Trim Var.612.234044758065
V(Y[t],d=2,D=2)3373.50416336135Range289.22Trim Var.1636.65642494624
V(Y[t],d=3,D=2)10795.5768178253Range526.48Trim Var.5378.66290678161



Parameters (Session):
par1 = 60 ; par2 = 1 ; par3 = 1 ; par4 = 1 ; par5 = 12 ; par6 = White Noise ; par7 = 0.95 ;
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(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')