Free Statistics

of Irreproducible Research!

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_sdplot.wasp
Title produced by softwareStandard Deviation Plot
Date of computationWed, 11 Nov 2009 14:51:07 -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/2009/Nov/11/t1257976313dl9il9aqdn3tn89.htm/, Retrieved Fri, 29 Mar 2024 09:30:44 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=55856, Retrieved Fri, 29 Mar 2024 09:30:44 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact153
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Standard Deviation Plot] [3/11/2009] [2009-11-02 22:09:58] [b98453cac15ba1066b407e146608df68]
-   PD    [Standard Deviation Plot] [IndexBouwnijverhe...] [2009-11-11 21:51:07] [a4292616308a56e4faddaa97386e0403] [Current]
Feedback Forum

Post a new message
Dataseries X:
100
99.94940551
102.0743739
102.0237794
102.6309132
102.8838857
103.0103719
104.3258285
105.211232
104.7305844
104.174045
103.6933974
104.3258285
105.211232
105.3124209
105.5400961
106.7037693
106.6025803
105.008854
104.3258285
104.174045
103.4151278
102.3273463
101.8214015
103.5669112
103.8957754
104.5788009
105.0341513
105.6665823
105.6918796
105.7171768
105.843663
105.7930686
105.3883127
105.9701493
106.5013913
107.1338224
109.3599798
109.3599798
108.4239818
107.9433342
108.0951176
108.3986845
110.5489502
111.8138123
112.5474323
111.6620288
111.3837592
113.1039717
115.1783456
121.0726031
123.0710853
123.3999494
122.9193018
122.3880597
123.5517329
124.9683784
124.8671895
123.2734632
121.9580066
122.4892487
125.6767012
126.7644827
126.4356185
125.3478371
126.0055654
127.4475082
130.5843663
133.0887933
133.3417657
132.8358209
133.7971161
136.6304073
138.8818619
140.6526689
143.9160132
149.0260562
149.3296231
152.6182646
161.5987857
162.6612699
166.1269922
165.9752087
165.4692639
166.9365039
167.4930433
169.2891475
170.2757399
171.5911966
171.3635214
171.641791
170.0986592
168.302555
168.9602833
170.0227675
167.5942322
172.6030863
175.082216
177.7637238
182.4943081
185.4793827
189.1474829
186.5924614
184.4927903
178.62383
170.5793069
166.7088287
163.3442955




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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=55856&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



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