Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_variancereduction.wasp
Title produced by softwareVariance Reduction Matrix
Date of computationTue, 29 Dec 2009 08:19:26 -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/29/t126210002668jqgdpgmbh4ggn.htm/, Retrieved Fri, 03 May 2024 09:18:24 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=71139, Retrieved Fri, 03 May 2024 09:18:24 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact131
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Variance Reduction Matrix] [vrm: olieprijs] [2009-12-29 15:19:26] [dbd46bd47d5f87b1007a5a1708bef00e] [Current]
Feedback Forum

Post a new message
Dataseries X:
32,68
31,54
32,43
26,54
25,85
27,6
25,71
25,38
28,57
27,64
25,36
25,9
26,29
21,74
19,2
19,32
19,82
20,36
24,31
25,97
25,61
24,67
25,59
26,09
28,37
27,34
24,46
27,46
30,23
32,33
29,87
24,87
25,48
27,28
28,24
29,58
26,95
29,08
28,76
29,59
30,7
30,52
32,67
33,19
37,13
35,54
37,75
41,84
42,94
49,14
44,61
40,22
44,23
45,85
53,38
53,26
51,8
55,3
57,81
63,96
63,77
59,15
56,12
57,42
63,52
61,71
63,01
68,18
72,03
69,75
74,41
74,33
64,24
60,03
59,44
62,5
55,04
58,34
61,92
67,65
67,68
70,3
75,26
71,44
76,36
81,71
92,6
90,6
92,23
94,09
102,79
109,65
124,05
132,69
135,81
116,07
101,42
75,73
55,48
43,80
45,29
44,01
47,48
51,07
57,84
69,04
65,61
72,87
68,41
73,25
77,43




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 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 & 1 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=71139&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]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=71139&T=0

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







Variance Reduction Matrix
V(Y[t],d=0,D=0)693.720836871417Range116.61Trim Var.417.237441517213
V(Y[t],d=1,D=0)33.0055796914095Range40.09Trim Var.10.0025212392173
V(Y[t],d=2,D=0)33.5214476384642Range36.03Trim Var.17.6619791881443
V(Y[t],d=3,D=0)90.8775747577017Range50.3600000000001Trim Var.50.5423649890351
V(Y[t],d=0,D=1)616.263325190682Range132.59Trim Var.289.160598237998
V(Y[t],d=1,D=1)75.0933890174627Range61.67Trim Var.21.1495315438872
V(Y[t],d=2,D=1)71.5439128436427Range50.3600000000001Trim Var.35.5145239775461
V(Y[t],d=3,D=1)201.900651929825Range93.02Trim Var.99.5338814363886
V(Y[t],d=0,D=2)1587.00939991981Range192.59Trim Var.714.737446753247
V(Y[t],d=1,D=2)232.184194965800Range104.19Trim Var.55.3030887017544
V(Y[t],d=2,D=2)184.571329467787Range76.83Trim Var.93.6433496216217
V(Y[t],d=3,D=2)505.665034423408Range109.67Trim Var.281.755548556090

\begin{tabular}{lllllllll}
\hline
Variance Reduction Matrix \tabularnewline
V(Y[t],d=0,D=0) & 693.720836871417 & Range & 116.61 & Trim Var. & 417.237441517213 \tabularnewline
V(Y[t],d=1,D=0) & 33.0055796914095 & Range & 40.09 & Trim Var. & 10.0025212392173 \tabularnewline
V(Y[t],d=2,D=0) & 33.5214476384642 & Range & 36.03 & Trim Var. & 17.6619791881443 \tabularnewline
V(Y[t],d=3,D=0) & 90.8775747577017 & Range & 50.3600000000001 & Trim Var. & 50.5423649890351 \tabularnewline
V(Y[t],d=0,D=1) & 616.263325190682 & Range & 132.59 & Trim Var. & 289.160598237998 \tabularnewline
V(Y[t],d=1,D=1) & 75.0933890174627 & Range & 61.67 & Trim Var. & 21.1495315438872 \tabularnewline
V(Y[t],d=2,D=1) & 71.5439128436427 & Range & 50.3600000000001 & Trim Var. & 35.5145239775461 \tabularnewline
V(Y[t],d=3,D=1) & 201.900651929825 & Range & 93.02 & Trim Var. & 99.5338814363886 \tabularnewline
V(Y[t],d=0,D=2) & 1587.00939991981 & Range & 192.59 & Trim Var. & 714.737446753247 \tabularnewline
V(Y[t],d=1,D=2) & 232.184194965800 & Range & 104.19 & Trim Var. & 55.3030887017544 \tabularnewline
V(Y[t],d=2,D=2) & 184.571329467787 & Range & 76.83 & Trim Var. & 93.6433496216217 \tabularnewline
V(Y[t],d=3,D=2) & 505.665034423408 & Range & 109.67 & Trim Var. & 281.755548556090 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=71139&T=1

