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 computationTue, 07 Dec 2010 00:10:56 +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/07/t1291680932bceb3ctnb0qn87h.htm/, Retrieved Sat, 04 May 2024 03:19:16 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=105991, Retrieved Sat, 04 May 2024 03:19:16 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Univariate Data Series] [data set] [2008-12-01 19:54:57] [b98453cac15ba1066b407e146608df68]
- RMP   [Standard Deviation-Mean Plot] [Unemployment] [2010-11-29 10:34:47] [b98453cac15ba1066b407e146608df68]
-   PD    [Standard Deviation-Mean Plot] [SMP 1] [2010-12-07 00:05:33] [b8e188bcc949964bed729335b3416734]
-    D        [Standard Deviation-Mean Plot] [SMP 2] [2010-12-07 00:10:56] [278a0539dc236556c5f30b5bc56ff9eb] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
431
465
511
540
552
512
413
542
544
491
458
529
525
483
528
502
563
537
465
528
505
493
456
488
488
468
542
499
477
534
528
598
474
537
376
447
545
425
458
510
472
541
507
472
540
496
432
452
420
435
509
441
416
490
396
463
403
448
398
387
426
428
510
437
453
451
434

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'Gwilym Jenkins' @ 72.249.127.135

\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 & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=105991&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]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=105991&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=105991&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'Gwilym Jenkins' @ 72.249.127.135






Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
149947.1920640477921139
2506.08333333333331.2743238669253107
3497.33333333333356.5288877611717222
4487.542.0832508250016120
5433.83333333333338.4680111103743122

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 499 & 47.1920640477921 & 139 \tabularnewline
2 & 506.083333333333 & 31.2743238669253 & 107 \tabularnewline
3 & 497.333333333333 & 56.5288877611717 & 222 \tabularnewline
4 & 487.5 & 42.0832508250016 & 120 \tabularnewline
5 & 433.833333333333 & 38.4680111103743 & 122 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=105991&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]499[/C][C]47.1920640477921[/C][C]139[/C][/ROW]
[ROW][C]2[/C][C]506.083333333333[/C][C]31.2743238669253[/C][C]107[/C][/ROW]
[ROW][C]3[/C][C]497.333333333333[/C][C]56.5288877611717[/C][C]222[/C][/ROW]
[ROW][C]4[/C][C]487.5[/C][C]42.0832508250016[/C][C]120[/C][/ROW]
[ROW][C]5[/C][C]433.833333333333[/C][C]38.4680111103743[/C][C]122[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=105991&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=105991&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
149947.1920640477921139
2506.08333333333331.2743238669253107
3497.33333333333356.5288877611717222
4487.542.0832508250016120
5433.83333333333338.4680111103743122






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha13.5893886395686
beta0.060897202439782
S.D.0.183928262952363
T-STAT0.331092141372282
p-value0.762353043788209

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 13.5893886395686 \tabularnewline
beta & 0.060897202439782 \tabularnewline
S.D. & 0.183928262952363 \tabularnewline
T-STAT & 0.331092141372282 \tabularnewline
p-value & 0.762353043788209 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=105991&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]13.5893886395686[/C][/ROW]
[ROW][C]beta[/C][C]0.060897202439782[/C][/ROW]
[ROW][C]S.D.[/C][C]0.183928262952363[/C][/ROW]
[ROW][C]T-STAT[/C][C]0.331092141372282[/C][/ROW]
[ROW][C]p-value[/C][C]0.762353043788209[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=105991&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=105991&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)
alpha13.5893886395686
beta0.060897202439782
S.D.0.183928262952363
T-STAT0.331092141372282
p-value0.762353043788209






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha0.551798624090398
beta0.516402162834365
S.D.2.02049483021219
T-STAT0.255582026300029
p-value0.814795233554622
Lambda0.483597837165635

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & 0.551798624090398 \tabularnewline
beta & 0.516402162834365 \tabularnewline
S.D. & 2.02049483021219 \tabularnewline
T-STAT & 0.255582026300029 \tabularnewline
p-value & 0.814795233554622 \tabularnewline
Lambda & 0.483597837165635 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=105991&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]0.551798624090398[/C][/ROW]
[ROW][C]beta[/C][C]0.516402162834365[/C][/ROW]
[ROW][C]S.D.[/C][C]2.02049483021219[/C][/ROW]
[ROW][C]T-STAT[/C][C]0.255582026300029[/C][/ROW]
[ROW][C]p-value[/C][C]0.814795233554622[/C][/ROW]
[ROW][C]Lambda[/C][C]0.483597837165635[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=105991&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=105991&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)
alpha0.551798624090398
beta0.516402162834365
S.D.2.02049483021219
T-STAT0.255582026300029
p-value0.814795233554622
Lambda0.483597837165635


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