FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*Unverified author*
R Software Modulerwasp_smp.wasp
Title produced by softwareStandard Deviation-Mean Plot
Date of computationSun, 23 Nov 2014 16:09:59 +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/2014/Nov/23/t1416759027rmvnhr688rqy1gz.htm/, Retrieved Tue, 28 May 2024 16:26:44 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=258039, Retrieved Tue, 28 May 2024 16:26:44 +0000
QR Codes:

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

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] [Verkopen BMW] [2014-11-23 16:09:59] [76c30f62b7052b57088120e90a652e05] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
2135
1157
1290
1071
1169
1431
945
1034
1100
1297
921
236
1990
966
1326
908
1206
1861
929
1296
1332
1352
1040
148
2090
1435
1124
1319
1436
1774
1566
1385
1147
1274
625
52
1990
1154
954
887
825
966
954
770
1838
1371
589
116
1898
712
1175
1240
1329
1550
1201
938
1030
1060
1035
635
2565
910
1304
1331
1681
1983
1021
1061
1292
1274
1024
568
2570
1125
1600
1492
2492
3523
990
869
1310
979
1244
442
2956
1055
2004
1462
1144
1454
4060
1538
1388
1547
4473
1570
1535
1352

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Gwilym Jenkins' @ jenkins.wessa.net

\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' @ jenkins.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=258039&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' @ jenkins.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=258039&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=258039&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 Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'Gwilym Jenkins' @ jenkins.wessa.net






Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
11148.83333333333430.9448041715141899
21196.16666666667474.0795547610381842
31268.91666666667524.0085804872312038
41034.5511.8129097096461874
51150.25343.8562397810511263
61334.5530.7260036111771997
71553878.4891989815663081
82054.251145.321798765433418

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 1148.83333333333 & 430.944804171514 & 1899 \tabularnewline
2 & 1196.16666666667 & 474.079554761038 & 1842 \tabularnewline
3 & 1268.91666666667 & 524.008580487231 & 2038 \tabularnewline
4 & 1034.5 & 511.812909709646 & 1874 \tabularnewline
5 & 1150.25 & 343.856239781051 & 1263 \tabularnewline
6 & 1334.5 & 530.726003611177 & 1997 \tabularnewline
7 & 1553 & 878.489198981566 & 3081 \tabularnewline
8 & 2054.25 & 1145.32179876543 & 3418 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=258039&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]1148.83333333333[/C][C]430.944804171514[/C][C]1899[/C][/ROW]
[ROW][C]2[/C][C]1196.16666666667[/C][C]474.079554761038[/C][C]1842[/C][/ROW]
[ROW][C]3[/C][C]1268.91666666667[/C][C]524.008580487231[/C][C]2038[/C][/ROW]
[ROW][C]4[/C][C]1034.5[/C][C]511.812909709646[/C][C]1874[/C][/ROW]
[ROW][C]5[/C][C]1150.25[/C][C]343.856239781051[/C][C]1263[/C][/ROW]
[ROW][C]6[/C][C]1334.5[/C][C]530.726003611177[/C][C]1997[/C][/ROW]
[ROW][C]7[/C][C]1553[/C][C]878.489198981566[/C][C]3081[/C][/ROW]
[ROW][C]8[/C][C]2054.25[/C][C]1145.32179876543[/C][C]3418[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=258039&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=258039&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
11148.83333333333430.9448041715141899
21196.16666666667474.0795547610381842
31268.91666666667524.0085804872312038
41034.5511.8129097096461874
51150.25343.8562397810511263
61334.5530.7260036111771997
71553878.4891989815663081
82054.251145.321798765433418






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha-438.181733389776
beta0.776943131385672
S.D.0.107851448951459
T-STAT7.20382654975133
p-value0.000362263140876121

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & -438.181733389776 \tabularnewline
beta & 0.776943131385672 \tabularnewline
S.D. & 0.107851448951459 \tabularnewline
T-STAT & 7.20382654975133 \tabularnewline
p-value & 0.000362263140876121 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=258039&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-438.181733389776[/C][/ROW]
[ROW][C]beta[/C][C]0.776943131385672[/C][/ROW]
[ROW][C]S.D.[/C][C]0.107851448951459[/C][/ROW]
[ROW][C]T-STAT[/C][C]7.20382654975133[/C][/ROW]
[ROW][C]p-value[/C][C]0.000362263140876121[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=258039&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=258039&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-438.181733389776
beta0.776943131385672
S.D.0.107851448951459
T-STAT7.20382654975133
p-value0.000362263140876121






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha-5.25129360493126
beta1.61337520314422
S.D.0.319781658311324
T-STAT5.04523996674481
p-value0.00234464875907615
Lambda-0.613375203144222

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & -5.25129360493126 \tabularnewline
beta & 1.61337520314422 \tabularnewline
S.D. & 0.319781658311324 \tabularnewline
T-STAT & 5.04523996674481 \tabularnewline
p-value & 0.00234464875907615 \tabularnewline
Lambda & -0.613375203144222 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=258039&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-5.25129360493126[/C][/ROW]
[ROW][C]beta[/C][C]1.61337520314422[/C][/ROW]
[ROW][C]S.D.[/C][C]0.319781658311324[/C][/ROW]
[ROW][C]T-STAT[/C][C]5.04523996674481[/C][/ROW]
[ROW][C]p-value[/C][C]0.00234464875907615[/C][/ROW]
[ROW][C]Lambda[/C][C]-0.613375203144222[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=258039&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=258039&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.25129360493126
beta1.61337520314422
S.D.0.319781658311324
T-STAT5.04523996674481
p-value0.00234464875907615
Lambda-0.613375203144222


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