FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_smp.wasp
Title produced by softwareStandard Deviation-Mean Plot
Date of computationSun, 26 Dec 2010 10:52:00 +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/2010/Dec/26/t1293360574l20iovvk3pdf201.htm/, Retrieved Mon, 06 May 2024 16:55:46 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=115507, Retrieved Mon, 06 May 2024 16:55:46 +0000
QR Codes:

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

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] [tijdreeks bevolki...] [2010-12-26 10:52:00] [531024149246456e4f6d79ace2e85c12] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
5140
4749
3635
4305
5805
4260
3869
7325
9280
6222
3272
7598
1345
1900
1480
1472
3823
4454
3357
5393
8329
4152
4042
7747
1451
911
-406
1387
2150
1577
2642
4273
8064
3243
1112
2280
505
744
-1369
-531
1041
2076
577
5080
6584
3761
294
5020
1141
3805
2127
2531
3682
3263
2798
5936
10568
5296
1870
4390
3707
5201
3748
5282
5349
6249
5517
8640
15767
8850
5582
6496
3255
6189
6452
5099
6833
7046
7739
10142
16054
7721
6182
6490
3704
6235
4655
5072
3640
5147
5703
11889
15603
9589

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'Sir Ronald Aylmer Fisher' @ 193.190.124.24

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

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






Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
154551843.224299476826008
23957.833333333332327.217098648606984
32390.333333333332148.737825427818470
41981.833333333332525.263652601757953
53950.583333333332508.043875878849427
666993264.7746128526712060
77433.53165.3529890763412799

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 5455 & 1843.22429947682 & 6008 \tabularnewline
2 & 3957.83333333333 & 2327.21709864860 & 6984 \tabularnewline
3 & 2390.33333333333 & 2148.73782542781 & 8470 \tabularnewline
4 & 1981.83333333333 & 2525.26365260175 & 7953 \tabularnewline
5 & 3950.58333333333 & 2508.04387587884 & 9427 \tabularnewline
6 & 6699 & 3264.77461285267 & 12060 \tabularnewline
7 & 7433.5 & 3165.35298907634 & 12799 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=115507&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]5455[/C][C]1843.22429947682[/C][C]6008[/C][/ROW]
[ROW][C]2[/C][C]3957.83333333333[/C][C]2327.21709864860[/C][C]6984[/C][/ROW]
[ROW][C]3[/C][C]2390.33333333333[/C][C]2148.73782542781[/C][C]8470[/C][/ROW]
[ROW][C]4[/C][C]1981.83333333333[/C][C]2525.26365260175[/C][C]7953[/C][/ROW]
[ROW][C]5[/C][C]3950.58333333333[/C][C]2508.04387587884[/C][C]9427[/C][/ROW]
[ROW][C]6[/C][C]6699[/C][C]3264.77461285267[/C][C]12060[/C][/ROW]
[ROW][C]7[/C][C]7433.5[/C][C]3165.35298907634[/C][C]12799[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=115507&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=115507&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
154551843.224299476826008
23957.833333333332327.217098648606984
32390.333333333332148.737825427818470
41981.833333333332525.263652601757953
53950.583333333332508.043875878849427
666993264.7746128526712060
77433.53165.3529890763412799






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha1875.19384152553
beta0.146110370510226
S.D.0.0904584831654561
T-STAT1.61522021370818
p-value0.167184417009976

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 1875.19384152553 \tabularnewline
beta & 0.146110370510226 \tabularnewline
S.D. & 0.0904584831654561 \tabularnewline
T-STAT & 1.61522021370818 \tabularnewline
p-value & 0.167184417009976 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=115507&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]1875.19384152553[/C][/ROW]
[ROW][C]beta[/C][C]0.146110370510226[/C][/ROW]
[ROW][C]S.D.[/C][C]0.0904584831654561[/C][/ROW]
[ROW][C]T-STAT[/C][C]1.61522021370818[/C][/ROW]
[ROW][C]p-value[/C][C]0.167184417009976[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=115507&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=115507&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)
alpha1875.19384152553
beta0.146110370510226
S.D.0.0904584831654561
T-STAT1.61522021370818
p-value0.167184417009976






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha6.37136397647843
beta0.174329771772032
S.D.0.164350256486391
T-STAT1.06072102045346
p-value0.337343181057719
Lambda0.825670228227968

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & 6.37136397647843 \tabularnewline
beta & 0.174329771772032 \tabularnewline
S.D. & 0.164350256486391 \tabularnewline
T-STAT & 1.06072102045346 \tabularnewline
p-value & 0.337343181057719 \tabularnewline
Lambda & 0.825670228227968 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=115507&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]6.37136397647843[/C][/ROW]
[ROW][C]beta[/C][C]0.174329771772032[/C][/ROW]
[ROW][C]S.D.[/C][C]0.164350256486391[/C][/ROW]
[ROW][C]T-STAT[/C][C]1.06072102045346[/C][/ROW]
[ROW][C]p-value[/C][C]0.337343181057719[/C][/ROW]
[ROW][C]Lambda[/C][C]0.825670228227968[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=115507&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=115507&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)
alpha6.37136397647843
beta0.174329771772032
S.D.0.164350256486391
T-STAT1.06072102045346
p-value0.337343181057719
Lambda0.825670228227968


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