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, 15 Aug 2017 21:38:14 +0200
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2017/Aug/15/t1502826042eq3k6wbk5ddgk6e.htm/, Retrieved Sun, 19 May 2024 23:28:04 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=307324, Retrieved Sun, 19 May 2024 23:28:04 +0000
QR Codes:

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

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] [] [2017-08-15 19:38:14] [f8975010d6e80ebfdd11eb899305ce74] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
1474200
1419600
1501500
1201200
1556100
1528800
1638000
1692600
1883700
1638000
1556100
1938300
1638000
1228500
1446900
1092000
1528800
1255800
1665300
1501500
1583400
1774500
1747200
2074800
1501500
1255800
1392300
1010100
1446900
1119300
1583400
1501500
1337700
1911000
1719900
1965600
1474200
1365000
1228500
1010100
1337700
1201200
1638000
1583400
1365000
1829100
1692600
2184000
1747200
1064700
1064700
1064700
1255800
1255800
1692600
1556100
1392300
1747200
1610700
2320500
1829100
1064700
1119300
928200
1283100
1474200
1856400
1829100
1474200
1719900
1528800
2184000
1665300
1337700
1201200
900900
1337700
1610700
1883700
1774500
1310400
1883700
1474200
2265900
1883700
1365000
1255800
846300
1337700
1283100
1938300
1938300
1474200
1911000
1419600
2211300
1883700
1392300
1064700
737100
1446900
1392300
1829100
2102100
1556100
1747200
1310400
2265900

Tables (Output of Computation)





Summary of computational transaction
Raw Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R ServerBig Analytics Cloud Computing Center

\begin{tabular}{lllllllll}
\hline
Summary of computational transaction \tabularnewline
Raw Input view raw input (R code)  \tabularnewline
Raw Outputview raw output of R engine  \tabularnewline
Computing time1 seconds \tabularnewline
R ServerBig Analytics Cloud Computing Center \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=307324&T=0

[TABLE]
[ROW]
Summary of computational transaction[/C][/ROW] [ROW]Raw Input[/C] view raw input (R code) [/C][/ROW] [ROW]Raw Output[/C]view raw output of R engine [/C][/ROW] [ROW]Computing time[/C]1 seconds[/C][/ROW] [ROW]R Server[/C]Big Analytics Cloud Computing Center[/C][/ROW] [/TABLE] Source: https://freestatistics.org/blog/index.php?pk=307324&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=307324&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 Input view raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R ServerBig Analytics Cloud Computing Center






Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
11585675197707.364671766737100
21544725269618.259177338982800
31478750289345.658339641955500
41492400316948.0370373321173900
51481025374522.3166672691255800
61524250376994.0908627811255800
71553825369116.3997859061365000
81572025396548.2869256111365000
91560650432565.5546745591528800

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 1585675 & 197707.364671766 & 737100 \tabularnewline
2 & 1544725 & 269618.259177338 & 982800 \tabularnewline
3 & 1478750 & 289345.658339641 & 955500 \tabularnewline
4 & 1492400 & 316948.037037332 & 1173900 \tabularnewline
5 & 1481025 & 374522.316667269 & 1255800 \tabularnewline
6 & 1524250 & 376994.090862781 & 1255800 \tabularnewline
7 & 1553825 & 369116.399785906 & 1365000 \tabularnewline
8 & 1572025 & 396548.286925611 & 1365000 \tabularnewline
9 & 1560650 & 432565.554674559 & 1528800 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=307324&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]1585675[/C][C]197707.364671766[/C][C]737100[/C][/ROW]
[ROW][C]2[/C][C]1544725[/C][C]269618.259177338[/C][C]982800[/C][/ROW]
[ROW][C]3[/C][C]1478750[/C][C]289345.658339641[/C][C]955500[/C][/ROW]
[ROW][C]4[/C][C]1492400[/C][C]316948.037037332[/C][C]1173900[/C][/ROW]
[ROW][C]5[/C][C]1481025[/C][C]374522.316667269[/C][C]1255800[/C][/ROW]
[ROW][C]6[/C][C]1524250[/C][C]376994.090862781[/C][C]1255800[/C][/ROW]
[ROW][C]7[/C][C]1553825[/C][C]369116.399785906[/C][C]1365000[/C][/ROW]
[ROW][C]8[/C][C]1572025[/C][C]396548.286925611[/C][C]1365000[/C][/ROW]
[ROW][C]9[/C][C]1560650[/C][C]432565.554674559[/C][C]1528800[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=307324&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=307324&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
11585675197707.364671766737100
21544725269618.259177338982800
31478750289345.658339641955500
41492400316948.0370373321173900
51481025374522.3166672691255800
61524250376994.0908627811255800
71553825369116.3997859061365000
81572025396548.2869256111365000
91560650432565.5546745591528800






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha500603.947339108
beta-0.107448317059866
S.D.0.68875700161761
T-STAT-0.156003230177716
p-value0.88043412050898

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 500603.947339108 \tabularnewline
beta & -0.107448317059866 \tabularnewline
S.D. & 0.68875700161761 \tabularnewline
T-STAT & -0.156003230177716 \tabularnewline
p-value & 0.88043412050898 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=307324&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]500603.947339108[/C][/ROW]
[ROW][C]beta[/C][C]-0.107448317059866[/C][/ROW]
[ROW][C]S.D.[/C][C]0.68875700161761[/C][/ROW]
[ROW][C]T-STAT[/C][C]-0.156003230177716[/C][/ROW]
[ROW][C]p-value[/C][C]0.88043412050898[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=307324&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=307324&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)
alpha500603.947339108
beta-0.107448317059866
S.D.0.68875700161761
T-STAT-0.156003230177716
p-value0.88043412050898






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha31.1627843287377
beta-1.29633550728117
S.D.3.45582772786278
T-STAT-0.375115778147562
p-value0.718681736243621
Lambda2.29633550728117

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & 31.1627843287377 \tabularnewline
beta & -1.29633550728117 \tabularnewline
S.D. & 3.45582772786278 \tabularnewline
T-STAT & -0.375115778147562 \tabularnewline
p-value & 0.718681736243621 \tabularnewline
Lambda & 2.29633550728117 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=307324&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]31.1627843287377[/C][/ROW]
[ROW][C]beta[/C][C]-1.29633550728117[/C][/ROW]
[ROW][C]S.D.[/C][C]3.45582772786278[/C][/ROW]
[ROW][C]T-STAT[/C][C]-0.375115778147562[/C][/ROW]
[ROW][C]p-value[/C][C]0.718681736243621[/C][/ROW]
[ROW][C]Lambda[/C][C]2.29633550728117[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=307324&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=307324&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)
alpha31.1627843287377
beta-1.29633550728117
S.D.3.45582772786278
T-STAT-0.375115778147562
p-value0.718681736243621
Lambda2.29633550728117


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