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, 15 Nov 2011 12:27:26 -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/2011/Nov/15/t1321378055pytaat62kjja1fe.htm/, Retrieved Fri, 29 Mar 2024 15:17:32 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=143244, Retrieved Fri, 29 Mar 2024 15:17:32 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact85
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Linear Regression Graphical Model Validation] [Colombia Coffee -...] [2008-02-26 10:22:06] [74be16979710d4c4e7c6647856088456]
-  M D  [Linear Regression Graphical Model Validation] [Regression Model 1] [2010-11-16 09:46:35] [1429a1a14191a86916b95357f6de790b]
-    D    [Linear Regression Graphical Model Validation] [Regression Model 2] [2010-11-16 17:18:05] [1429a1a14191a86916b95357f6de790b]
- RMPD        [Mean Plot] [] [2011-11-15 17:27:26] [7dc03dd48c8acabd98b217fada4a6bc0] [Current]
Feedback Forum

Post a new message
Dataseries X:
61.04
64.21
63.92
68.48
67.76
65.37
61.84
59.88
62.91
63.67
59.96
62.80
64.26
66.35
67.43
69.45
72.98
71.58
69.97
71.93
68.53
71.29
72.75
71.64
69.97
70.78
73.52
73.86
76.80
74.57
73.36
74.77
74.35
71.14
72.51
69.19
66.25
63.33
60.55
62.91
59.76
58.57
59.33
60.75
59.05
59.59
58.97
59.24
63.28
62.03
63.43
62.41
61.05
56.31
53.30
51.99
55.42
59.02
59.89
59.39
61.44
61.64
59.69
57.11
62.28
65.87
64.28
63.63
64.11
66.46
61.93
62.37
64.94
65.20
65.08
64.76
68.00
69.14
70.68
72.81
73.93
75.79
77.02
75.48
71.47
71.98
71.09
74.04
76.70
79.10
81.62
81.66
81.22
83.69
88.60
91.86
95.93
96.32
93.84
98.18
88.71
88.28
91.27
93.31
96.00
95.98
97.91
92.69
90.57
90.71
88.96
91.77
95.50
98.81
01.84
05.15
10.21
01.84
05.62
06.23
10.14
16.69
18.52
16.32
25.96
26.29
32.19
27.35
38.54
34.86
35.36
40.21
45.29
45.08
28.88
23.26
25.10
15.20
13.77
14.59
15.46
06.23
01.18
04.55
06.89
93.88
77.70
71.85
64.15
67.81
61.04
57.04
49.93
54.43
40.81
46.28
42.36
37.71
44.60
46.34
40.83
36.51
46.47
41.68
40.17
37.51
40.03
44.76
45.52
46.25
52.07
52.38
52.51
52.24
50.33
51.55
53.20
58.63
56.34
61.67
66.31
68.44
72.04
69.55
69.16
65.63
59.89
63.56
68.05
69.45
70.93
67.51
73.89
72.74
68.02
69.29
72.04
66.02
69.95
71.77
78.53
80.50
77.00
77.43
76.35
77.47
76.05
75.47
69.87
73.36
78.05
79.36
82.75
78.00
74.54
72.89
71.19
74.13
79.81
79.66
81.50
81.24
80.68
80.00
84.87
84.92
83.24
85.12
86.15
75.11
71.61
70.04
73.97
71.51
73.78
77.18
78.86
72.14
76.09
76.01
78.98
78.95
80.70
75.39
73.82
75.17
74.60
76.45
73.66
76.49
81.58
82.66
81.25
81.69
81.43
86.85
84.88
81.98
83.76
89.19
87.79
88.02
91.51
91.38
88.03
91.54
89.11
89.34
89.03
85.58
86.20
97.88
04.42
01.16
01.07
05.40
07.94
12.79
09.66
12.29
13.93
97.18
99.65
00.10
00.59
00.22
99.29
93.01
91.16
94.94
96.20
97.24
99.87
95.70
86.88
85.38
82.26
85.37
86.45
87.24
87.96
79.85
79.20
82.98
86.80
87.40
93.32
94.26
98.99




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

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