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 computationTue, 31 May 2016 11:25:48 +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/2016/May/31/t146469036085mtaj632h1h9ff.htm/, Retrieved Tue, 07 May 2024 00:38:08 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=295743, Retrieved Tue, 07 May 2024 00:38:08 +0000
QR Codes:

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

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] [] [2016-05-31 10:25:48] [9b4dafad127b39cd929ee42874de7246] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
113149
112534
108783
106640
102617
102191
117359
116083
108666
105017
100918
103907
105732
103409
100255
97036
94055
92523
106380
104846
101411
98072
95678
99148
106813
106782
103496
100854
99592
98923
110497
114783
113551
112376
111683
113467
117277
117442
115640
114872
111628
111098
124301
125847
125323
122394
121164
123963
130549
128563
125418
121982
117708
116905
128862
129655
128649
126084
123725
123974

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=295743&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'Sir Maurice George Kendall' @ kendall.wessa.net






Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
1108155.3333333335568.8962727788116441
299878.754609.2926912527913857
3107734.755816.4780053193415860
4119245.755270.7183247723414749
5125172.8333333334521.3219875553613644

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 108155.333333333 & 5568.89627277881 & 16441 \tabularnewline
2 & 99878.75 & 4609.29269125279 & 13857 \tabularnewline
3 & 107734.75 & 5816.47800531934 & 15860 \tabularnewline
4 & 119245.75 & 5270.71832477234 & 14749 \tabularnewline
5 & 125172.833333333 & 4521.32198755536 & 13644 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=295743&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]108155.333333333[/C][C]5568.89627277881[/C][C]16441[/C][/ROW]
[ROW][C]2[/C][C]99878.75[/C][C]4609.29269125279[/C][C]13857[/C][/ROW]
[ROW][C]3[/C][C]107734.75[/C][C]5816.47800531934[/C][C]15860[/C][/ROW]
[ROW][C]4[/C][C]119245.75[/C][C]5270.71832477234[/C][C]14749[/C][/ROW]
[ROW][C]5[/C][C]125172.833333333[/C][C]4521.32198755536[/C][C]13644[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=295743&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=295743&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
1108155.3333333335568.8962727788116441
299878.754609.2926912527913857
3107734.755816.4780053193415860
4119245.755270.7183247723414749
5125172.8333333334521.3219875553613644






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha6622.22173185903
beta-0.0130749123591522
S.D.0.032066709439728
T-STAT-0.407741005784442
p-value0.710816456346865

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 6622.22173185903 \tabularnewline
beta & -0.0130749123591522 \tabularnewline
S.D. & 0.032066709439728 \tabularnewline
T-STAT & -0.407741005784442 \tabularnewline
p-value & 0.710816456346865 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=295743&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]6622.22173185903[/C][/ROW]
[ROW][C]beta[/C][C]-0.0130749123591522[/C][/ROW]
[ROW][C]S.D.[/C][C]0.032066709439728[/C][/ROW]
[ROW][C]T-STAT[/C][C]-0.407741005784442[/C][/ROW]
[ROW][C]p-value[/C][C]0.710816456346865[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=295743&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=295743&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)
alpha6622.22173185903
beta-0.0130749123591522
S.D.0.032066709439728
T-STAT-0.407741005784442
p-value0.710816456346865






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha11.3776001694507
beta-0.243858462754407
S.D.0.711438740124372
T-STAT-0.342768040311913
p-value0.754384554259062
Lambda1.24385846275441

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & 11.3776001694507 \tabularnewline
beta & -0.243858462754407 \tabularnewline
S.D. & 0.711438740124372 \tabularnewline
T-STAT & -0.342768040311913 \tabularnewline
p-value & 0.754384554259062 \tabularnewline
Lambda & 1.24385846275441 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=295743&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]11.3776001694507[/C][/ROW]
[ROW][C]beta[/C][C]-0.243858462754407[/C][/ROW]
[ROW][C]S.D.[/C][C]0.711438740124372[/C][/ROW]
[ROW][C]T-STAT[/C][C]-0.342768040311913[/C][/ROW]
[ROW][C]p-value[/C][C]0.754384554259062[/C][/ROW]
[ROW][C]Lambda[/C][C]1.24385846275441[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=295743&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=295743&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)
alpha11.3776001694507
beta-0.243858462754407
S.D.0.711438740124372
T-STAT-0.342768040311913
p-value0.754384554259062
Lambda1.24385846275441


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