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 computationFri, 13 Nov 2009 09:08:29 -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/13/t1258128583k8w2tc34c6lcl6p.htm/, Retrieved Sun, 28 Apr 2024 09:52:01 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=56852, Retrieved Sun, 28 Apr 2024 09:52:01 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact187
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Box-Cox Linearity Plot] [3/11/2009] [2009-11-02 21:47:57] [b98453cac15ba1066b407e146608df68]
- RMPD  [Mean Plot] [] [2009-11-09 13:53:17] [023d83ebdf42a2acf423907b4076e8a1]
- RM      [Standard Deviation Plot] [] [2009-11-09 19:25:16] [023d83ebdf42a2acf423907b4076e8a1]
-   PD        [Standard Deviation Plot] [W6 Mean Plot] [2009-11-13 16:08:29] [852eae237d08746109043531619a60c9] [Current]
Feedback Forum

Post a new message
Dataseries X:
112,1294759
123,1474443
119,3806175
121,5465429
124,3716631
85,47917603
88,86932018
128,3268312
141,0398718
117,9680574
101,9590434
103,0890915
111,9411345
120,6990069
111,9411345
111,9411345
117,0263507
89,99936823
68,90513797
113,4478653
105,9142116
119,2864468
114,954596
106,7617476
110,2460624
143,0174559
125,5958818
106,7617476
147,2551361
79,16974108
92,63614701
137,6497276
141,0398718
123,4299564
112,5061585
127,5734659
130,5869273
152,1520109
131,3402927
106,1967236
144,6183573
88,02178414
99,60477665
142,8291145
116,8380094
125,0308578
111,2819398
115,0487667
134,2595835
152,3403523
107,8917957
139,3447997
147,7259894
89,24600286
97,43885122
145,6542347
144,9008693
144,4300159
106,7617476
130,7752687
125,0308578
139,250629
115,0487667
121,9232256
150,7394509
104,1249688
96,59131518
136,8021916
144,0533333
133,8829008
118,3447401
127,1967832
127,0084418
128,3268312
136,3313383
124,7483457
137,1788743
104,9725049
92,54197634
136,4255089
134,1654128
97,8155339
86,42088274
77,28632767
86,89173609




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

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



Parameters (Session):
par1 = 50 ; par2 = 50 ; par3 = 0 ; par4 = 0 ; par5 = 0 ; par6 = Y ; par7 = Y ;
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()