FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

Author*The author of this computation has been verified*
R Software Modulerwasp_smp.wasp
Title produced by softwareStandard Deviation-Mean Plot
Date of computationMon, 15 Dec 2008 15:25: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/2008/Dec/15/t1229379970cfg5erjmmzxwtcs.htm/, Retrieved Tue, 14 May 2024 11:18:26 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=33836, Retrieved Tue, 14 May 2024 11:18:26 +0000
QR Codes:

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

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] [SMP prof bach] [2008-12-15 22:25:20] [21d7d81e7693ad6dde5aadefb1046611] [Current]
- RM      [Variance Reduction Matrix] [VRM prof bach] [2008-12-15 22:31:00] [bc937651ef42bf891200cf0e0edc7238]
- RMP       [(Partial) Autocorrelation Function] [ARIMA Prof bach A...] [2008-12-15 22:38:57] [bc937651ef42bf891200cf0e0edc7238]
-   P         [(Partial) Autocorrelation Function] [ARIMA ACF prof ba...] [2008-12-15 22:41:53] [bc937651ef42bf891200cf0e0edc7238]
-   P           [(Partial) Autocorrelation Function] [ARIMA ACF prof ba...] [2008-12-15 22:44:07] [bc937651ef42bf891200cf0e0edc7238]
-   P             [(Partial) Autocorrelation Function] [ARIMA ACF prof ba...] [2008-12-15 22:46:08] [bc937651ef42bf891200cf0e0edc7238]
-   P           [(Partial) Autocorrelation Function] [acf prof bach L =...] [2008-12-19 15:35:04] [bc937651ef42bf891200cf0e0edc7238]
-   P             [(Partial) Autocorrelation Function] [acf prof bach lam...] [2008-12-19 15:40:07] [bc937651ef42bf891200cf0e0edc7238]
-   P             [(Partial) Autocorrelation Function] [acf prof bach lam...] [2008-12-19 15:40:07] [bc937651ef42bf891200cf0e0edc7238]
-   P               [(Partial) Autocorrelation Function] [acf lambda = 1,1,1] [2008-12-19 15:45:45] [bc937651ef42bf891200cf0e0edc7238]
-  MPD                [(Partial) Autocorrelation Function] [ACF bij d=0 en D=0] [2010-12-17 13:29:13] [616fb52b46273b7e6805de1e68b3a688]
-   P                   [(Partial) Autocorrelation Function] [ACF bij d=0 en D=1] [2010-12-17 13:34:11] [616fb52b46273b7e6805de1e68b3a688]
-   P                     [(Partial) Autocorrelation Function] [ACF bij d=1 en D=1] [2010-12-17 13:51:58] [616fb52b46273b7e6805de1e68b3a688]
-  MPD            [(Partial) Autocorrelation Function] [] [2010-12-24 12:03:15] [4dfa50539945b119a90a7606969443b9]
-  MPD            [(Partial) Autocorrelation Function] [] [2010-12-24 12:12:12] [4dfa50539945b119a90a7606969443b9]
- RMP         [ARIMA Backward Selection] [Arima backward se...] [2008-12-19 17:26:16] [bc937651ef42bf891200cf0e0edc7238]
- RMP           [ARIMA Forecasting] [ARIMA forecast pr...] [2008-12-20 11:34:44] [bc937651ef42bf891200cf0e0edc7238]
-  MPD            [ARIMA Forecasting] [] [2010-12-21 19:37:30] [94f4aa1c01e87d8321fffb341ed4df07]
-    D              [ARIMA Forecasting] [] [2010-12-22 16:40:18] [94f4aa1c01e87d8321fffb341ed4df07]
-  MPD            [ARIMA Forecasting] [ARIMA Forecasting] [2010-12-22 13:25:50] [616fb52b46273b7e6805de1e68b3a688]
-    D              [ARIMA Forecasting] [ARIMA Forecasting] [2010-12-24 15:04:34] [616fb52b46273b7e6805de1e68b3a688]
- RMPD            [ARIMA Forecasting] [] [2010-12-24 13:59:54] [4dfa50539945b119a90a7606969443b9]
-   PD              [ARIMA Forecasting] [] [2010-12-26 10:14:56] [4dfa50539945b119a90a7606969443b9]
- RMP           [ARIMA Forecasting] [ARIMA forecast pr...] [2008-12-20 11:40:03] [bc937651ef42bf891200cf0e0edc7238]
-   P             [ARIMA Forecasting] [ARIMA forecast pr...] [2008-12-20 13:00:09] [bc937651ef42bf891200cf0e0edc7238]
-   P               [ARIMA Forecasting] [ARIMA voorspellin...] [2008-12-20 13:17:10] [bc937651ef42bf891200cf0e0edc7238]
-  M D          [ARIMA Backward Selection] [] [2010-12-21 17:53:23] [94f4aa1c01e87d8321fffb341ed4df07]
-  MPD          [ARIMA Backward Selection] [] [2010-12-21 19:17:29] [94f4aa1c01e87d8321fffb341ed4df07]
-  MPD          [ARIMA Backward Selection] [] [2010-12-24 13:46:46] [4dfa50539945b119a90a7606969443b9]
-  MPD          [ARIMA Backward Selection] [Paper Statistiek] [2010-12-28 15:46:48] [82c18f3ebe9df70882495121eb816e07]
-  MP           [ARIMA Backward Selection] [Paper Statistiek] [2010-12-28 16:09:56] [82c18f3ebe9df70882495121eb816e07]
- RMPD          [ARIMA Backward Selection] [ARIMA Backward Se...] [2010-12-28 18:28:51] [74be16979710d4c4e7c6647856088456]
-  MPD        [(Partial) Autocorrelation Function] [] [2010-12-20 16:16:49] [94f4aa1c01e87d8321fffb341ed4df07]
-  MPD        [(Partial) Autocorrelation Function] [b-r0245787] [2010-12-22 14:53:44] [ebb35fb07def4d07c0eb7ec8d2fd3b0e]
-               [(Partial) Autocorrelation Function] [b-r0245095] [2010-12-22 14:58:29] [ec8d68d52c1e9c5e97bb689b42436a8c]
-   PD            [(Partial) Autocorrelation Function] [b-r0245095] [2010-12-22 15:32:51] [ec8d68d52c1e9c5e97bb689b42436a8c]
-   PD          [(Partial) Autocorrelation Function] [b-r0245787] [2010-12-22 15:30:39] [ebb35fb07def4d07c0eb7ec8d2fd3b0e]
-  MPD        [(Partial) Autocorrelation Function] [] [2010-12-24 12:53:41] [4dfa50539945b119a90a7606969443b9]
-  MPD        [(Partial) Autocorrelation Function] [] [2010-12-28 12:28:43] [c6813a60da787bb62b5d86150b8926dd]
- R P           [(Partial) Autocorrelation Function] [Deel 4: ARIMA bac...] [2012-12-13 15:17:46] [b4e5b8b5af0253f45dc68b47bb41cf13]
- RMPD        [(Partial) Autocorrelation Function] [Autocorrelation] [2010-12-28 16:20:44] [74be16979710d4c4e7c6647856088456]
- RMPD        [ARIMA Backward Selection] [arima-model] [2010-12-29 22:10:14] [5a05da414fd67612c3b80d44effe0727]
- RMPD    [] [Standard Deviatio...] [-0001-11-30 00:00:00] [616fb52b46273b7e6805de1e68b3a688]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
13363
12530
11420
10948
10173
10602
16094
19631
17140
14345
12632
12894
11808
10673
9939
9890
9283
10131
15864
19283
16203
13919
11937
11795
11268
10522
9929
9725
9372
10068
16230
19115
18351
16265
14103
14115
13327
12618
12129
11775
11493
12470
20792
22337
21325
18581
16475
16581
15745
14453
13712
13766
13336
15346
24446
26178
24628
21282
18850
18822
18060
17536
16417
15842
15188
16905
25430
27962
26607
23364
20827
20506
19181
18016
17354
16256
15770
17538
26899
28915
25247
22856
19980
19856
16994
16839
15618
15883
15513
17106
25272
26731
22891
19583
16939
16757
15435
14786
13680
13208
12707
14277
22436
23229
18241
16145
13994
14780
13100
12329
12463
11532
10784
13106
19491
20418
16094
14491
13067

