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 computationWed, 28 Nov 2012 13:33:03 -0500
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2012/Nov/28/t1354127704ptnn7omjc04rouf.htm/, Retrieved Wed, 24 Apr 2024 22:22:19 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=194211, Retrieved Wed, 24 Apr 2024 22:22:19 +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] [inschrijving nieu...] [2012-11-28 18:22:24] [6e17e56f9bbdc65ea77b7e28bc808fc4]
-    D    [Standard Deviation-Mean Plot] [verkoop aantal pr...] [2012-11-28 18:33:03] [32426a9b6b43a24d932e8bdfe421cad9] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
124275
58605
21828
42811
65963
128457
26867
143540
107458
258558
22475
294584
159848
289458
117950
78351
84589
44324
52285
185486
23846
13257
166999
148488
74747
55700
218584
187888
44788
18840
48787
69100
41892
90588
148574
50201
86828
102785
118844
145288
56790
287525
187880
87740
55258
58769
43366
77051
91574
15533
18425
65192
81059
73322
91261
86166
61842
25192
21059
15855
12618
101667
224275
55700
60748
41848
61781
120077
42032
46485
36861
55027
48999
68352
126987
86526
125340
69029
153287
135724
92108
119906
79798
97206

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=194211&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'Gertrude Mary Cox' @ cox.wessa.net






Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
1107951.7590054.178810844272756
2113740.08333333380389.4752931098276201
387474.083333333363260.8436604254199744
4109010.33333333370035.6818243819244159
553873.333333333331969.126721468876041
671593.2555906.7713318342211657
7100271.83333333331770.6908324276104288

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 107951.75 & 90054.178810844 & 272756 \tabularnewline
2 & 113740.083333333 & 80389.4752931098 & 276201 \tabularnewline
3 & 87474.0833333333 & 63260.8436604254 & 199744 \tabularnewline
4 & 109010.333333333 & 70035.6818243819 & 244159 \tabularnewline
5 & 53873.3333333333 & 31969.1267214688 & 76041 \tabularnewline
6 & 71593.25 & 55906.7713318342 & 211657 \tabularnewline
7 & 100271.833333333 & 31770.6908324276 & 104288 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=194211&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]107951.75[/C][C]90054.178810844[/C][C]272756[/C][/ROW]
[ROW][C]2[/C][C]113740.083333333[/C][C]80389.4752931098[/C][C]276201[/C][/ROW]
[ROW][C]3[/C][C]87474.0833333333[/C][C]63260.8436604254[/C][C]199744[/C][/ROW]
[ROW][C]4[/C][C]109010.333333333[/C][C]70035.6818243819[/C][C]244159[/C][/ROW]
[ROW][C]5[/C][C]53873.3333333333[/C][C]31969.1267214688[/C][C]76041[/C][/ROW]
[ROW][C]6[/C][C]71593.25[/C][C]55906.7713318342[/C][C]211657[/C][/ROW]
[ROW][C]7[/C][C]100271.833333333[/C][C]31770.6908324276[/C][C]104288[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=194211&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=194211&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
1107951.7590054.178810844272756
2113740.08333333380389.4752931098276201
387474.083333333363260.8436604254199744
4109010.33333333370035.6818243819244159
553873.333333333331969.126721468876041
671593.2555906.7713318342211657
7100271.83333333331770.6908324276104288






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha-1268.24331959565
beta0.671307075438041
S.D.0.336761245684702
T-STAT1.99342140474965
p-value0.102799650948082

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & -1268.24331959565 \tabularnewline
beta & 0.671307075438041 \tabularnewline
S.D. & 0.336761245684702 \tabularnewline
T-STAT & 1.99342140474965 \tabularnewline
p-value & 0.102799650948082 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=194211&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-1268.24331959565[/C][/ROW]
[ROW][C]beta[/C][C]0.671307075438041[/C][/ROW]
[ROW][C]S.D.[/C][C]0.336761245684702[/C][/ROW]
[ROW][C]T-STAT[/C][C]1.99342140474965[/C][/ROW]
[ROW][C]p-value[/C][C]0.102799650948082[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=194211&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=194211&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-1268.24331959565
beta0.671307075438041
S.D.0.336761245684702
T-STAT1.99342140474965
p-value0.102799650948082






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha-0.087871432974524
beta0.967439862011448
S.D.0.529017361396222
T-STAT1.82874879466736
p-value0.126971873783101
Lambda0.0325601379885515

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & -0.087871432974524 \tabularnewline
beta & 0.967439862011448 \tabularnewline
S.D. & 0.529017361396222 \tabularnewline
T-STAT & 1.82874879466736 \tabularnewline
p-value & 0.126971873783101 \tabularnewline
Lambda & 0.0325601379885515 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=194211&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-0.087871432974524[/C][/ROW]
[ROW][C]beta[/C][C]0.967439862011448[/C][/ROW]
[ROW][C]S.D.[/C][C]0.529017361396222[/C][/ROW]
[ROW][C]T-STAT[/C][C]1.82874879466736[/C][/ROW]
[ROW][C]p-value[/C][C]0.126971873783101[/C][/ROW]
[ROW][C]Lambda[/C][C]0.0325601379885515[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=194211&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=194211&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-0.087871432974524
beta0.967439862011448
S.D.0.529017361396222
T-STAT1.82874879466736
p-value0.126971873783101
Lambda0.0325601379885515


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