Repository of Reproducible Computations

Free Statistics

of Irreproducible Research!

Author's title

Author

*The author of this computation has been verified*

R Software Module

rwasp_smp.wasp

Title produced by software

Standard Deviation-Mean Plot

Date of computation

Fri, 27 Nov 2009 08:52:55 -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/2009/Nov/27/t12593372119nyydz01w5zlbhr.htm/, Retrieved Sun, 28 Apr 2024 20:37:01 +0000

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=60916, Retrieved Sun, 28 Apr 2024 20:37:01 +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

Estimated Impact

130

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

-     [Univariate Explorative Data Analysis] [Run Sequence gebo...] [2008-12-12 13:32:37] [76963dc1903f0f612b6153510a3818cf]
- R  D  [Univariate Explorative Data Analysis] [Run Sequence gebo...] [2008-12-17 12:14:40] [76963dc1903f0f612b6153510a3818cf]
-         [Univariate Explorative Data Analysis] [Run Sequence Plot...] [2008-12-22 18:19:51] [1ce0d16c8f4225c977b42c8fa93bc163]
- RMP       [Standard Deviation-Mean Plot] [Identifying Integ...] [2009-11-22 12:50:05] [b98453cac15ba1066b407e146608df68]
- R PD        [Standard Deviation-Mean Plot] [workshop 8] [2009-11-26 19:55:09] [3d8acb8ffdb376c5fec19e610f8198c2]
-    D          [Standard Deviation-Mean Plot] [workshop 8] [2009-11-26 20:14:23] [3d8acb8ffdb376c5fec19e610f8198c2]
-   PD              [Standard Deviation-Mean Plot] [ws 8: SPM transfo...] [2009-11-27 15:52:55] [ac86848d66148c9c4c9404e0c9a511eb] [Current]

Feedback Forum

Post a new message

Dataseries X:

Download CSV

Histogram

Boxplots

29.93995649
30.98309281
39.38169911
39.19321992
36.77331228
38.38246202
35.39357642
37.4029774
43.43958773
38.70735233
36.47106861
43.72703685
32.0161951
32.15891705
40.45379715
36.60854129
38.62620371
39.11236311
36.36098279
39.05843178
45.1820972
37.0748412
43.12537983
45.5437433
33.15293223
33.01146129
41.72960782
41.96756455
42.23150544
41.0931305
40.58723059
40.82709272
47.00773541
42.04679696
43.62257313
49.26778171
38.517929
35.44906073
43.85751539
45.33719023
43.64869582
43.20399353
41.51776177
45.23381191
47.28874003
48.23043482
45.33719023
50.52438459
41.96756455
36.91045893
46.26456675
45.41467937
38.97749448
40.69388631
34.11038175
36.30590383
38.78845191
35.44906073
32.24446257
39.97237997

Summary of computational 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 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=60916&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=60916&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=60916&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 computational 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	37.4829451641667	4.13799435118977	13.78708036
2	38.7767911258333	4.38688908062706	13.5275482
3	41.3787843625	4.6640630768072	16.25632042
4	44.0122256708333	4.11744792920394	15.07532386
5	38.9249409291667	4.26721183783709	14.02010418

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 37.4829451641667 & 4.13799435118977 & 13.78708036 \tabularnewline
2 & 38.7767911258333 & 4.38688908062706 & 13.5275482 \tabularnewline
3 & 41.3787843625 & 4.6640630768072 & 16.25632042 \tabularnewline
4 & 44.0122256708333 & 4.11744792920394 & 15.07532386 \tabularnewline
5 & 38.9249409291667 & 4.26721183783709 & 14.02010418 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=60916&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]37.4829451641667[/C][C]4.13799435118977[/C][C]13.78708036[/C][/ROW]
[ROW][C]2[/C][C]38.7767911258333[/C][C]4.38688908062706[/C][C]13.5275482[/C][/ROW]
[ROW][C]3[/C][C]41.3787843625[/C][C]4.6640630768072[/C][C]16.25632042[/C][/ROW]
[ROW][C]4[/C][C]44.0122256708333[/C][C]4.11744792920394[/C][C]15.07532386[/C][/ROW]
[ROW][C]5[/C][C]38.9249409291667[/C][C]4.26721183783709[/C][C]14.02010418[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=60916&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=60916&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	37.4829451641667	4.13799435118977	13.78708036
2	38.7767911258333	4.38688908062706	13.5275482
3	41.3787843625	4.6640630768072	16.25632042
4	44.0122256708333	4.11744792920394	15.07532386
5	38.9249409291667	4.26721183783709	14.02010418

Regression: S.E.(k) = alpha + beta * Mean(k)
alpha	4.16899638249528
beta	0.00363266542006865
S.D.	0.0496766325290915
T-STAT	0.0731262413558587
p-value	0.946308296663852

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 4.16899638249528 \tabularnewline
beta & 0.00363266542006865 \tabularnewline
S.D. & 0.0496766325290915 \tabularnewline
T-STAT & 0.0731262413558587 \tabularnewline
p-value & 0.946308296663852 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=60916&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]4.16899638249528[/C][/ROW]
[ROW][C]beta[/C][C]0.00363266542006865[/C][/ROW]
[ROW][C]S.D.[/C][C]0.0496766325290915[/C][/ROW]
[ROW][C]T-STAT[/C][C]0.0731262413558587[/C][/ROW]
[ROW][C]p-value[/C][C]0.946308296663852[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=60916&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=60916&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	4.16899638249528
beta	0.00363266542006865
S.D.	0.0496766325290915
T-STAT	0.0731262413558587
p-value	0.946308296663852

Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha	1.30707388423646
beta	0.0417083547743476
S.D.	0.461676618022885
T-STAT	0.090341059404226
p-value	0.93370996990255
Lambda	0.958291645225652

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & 1.30707388423646 \tabularnewline
beta & 0.0417083547743476 \tabularnewline
S.D. & 0.461676618022885 \tabularnewline
T-STAT & 0.090341059404226 \tabularnewline
p-value & 0.93370996990255 \tabularnewline
Lambda & 0.958291645225652 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=60916&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]1.30707388423646[/C][/ROW]
[ROW][C]beta[/C][C]0.0417083547743476[/C][/ROW]
[ROW][C]S.D.[/C][C]0.461676618022885[/C][/ROW]
[ROW][C]T-STAT[/C][C]0.090341059404226[/C][/ROW]
[ROW][C]p-value[/C][C]0.93370996990255[/C][/ROW]
[ROW][C]Lambda[/C][C]0.958291645225652[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=60916&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=60916&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	1.30707388423646
beta	0.0417083547743476
S.D.	0.461676618022885
T-STAT	0.090341059404226
p-value	0.93370996990255
Lambda	0.958291645225652

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