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 computationSun, 30 May 2010 07:25:06 +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/May/30/t12752043669v00fsyj5kkd28c.htm/, Retrieved Thu, 02 May 2024 22:13:13 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=76681, Retrieved Thu, 02 May 2024 22:13:13 +0000
QR Codes:

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

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] [Werkloosheid V ...] [2010-05-30 07:25:06] [4942f64bbdc4cce21b299a740a533758] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
43657
42811
45419
50846
54500
51035
38675
36214
38763
39486
40540
40719
40471
39947
42683
47090
51520
48823
36122
33812
36928
37737
40123
41713
42025
42169
46352
50939
56139
52713
38532
37860
40880
41988
44576
46728
46913
49357
54709
60819
63695
60109
45544
43596
44431
45575
47980
49211
51374
52954
57529
62960
64530
61008
44964
43480
45429
47616
49364
51010
53188
55317
60106
65845
67028
63617
47605
45844
47925
50156
52258
53476
54327
55214
59347
64718
66208
62744
45587
43684
45676
47088
48907
50964
51798

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=76681&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
143555.41666666675778.4739013677218286
241414.08333333335352.0553988560617708
345075.08333333335716.7360277177418279
450994.91666666677027.2166668840820099
552684.83333333337242.1256652667421050
655197.08333333337307.5796696846121184
753705.33333333337971.0647553092522524

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 43555.4166666667 & 5778.47390136772 & 18286 \tabularnewline
2 & 41414.0833333333 & 5352.05539885606 & 17708 \tabularnewline
3 & 45075.0833333333 & 5716.73602771774 & 18279 \tabularnewline
4 & 50994.9166666667 & 7027.21666688408 & 20099 \tabularnewline
5 & 52684.8333333333 & 7242.12566526674 & 21050 \tabularnewline
6 & 55197.0833333333 & 7307.57966968461 & 21184 \tabularnewline
7 & 53705.3333333333 & 7971.06475530925 & 22524 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=76681&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]43555.4166666667[/C][C]5778.47390136772[/C][C]18286[/C][/ROW]
[ROW][C]2[/C][C]41414.0833333333[/C][C]5352.05539885606[/C][C]17708[/C][/ROW]
[ROW][C]3[/C][C]45075.0833333333[/C][C]5716.73602771774[/C][C]18279[/C][/ROW]
[ROW][C]4[/C][C]50994.9166666667[/C][C]7027.21666688408[/C][C]20099[/C][/ROW]
[ROW][C]5[/C][C]52684.8333333333[/C][C]7242.12566526674[/C][C]21050[/C][/ROW]
[ROW][C]6[/C][C]55197.0833333333[/C][C]7307.57966968461[/C][C]21184[/C][/ROW]
[ROW][C]7[/C][C]53705.3333333333[/C][C]7971.06475530925[/C][C]22524[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=76681&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=76681&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
143555.41666666675778.4739013677218286
241414.08333333335352.0553988560617708
345075.08333333335716.7360277177418279
450994.91666666677027.2166668840820099
552684.83333333337242.1256652667421050
655197.08333333337307.5796696846121184
753705.33333333337971.0647553092522524






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha-1897.97584926868
beta0.174186875455483
S.D.0.0235353010884478
T-STAT7.40108974178291
p-value0.000708630393647484

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & -1897.97584926868 \tabularnewline
beta & 0.174186875455483 \tabularnewline
S.D. & 0.0235353010884478 \tabularnewline
T-STAT & 7.40108974178291 \tabularnewline
p-value & 0.000708630393647484 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=76681&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-1897.97584926868[/C][/ROW]
[ROW][C]beta[/C][C]0.174186875455483[/C][/ROW]
[ROW][C]S.D.[/C][C]0.0235353010884478[/C][/ROW]
[ROW][C]T-STAT[/C][C]7.40108974178291[/C][/ROW]
[ROW][C]p-value[/C][C]0.000708630393647484[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=76681&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=76681&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-1897.97584926868
beta0.174186875455483
S.D.0.0235353010884478
T-STAT7.40108974178291
p-value0.000708630393647484






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha-5.18121879446483
beta1.29439185627456
S.D.0.155204726237380
T-STAT8.33989974180834
p-value0.000405343757807107
Lambda-0.294391856274558

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & -5.18121879446483 \tabularnewline
beta & 1.29439185627456 \tabularnewline
S.D. & 0.155204726237380 \tabularnewline
T-STAT & 8.33989974180834 \tabularnewline
p-value & 0.000405343757807107 \tabularnewline
Lambda & -0.294391856274558 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=76681&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-5.18121879446483[/C][/ROW]
[ROW][C]beta[/C][C]1.29439185627456[/C][/ROW]
[ROW][C]S.D.[/C][C]0.155204726237380[/C][/ROW]
[ROW][C]T-STAT[/C][C]8.33989974180834[/C][/ROW]
[ROW][C]p-value[/C][C]0.000405343757807107[/C][/ROW]
[ROW][C]Lambda[/C][C]-0.294391856274558[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=76681&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=76681&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-5.18121879446483
beta1.29439185627456
S.D.0.155204726237380
T-STAT8.33989974180834
p-value0.000405343757807107
Lambda-0.294391856274558


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