FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_smp.wasp
Title produced by softwareStandard Deviation-Mean Plot
Date of computationSun, 09 Dec 2018 17:48:27 +0100
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2018/Dec/09/t1544374135w83s6dfoqi0m0va.htm/, Retrieved Sat, 27 Apr 2024 15:28:40 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=315818, Retrieved Sat, 27 Apr 2024 15:28:40 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact110

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Standard Deviation-Mean Plot] [bouwv: std mean plot] [2018-12-09 16:48:27] [ed325faff554f097403cb8f6ee07a063] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
2570
2669
2450
2842
3440
2678
2981
2260
2844
2546
2456
2295
2379
2471
2057
2280
2351
2276
2548
2311
2201
2725
2408
2139
1898
2539
2070
2063
2565
2443
2196
2799
2076
2628
2292
2155
2476
2138
1854
2081
1795
1756
2237
1960
1829
2524
2077
2366
2185
2098
1836
1863
2044
2136
2931
3263
3328
3570
2313
1623
1316
1507
1419
1660
1790
1733
2086
1814
2241
1943
1773
2143
2087
1805
1913
2296
2500
2210
2526
2249
2024
2091
2045
1882
1831
1964
1763
1688
2149
1823
2094
2145
1791
1996
2097
1796
1963
2042
1746
2210
2968
3126
3708
3015
1569
1518
1393
1615
1777
1648
1463
1779

Tables (Output of Computation)





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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=315818&T=0

As an alternative you can also use a QR Code:  
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






Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
12669.25327.2247281720511180
22345.5181.634247871925668
32310.33333333333278.935096176132901
42091.08333333333265.77073453135768
52432.5660.2866457421771947
61785.41666666667286.388629260481925
72135.66666666667229.128123480504721
81928.08333333333164.984549322509461
92239.41666666667769.3594098437682315

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 2669.25 & 327.224728172051 & 1180 \tabularnewline
2 & 2345.5 & 181.634247871925 & 668 \tabularnewline
3 & 2310.33333333333 & 278.935096176132 & 901 \tabularnewline
4 & 2091.08333333333 & 265.77073453135 & 768 \tabularnewline
5 & 2432.5 & 660.286645742177 & 1947 \tabularnewline
6 & 1785.41666666667 & 286.388629260481 & 925 \tabularnewline
7 & 2135.66666666667 & 229.128123480504 & 721 \tabularnewline
8 & 1928.08333333333 & 164.984549322509 & 461 \tabularnewline
9 & 2239.41666666667 & 769.359409843768 & 2315 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=315818&T=1

