FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_meanplot.wasp
Title produced by softwareMean Plot
Date of computationSat, 13 Dec 2014 15:52:04 +0000
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2014/Dec/13/t1418485938f3rfy33xyrbsrk8.htm/, Retrieved Thu, 16 May 2024 21:16:43 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=267169, Retrieved Thu, 16 May 2024 21:16:43 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact75

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Mean Plot] [] [2014-12-13 15:52:04] [f403c9f98aaaf69e3ef2a3935929401f] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
21
26
22
22
18
23
12
20
22
21
19
22
15
20
19
18
15
20
21
21
15
16
23
21
18
25
9
30
20
23
16
16
19
25
25
18
23
21
10
14
22
26
23
23
24
24
18
23
15
19
16
25
23
17
19
21
18
27
21
13
8
29
28
23
21
19
19
20
18
19
17
19
25
19
22
23
26
14
28
16
24
20
12
24
22
12
22
20
10
23
17
22
24
18
21
20
20
22
19
20
26
23
24
21
21
19
8
17
20
11
8
15
18
18
19
19
23
22
21
25
30
17
27
23
23
18
18
23
19
15
20
16
24
25
25
19
19
16
19
19
23
21
22
19
20
20
3
23
14
23
20
15
13
16
7
24
17
24
24
19
25
20
28
23
27
18
28
21
19
23
27
22
28
25
21
22
28
20
29
25
25
20
20
16
20
20
23
18
25
18
19
25
25
25
24
19
26
10
17
13
17
30
25
4
16
21
23
22
17
20
20
22
16
23
16
0
18
25
23
12
18
24
11
18
14
23
24
29
18
15
29
16
19
22
16
23
23
19
4
20
24
20
4
24
22
16
3
15
24
17
20
27
23
26
23
17
20
22
19
24
19
23
15
27
26
22
22
18
15
22
27
10
20
17
23
19
13
27
23
16
25
2
26
20
23
22
24
22
17
23
23
28
29
21
24
20
7
19
28
18
26
21
19
20
NA
23
24
16
19
24
21
16
16
21
NA
28
16
23
26
29
18
19
19
16
16
16
18
22
14
20
15
22
24
16
19
24
19
15
11
15
17
20
21
16
17
20
15
21
16
18
25
21
21
16
20
24
28
27
22
20
27
17
22
23
15
22
13
21
18
22
19
15
20
17
21
23
20
18
22
24
24
18
27
19
20
15
20
27
20
20
13
21
23
26
24
25
18
21
23
16
19
20
25
22
20
25
27
20
18
26
26
24
27
16
15
25
27
18
16
18
23
21
21
14
24
18
16
25
22
13
20
17
23
22
23
22
23
10
18
25
26
14
23
22
23
19
14
26
24
21
17
16
15
11
NA
19
21
20
16
19
16
11
22
20
26
26
20
24
20
15
23
25
27
23
20
25
24
22
27
20
17
22
26
19
19
24
22
16
22
23
19
20
16
19
20
15
22
22
12
15
21
26
27
23
21
22
26
24
27
18
18
25
12
19
24
17
22
15
20
24
17
24
15
20
20
17
11
21
28
14
13
12
21
13
19
23
27
25
22
27
16
20
18
19
17
10
11
16
13
14
12
15
19
15
14
14
10
13
21
11
14
20
7
22
24
16
22
25
5
19
23
13
10
12
21
22
20
17
20
13
9
22
15
12
25
14
14
17
9
10
15
15
15
14
21
13
18
20
16
28
12
20
26
18
21
23
13
22
14
23
16
14
22
19
23
16
20
8
16
11
16
10
17
16
17
10
15
13
19
14
18
25
10
22
15
18
22
18
15
20
18
6
17
12
12
19
23
26
28
19
16
3
11
15
22
12
21
25
12
14
24
12
13
15
17
12
28
25
14
21
18
23
16
15
5
19
22
DELETED
19
12
22
18
24
19
4
20
24
26
22
19
9
22
18
16
19
20
21
17
9
26
28
13
16
22
18
21
10
15
15
13
10
23
21
14
17
15
15
17
26
12
14
26
18
17
20
16
19
12
20
19
25
19
15
12

Tables (Output of Computation)





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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=267169&T=0

As an alternative you can also use a QR Code:  
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'Gertrude Mary Cox' @ cox.wessa.net


Figures (Output of Computation)

Input Parameters & R Code

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