[TABLE]
[ROW][C]Variance Reduction Matrix[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=0)[/C][C]693.720836871417[/C][C]Range[/C][C]116.61[/C][C]Trim Var.[/C][C]417.237441517213[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=0)[/C][C]33.0055796914095[/C][C]Range[/C][C]40.09[/C][C]Trim Var.[/C][C]10.0025212392173[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=0)[/C][C]33.5214476384642[/C][C]Range[/C][C]36.03[/C][C]Trim Var.[/C][C]17.6619791881443[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=0)[/C][C]90.8775747577017[/C][C]Range[/C][C]50.3600000000001[/C][C]Trim Var.[/C][C]50.5423649890351[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=1)[/C][C]616.263325190682[/C][C]Range[/C][C]132.59[/C][C]Trim Var.[/C][C]289.160598237998[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=1)[/C][C]75.0933890174627[/C][C]Range[/C][C]61.67[/C][C]Trim Var.[/C][C]21.1495315438872[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=1)[/C][C]71.5439128436427[/C][C]Range[/C][C]50.3600000000001[/C][C]Trim Var.[/C][C]35.5145239775461[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=1)[/C][C]201.900651929825[/C][C]Range[/C][C]93.02[/C][C]Trim Var.[/C][C]99.5338814363886[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=2)[/C][C]1587.00939991981[/C][C]Range[/C][C]192.59[/C][C]Trim Var.[/C][C]714.737446753247[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=2)[/C][C]232.184194965800[/C][C]Range[/C][C]104.19[/C][C]Trim Var.[/C][C]55.3030887017544[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=2)[/C][C]184.571329467787[/C][C]Range[/C][C]76.83[/C][C]Trim Var.[/C][C]93.6433496216217[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=2)[/C][C]505.665034423408[/C][C]Range[/C][C]109.67[/C][C]Trim Var.[/C][C]281.755548556090[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=71139&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=71139&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)693.720836871417Range116.61Trim Var.417.237441517213
V(Y[t],d=1,D=0)33.0055796914095Range40.09Trim Var.10.0025212392173
V(Y[t],d=2,D=0)33.5214476384642Range36.03Trim Var.17.6619791881443
V(Y[t],d=3,D=0)90.8775747577017Range50.3600000000001Trim Var.50.5423649890351
V(Y[t],d=0,D=1)616.263325190682Range132.59Trim Var.289.160598237998
V(Y[t],d=1,D=1)75.0933890174627Range61.67Trim Var.21.1495315438872
V(Y[t],d=2,D=1)71.5439128436427Range50.3600000000001Trim Var.35.5145239775461
V(Y[t],d=3,D=1)201.900651929825Range93.02Trim Var.99.5338814363886
V(Y[t],d=0,D=2)1587.00939991981Range192.59Trim Var.714.737446753247
V(Y[t],d=1,D=2)232.184194965800Range104.19Trim Var.55.3030887017544
V(Y[t],d=2,D=2)184.571329467787Range76.83Trim Var.93.6433496216217
V(Y[t],d=3,D=2)505.665034423408Range109.67Trim Var.281.755548556090



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(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')