Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_smp.wasp
Title produced by softwareStandard Deviation-Mean Plot
Date of computationWed, 01 Apr 2015 16:11:23 +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/2015/Apr/01/t1427902649lebpbl19focftbd.htm/, Retrieved Thu, 31 Oct 2024 23:18:16 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=278519, Retrieved Thu, 31 Oct 2024 23:18:16 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact181
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Standard Deviation-Mean Plot] [] [2015-04-01 15:11:23] [9baff654455058ed055e965df18e01ff] [Current]
Feedback Forum

Post a new message
Dataseries X:
2341
2115
2402
2180
2453
2507
2679
2622
2618
2648
2523
2473
2513
2466
2544
2537
2564
2582
2716
2904
2851
2932
2772
2811
2935
2783
3003
2995
3127
2985
3287
3236
3252
3228
2856
3176
3362
3036
3330
3251
3318
3238
3597
3708
3902
3745
3426
3526
3483
3458
3824
3696
3518
3814
3996
4136
4037
3915
3760
3955
4160
4115
4202
4018
4233
4029
4401
4645
4491
4379
4394
4472
4614
4160
4328
4202
4635
4542
4920
4774
4698
4916
4703
4616
4873
4375
4801
4427
4684
4648
5225
5174
5181
5266
4839
5032
5221
4658
5014
4980
4952
4946
5365
5456
5397
5436
4995
5019
5249
4799
5137
4979
4951
5265
5612
5572
5403
5373
5252
5437
5296
5011
5294
5335
5398
5396
5724
5898
5718
5625
5380
5488
5678
5224
5596
5184
5620
5531
5816
6086
6175
6112
5813
5740
5821
5294
5881
5589
5845
5706
6355
6404
6426
6375
5869
5994
6105
5792
6011
5968
6255
6208
6897
6814
6897
6596
6188
6406
6548
5842
6555
6424
6596
6645
7203
7128
7133
6778
6593
6591
6120
5612
6070
5983
6145
6303
6588
6640
6719
6575
6487
6510
6365
5844
5974
5880
6279
6342
6598
6801
6529
6369
6028
6187
6164
5866
6198
5898
6462
6063
6496
6678
6554
6513
6210
5928
6268
5582
5869
5764
6082
6062
6810
6727
6537
6175
6014
6109




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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=278519&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'Sir Ronald Aylmer Fisher' @ fisher.wessa.net







Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
12463.41666666667179.827927011722564
22682.66666666667166.636967050768466
33071.91666666667168.118712566899504
43453.25248.858530458207866
53799.33333333333224.598766827884678
64294.91666666667197.93040941123627
74592.33333333333249.193000534027760
84877.08333333333305.778427901147891
95119.91666666667250.147759364643798
105252.41666666667249.983075184677813
115463.58333333333241.581327669975887
125714.58333333333317.480266626084991
135963.25361.3550570736871132
146344.75377.1546885442481105
156669.66666666667370.7933345809991361
166312.66666666667331.6680523802331107
176266.33333333333296.525434547078957
186252.5281.532655945339812
196166.58333333333370.8311201840671228

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 2463.41666666667 & 179.827927011722 & 564 \tabularnewline
2 & 2682.66666666667 & 166.636967050768 & 466 \tabularnewline
3 & 3071.91666666667 & 168.118712566899 & 504 \tabularnewline
4 & 3453.25 & 248.858530458207 & 866 \tabularnewline
5 & 3799.33333333333 & 224.598766827884 & 678 \tabularnewline
6 & 4294.91666666667 & 197.93040941123 & 627 \tabularnewline
7 & 4592.33333333333 & 249.193000534027 & 760 \tabularnewline
8 & 4877.08333333333 & 305.778427901147 & 891 \tabularnewline
9 & 5119.91666666667 & 250.147759364643 & 798 \tabularnewline
10 & 5252.41666666667 & 249.983075184677 & 813 \tabularnewline
11 & 5463.58333333333 & 241.581327669975 & 887 \tabularnewline
12 & 5714.58333333333 & 317.480266626084 & 991 \tabularnewline
13 & 5963.25 & 361.355057073687 & 1132 \tabularnewline
14 & 6344.75 & 377.154688544248 & 1105 \tabularnewline
15 & 6669.66666666667 & 370.793334580999 & 1361 \tabularnewline
16 & 6312.66666666667 & 331.668052380233 & 1107 \tabularnewline
17 & 6266.33333333333 & 296.525434547078 & 957 \tabularnewline
18 & 6252.5 & 281.532655945339 & 812 \tabularnewline
19 & 6166.58333333333 & 370.831120184067 & 1228 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=278519&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]2463.41666666667[/C][C]179.827927011722[/C][C]564[/C][/ROW]
[ROW][C]2[/C][C]2682.66666666667[/C][C]166.636967050768[/C][C]466[/C][/ROW]
[ROW][C]3[/C][C]3071.91666666667[/C][C]168.118712566899[/C][C]504[/C][/ROW]
[ROW][C]4[/C][C]3453.25[/C][C]248.858530458207[/C][C]866[/C][/ROW]
[ROW][C]5[/C][C]3799.33333333333[/C][C]224.598766827884[/C][C]678[/C][/ROW]
[ROW][C]6[/C][C]4294.91666666667[/C][C]197.93040941123[/C][C]627[/C][/ROW]
[ROW][C]7[/C][C]4592.33333333333[/C][C]249.193000534027[/C][C]760[/C][/ROW]
[ROW][C]8[/C][C]4877.08333333333[/C][C]305.778427901147[/C][C]891[/C][/ROW]
[ROW][C]9[/C][C]5119.91666666667[/C][C]250.147759364643[/C][C]798[/C][/ROW]
[ROW][C]10[/C][C]5252.41666666667[/C][C]249.983075184677[/C][C]813[/C][/ROW]
[ROW][C]11[/C][C]5463.58333333333[/C][C]241.581327669975[/C][C]887[/C][/ROW]
[ROW][C]12[/C][C]5714.58333333333[/C][C]317.480266626084[/C][C]991[/C][/ROW]
[ROW][C]13[/C][C]5963.25[/C][C]361.355057073687[/C][C]1132[/C][/ROW]
[ROW][C]14[/C][C]6344.75[/C][C]377.154688544248[/C][C]1105[/C][/ROW]
[ROW][C]15[/C][C]6669.66666666667[/C][C]370.793334580999[/C][C]1361[/C][/ROW]
[ROW][C]16[/C][C]6312.66666666667[/C][C]331.668052380233[/C][C]1107[/C][/ROW]
[ROW][C]17[/C][C]6266.33333333333[/C][C]296.525434547078[/C][C]957[/C][/ROW]
[ROW][C]18[/C][C]6252.5[/C][C]281.532655945339[/C][C]812[/C][/ROW]
[ROW][C]19[/C][C]6166.58333333333[/C][C]370.831120184067[/C][C]1228[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=278519&T=1

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

As an alternative you can also use a QR Code:  

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

Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
12463.41666666667179.827927011722564
22682.66666666667166.636967050768466
33071.91666666667168.118712566899504
43453.25248.858530458207866
53799.33333333333224.598766827884678
64294.91666666667197.93040941123627
74592.33333333333249.193000534027760
84877.08333333333305.778427901147891
95119.91666666667250.147759364643798
105252.41666666667249.983075184677813
115463.58333333333241.581327669975887
125714.58333333333317.480266626084991
135963.25361.3550570736871132
146344.75377.1546885442481105
156669.66666666667370.7933345809991361
166312.66666666667331.6680523802331107
176266.33333333333296.525434547078957
186252.5281.532655945339812
196166.58333333333370.8311201840671228







Regression: S.E.(k) = alpha + beta * Mean(k)
alpha48.6944148740964
beta0.0450057948975767
S.D.0.00615542599026836
T-STAT7.31156462099134
p-value1.21515460206128e-06

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 48.6944148740964 \tabularnewline
beta & 0.0450057948975767 \tabularnewline
S.D. & 0.00615542599026836 \tabularnewline
T-STAT & 7.31156462099134 \tabularnewline
p-value & 1.21515460206128e-06 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=278519&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]48.6944148740964[/C][/ROW]
[ROW][C]beta[/C][C]0.0450057948975767[/C][/ROW]
[ROW][C]S.D.[/C][C]0.00615542599026836[/C][/ROW]
[ROW][C]T-STAT[/C][C]7.31156462099134[/C][/ROW]
[ROW][C]p-value[/C][C]1.21515460206128e-06[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=278519&T=2

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

As an alternative you can also use a QR Code:  

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

Regression: S.E.(k) = alpha + beta * Mean(k)
alpha48.6944148740964
beta0.0450057948975767
S.D.0.00615542599026836
T-STAT7.31156462099134
p-value1.21515460206128e-06







Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha-0.817167886406544
beta0.75469743710828
S.D.0.0989364695654391
T-STAT7.62810155267471
p-value6.93558731192016e-07
Lambda0.24530256289172

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & -0.817167886406544 \tabularnewline
beta & 0.75469743710828 \tabularnewline
S.D. & 0.0989364695654391 \tabularnewline
T-STAT & 7.62810155267471 \tabularnewline
p-value & 6.93558731192016e-07 \tabularnewline
Lambda & 0.24530256289172 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=278519&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-0.817167886406544[/C][/ROW]
[ROW][C]beta[/C][C]0.75469743710828[/C][/ROW]
[ROW][C]S.D.[/C][C]0.0989364695654391[/C][/ROW]
[ROW][C]T-STAT[/C][C]7.62810155267471[/C][/ROW]
[ROW][C]p-value[/C][C]6.93558731192016e-07[/C][/ROW]
[ROW][C]Lambda[/C][C]0.24530256289172[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=278519&T=3

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

As an alternative you can also use a 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.817167886406544
beta0.75469743710828
S.D.0.0989364695654391
T-STAT7.62810155267471
p-value6.93558731192016e-07
Lambda0.24530256289172



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