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, 12 Mar 2015 20:24:01 +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/2015/Mar/12/t1426191872ayijkqh870r3pqw.htm/, Retrieved Sun, 19 May 2024 13:22:16 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=278346, Retrieved Sun, 19 May 2024 13:22:16 +0000
QR Codes:

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

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] [] [2015-03-12 20:24:01] [bad5dfd772bc354c7f8aa9414b1d4071] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
1329
1385
1681
1591
1598
1557
1190
932
1664
1717
1567
1355
1430
1863
1868
1711
1873
2095
1379
1021
1999
2094
2026
1390
1744
2117
1823
1963
1816
1966
1309
1250
2184
2295
1870
1222
1640
2194
2179
1976
1850
2077
1658
1156
2400
2218
1802
1444
1804
1541
2206
1972
1815
1749
1492
1307
1916
2035
1855
1086
1951
1733
1868
1532
1894
1586
1247
1212
2119
1931
1649
1296
1625
1454
1562
1612
1648
1412
1219
1207
1614
1537
1497
1141
1135
1368
1203
1201
1190
1347
607
914
1606
1518
1120
910

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=278346&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=278346&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=278346&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
11463.83333333333233.333095237974785
21729.08333333333345.7157156691121074
31796.58333333333360.9886571709151073
41882.83333333333362.9061569991321244
51731.5320.1492549648451120
61668.16666666667301.074087328362907
71460.66666666667179.209543847588507
81176.58333333333275.927018787522999

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 1463.83333333333 & 233.333095237974 & 785 \tabularnewline
2 & 1729.08333333333 & 345.715715669112 & 1074 \tabularnewline
3 & 1796.58333333333 & 360.988657170915 & 1073 \tabularnewline
4 & 1882.83333333333 & 362.906156999132 & 1244 \tabularnewline
5 & 1731.5 & 320.149254964845 & 1120 \tabularnewline
6 & 1668.16666666667 & 301.074087328362 & 907 \tabularnewline
7 & 1460.66666666667 & 179.209543847588 & 507 \tabularnewline
8 & 1176.58333333333 & 275.927018787522 & 999 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=278346&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]1463.83333333333[/C][C]233.333095237974[/C][C]785[/C][/ROW]
[ROW][C]2[/C][C]1729.08333333333[/C][C]345.715715669112[/C][C]1074[/C][/ROW]
[ROW][C]3[/C][C]1796.58333333333[/C][C]360.988657170915[/C][C]1073[/C][/ROW]
[ROW][C]4[/C][C]1882.83333333333[/C][C]362.906156999132[/C][C]1244[/C][/ROW]
[ROW][C]5[/C][C]1731.5[/C][C]320.149254964845[/C][C]1120[/C][/ROW]
[ROW][C]6[/C][C]1668.16666666667[/C][C]301.074087328362[/C][C]907[/C][/ROW]
[ROW][C]7[/C][C]1460.66666666667[/C][C]179.209543847588[/C][C]507[/C][/ROW]
[ROW][C]8[/C][C]1176.58333333333[/C][C]275.927018787522[/C][C]999[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=278346&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=278346&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
11463.83333333333233.333095237974785
21729.08333333333345.7157156691121074
31796.58333333333360.9886571709151073
41882.83333333333362.9061569991321244
51731.5320.1492549648451120
61668.16666666667301.074087328362907
71460.66666666667179.209543847588507
81176.58333333333275.927018787522999






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha-26.1471156725249
beta0.200513620495819
S.D.0.0811529781964679
T-STAT2.47081037507193
p-value0.0484044975548931

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & -26.1471156725249 \tabularnewline
beta & 0.200513620495819 \tabularnewline
S.D. & 0.0811529781964679 \tabularnewline
T-STAT & 2.47081037507193 \tabularnewline
p-value & 0.0484044975548931 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=278346&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-26.1471156725249[/C][/ROW]
[ROW][C]beta[/C][C]0.200513620495819[/C][/ROW]
[ROW][C]S.D.[/C][C]0.0811529781964679[/C][/ROW]
[ROW][C]T-STAT[/C][C]2.47081037507193[/C][/ROW]
[ROW][C]p-value[/C][C]0.0484044975548931[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=278346&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=278346&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-26.1471156725249
beta0.200513620495819
S.D.0.0811529781964679
T-STAT2.47081037507193
p-value0.0484044975548931






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha-1.65063369072144
beta0.992533708456184
S.D.0.512643515348369
T-STAT1.93610896995684
p-value0.100985226921614
Lambda0.00746629154381551

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & -1.65063369072144 \tabularnewline
beta & 0.992533708456184 \tabularnewline
S.D. & 0.512643515348369 \tabularnewline
T-STAT & 1.93610896995684 \tabularnewline
p-value & 0.100985226921614 \tabularnewline
Lambda & 0.00746629154381551 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=278346&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-1.65063369072144[/C][/ROW]
[ROW][C]beta[/C][C]0.992533708456184[/C][/ROW]
[ROW][C]S.D.[/C][C]0.512643515348369[/C][/ROW]
[ROW][C]T-STAT[/C][C]1.93610896995684[/C][/ROW]
[ROW][C]p-value[/C][C]0.100985226921614[/C][/ROW]
[ROW][C]Lambda[/C][C]0.00746629154381551[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=278346&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=278346&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-1.65063369072144
beta0.992533708456184
S.D.0.512643515348369
T-STAT1.93610896995684
p-value0.100985226921614
Lambda0.00746629154381551


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