FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Spreidings- en gemiddeldegrafieken Eigen Reeks verbetering - Shari Van Else...

Author*Unverified author*
R Software Modulerwasp_smp.wasp
Title produced by softwareStandard Deviation-Mean Plot
Date of computationTue, 27 May 2008 13:57:52 -0600
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2008/May/27/t1211918307ag3h23giwzfdmfm.htm/, Retrieved Sun, 19 May 2024 22:06:35 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=13392, Retrieved Sun, 19 May 2024 22:06:35 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Mean Plot] [Shari Van Elsen-n...] [2008-04-27 15:50:31] [88b1f635a304a21eaec9f0fa0b832dcd]
- RM D    [Standard Deviation-Mean Plot] [Spreidings- en ge...] [2008-05-27 19:57:52] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
516.922
514.258
509.846
527.070
541.657
564.591
555.362
498.662
511.038
525.919
531.673
548.854
560.576
557.274
565.742
587.625
619.916
625.809
619.567
572.942
572.775
574.205
579.799
590.072
593.408
597.141
595.404
612.117
628.232
628.884
620.735
569.028
567.456
573.100
584.428
589.379
590.865
595.454
594.167
611.324
612.613
610.763
593.530
542.722
536.662
543.599
555.332
560.854
562.325
554.788
547.344
565.464
577.992
579.714
569.323
506.971
500.857
509.127
509.933
517.009
519.164
512.238
509.239
518.585
522.975
525.192
516.847
455.626
454.724
461.251
470.439
474.605
476.049
471.067
470.984
502.831
512.927
509.673
484.015
431.328
436.087
442.867
447.988
460.070
467.037
460.170
464.196
485.025
501.492
520.564
488.180
439.148
441.977
456.608
461.935
480.961
492.865

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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=13392&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
1528.82120.246528860747765.929
2585.52516666666723.899880307168368.535
3596.60933333333321.826099586309861.428
4578.99041666666729.057176880508475.951
5541.7372530.559179439617378.857
6495.0737528.828796074647970.468
7470.490528.109731399706381.599
8472.27441666666723.900197905070181.416

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 528.821 & 20.2465288607477 & 65.929 \tabularnewline
2 & 585.525166666667 & 23.8998803071683 & 68.535 \tabularnewline
3 & 596.609333333333 & 21.8260995863098 & 61.428 \tabularnewline
4 & 578.990416666667 & 29.0571768805084 & 75.951 \tabularnewline
5 & 541.73725 & 30.5591794396173 & 78.857 \tabularnewline
6 & 495.07375 & 28.8287960746479 & 70.468 \tabularnewline
7 & 470.4905 & 28.1097313997063 & 81.599 \tabularnewline
8 & 472.274416666667 & 23.9001979050701 & 81.416 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=13392&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]528.821[/C][C]20.2465288607477[/C][C]65.929[/C][/ROW]
[ROW][C]2[/C][C]585.525166666667[/C][C]23.8998803071683[/C][C]68.535[/C][/ROW]
[ROW][C]3[/C][C]596.609333333333[/C][C]21.8260995863098[/C][C]61.428[/C][/ROW]
[ROW][C]4[/C][C]578.990416666667[/C][C]29.0571768805084[/C][C]75.951[/C][/ROW]
[ROW][C]5[/C][C]541.73725[/C][C]30.5591794396173[/C][C]78.857[/C][/ROW]
[ROW][C]6[/C][C]495.07375[/C][C]28.8287960746479[/C][C]70.468[/C][/ROW]
[ROW][C]7[/C][C]470.4905[/C][C]28.1097313997063[/C][C]81.599[/C][/ROW]
[ROW][C]8[/C][C]472.274416666667[/C][C]23.9001979050701[/C][C]81.416[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=13392&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=13392&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
1528.82120.246528860747765.929
2585.52516666666723.899880307168368.535
3596.60933333333321.826099586309861.428
4578.99041666666729.057176880508475.951
5541.7372530.559179439617378.857
6495.0737528.828796074647970.468
7470.490528.109731399706381.599
8472.27441666666723.900197905070181.416






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha34.1413145903349
beta-0.0156230437207594
S.D.0.0299661981585141
T-STAT-0.521355549947215
p-value0.620794975669954

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 34.1413145903349 \tabularnewline
beta & -0.0156230437207594 \tabularnewline
S.D. & 0.0299661981585141 \tabularnewline
T-STAT & -0.521355549947215 \tabularnewline
p-value & 0.620794975669954 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=13392&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]34.1413145903349[/C][/ROW]
[ROW][C]beta[/C][C]-0.0156230437207594[/C][/ROW]
[ROW][C]S.D.[/C][C]0.0299661981585141[/C][/ROW]
[ROW][C]T-STAT[/C][C]-0.521355549947215[/C][/ROW]
[ROW][C]p-value[/C][C]0.620794975669954[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=13392&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=13392&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)
alpha34.1413145903349
beta-0.0156230437207594
S.D.0.0299661981585141
T-STAT-0.521355549947215
p-value0.620794975669954






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha5.3287463611373
beta-0.332719651738732
S.D.0.629784535033919
T-STAT-0.528307116529643
p-value0.616247393407135
Lambda1.33271965173873

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & 5.3287463611373 \tabularnewline
beta & -0.332719651738732 \tabularnewline
S.D. & 0.629784535033919 \tabularnewline
T-STAT & -0.528307116529643 \tabularnewline
p-value & 0.616247393407135 \tabularnewline
Lambda & 1.33271965173873 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=13392&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]5.3287463611373[/C][/ROW]
[ROW][C]beta[/C][C]-0.332719651738732[/C][/ROW]
[ROW][C]S.D.[/C][C]0.629784535033919[/C][/ROW]
[ROW][C]T-STAT[/C][C]-0.528307116529643[/C][/ROW]
[ROW][C]p-value[/C][C]0.616247393407135[/C][/ROW]
[ROW][C]Lambda[/C][C]1.33271965173873[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=13392&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=13392&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)
alpha5.3287463611373
beta-0.332719651738732
S.D.0.629784535033919
T-STAT-0.528307116529643
p-value0.616247393407135
Lambda1.33271965173873


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