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 computationThu, 19 May 2011 14:30:20 +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/2011/May/19/t1305815200kskf9nzpgclqsh0.htm/, Retrieved Sun, 12 May 2024 00:52:49 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=122041, Retrieved Sun, 12 May 2024 00:52:49 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordsKEYWORD: KDGP2W83
Estimated Impact84

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] [Spreidings-en gem...] [2011-05-19 14:30:20] [be417f314f65e9d8a38b0902dfa3287c] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
32819
32700
32242
32810
33865
32226
31077
31293
30236
30160
32436
30695
27525
26434
25739
25204
24977
24320
22680
22052
21467
21383
21777
21928
21814
22937
23595
20830
19650
19195
19644
18483
18079
19178
18391
18441
18584
20108
20148
19394
17745
17696
17032
16438
15683
15594
15713
15937
16171
15928
16348
15579
15305
15648
14954
15137
15839
16050
15168
17064
16005
14886
14931
14544
13812
13031
12574
11964
11451
11346
11353
10702
10646
10556
10463
10407
10625
10872
10805
10653
10574
10431
10383
10296
10872
10635
10297
10570
10662
10709
10413
10846
10371
9924
9828

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'George Udny Yule' @ 216.218.223.82

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=122041&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'George Udny Yule' @ 216.218.223.82






Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
131879.91666666671166.731289335793705
223790.52161.454876192936142
320019.751861.2958889885516
4175061721.294438919314554
515765.9166666667601.0434046243452110
613049.91666666671750.673244956375303
710559.25172.765803118767576

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 31879.9166666667 & 1166.73128933579 & 3705 \tabularnewline
2 & 23790.5 & 2161.45487619293 & 6142 \tabularnewline
3 & 20019.75 & 1861.295888988 & 5516 \tabularnewline
4 & 17506 & 1721.29443891931 & 4554 \tabularnewline
5 & 15765.9166666667 & 601.043404624345 & 2110 \tabularnewline
6 & 13049.9166666667 & 1750.67324495637 & 5303 \tabularnewline
7 & 10559.25 & 172.765803118767 & 576 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=122041&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]31879.9166666667[/C][C]1166.73128933579[/C][C]3705[/C][/ROW]
[ROW][C]2[/C][C]23790.5[/C][C]2161.45487619293[/C][C]6142[/C][/ROW]
[ROW][C]3[/C][C]20019.75[/C][C]1861.295888988[/C][C]5516[/C][/ROW]
[ROW][C]4[/C][C]17506[/C][C]1721.29443891931[/C][C]4554[/C][/ROW]
[ROW][C]5[/C][C]15765.9166666667[/C][C]601.043404624345[/C][C]2110[/C][/ROW]
[ROW][C]6[/C][C]13049.9166666667[/C][C]1750.67324495637[/C][C]5303[/C][/ROW]
[ROW][C]7[/C][C]10559.25[/C][C]172.765803118767[/C][C]576[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=122041&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=122041&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
131879.91666666671166.731289335793705
223790.52161.454876192936142
320019.751861.2958889885516
4175061721.294438919314554
515765.9166666667601.0434046243452110
613049.91666666671750.673244956375303
710559.25172.765803118767576






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha645.4938638515
beta0.0370879953170466
S.D.0.0423545841802254
T-STAT0.875654808915873
p-value0.421284965086895

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 645.4938638515 \tabularnewline
beta & 0.0370879953170466 \tabularnewline
S.D. & 0.0423545841802254 \tabularnewline
T-STAT & 0.875654808915873 \tabularnewline
p-value & 0.421284965086895 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=122041&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]645.4938638515[/C][/ROW]
[ROW][C]beta[/C][C]0.0370879953170466[/C][/ROW]
[ROW][C]S.D.[/C][C]0.0423545841802254[/C][/ROW]
[ROW][C]T-STAT[/C][C]0.875654808915873[/C][/ROW]
[ROW][C]p-value[/C][C]0.421284965086895[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=122041&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=122041&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)
alpha645.4938638515
beta0.0370879953170466
S.D.0.0423545841802254
T-STAT0.875654808915873
p-value0.421284965086895






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha-7.11178585171747
beta1.4376689308096
S.D.0.887442993530035
T-STAT1.62001271213028
p-value0.16615512459958
Lambda-0.437668930809604

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & -7.11178585171747 \tabularnewline
beta & 1.4376689308096 \tabularnewline
S.D. & 0.887442993530035 \tabularnewline
T-STAT & 1.62001271213028 \tabularnewline
p-value & 0.16615512459958 \tabularnewline
Lambda & -0.437668930809604 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=122041&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-7.11178585171747[/C][/ROW]
[ROW][C]beta[/C][C]1.4376689308096[/C][/ROW]
[ROW][C]S.D.[/C][C]0.887442993530035[/C][/ROW]
[ROW][C]T-STAT[/C][C]1.62001271213028[/C][/ROW]
[ROW][C]p-value[/C][C]0.16615512459958[/C][/ROW]
[ROW][C]Lambda[/C][C]-0.437668930809604[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=122041&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=122041&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-7.11178585171747
beta1.4376689308096
S.D.0.887442993530035
T-STAT1.62001271213028
p-value0.16615512459958
Lambda-0.437668930809604


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