Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*Unverified author*

R Software Module

rwasp_smp.wasp

Title produced by software

Standard Deviation-Mean Plot

Date of computation

Fri, 14 Dec 2007 02:54:56 -0700

Cite this page as follows

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2007/Dec/14/t1197625168sls7rfwqw79h3ea.htm/, Retrieved Fri, 03 May 2024 00:10:18 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=3780, Retrieved Fri, 03 May 2024 00:10:18 +0000

QR Codes:

Paste this QR Code to cite your computation.

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords

s0650550 s0650062

Estimated Impact

222

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)

-     [Standard Deviation-Mean Plot] [Inducing Stationa...] [2007-11-29 14:28:06] [68a1fecd8f1c75119cd425050381cede]
-    D    [Standard Deviation-Mean Plot] [Time seriues] [2007-12-14 09:54:56] [d41d8cd98f00b204e9800998ecf8427e] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

263418000000
262752000000
266433000000
267722000000
266003000000
262971000000
265521000000
264676000000
270223000000
269508000000
268457000000
265814000000
266680000000
263018000000
269285000000
269829000000
270911000000
266844000000
271244000000
269907000000
271296000000
270157000000
271322000000
267179000000
264101000000
265518000000
269419000000
268714000000
272482000000
268351000000
268175000000
270674000000
272764000000
272599000000
270333000000
270846000000
270491000000
269160000000
274027000000
273784000000
276663000000
274525000000
271344000000
271115000000
270798000000
273911000000
273985000000
271917000000
273338000000
270601000000
273547000000
275363000000
281229000000
277793000000
279913000000
282500000000
280041000000
282166000000
290304000000
283519000000

Summary of compuational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 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=3780&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=3780&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=3780&T=0

As an alternative you can also use a QR Code:

The GUIDs for individual cells are displayed in the table below:

Summary of compuational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 seconds
R Server	'Gwilym Jenkins' @ 72.249.127.135

Standard Deviation-Mean Plot
Section	Mean	Standard Deviation	Range
1	266124833333.333	2478095046.72562	6.805e+09
2	268972666666.667	2547773656.75615	8.57e+09
3	2.69498e+11	2725958747.10332	6.331e+09
4	272643333333.333	2154153000.95308	8.941e+09
5	279192833333.333	5426668823.05817	2.4301e+10

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 266124833333.333 & 2478095046.72562 & 6.805e+09 \tabularnewline
2 & 268972666666.667 & 2547773656.75615 & 8.57e+09 \tabularnewline
3 & 2.69498e+11 & 2725958747.10332 & 6.331e+09 \tabularnewline
4 & 272643333333.333 & 2154153000.95308 & 8.941e+09 \tabularnewline
5 & 279192833333.333 & 5426668823.05817 & 2.4301e+10 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=3780&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]266124833333.333[/C][C]2478095046.72562[/C][C]6.805e+09[/C][/ROW]
[ROW][C]2[/C][C]268972666666.667[/C][C]2547773656.75615[/C][C]8.57e+09[/C][/ROW]
[ROW][C]3[/C][C]2.69498e+11[/C][C]2725958747.10332[/C][C]6.331e+09[/C][/ROW]
[ROW][C]4[/C][C]272643333333.333[/C][C]2154153000.95308[/C][C]8.941e+09[/C][/ROW]
[ROW][C]5[/C][C]279192833333.333[/C][C]5426668823.05817[/C][C]2.4301e+10[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=3780&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=3780&T=1

As an alternative you can also use a QR Code:

The GUIDs for individual cells are displayed in the table below:

Standard Deviation-Mean Plot
Section	Mean	Standard Deviation	Range
1	266124833333.333	2478095046.72562	6.805e+09
2	268972666666.667	2547773656.75615	8.57e+09
3	2.69498e+11	2725958747.10332	6.331e+09
4	272643333333.333	2154153000.95308	8.941e+09
5	279192833333.333	5426668823.05817	2.4301e+10

Regression: S.E.(k) = alpha + beta * Mean(k)
alpha	-57621011662.2493
beta	0.223702907446506
S.D.	0.0848968689441447
T-STAT	2.63499596897601
p-value	0.0779899249818845

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & -57621011662.2493 \tabularnewline
beta & 0.223702907446506 \tabularnewline
S.D. & 0.0848968689441447 \tabularnewline
T-STAT & 2.63499596897601 \tabularnewline
p-value & 0.0779899249818845 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=3780&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-57621011662.2493[/C][/ROW]
[ROW][C]beta[/C][C]0.223702907446506[/C][/ROW]
[ROW][C]S.D.[/C][C]0.0848968689441447[/C][/ROW]
[ROW][C]T-STAT[/C][C]2.63499596897601[/C][/ROW]
[ROW][C]p-value[/C][C]0.0779899249818845[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=3780&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=3780&T=2

As an alternative you can also use a QR Code:

The GUIDs for individual cells are displayed in the table below:

Regression: S.E.(k) = alpha + beta * Mean(k)
alpha	-57621011662.2493
beta	0.223702907446506
S.D.	0.0848968689441447
T-STAT	2.63499596897601
p-value	0.0779899249818845

Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha	-393.025292899733
beta	15.7564605247205
S.D.	6.96658284858452
T-STAT	2.26172010972666
p-value	0.108753016411034
Lambda	-14.7564605247205

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & -393.025292899733 \tabularnewline
beta & 15.7564605247205 \tabularnewline
S.D. & 6.96658284858452 \tabularnewline
T-STAT & 2.26172010972666 \tabularnewline
p-value & 0.108753016411034 \tabularnewline
Lambda & -14.7564605247205 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=3780&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-393.025292899733[/C][/ROW]
[ROW][C]beta[/C][C]15.7564605247205[/C][/ROW]
[ROW][C]S.D.[/C][C]6.96658284858452[/C][/ROW]
[ROW][C]T-STAT[/C][C]2.26172010972666[/C][/ROW]
[ROW][C]p-value[/C][C]0.108753016411034[/C][/ROW]
[ROW][C]Lambda[/C][C]-14.7564605247205[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=3780&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=3780&T=3

As an alternative you can also use a QR Code:

The GUIDs for individual cells are displayed in the table below:

Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha	-393.025292899733
beta	15.7564605247205
S.D.	6.96658284858452
T-STAT	2.26172010972666
p-value	0.108753016411034
Lambda	-14.7564605247205

Figure 1

PNG link

Postscript link

PDF link

Figure 2

PNG link

Postscript link

PDF link

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

Free Statistics

Description of Statistical Computation

Tree of Dependent Computations

Dataset

Tables (Output of Computation)

Figures (Output of Computation)

Input Parameters & R Code