Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_meanplot.wasp
Title produced by softwareMean Plot
Date of computationTue, 18 Dec 2012 14:33:01 -0500
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2012/Dec/18/t1355859236jn7oa7fiwm995ya.htm/, Retrieved Fri, 29 Mar 2024 10:22:30 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=201604, Retrieved Fri, 29 Mar 2024 10:22:30 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact145
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Univariate Explorative Data Analysis] [time effect in su...] [2010-11-17 08:55:33] [b98453cac15ba1066b407e146608df68]
- RMPD    [Mean Plot] [Paper Statistiek ...] [2012-12-18 19:33:01] [3353489d44052879174bf0d9e8b7362f] [Current]
Feedback Forum

Post a new message
Dataseries X:
210907
120982
176508
179321
123185
52746
385534
33170
101645
149061
165446
237213
173326
133131
258873
180083
324799
230964
236785
135473
202925
215147
344297
153935
132943
174724
174415
225548
223632
124817
221698
210767
170266
260561
84853
294424
101011
215641
325107
7176
167542
106408
96560
265769
269651
149112
175824
152871
111665
116408
362301
78800
183167
277965
150629
168809
24188
329267
65029
101097
218946
244052
341570
103597
233328
256462
206161
311473
235800
177939
207176
196553
174184
143246
187559
187681
119016
182192
73566
194979
167488
143756
275541
243199
182999
135649
152299
120221
346485
145790
193339
80953
122774
130585
112611
286468
241066
148446
204713
182079
140344
220516
243060
162765
182613
232138
265318
85574
310839
225060
232317
144966
43287
155754
164709
201940
235454
220801
99466
92661
133328
61361
125930
100750
224549
82316
102010
101523
243511
22938
41566
152474
61857
99923
132487
317394
21054
209641
22648
31414
46698
131698
91735
244749
184510
79863
128423
97839
38214
151101
272458
172494
108043
328107
250579
351067
158015
98866
85439
229242
351619
84207
120445
324598
131069
204271
165543
141722
116048
250047
299775
195838
173260
254488
104389
136084
199476
92499
224330
135781
74408
81240
14688
181633
271856
7199
46660
17547
133368
95227
152601
98146
79619
59194
139942
118612
72880
65475
99643
71965
77272
49289
135131
108446
89746
44296
77648
181528
134019
124064
92630
121848
52915
81872
58981
53515
60812
56375
65490
80949
76302
104011
98104
67989
30989
135458
73504
63123
61254
74914
31774
81437
87186
50090
65745
56653
158399
46455
73624
38395
91899
139526
52164
51567
70551
84856
102538
86678
85709
34662
150580
99611
19349
99373
86230
30837
31706
89806
62088
40151
27634
76990
37460
54157
49862
84337
64175
59382
119308
76702
103425
70344
43410
104838
62215
69304
53117
19764
86680
84105
77945
89113
91005
40248
64187
50857
56613
62792
72535




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time6 seconds
R Server'Sir Maurice George Kendall' @ kendall.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 & 6 seconds \tabularnewline
R Server & 'Sir Maurice George Kendall' @ kendall.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=201604&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]6 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Maurice George Kendall' @ kendall.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=201604&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=201604&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 time6 seconds
R Server'Sir Maurice George Kendall' @ kendall.wessa.net



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))
(np <- floor(n / par1))
arr <- array(NA,dim=c(par1,np+1))
darr <- array(NA,dim=c(par1,np+1))
ari <- array(0,dim=par1)
dx <- diff(x)
j <- 0
for (i in 1:n)
{
j = j + 1
ari[j] = ari[j] + 1
arr[j,ari[j]] <- x[i]
darr[j,ari[j]] <- dx[i]
if (j == par1) j = 0
}
ari
arr
darr
arr.mean <- array(NA,dim=par1)
arr.median <- array(NA,dim=par1)
arr.midrange <- array(NA,dim=par1)
for (j in 1:par1)
{
arr.mean[j] <- mean(arr[j,],na.rm=TRUE)
arr.median[j] <- median(arr[j,],na.rm=TRUE)
arr.midrange[j] <- (quantile(arr[j,],0.75,na.rm=TRUE) + quantile(arr[j,],0.25,na.rm=TRUE)) / 2
}
overall.mean <- mean(x)
overall.median <- median(x)
overall.midrange <- (quantile(x,0.75) + quantile(x,0.25)) / 2
bitmap(file='plot1.png')
plot(arr.mean,type='b',ylab='mean',main='Mean Plot',xlab='Periodic Index')
mtext(paste('#blocks = ',np))
abline(overall.mean,0)
dev.off()
bitmap(file='plot2.png')
plot(arr.median,type='b',ylab='median',main='Median Plot',xlab='Periodic Index')
mtext(paste('#blocks = ',np))
abline(overall.median,0)
dev.off()
bitmap(file='plot3.png')
plot(arr.midrange,type='b',ylab='midrange',main='Midrange Plot',xlab='Periodic Index')
mtext(paste('#blocks = ',np))
abline(overall.midrange,0)
dev.off()
bitmap(file='plot4.png')
z <- data.frame(t(arr))
names(z) <- c(1:par1)
(boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Periodic Subseries'))
dev.off()
bitmap(file='plot4b.png')
z <- data.frame(t(darr))
names(z) <- c(1:par1)
(boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Differenced Periodic Subseries'))
dev.off()
bitmap(file='plot5.png')
z <- data.frame(arr)
names(z) <- c(1:np)
(boxplot(z,notch=TRUE,col='grey',xlab='Block Index',ylab='Value',main='Notched Box Plots - Sequential Blocks'))
dev.off()
bitmap(file='plot6.png')
z <- data.frame(cbind(arr.mean,arr.median,arr.midrange))
names(z) <- list('mean','median','midrange')
(boxplot(z,notch=TRUE,col='grey',ylab='Overall Central Tendency',main='Notched Box Plots'))
dev.off()