Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_smp.wasp
Title produced by softwareStandard Deviation-Mean Plot
Date of computationThu, 18 Dec 2008 06:23:55 -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/2008/Dec/18/t12296074637llozetxlorwv9i.htm/, Retrieved Sun, 12 May 2024 12:08:54 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=34750, Retrieved Sun, 12 May 2024 12:08:54 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywordspaper sdmp werk 20
Estimated Impact200
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
F     [(Partial) Autocorrelation Function] [Paper: autocorrel...] [2008-12-05 10:20:24] [27f46dbe13ae2811dfd3a6f3c54d4d50]
F   PD  [(Partial) Autocorrelation Function] [Paper: autocorrel...] [2008-12-05 10:26:40] [27f46dbe13ae2811dfd3a6f3c54d4d50]
F   PD    [(Partial) Autocorrelation Function] [Paper: werklooshe...] [2008-12-05 10:43:32] [74be16979710d4c4e7c6647856088456]
F RMP       [Spectral Analysis] [paper: spectral a...] [2008-12-05 11:18:29] [74be16979710d4c4e7c6647856088456]
-   P         [Spectral Analysis] [Paper spectral an...] [2008-12-05 11:38:42] [74be16979710d4c4e7c6647856088456]
-               [Spectral Analysis] [Paper spectral an...] [2008-12-05 11:44:28] [74be16979710d4c4e7c6647856088456]
F                 [Spectral Analysis] [Paper: spectral a...] [2008-12-05 11:46:39] [74be16979710d4c4e7c6647856088456]
F   P               [Spectral Analysis] [Paper: spectral a...] [2008-12-05 11:50:13] [74be16979710d4c4e7c6647856088456]
F RMP                 [Standard Deviation-Mean Plot] [Paper: SD mean-pl...] [2008-12-05 13:05:29] [74be16979710d4c4e7c6647856088456]
-                       [Standard Deviation-Mean Plot] [paper stan dev me...] [2008-12-16 17:54:29] [5de5fb433ddcb9578e0fa830f795b7e9]
-    D                      [Standard Deviation-Mean Plot] [paper sdmp werk 20] [2008-12-18 13:23:55] [d8c5724db236abb5950452133b88474d] [Current]
Feedback Forum

Post a new message
Dataseries X:
106,2
101,6
99
96,4
94
93,2
103
103,6
103,2
102,2
100
99,6
98,8
95,2
91,6
88,6
86
84,8
95,2
96,2
94
92
90,2
90
88,8
85,8
84,2
80
77,8
76,8
86,4
89,2
86,2
84,6
83,2
83,2
82,6
79,8
77,2
74,8
73
73
83,6
85,6
84,8
84,2
83,4
84,6
84,6
83,8
81,2
79,6
78
78,2
88,8
92
91
91,2
90,4
91,8
92,2
90,2
88,6
87,8
86
87,2
97,6
101,2
100,4
100,2
100,2
103
104,2
104
102,4
101,8
101
102,2
114
118,4
118,8
117,2
117,2
118,4
118,8
117,2
114,4
112,6
111
110,8
120,2
124,4
123,4
121,2
119
119,8
120
118,4
115
113,4
111
111
121,6
126,2
125,8
124,8
122
123,2
124,2
120,8
116,8
114,8
111
109
119,8
124
121,6
118
115,8
116
115,8
114,4
112
110,2
107,4
108,2
117,6
121,4
119,8
115,6
112,6
113,2
112,2
110,8
108
105,2
102,4
101
110,8
116,8
113,8
108
104,4
105,2
105,4
103,2
100,6
97,8
95,8
95
104,8
110,4
106,4
102,2
98,4
98,4
98,6
96,2
92,4
91,4
88,4
87,8
97,6
104,2
100,2
97
92,8
92
93,4
92
89,6
88,6
87,2
86,2
96,8
102
102,6
100,6
94,2
94,2
95,2
95
94
92,2
91
91,2
103,4
105
104,6
103,8
101,8
102,4
103,8
103,4
102
101,8
100,2
101,4
113,8
116
115,6
113
109,4
111
112,4
112,2
111
108,8
107,4
108,6
118,8
122,2
122,6
122,2
118,8
119
118,2
117,8
116,8
114,6
113,4
113,8
124,2
125,8
125,6
122,4
119
119,4
118,6
118
116
114,8
114,6
114,6
124
125,2
124
117,6
113,2
111,4
112,2
109,8
106,4
105,2
102,2
99,8
111
113
108,4
105,4
102
102,8




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'George Udny Yule' @ 72.249.76.132

