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, 16 Dec 2007 02:50:38 -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/Dec/16/t1197797796ab4uwvosmy7n2dl.htm/, Retrieved Thu, 02 May 2024 11:53:19 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=4117, Retrieved Thu, 02 May 2024 11:53:19 +0000
QR Codes:

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

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] [standard deviatio...] [2007-12-16 09:50:38] [bd02e85be52eb1cb060a2c60779eb820] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
51978
50330
48538
46636
44574
43028
56924
65193
67132
64398
57591
56279
56331
55015
53405
51200
48695
48057
65964
71371
72184
69400
65605
64727
64506
62751
59938
57870
55263
55777
73944
78201
78715
75009
69705
69979
69755
68155
64211
61110
58373
58288
74822
79606
80232
75631
70996
69144
67141
65413
63391
60964
58412
57539
73377
77413
78932
74789
70076
67944
64076
63136
60198
59057
57388
56708
70019
72263
74152
67057
61941
58331

Tables (Output of Computation)





Summary of compuational 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 compuational 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=4117&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]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=4117&T=0

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






Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
154383.41666666678210.7438424889215154
260162.83333333338968.6588420032121854
366804.83333333338591.84024984130177
469193.58333333337593.1590438467433596
567949.257179.12846223534358
663693.83333333335944.4290318878831124

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 54383.4166666667 & 8210.74384248892 & 15154 \tabularnewline
2 & 60162.8333333333 & 8968.65884200321 & 21854 \tabularnewline
3 & 66804.8333333333 & 8591.840249841 & 30177 \tabularnewline
4 & 69193.5833333333 & 7593.15904384674 & 33596 \tabularnewline
5 & 67949.25 & 7179.128462235 & 34358 \tabularnewline
6 & 63693.8333333333 & 5944.42903188788 & 31124 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=4117&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]54383.4166666667[/C][C]8210.74384248892[/C][C]15154[/C][/ROW]
[ROW][C]2[/C][C]60162.8333333333[/C][C]8968.65884200321[/C][C]21854[/C][/ROW]
[ROW][C]3[/C][C]66804.8333333333[/C][C]8591.840249841[/C][C]30177[/C][/ROW]
[ROW][C]4[/C][C]69193.5833333333[/C][C]7593.15904384674[/C][C]33596[/C][/ROW]
[ROW][C]5[/C][C]67949.25[/C][C]7179.128462235[/C][C]34358[/C][/ROW]
[ROW][C]6[/C][C]63693.8333333333[/C][C]5944.42903188788[/C][C]31124[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=4117&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=4117&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
154383.41666666678210.7438424889215154
260162.83333333338968.6588420032121854
366804.83333333338591.84024984130177
469193.58333333337593.1590438467433596
567949.257179.12846223534358
663693.83333333335944.4290318878831124






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha11502.8106880387
beta-0.0589472181040068
S.D.0.0932773005060287
T-STAT-0.631956733141059
p-value0.561732881846253

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 11502.8106880387 \tabularnewline
beta & -0.0589472181040068 \tabularnewline
S.D. & 0.0932773005060287 \tabularnewline
T-STAT & -0.631956733141059 \tabularnewline
p-value & 0.561732881846253 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=4117&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]11502.8106880387[/C][/ROW]
[ROW][C]beta[/C][C]-0.0589472181040068[/C][/ROW]
[ROW][C]S.D.[/C][C]0.0932773005060287[/C][/ROW]
[ROW][C]T-STAT[/C][C]-0.631956733141059[/C][/ROW]
[ROW][C]p-value[/C][C]0.561732881846253[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=4117&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=4117&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)
alpha11502.8106880387
beta-0.0589472181040068
S.D.0.0932773005060287
T-STAT-0.631956733141059
p-value0.561732881846253






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha14.1088794243720
beta-0.466847355796923
S.D.0.785532265794973
T-STAT-0.594307040111796
p-value0.584290067045958
Lambda1.46684735579692

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & 14.1088794243720 \tabularnewline
beta & -0.466847355796923 \tabularnewline
S.D. & 0.785532265794973 \tabularnewline
T-STAT & -0.594307040111796 \tabularnewline
p-value & 0.584290067045958 \tabularnewline
Lambda & 1.46684735579692 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=4117&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]14.1088794243720[/C][/ROW]
[ROW][C]beta[/C][C]-0.466847355796923[/C][/ROW]
[ROW][C]S.D.[/C][C]0.785532265794973[/C][/ROW]
[ROW][C]T-STAT[/C][C]-0.594307040111796[/C][/ROW]
[ROW][C]p-value[/C][C]0.584290067045958[/C][/ROW]
[ROW][C]Lambda[/C][C]1.46684735579692[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=4117&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=4117&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)
alpha14.1088794243720
beta-0.466847355796923
S.D.0.785532265794973
T-STAT-0.594307040111796
p-value0.584290067045958
Lambda1.46684735579692


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