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, 29 Dec 2010 12:28:03 +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/2010/Dec/29/t12936255449ruynhiy2qyl8vm.htm/, Retrieved Fri, 03 May 2024 09:56:22 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=116751, Retrieved Fri, 03 May 2024 09:56:22 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact117
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]
- R PD  [Standard Deviation Plot] [] [2009-12-19 15:17:06] [69c775ce4d55db2aa75a88e773e8d700]
-         [Standard Deviation Plot] [] [2010-12-26 18:03:29] [69c775ce4d55db2aa75a88e773e8d700]
-    D        [Standard Deviation Plot] [] [2010-12-29 12:28:03] [e180d4cd19004beeddc12e67012247dc] [Current]
Feedback Forum

Post a new message
Dataseries X:
09,166456
07,970589
07,104091
06,621064
07,529215
08,170938
08,157450
07,378962
07,921496
08,156740
08,856365
08,817177
08,734347
09,345927
08,992970
10,785120
08,886867
08,818847
08,823744
09,165298
08,652657
08,173054
07,563416
07,595809
08,381467
07,216432
06,540178
06,238914
05,487288
05,759462
05,993215
07,474726
07,348907
07,303379
07,119314
06,993780
06,958153
07,595706
08,088153
07,555753
07,315433
07,893427
08,858794
08,839367
08,014733
07,873465
08,930377
10,500550
12,611440
11,417870
11,872490
11,060820
12,043310
09,776299
09,557194
09,202590
10,224020
09,350807
08,300913
08,365779
08,133595
07,660470
08,074839
07,848597
07,998220
07,396895
07,900419
08,100500
07,899453
07,599783
08,100929
09,002175
10,298900
10,101520
10,699150
09,698140
09,800951
10,900470
10,697850
09,297252
10,397440
10,900720
12,901270
13,099060
11,698280
11,099870
11,301570
10,702110
10,099310
09,591119




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24

\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 & 'Sir Ronald Aylmer Fisher' @ 193.190.124.24 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=116751&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]'Sir Ronald Aylmer Fisher' @ 193.190.124.24[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=116751&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=116751&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'Sir Ronald Aylmer Fisher' @ 193.190.124.24



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