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, 27 Jul 2017 23:07:55 +0200
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2017/Jul/27/t1501189743u304j8ynhf5cvr1.htm/, Retrieved Thu, 31 Oct 2024 23:51:37 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=306771, Retrieved Thu, 31 Oct 2024 23:51:37 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact146
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Standard Deviation Plot] [Reeks A stap 24] [2017-07-27 21:07:55] [5e513ceaaef205c0c6f269c0b513af8d] [Current]
Feedback Forum

Post a new message
Dataseries X:
6195800
6172725
6149325
6100900
6579950
6554600
6195800
5957250
5980325
5980325
6006000
6052150
6123975
6123975
6077825
5957250
6579950
6674850
6531525
6195800
6339450
6123975
6221150
6267625
6316050
6195800
6221150
6052150
6579950
6746675
6603350
6339450
6626425
6316050
6603350
6579950
6651775
6387875
6674850
6651775
7082400
6985225
6603350
6410950
6674850
6316050
6579950
6626425
6723600
6508450
6626425
6698250
6962150
6746675
6459700
6149325
6436625
5646875
6029075
6244225
6459700
6149325
6149325
6149325
6316050
6077825
5765175
5503550
5693350
4952350
5406375
5670275
5718700
5454800
5477875
5406375
5646875
5477875
5144750
4903925
5311150
4426825
5001100
5262725
5262725
4952350
4665375
4642300
4903925
4665375
4211675
3899025
4234750
3445325
4162925
4544800
4665375
4401475
4068025
4306575
4401475
4329650
3611725
3278600
3516825
2799225
3540225
3804125
4019275
3660475
3324750
3516825
3611725
3421925
2704325
2391675
2678650
1889225
2750475
3278600




Summary of computational transaction
Raw Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R ServerBig Analytics Cloud Computing Center

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input view raw input (R code)  \tabularnewline
Raw Outputview raw output of R engine  \tabularnewline
Computing time2 seconds \tabularnewline
R ServerBig Analytics Cloud Computing Center \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=306771&T=0

[TABLE]
[ROW]
Summary of computational transaction[/C][/ROW] [ROW]Raw Input[/C] view raw input (R code) [/C][/ROW] [ROW]Raw Output[/C]view raw output of R engine [/C][/ROW] [ROW]Computing time[/C]2 seconds[/C][/ROW] [ROW]R Server[/C]Big Analytics Cloud Computing Center[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=306771&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=306771&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 Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R ServerBig Analytics Cloud Computing Center



Parameters (Session):
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()