Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_variancereduction.wasp
Title produced by softwareVariance Reduction Matrix
Date of computationThu, 04 Dec 2008 08:53:33 -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/04/t12284060379rzy4pk0q2rrsm4.htm/, Retrieved Fri, 17 May 2024 07:00:40 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=28940, Retrieved Fri, 17 May 2024 07:00:40 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact231
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
F     [Univariate Data Series] [Airline data] [2007-10-18 09:58:47] [42daae401fd3def69a25014f2252b4c2]
F RMPD  [Cross Correlation Function] [] [2008-12-02 21:28:41] [74be16979710d4c4e7c6647856088456]
F   PD    [Cross Correlation Function] [] [2008-12-02 21:39:53] [74be16979710d4c4e7c6647856088456]
- RMPD        [Variance Reduction Matrix] [] [2008-12-04 15:53:33] [e02910eed3830f1815f587e12f46cbdb] [Current]
Feedback Forum

Post a new message
Dataseries X:
100.8
100.7
86.2
83.2
71.7
77.5
89.8
80.3
78.7
93.8
57.6
60.6
91
85.3
77.4
77.3
68.3
69.9
81.7
75.1
69.9
84
54.3
60
89.9
77
85.3
77.6
69.2
75.5
85.7
72.2
79.9
85.3
52.2
61.2
82.4
85.4
78.2
70.2
70.2
69.3
77.5
66.1
69
79.2
56.2
63.3
77.8
92
78.1
65.1
71.1
70.9
72
81.9
70.6
72.5
65.1
61.1




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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=28940&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)117.287827683616Range48.6Trim Var.74.4323864430468
V(Y[t],d=1,D=0)172.441665692577Range66.6Trim Var.85.8323294629898
V(Y[t],d=2,D=0)440.207153660012Range93.4Trim Var.277.296753393665
V(Y[t],d=3,D=0)1286.0096679198Range158.5Trim Var.757.898666666667
V(Y[t],d=0,D=1)40.9986125886525Range31.2Trim Var.24.1392195121951
V(Y[t],d=1,D=1)68.6276318223867Range36.8Trim Var.40.5862560975610
V(Y[t],d=2,D=1)211.282942028986Range63.9Trim Var.117.005743589744
V(Y[t],d=3,D=1)668.960525252525Range110.5Trim Var.380.578515519568
V(Y[t],d=0,D=2)92.794753968254Range42.8Trim Var.56.3976512096774
V(Y[t],d=1,D=2)130.026218487395Range54.8Trim Var.73.4722365591397
V(Y[t],d=2,D=2)401.16367201426Range93.8Trim Var.238.531091954023
V(Y[t],d=3,D=2)1342.35257575758Range179.9Trim Var.734.348645320197

\begin{tabular}{lllllllll}
\hline
Variance Reduction Matrix \tabularnewline
V(Y[t],d=0,D=0) & 117.287827683616 & Range & 48.6 & Trim Var. & 74.4323864430468 \tabularnewline
V(Y[t],d=1,D=0) & 172.441665692577 & Range & 66.6 & Trim Var. & 85.8323294629898 \tabularnewline
V(Y[t],d=2,D=0) & 440.207153660012 & Range & 93.4 & Trim Var. & 277.296753393665 \tabularnewline
V(Y[t],d=3,D=0) & 1286.0096679198 & Range & 158.5 & Trim Var. & 757.898666666667 \tabularnewline
V(Y[t],d=0,D=1) & 40.9986125886525 & Range & 31.2 & Trim Var. & 24.1392195121951 \tabularnewline
V(Y[t],d=1,D=1) & 68.6276318223867 & Range & 36.8 & Trim Var. & 40.5862560975610 \tabularnewline
V(Y[t],d=2,D=1) & 211.282942028986 & Range & 63.9 & Trim Var. & 117.005743589744 \tabularnewline
V(Y[t],d=3,D=1) & 668.960525252525 & Range & 110.5 & Trim Var. & 380.578515519568 \tabularnewline
V(Y[t],d=0,D=2) & 92.794753968254 & Range & 42.8 & Trim Var. & 56.3976512096774 \tabularnewline
V(Y[t],d=1,D=2) & 130.026218487395 & Range & 54.8 & Trim Var. & 73.4722365591397 \tabularnewline
V(Y[t],d=2,D=2) & 401.16367201426 & Range & 93.8 & Trim Var. & 238.531091954023 \tabularnewline
V(Y[t],d=3,D=2) & 1342.35257575758 & Range & 179.9 & Trim Var. & 734.348645320197 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=28940&T=1

