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, 14 Jul 2015 13:54:19 +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/Jul/14/t14368800010z92uoig5k4irg8.htm/, Retrieved Wed, 15 May 2024 18:29:31 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=279565, Retrieved Wed, 15 May 2024 18:29:31 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Mean Plot] [meanplot consumpt...] [2015-02-25 21:22:56] [6b666be92a2717ec998be34b30849c21]
- R PD  [Mean Plot] [REEKS A: oefening 7] [2015-07-10 22:53:58] [6b666be92a2717ec998be34b30849c21]
- RM D    [(Partial) Autocorrelation Function] [REEKS B: oefening 9] [2015-07-13 19:34:17] [6b666be92a2717ec998be34b30849c21]
- RM        [Variability] [REEKS B: oefening 10] [2015-07-14 12:47:52] [6b666be92a2717ec998be34b30849c21]
- RMP           [Standard Deviation-Mean Plot] [REEKS B: oefening 11] [2015-07-14 12:54:19] [b81b5adcb18a6dd731e9cb79a54989dd] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
1263600
1216800
1287000
1029600
1333800
1310400
1404000
1450800
1614600
1404000
1333800
1661400
1404000
1053000
1240200
936000
1310400
1076400
1427400
1287000
1357200
1521000
1497600
1778400
1287000
1076400
1193400
865800
1240200
959400
1357200
1287000
1146600
1638000
1474200
1684800
1263600
1170000
1053000
865800
1146600
1029600
1404000
1357200
1170000
1567800
1450800
1872000
1497600
912600
912600
912600
1076400
1076400
1450800
1333800
1193400
1497600
1380600
1989000
1567800
912600
959400
795600
1099800
1263600
1591200
1567800
1263600
1474200
1310400
1872000
1427400
1146600
1029600
772200
1146600
1380600
1614600
1521000
1123200
1614600
1263600
1942200
1614600
1170000
1076400
725400
1146600
1099800
1661400
1661400
1263600
1638000
1216800
1895400
1614600
1193400
912600
631800
1240200
1193400
1567800
1801800
1333800
1497600
1123200
1942200

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 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 & 1 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ fisher.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=279565&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 Ronald Aylmer Fisher' @ fisher.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=279565&T=0

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






Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
11359150169463.455432942631800
21324050231101.365009147842400
31267500248010.56429112819000
41279200271669.7460319991006200
51269450321019.1285719451076400
61306500323137.7921680981076400
71331850316385.4855307771170000
81347450339898.5316505231170000
91337700370770.4754353361310400

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 1359150 & 169463.455432942 & 631800 \tabularnewline
2 & 1324050 & 231101.365009147 & 842400 \tabularnewline
3 & 1267500 & 248010.56429112 & 819000 \tabularnewline
4 & 1279200 & 271669.746031999 & 1006200 \tabularnewline
5 & 1269450 & 321019.128571945 & 1076400 \tabularnewline
6 & 1306500 & 323137.792168098 & 1076400 \tabularnewline
7 & 1331850 & 316385.485530777 & 1170000 \tabularnewline
8 & 1347450 & 339898.531650523 & 1170000 \tabularnewline
9 & 1337700 & 370770.475435336 & 1310400 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=279565&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]1359150[/C][C]169463.455432942[/C][C]631800[/C][/ROW]
[ROW][C]2[/C][C]1324050[/C][C]231101.365009147[/C][C]842400[/C][/ROW]
[ROW][C]3[/C][C]1267500[/C][C]248010.56429112[/C][C]819000[/C][/ROW]
[ROW][C]4[/C][C]1279200[/C][C]271669.746031999[/C][C]1006200[/C][/ROW]
[ROW][C]5[/C][C]1269450[/C][C]321019.128571945[/C][C]1076400[/C][/ROW]
[ROW][C]6[/C][C]1306500[/C][C]323137.792168098[/C][C]1076400[/C][/ROW]
[ROW][C]7[/C][C]1331850[/C][C]316385.485530777[/C][C]1170000[/C][/ROW]
[ROW][C]8[/C][C]1347450[/C][C]339898.531650523[/C][C]1170000[/C][/ROW]
[ROW][C]9[/C][C]1337700[/C][C]370770.475435336[/C][C]1310400[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=279565&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=279565&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
11359150169463.455432942631800
21324050231101.365009147842400
31267500248010.56429112819000
41279200271669.7460319991006200
51269450321019.1285719451076400
61306500323137.7921680981076400
71331850316385.4855307771170000
81347450339898.5316505231170000
91337700370770.4754353361310400






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha429089.097719235
beta-0.107448317059865
S.D.0.688757001617609
T-STAT-0.156003230177716
p-value0.88043412050898

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

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]429089.097719235[/C][/ROW]
[ROW][C]beta[/C][C]-0.107448317059865[/C][/ROW]
[ROW][C]S.D.[/C][C]0.688757001617609[/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=279565&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=279565&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)
alpha429089.097719235
beta-0.107448317059865
S.D.0.688757001617609
T-STAT-0.156003230177716
p-value0.88043412050898






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha30.8088026491764
beta-1.296335507281
S.D.3.45582772786276
T-STAT-0.375115778147515
p-value0.718681736243654
Lambda2.296335507281

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & 30.8088026491764 \tabularnewline
beta & -1.296335507281 \tabularnewline
S.D. & 3.45582772786276 \tabularnewline
T-STAT & -0.375115778147515 \tabularnewline
p-value & 0.718681736243654 \tabularnewline
Lambda & 2.296335507281 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=279565&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]30.8088026491764[/C][/ROW]
[ROW][C]beta[/C][C]-1.296335507281[/C][/ROW]
[ROW][C]S.D.[/C][C]3.45582772786276[/C][/ROW]
[ROW][C]T-STAT[/C][C]-0.375115778147515[/C][/ROW]
[ROW][C]p-value[/C][C]0.718681736243654[/C][/ROW]
[ROW][C]Lambda[/C][C]2.296335507281[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=279565&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=279565&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)
alpha30.8088026491764
beta-1.296335507281
S.D.3.45582772786276
T-STAT-0.375115778147515
p-value0.718681736243654
Lambda2.296335507281


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