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, 27 Nov 2007 02:37:20 -0700
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2007/Nov/27/t1196155741hjaadcs7pa59l35.htm/, Retrieved Sun, 05 May 2024 20:18:09 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=6754, Retrieved Sun, 05 May 2024 20:18:09 +0000
QR Codes:

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

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] [statistiek opdrac...] [2007-11-27 09:37:20] [0d8f90662bd29255b957138f6acc6b00] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
3804
3491
4151
4254
4717
4866
4001
3758
4780
5016
4296
4467
3891
3872
3867
3973
4640
4538
3836
3770
4374
4497
3945
3862
3608
3301
3882
3605
4305
4216
3971
3988
4317
4484
4247
3520
3686
3403
3990
4053
4548
4559
3922
4209
4517
4386
3221
3127
3777
3322
3899
4033
4463
4819
4246
4255
4760
4581
4309
4016
3601
3257
3823
3940
4534
4575
3953
4206
4649
4353
3835
3944

Tables (Output of Computation)





Summary of compuational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135

\begin{tabular}{lllllllll}
\hline
Summary of compuational 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 & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=6754&T=0

[TABLE]
[ROW][C]Summary of compuational 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]'Gwilym Jenkins' @ 72.249.127.135[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=6754&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=6754&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 compuational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Gwilym Jenkins' @ 72.249.127.135






Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
14300.08333333333483.6588689511971415
24088.75321.9201270219341383
33953.66666666667375.774674586521661
43968.41666666667512.901631244563954
54206.66666666667429.559887524076514
64055.83333333333420.537282604772433

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 4300.08333333333 & 483.658868951197 & 1415 \tabularnewline
2 & 4088.75 & 321.920127021934 & 1383 \tabularnewline
3 & 3953.66666666667 & 375.774674586521 & 661 \tabularnewline
4 & 3968.41666666667 & 512.901631244563 & 954 \tabularnewline
5 & 4206.66666666667 & 429.559887524076 & 514 \tabularnewline
6 & 4055.83333333333 & 420.537282604772 & 433 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=6754&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]4300.08333333333[/C][C]483.658868951197[/C][C]1415[/C][/ROW]
[ROW][C]2[/C][C]4088.75[/C][C]321.920127021934[/C][C]1383[/C][/ROW]
[ROW][C]3[/C][C]3953.66666666667[/C][C]375.774674586521[/C][C]661[/C][/ROW]
[ROW][C]4[/C][C]3968.41666666667[/C][C]512.901631244563[/C][C]954[/C][/ROW]
[ROW][C]5[/C][C]4206.66666666667[/C][C]429.559887524076[/C][C]514[/C][/ROW]
[ROW][C]6[/C][C]4055.83333333333[/C][C]420.537282604772[/C][C]433[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=6754&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=6754&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
14300.08333333333483.6588689511971415
24088.75321.9201270219341383
33953.66666666667375.774674586521661
43968.41666666667512.901631244563954
54206.66666666667429.559887524076514
64055.83333333333420.537282604772433






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha15.3078048648578
beta0.0998032009960866
S.D.0.251675368631013
T-STAT0.396555298752379
p-value0.711941426510526

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 15.3078048648578 \tabularnewline
beta & 0.0998032009960866 \tabularnewline
S.D. & 0.251675368631013 \tabularnewline
T-STAT & 0.396555298752379 \tabularnewline
p-value & 0.711941426510526 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=6754&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]15.3078048648578[/C][/ROW]
[ROW][C]beta[/C][C]0.0998032009960866[/C][/ROW]
[ROW][C]S.D.[/C][C]0.251675368631013[/C][/ROW]
[ROW][C]T-STAT[/C][C]0.396555298752379[/C][/ROW]
[ROW][C]p-value[/C][C]0.711941426510526[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=6754&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=6754&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)
alpha15.3078048648578
beta0.0998032009960866
S.D.0.251675368631013
T-STAT0.396555298752379
p-value0.711941426510526






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha-2.43213656892838
beta1.01840818379384
S.D.2.51995998698892
T-STAT0.404136648618268
p-value0.706787345092888
Lambda-0.0184081837938355

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & -2.43213656892838 \tabularnewline
beta & 1.01840818379384 \tabularnewline
S.D. & 2.51995998698892 \tabularnewline
T-STAT & 0.404136648618268 \tabularnewline
p-value & 0.706787345092888 \tabularnewline
Lambda & -0.0184081837938355 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=6754&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-2.43213656892838[/C][/ROW]
[ROW][C]beta[/C][C]1.01840818379384[/C][/ROW]
[ROW][C]S.D.[/C][C]2.51995998698892[/C][/ROW]
[ROW][C]T-STAT[/C][C]0.404136648618268[/C][/ROW]
[ROW][C]p-value[/C][C]0.706787345092888[/C][/ROW]
[ROW][C]Lambda[/C][C]-0.0184081837938355[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=6754&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=6754&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-2.43213656892838
beta1.01840818379384
S.D.2.51995998698892
T-STAT0.404136648618268
p-value0.706787345092888
Lambda-0.0184081837938355


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