[TABLE]
[ROW][C]Standard Deviation-Mean Plot[/C][/ROW]
[ROW][C]Section[/C][C]Mean[/C][C]Standard Deviation[/C][C]Range[/C][/ROW]
[ROW][C]1[/C][C]2669.25[/C][C]327.224728172051[/C][C]1180[/C][/ROW]
[ROW][C]2[/C][C]2345.5[/C][C]181.634247871925[/C][C]668[/C][/ROW]
[ROW][C]3[/C][C]2310.33333333333[/C][C]278.935096176132[/C][C]901[/C][/ROW]
[ROW][C]4[/C][C]2091.08333333333[/C][C]265.77073453135[/C][C]768[/C][/ROW]
[ROW][C]5[/C][C]2432.5[/C][C]660.286645742177[/C][C]1947[/C][/ROW]
[ROW][C]6[/C][C]1785.41666666667[/C][C]286.388629260481[/C][C]925[/C][/ROW]
[ROW][C]7[/C][C]2135.66666666667[/C][C]229.128123480504[/C][C]721[/C][/ROW]
[ROW][C]8[/C][C]1928.08333333333[/C][C]164.984549322509[/C][C]461[/C][/ROW]
[ROW][C]9[/C][C]2239.41666666667[/C][C]769.359409843768[/C][C]2315[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=315818&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=315818&T=1

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
12669.25327.2247281720511180
22345.5181.634247871925668
32310.33333333333278.935096176132901
42091.08333333333265.77073453135768
52432.5660.2866457421771947
61785.41666666667286.388629260481925
72135.66666666667229.128123480504721
81928.08333333333164.984549322509461
92239.41666666667769.3594098437682315






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha-190.018926362714
beta0.24446112185308
S.D.0.288751668455603
T-STAT0.846613711915804
p-value0.425204096620899

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & -190.018926362714 \tabularnewline
beta & 0.24446112185308 \tabularnewline
S.D. & 0.288751668455603 \tabularnewline
T-STAT & 0.846613711915804 \tabularnewline
p-value & 0.425204096620899 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=315818&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-190.018926362714[/C][/ROW]
[ROW][C]beta[/C][C]0.24446112185308[/C][/ROW]
[ROW][C]S.D.[/C][C]0.288751668455603[/C][/ROW]
[ROW][C]T-STAT[/C][C]0.846613711915804[/C][/ROW]
[ROW][C]p-value[/C][C]0.425204096620899[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=315818&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=315818&T=2

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Regression: S.E.(k) = alpha + beta * Mean(k)
alpha-190.018926362714
beta0.24446112185308
S.D.0.288751668455603
T-STAT0.846613711915804
p-value0.425204096620899






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha-5.92213324419226
beta1.5135309675672
S.D.1.52931306213866
T-STAT0.989680272167822
p-value0.355298524251314
Lambda-0.5135309675672

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & -5.92213324419226 \tabularnewline
beta & 1.5135309675672 \tabularnewline
S.D. & 1.52931306213866 \tabularnewline
T-STAT & 0.989680272167822 \tabularnewline
p-value & 0.355298524251314 \tabularnewline
Lambda & -0.5135309675672 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=315818&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-5.92213324419226[/C][/ROW]
[ROW][C]beta[/C][C]1.5135309675672[/C][/ROW]
[ROW][C]S.D.[/C][C]1.52931306213866[/C][/ROW]
[ROW][C]T-STAT[/C][C]0.989680272167822[/C][/ROW]
[ROW][C]p-value[/C][C]0.355298524251314[/C][/ROW]
[ROW][C]Lambda[/C][C]-0.5135309675672[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=315818&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=315818&T=3

As an alternative you can also use a QR Code:  
QR Code


The GUIDs for individual cells are displayed in the table below:

Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha-5.92213324419226
beta1.5135309675672
S.D.1.52931306213866
T-STAT0.989680272167822
p-value0.355298524251314
Lambda-0.5135309675672


Figures (Output of Computation)

Input Parameters & R Code

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))
j <- 0
k <- 1
for (i in 1:(np*par1))
{
j = j + 1
arr[j,k] <- x[i]
if (j == par1) {
j = 0
k=k+1
}
}
arr
arr.mean <- array(NA,dim=np)
arr.sd <- array(NA,dim=np)
arr.range <- array(NA,dim=np)
for (j in 1:np)
{
arr.mean[j] <- mean(arr[,j],na.rm=TRUE)
arr.sd[j] <- sd(arr[,j],na.rm=TRUE)
arr.range[j] <- max(arr[,j],na.rm=TRUE) - min(arr[,j],na.rm=TRUE)
}
arr.mean
arr.sd
arr.range
(lm1 <- lm(arr.sd~arr.mean))
(lnlm1 <- lm(log(arr.sd)~log(arr.mean)))
(lm2 <- lm(arr.range~arr.mean))
bitmap(file='test1.png')
plot(arr.mean,arr.sd,main='Standard Deviation-Mean Plot',xlab='mean',ylab='standard deviation')
dev.off()
bitmap(file='test2.png')
plot(arr.mean,arr.range,main='Range-Mean Plot',xlab='mean',ylab='range')
dev.off()
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Standard Deviation-Mean Plot',4,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Section',header=TRUE)
a<-table.element(a,'Mean',header=TRUE)
a<-table.element(a,'Standard Deviation',header=TRUE)
a<-table.element(a,'Range',header=TRUE)
a<-table.row.end(a)
for (j in 1:np) {
a<-table.row.start(a)
a<-table.element(a,j,header=TRUE)
a<-table.element(a,arr.mean[j])
a<-table.element(a,arr.sd[j] )
a<-table.element(a,arr.range[j] )
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Regression: S.E.(k) = alpha + beta * Mean(k)',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'alpha',header=TRUE)
a<-table.element(a,lm1$coefficients[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'beta',header=TRUE)
a<-table.element(a,lm1$coefficients[[2]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'S.D.',header=TRUE)
a<-table.element(a,summary(lm1)$coefficients[2,2])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'T-STAT',header=TRUE)
a<-table.element(a,summary(lm1)$coefficients[2,3])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-value',header=TRUE)
a<-table.element(a,summary(lm1)$coefficients[2,4])
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable1.tab')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Regression: ln S.E.(k) = alpha + beta * ln Mean(k)',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'alpha',header=TRUE)
a<-table.element(a,lnlm1$coefficients[[1]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'beta',header=TRUE)
a<-table.element(a,lnlm1$coefficients[[2]])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'S.D.',header=TRUE)
a<-table.element(a,summary(lnlm1)$coefficients[2,2])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'T-STAT',header=TRUE)
a<-table.element(a,summary(lnlm1)$coefficients[2,3])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'p-value',header=TRUE)
a<-table.element(a,summary(lnlm1)$coefficients[2,4])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Lambda',header=TRUE)
a<-table.element(a,1-lnlm1$coefficients[[2]])
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable2.tab')