\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 & 'George Udny Yule' @ 72.249.76.132 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=34750&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]'George Udny Yule' @ 72.249.76.132[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=34750&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=34750&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 Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time1 seconds
R Server'George Udny Yule' @ 72.249.76.132







Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
1100.1666666666673.9779696359125713
291.88333333333334.2103408053232114
383.853.9497986139687012.4
480.554.8031239834091312.6
585.88333333333335.5951006706896314
694.556.4352014872743717
7109.9666666666677.831792240411517.8
8117.7333333333334.6006587143647813.6
9119.3666666666675.5614473808260515.2
10117.654.7437039613288815.2
11114.0166666666674.3228427795803214
12108.2166666666674.7942832118212915.8
13101.5333333333334.6821388533915.4
1494.88333333333334.9141416217675216.4
1593.955.6356817769378316.4
1698.35.6391488719486714
17107.6166666666676.084530807836815.8
18115.3333333333335.821173473283615.2
19119.254.3951419462029612.4
20117.6666666666674.5425930386896313.8
21106.5166666666674.3663659392166813.2

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 100.166666666667 & 3.97796963591257 & 13 \tabularnewline
2 & 91.8833333333333 & 4.21034080532321 & 14 \tabularnewline
3 & 83.85 & 3.94979861396870 & 12.4 \tabularnewline
4 & 80.55 & 4.80312398340913 & 12.6 \tabularnewline
5 & 85.8833333333333 & 5.59510067068963 & 14 \tabularnewline
6 & 94.55 & 6.43520148727437 & 17 \tabularnewline
7 & 109.966666666667 & 7.8317922404115 & 17.8 \tabularnewline
8 & 117.733333333333 & 4.60065871436478 & 13.6 \tabularnewline
9 & 119.366666666667 & 5.56144738082605 & 15.2 \tabularnewline
10 & 117.65 & 4.74370396132888 & 15.2 \tabularnewline
11 & 114.016666666667 & 4.32284277958032 & 14 \tabularnewline
12 & 108.216666666667 & 4.79428321182129 & 15.8 \tabularnewline
13 & 101.533333333333 & 4.68213885339 & 15.4 \tabularnewline
14 & 94.8833333333333 & 4.91414162176752 & 16.4 \tabularnewline
15 & 93.95 & 5.63568177693783 & 16.4 \tabularnewline
16 & 98.3 & 5.63914887194867 & 14 \tabularnewline
17 & 107.616666666667 & 6.0845308078368 & 15.8 \tabularnewline
18 & 115.333333333333 & 5.8211734732836 & 15.2 \tabularnewline
19 & 119.25 & 4.39514194620296 & 12.4 \tabularnewline
20 & 117.666666666667 & 4.54259303868963 & 13.8 \tabularnewline
21 & 106.516666666667 & 4.36636593921668 & 13.2 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=34750&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]100.166666666667[/C][C]3.97796963591257[/C][C]13[/C][/ROW]
[ROW][C]2[/C][C]91.8833333333333[/C][C]4.21034080532321[/C][C]14[/C][/ROW]
[ROW][C]3[/C][C]83.85[/C][C]3.94979861396870[/C][C]12.4[/C][/ROW]
[ROW][C]4[/C][C]80.55[/C][C]4.80312398340913[/C][C]12.6[/C][/ROW]
[ROW][C]5[/C][C]85.8833333333333[/C][C]5.59510067068963[/C][C]14[/C][/ROW]
[ROW][C]6[/C][C]94.55[/C][C]6.43520148727437[/C][C]17[/C][/ROW]
[ROW][C]7[/C][C]109.966666666667[/C][C]7.8317922404115[/C][C]17.8[/C][/ROW]
[ROW][C]8[/C][C]117.733333333333[/C][C]4.60065871436478[/C][C]13.6[/C][/ROW]
[ROW][C]9[/C][C]119.366666666667[/C][C]5.56144738082605[/C][C]15.2[/C][/ROW]
[ROW][C]10[/C][C]117.65[/C][C]4.74370396132888[/C][C]15.2[/C][/ROW]
[ROW][C]11[/C][C]114.016666666667[/C][C]4.32284277958032[/C][C]14[/C][/ROW]
[ROW][C]12[/C][C]108.216666666667[/C][C]4.79428321182129[/C][C]15.8[/C][/ROW]
[ROW][C]13[/C][C]101.533333333333[/C][C]4.68213885339[/C][C]15.4[/C][/ROW]
[ROW][C]14[/C][C]94.8833333333333[/C][C]4.91414162176752[/C][C]16.4[/C][/ROW]
[ROW][C]15[/C][C]93.95[/C][C]5.63568177693783[/C][C]16.4[/C][/ROW]
[ROW][C]16[/C][C]98.3[/C][C]5.63914887194867[/C][C]14[/C][/ROW]
[ROW][C]17[/C][C]107.616666666667[/C][C]6.0845308078368[/C][C]15.8[/C][/ROW]
[ROW][C]18[/C][C]115.333333333333[/C][C]5.8211734732836[/C][C]15.2[/C][/ROW]
[ROW][C]19[/C][C]119.25[/C][C]4.39514194620296[/C][C]12.4[/C][/ROW]
[ROW][C]20[/C][C]117.666666666667[/C][C]4.54259303868963[/C][C]13.8[/C][/ROW]
[ROW][C]21[/C][C]106.516666666667[/C][C]4.36636593921668[/C][C]13.2[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=34750&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=34750&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
SectionMeanStandard DeviationRange
1100.1666666666673.9779696359125713
291.88333333333334.2103408053232114
383.853.9497986139687012.4
480.554.8031239834091312.6
585.88333333333335.5951006706896314
694.556.4352014872743717
7109.9666666666677.831792240411517.8
8117.7333333333334.6006587143647813.6
9119.3666666666675.5614473808260515.2
10117.654.7437039613288815.2
11114.0166666666674.3228427795803214
12108.2166666666674.7942832118212915.8
13101.5333333333334.6821388533915.4
1494.88333333333334.9141416217675216.4
1593.955.6356817769378316.4
1698.35.6391488719486714
17107.6166666666676.084530807836815.8
18115.3333333333335.821173473283615.2
19119.254.3951419462029612.4
20117.6666666666674.5425930386896313.8
21106.5166666666674.3663659392166813.2







Regression: S.E.(k) = alpha + beta * Mean(k)
alpha4.68867808216129
beta0.00387581104486105
S.D.0.0174229472163172
T-STAT0.222454387121785
p-value0.82633399429414

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 4.68867808216129 \tabularnewline
beta & 0.00387581104486105 \tabularnewline
S.D. & 0.0174229472163172 \tabularnewline
T-STAT & 0.222454387121785 \tabularnewline
p-value & 0.82633399429414 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=34750&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]4.68867808216129[/C][/ROW]
[ROW][C]beta[/C][C]0.00387581104486105[/C][/ROW]
[ROW][C]S.D.[/C][C]0.0174229472163172[/C][/ROW]
[ROW][C]T-STAT[/C][C]0.222454387121785[/C][/ROW]
[ROW][C]p-value[/C][C]0.82633399429414[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=34750&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=34750&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)
alpha4.68867808216129
beta0.00387581104486105
S.D.0.0174229472163172
T-STAT0.222454387121785
p-value0.82633399429414







Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha1.18942548540492
beta0.0912866655123956
S.D.0.321613618733699
T-STAT0.283839552167666
p-value0.779603303816735
Lambda0.908713334487604

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & 1.18942548540492 \tabularnewline
beta & 0.0912866655123956 \tabularnewline
S.D. & 0.321613618733699 \tabularnewline
T-STAT & 0.283839552167666 \tabularnewline
p-value & 0.779603303816735 \tabularnewline
Lambda & 0.908713334487604 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=34750&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]1.18942548540492[/C][/ROW]
[ROW][C]beta[/C][C]0.0912866655123956[/C][/ROW]
[ROW][C]S.D.[/C][C]0.321613618733699[/C][/ROW]
[ROW][C]T-STAT[/C][C]0.283839552167666[/C][/ROW]
[ROW][C]p-value[/C][C]0.779603303816735[/C][/ROW]
[ROW][C]Lambda[/C][C]0.908713334487604[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=34750&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=34750&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)
alpha1.18942548540492
beta0.0912866655123956
S.D.0.321613618733699
T-STAT0.283839552167666
p-value0.779603303816735
Lambda0.908713334487604



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