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 computationFri, 08 Jan 2016 20:20:34 +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/2016/Jan/08/t1452284468awkjb738u4rd4mr.htm/, Retrieved Sun, 28 Apr 2024 09:18:53 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=287438, Retrieved Sun, 28 Apr 2024 09:18:53 +0000
QR Codes:

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

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] [Aantal niet-werke...] [2016-01-08 20:20:34] [f442d180d44854b5d66611a6a05f7502] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
64076
63136
60198
59057
57388
56708
70019
72263
74152
67057
61941
58331
59252
56568
53031
51840
48290
45817
59421
61621
60976
57497
53037
53088
53119
51644
47866
47691
42401
43069
55797
57170
58335
55439
54399
56316
58381
58468
59025
58298
54255
55670
67816
70485
71361
66953
64505
66770
66418
65277
62008
59096
55106
54954
67943
69411
69951
63966
60410
59440
59445
57614
55396
53030
50090
48764
61658
63943
64878
60634
57905
57224
60953
60621
57258
54903
53278
53042
63753
69210
71446
68408
65427
64630
66086
65058
62689
60841
57346
56222
68202
70745
73690
68992
65925
65546
67221
65315
62038
58774
55320
53900
65544
67906
70911
66544
63657
61720

Tables (Output of Computation)





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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=287438&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 time2 seconds
R Server'Herman Ole Andreas Wold' @ wold.wessa.net






Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
163693.83333333335944.4290318878817444
255036.55014.3056348810715804
351937.16666666675447.1062515019415934
462665.58333333335953.4181237387217106
562831.66666666675169.5918483752514997
657548.41666666675069.5136660833716114
761910.756296.3409449818218404
865111.83333333335180.13176093117468
963237.55141.9026280515717011

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 63693.8333333333 & 5944.42903188788 & 17444 \tabularnewline
2 & 55036.5 & 5014.30563488107 & 15804 \tabularnewline
3 & 51937.1666666667 & 5447.10625150194 & 15934 \tabularnewline
4 & 62665.5833333333 & 5953.41812373872 & 17106 \tabularnewline
5 & 62831.6666666667 & 5169.59184837525 & 14997 \tabularnewline
6 & 57548.4166666667 & 5069.51366608337 & 16114 \tabularnewline
7 & 61910.75 & 6296.34094498182 & 18404 \tabularnewline
8 & 65111.8333333333 & 5180.131760931 & 17468 \tabularnewline
9 & 63237.5 & 5141.90262805157 & 17011 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=287438&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]63693.8333333333[/C][C]5944.42903188788[/C][C]17444[/C][/ROW]
[ROW][C]2[/C][C]55036.5[/C][C]5014.30563488107[/C][C]15804[/C][/ROW]
[ROW][C]3[/C][C]51937.1666666667[/C][C]5447.10625150194[/C][C]15934[/C][/ROW]
[ROW][C]4[/C][C]62665.5833333333[/C][C]5953.41812373872[/C][C]17106[/C][/ROW]
[ROW][C]5[/C][C]62831.6666666667[/C][C]5169.59184837525[/C][C]14997[/C][/ROW]
[ROW][C]6[/C][C]57548.4166666667[/C][C]5069.51366608337[/C][C]16114[/C][/ROW]
[ROW][C]7[/C][C]61910.75[/C][C]6296.34094498182[/C][C]18404[/C][/ROW]
[ROW][C]8[/C][C]65111.8333333333[/C][C]5180.131760931[/C][C]17468[/C][/ROW]
[ROW][C]9[/C][C]63237.5[/C][C]5141.90262805157[/C][C]17011[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=287438&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=287438&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
163693.83333333335944.4290318878817444
255036.55014.3056348810715804
351937.16666666675447.1062515019415934
462665.58333333335953.4181237387217106
562831.66666666675169.5918483752514997
657548.41666666675069.5136660833716114
761910.756296.3409449818218404
865111.83333333335180.13176093117468
963237.55141.9026280515717011






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha3740.05790520222
beta0.0285973965514161
S.D.0.0381371199349498
T-STAT0.749857267675024
p-value0.477780588556407

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 3740.05790520222 \tabularnewline
beta & 0.0285973965514161 \tabularnewline
S.D. & 0.0381371199349498 \tabularnewline
T-STAT & 0.749857267675024 \tabularnewline
p-value & 0.477780588556407 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=287438&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]3740.05790520222[/C][/ROW]
[ROW][C]beta[/C][C]0.0285973965514161[/C][/ROW]
[ROW][C]S.D.[/C][C]0.0381371199349498[/C][/ROW]
[ROW][C]T-STAT[/C][C]0.749857267675024[/C][/ROW]
[ROW][C]p-value[/C][C]0.477780588556407[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=287438&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=287438&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)
alpha3740.05790520222
beta0.0285973965514161
S.D.0.0381371199349498
T-STAT0.749857267675024
p-value0.477780588556407






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha5.33826099058315
beta0.296659200228142
S.D.0.397510886876874
T-STAT0.746292013682709
p-value0.479798791276259
Lambda0.703340799771858

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & 5.33826099058315 \tabularnewline
beta & 0.296659200228142 \tabularnewline
S.D. & 0.397510886876874 \tabularnewline
T-STAT & 0.746292013682709 \tabularnewline
p-value & 0.479798791276259 \tabularnewline
Lambda & 0.703340799771858 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=287438&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]5.33826099058315[/C][/ROW]
[ROW][C]beta[/C][C]0.296659200228142[/C][/ROW]
[ROW][C]S.D.[/C][C]0.397510886876874[/C][/ROW]
[ROW][C]T-STAT[/C][C]0.746292013682709[/C][/ROW]
[ROW][C]p-value[/C][C]0.479798791276259[/C][/ROW]
[ROW][C]Lambda[/C][C]0.703340799771858[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=287438&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=287438&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)
alpha5.33826099058315
beta0.296659200228142
S.D.0.397510886876874
T-STAT0.746292013682709
p-value0.479798791276259
Lambda0.703340799771858


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