Tables (Output of Computation)





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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=33836&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 computational 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
1134812870.920916551159458
212560.41666666673118.426335606810000
313255.253565.663730767269743
415825.254073.4932769396710844
518380.33333333334724.8991972572712842
6203874468.9448012213812774
720655.66666666674340.1270893275913145
818843.83333333333921.5744016757811218
916076.53476.5332131044310522

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 13481 & 2870.92091655115 & 9458 \tabularnewline
2 & 12560.4166666667 & 3118.4263356068 & 10000 \tabularnewline
3 & 13255.25 & 3565.66373076726 & 9743 \tabularnewline
4 & 15825.25 & 4073.49327693967 & 10844 \tabularnewline
5 & 18380.3333333333 & 4724.89919725727 & 12842 \tabularnewline
6 & 20387 & 4468.94480122138 & 12774 \tabularnewline
7 & 20655.6666666667 & 4340.12708932759 & 13145 \tabularnewline
8 & 18843.8333333333 & 3921.57440167578 & 11218 \tabularnewline
9 & 16076.5 & 3476.53321310443 & 10522 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=33836&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]13481[/C][C]2870.92091655115[/C][C]9458[/C][/ROW]
[ROW][C]2[/C][C]12560.4166666667[/C][C]3118.4263356068[/C][C]10000[/C][/ROW]
[ROW][C]3[/C][C]13255.25[/C][C]3565.66373076726[/C][C]9743[/C][/ROW]
[ROW][C]4[/C][C]15825.25[/C][C]4073.49327693967[/C][C]10844[/C][/ROW]
[ROW][C]5[/C][C]18380.3333333333[/C][C]4724.89919725727[/C][C]12842[/C][/ROW]
[ROW][C]6[/C][C]20387[/C][C]4468.94480122138[/C][C]12774[/C][/ROW]
[ROW][C]7[/C][C]20655.6666666667[/C][C]4340.12708932759[/C][C]13145[/C][/ROW]
[ROW][C]8[/C][C]18843.8333333333[/C][C]3921.57440167578[/C][C]11218[/C][/ROW]
[ROW][C]9[/C][C]16076.5[/C][C]3476.53321310443[/C][C]10522[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=33836&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=33836&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
1134812870.920916551159458
212560.41666666673118.426335606810000
313255.253565.663730767269743
415825.254073.4932769396710844
518380.33333333334724.8991972572712842
6203874468.9448012213812774
720655.66666666674340.1270893275913145
818843.83333333333921.5744016757811218
916076.53476.5332131044310522






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha1032.83724379402
beta0.169036266077266
S.D.0.0420434255316822
T-STAT4.02051602455388
p-value0.00505629204279813

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 1032.83724379402 \tabularnewline
beta & 0.169036266077266 \tabularnewline
S.D. & 0.0420434255316822 \tabularnewline
T-STAT & 4.02051602455388 \tabularnewline
p-value & 0.00505629204279813 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=33836&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]1032.83724379402[/C][/ROW]
[ROW][C]beta[/C][C]0.169036266077266[/C][/ROW]
[ROW][C]S.D.[/C][C]0.0420434255316822[/C][/ROW]
[ROW][C]T-STAT[/C][C]4.02051602455388[/C][/ROW]
[ROW][C]p-value[/C][C]0.00505629204279813[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=33836&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=33836&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)
alpha1032.83724379402
beta0.169036266077266
S.D.0.0420434255316822
T-STAT4.02051602455388
p-value0.00505629204279813






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha1.01061763286943
beta0.745258285492207
S.D.0.180734952114381
T-STAT4.12348733199409
p-value0.00444009644399991
Lambda0.254741714507793

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & 1.01061763286943 \tabularnewline
beta & 0.745258285492207 \tabularnewline
S.D. & 0.180734952114381 \tabularnewline
T-STAT & 4.12348733199409 \tabularnewline
p-value & 0.00444009644399991 \tabularnewline
Lambda & 0.254741714507793 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=33836&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]1.01061763286943[/C][/ROW]
[ROW][C]beta[/C][C]0.745258285492207[/C][/ROW]
[ROW][C]S.D.[/C][C]0.180734952114381[/C][/ROW]
[ROW][C]T-STAT[/C][C]4.12348733199409[/C][/ROW]
[ROW][C]p-value[/C][C]0.00444009644399991[/C][/ROW]
[ROW][C]Lambda[/C][C]0.254741714507793[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=33836&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=33836&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)
alpha1.01061763286943
beta0.745258285492207
S.D.0.180734952114381
T-STAT4.12348733199409
p-value0.00444009644399991
Lambda0.254741714507793


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