[TABLE]
[ROW][C]Variance Reduction Matrix[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=0)[/C][C]117.287827683616[/C][C]Range[/C][C]48.6[/C][C]Trim Var.[/C][C]74.4323864430468[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=0)[/C][C]172.441665692577[/C][C]Range[/C][C]66.6[/C][C]Trim Var.[/C][C]85.8323294629898[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=0)[/C][C]440.207153660012[/C][C]Range[/C][C]93.4[/C][C]Trim Var.[/C][C]277.296753393665[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=0)[/C][C]1286.0096679198[/C][C]Range[/C][C]158.5[/C][C]Trim Var.[/C][C]757.898666666667[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=1)[/C][C]40.9986125886525[/C][C]Range[/C][C]31.2[/C][C]Trim Var.[/C][C]24.1392195121951[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=1)[/C][C]68.6276318223867[/C][C]Range[/C][C]36.8[/C][C]Trim Var.[/C][C]40.5862560975610[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=1)[/C][C]211.282942028986[/C][C]Range[/C][C]63.9[/C][C]Trim Var.[/C][C]117.005743589744[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=1)[/C][C]668.960525252525[/C][C]Range[/C][C]110.5[/C][C]Trim Var.[/C][C]380.578515519568[/C][/ROW]
[ROW][C]V(Y[t],d=0,D=2)[/C][C]92.794753968254[/C][C]Range[/C][C]42.8[/C][C]Trim Var.[/C][C]56.3976512096774[/C][/ROW]
[ROW][C]V(Y[t],d=1,D=2)[/C][C]130.026218487395[/C][C]Range[/C][C]54.8[/C][C]Trim Var.[/C][C]73.4722365591397[/C][/ROW]
[ROW][C]V(Y[t],d=2,D=2)[/C][C]401.16367201426[/C][C]Range[/C][C]93.8[/C][C]Trim Var.[/C][C]238.531091954023[/C][/ROW]
[ROW][C]V(Y[t],d=3,D=2)[/C][C]1342.35257575758[/C][C]Range[/C][C]179.9[/C][C]Trim Var.[/C][C]734.348645320197[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=28940&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=28940&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)117.287827683616Range48.6Trim Var.74.4323864430468
V(Y[t],d=1,D=0)172.441665692577Range66.6Trim Var.85.8323294629898
V(Y[t],d=2,D=0)440.207153660012Range93.4Trim Var.277.296753393665
V(Y[t],d=3,D=0)1286.0096679198Range158.5Trim Var.757.898666666667
V(Y[t],d=0,D=1)40.9986125886525Range31.2Trim Var.24.1392195121951
V(Y[t],d=1,D=1)68.6276318223867Range36.8Trim Var.40.5862560975610
V(Y[t],d=2,D=1)211.282942028986Range63.9Trim Var.117.005743589744
V(Y[t],d=3,D=1)668.960525252525Range110.5Trim Var.380.578515519568
V(Y[t],d=0,D=2)92.794753968254Range42.8Trim Var.56.3976512096774
V(Y[t],d=1,D=2)130.026218487395Range54.8Trim Var.73.4722365591397
V(Y[t],d=2,D=2)401.16367201426Range93.8Trim Var.238.531091954023
V(Y[t],d=3,D=2)1342.35257575758Range179.9Trim Var.734.348645320197



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