Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_sdplot.wasp
Title produced by softwareStandard Deviation Plot
Date of computationThu, 28 May 2009 08:32:24 -0600
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/May/28/t1243521171bst9u56jc7v8fqv.htm/, Retrieved Mon, 06 May 2024 03:19:17 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=40630, Retrieved Mon, 06 May 2024 03:19:17 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact122
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Standard Deviation Plot] [standard deviatio...] [2009-05-28 14:32:24] [1a4b490ae86a78f004f9a6b70ce3539b] [Current]
-    D    [Standard Deviation Plot] [spreidingsgrafiek...] [2009-06-02 15:30:48] [74be16979710d4c4e7c6647856088456]
Feedback Forum

Post a new message
Dataseries X:
18,09
18,13
18
17,72
17,62
17,13
17,39
17,09
17,14
17,38
16,8
16,51
16,01
15,05
13,56
15,22
14,91
15,13
15,25
14,61
14,87
15,1
15,22
15,46
14,96
14
14,2
13,9
13,63
13,32
13,8
14,5
14,12
13,88
14,11
14,26
14,71
14,52
14,32
14,69
15,25
15,04
14,82
14,5
14,72
14,6
14,58
14
14,75
14,41
15,19
14,96
14,83
14,25
14,32
14,93
14,65
15,65
15,65
15,61
15,95
15,83
15,77
16,7
16,69
16,4
16,77
16,78
16,84
16,68
16,67
16,3
16,37
16,6
16,72
16,82
17,5
17,2
17,29
17,2
17,2
17,32
17,16
17,41
17,31
17,3
17,34
17,19
17,05
17,07
17,07
16,81
16,81
16,96
17,05
17
16,77
16,66
16,2
16,26
15,84
15,85
15,71
15,84
15,73
15,77
15,3
15,41
15,4
15,61
15
14,12
14,01
13,46
13,85
13,92
13,59
13,67
13,05
12,87
12,28
11,88
12,49
11,9
10,8
10,99
10,15
10,07
10,05
10,31
9,94
9,65
9,74
9,85
9,96
9,63
9,43
8,77
9,53
9,5
9,78
9,9
9,93
10,35
9,79
9,63
9,02
9,25
9,11
8,95
9,3
9,13
9,75
9,65
9,27
9,59
9,58
9,98
9,57
9,6
9,64
9,46
9,19
9,02
8,9
9,12
8,86
8,94
9
9,23
9,39
9,62
9,9
9,8
9,2
9,87
9,6
9,37
9,21
9,15
8,7
8,2
8,1
6,68
7,7
8,2
7,55
7,53
7,02
6,6
6
3,95
4,91
5,15
5,7
1,93
1,36
1,1
0,98
1
1,1
1,06
1,01
0,93
0,89
0,9
0,88
0,85
0,84
0,94
1
1,1
1,15
1,05
1,06
0,99
0,93
0,84
0,9
0,86
0,78
0,77
0,6
0,57
0,62
0,62
0,58
0,6
0,73
0,75
0,63
0,71
0,68
0,64
0,66
0,69
0,72
0,92
0,85
0,95
1
1,15
1,07
1,01
0,99
0,95
0,92
0,94
0,96
1,05
1,04
1,1
1,14
1,12
1,19
1,35
1,62
1,43
1,45
1,47
1,35
1,15
1,46
1,3
1,3
1,5
1,52
1,63
1,9
1,65
1,5
1,38
1,39




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'George Udny Yule' @ 72.249.76.132

\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 & 2 seconds \tabularnewline
R Server & 'George Udny Yule' @ 72.249.76.132 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=40630&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]2 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'George Udny Yule' @ 72.249.76.132[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=40630&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=40630&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 time2 seconds
R Server'George Udny Yule' @ 72.249.76.132



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))
ari <- array(0,dim=par1)
j <- 0
for (i in 1:n)
{
j = j + 1
ari[j] = ari[j] + 1
arr[j,ari[j]] <- x[i]
if (j == par1) j = 0
}
ari
arr
arr.sd <- array(NA,dim=par1)
arr.range <- array(NA,dim=par1)
arr.iqr <- array(NA,dim=par1)
for (j in 1:par1)
{
arr.sd[j] <- sqrt(var(arr[j,],na.rm=TRUE))
arr.range[j] <- max(arr[j,],na.rm=TRUE) - min(arr[j,],na.rm=TRUE)
arr.iqr[j] <- quantile(arr[j,],0.75,na.rm=TRUE) - quantile(arr[j,],0.25,na.rm=TRUE)
}
overall.sd <- sqrt(var(x))
overall.range <- max(x) - min(x)
overall.iqr <- quantile(x,0.75) - quantile(x,0.25)
bitmap(file='plot1.png')
plot(arr.sd,type='b',ylab='S.D.',main='Standard Deviation Plot',xlab='Periodic Index')
mtext(paste('# blocks = ',np))
abline(overall.sd,0)
dev.off()
bitmap(file='plot2.png')
plot(arr.range,type='b',ylab='range',main='Range Plot',xlab='Periodic Index')
mtext(paste('# blocks = ',np))
abline(overall.range,0)
dev.off()
bitmap(file='plot3.png')
plot(arr.iqr,type='b',ylab='IQR',main='Interquartile Range Plot',xlab='Periodic Index')
mtext(paste('# blocks = ',np))
abline(overall.iqr,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='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.sd,arr.range,arr.iqr))
names(z) <- list('S.D.','Range','IQR')
(boxplot(z,notch=TRUE,col='grey',ylab='Overall Variability',main='Notched Box Plots'))
dev.off()