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 computationFri, 04 Dec 2009 07:44:24 -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/2009/Dec/04/t12599379235t6o2glmkoxuri4.htm/, Retrieved Fri, 03 May 2024 16:39:35 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=63652, Retrieved Fri, 03 May 2024 16:39:35 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-     [Univariate Data Series] [data set] [2008-12-01 19:54:57] [b98453cac15ba1066b407e146608df68]
- RMP   [Standard Deviation-Mean Plot] [] [2009-11-27 14:40:44] [b98453cac15ba1066b407e146608df68]
F RMPD    [Spectral Analysis] [ws 9 4] [2009-12-04 14:39:57] [55b7a497389226c9339ee8d75ebc3b97]
F RMP         [Standard Deviation-Mean Plot] [ws 9 5] [2009-12-04 14:44:24] [84778c3520b84fd5786bccf2e25a5aef] [Current]
Feedback Forum
2009-12-10 17:11:05 [Brecht Thijs] [reply] 
Je hebt hier blijkbaar een fout gemaakt bij het instellen
van de seasonal period. Die moet namelijk op 12 staan aangezien
het de bedoeling is de maandcijfers jaar op jaar te vergelijken.
Dan krijg je volgende SMP
http://www.freestatistics.org/blog/index.php?v=date/2009/Dec/10/t1260463665nzm0qe4cy3q1r0t.htm/

Zoals eerder gezegd geeft de p-waarde van tabel aan
dat we Ho NIET verwerpen en we Lambda dus op 1 laten staan.

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
29.837
29.571
30.167
30.524
30.996
31.033
31.198
30.937
31.649
33.115
34.106
33.926
33.382
32.851
32.948
36.112
36.113
35.210
35.193
34.383
35.349
37.058
38.076
36.630
36.045
35.638
35.114
35.465
35.254
35.299
35.916
36.683
37.288
38.536
38.977
36.407
34.955
34.951
32.680
34.791
34.178
35.213
34.871
35.299
35.443
37.108
36.419
34.471
33.868
34.385
33.643
34.627
32.919
35.500
36.110
37.086
37.711
40.427
39.884
38.512
38.767

Tables (Output of Computation)





Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time4 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 & 4 seconds \tabularnewline
R Server & 'Gwilym Jenkins' @ 72.249.127.135 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=63652&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]4 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=63652&T=0

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






Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
130.024750.4125630254882280.953
231.0410.1118838683635850.260999999999999
333.1991.119632975577272.457
433.823251.543201299247773.261
535.224750.7068136836064591.73000000000000
636.778251.129476095364572.72700000000000
735.56550.3869233343527730.931000000000004
835.7880.6687515732068731.42900000000000
937.8021.173257857420952.57
1034.344251.112126304277832.275
1134.890250.5095471028276011.12100000000000
1235.860251.150825030141422.637
1334.130750.4538152156990770.984000000000002
1435.403751.780636660485984.167
1539.13351.244475927181132.716

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 30.02475 & 0.412563025488228 & 0.953 \tabularnewline
2 & 31.041 & 0.111883868363585 & 0.260999999999999 \tabularnewline
3 & 33.199 & 1.11963297557727 & 2.457 \tabularnewline
4 & 33.82325 & 1.54320129924777 & 3.261 \tabularnewline
5 & 35.22475 & 0.706813683606459 & 1.73000000000000 \tabularnewline
6 & 36.77825 & 1.12947609536457 & 2.72700000000000 \tabularnewline
7 & 35.5655 & 0.386923334352773 & 0.931000000000004 \tabularnewline
8 & 35.788 & 0.668751573206873 & 1.42900000000000 \tabularnewline
9 & 37.802 & 1.17325785742095 & 2.57 \tabularnewline
10 & 34.34425 & 1.11212630427783 & 2.275 \tabularnewline
11 & 34.89025 & 0.509547102827601 & 1.12100000000000 \tabularnewline
12 & 35.86025 & 1.15082503014142 & 2.637 \tabularnewline
13 & 34.13075 & 0.453815215699077 & 0.984000000000002 \tabularnewline
14 & 35.40375 & 1.78063666048598 & 4.167 \tabularnewline
15 & 39.1335 & 1.24447592718113 & 2.716 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=63652&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]30.02475[/C][C]0.412563025488228[/C][C]0.953[/C][/ROW]
[ROW][C]2[/C][C]31.041[/C][C]0.111883868363585[/C][C]0.260999999999999[/C][/ROW]
[ROW][C]3[/C][C]33.199[/C][C]1.11963297557727[/C][C]2.457[/C][/ROW]
[ROW][C]4[/C][C]33.82325[/C][C]1.54320129924777[/C][C]3.261[/C][/ROW]
[ROW][C]5[/C][C]35.22475[/C][C]0.706813683606459[/C][C]1.73000000000000[/C][/ROW]
[ROW][C]6[/C][C]36.77825[/C][C]1.12947609536457[/C][C]2.72700000000000[/C][/ROW]
[ROW][C]7[/C][C]35.5655[/C][C]0.386923334352773[/C][C]0.931000000000004[/C][/ROW]
[ROW][C]8[/C][C]35.788[/C][C]0.668751573206873[/C][C]1.42900000000000[/C][/ROW]
[ROW][C]9[/C][C]37.802[/C][C]1.17325785742095[/C][C]2.57[/C][/ROW]
[ROW][C]10[/C][C]34.34425[/C][C]1.11212630427783[/C][C]2.275[/C][/ROW]
[ROW][C]11[/C][C]34.89025[/C][C]0.509547102827601[/C][C]1.12100000000000[/C][/ROW]
[ROW][C]12[/C][C]35.86025[/C][C]1.15082503014142[/C][C]2.637[/C][/ROW]
[ROW][C]13[/C][C]34.13075[/C][C]0.453815215699077[/C][C]0.984000000000002[/C][/ROW]
[ROW][C]14[/C][C]35.40375[/C][C]1.78063666048598[/C][C]4.167[/C][/ROW]
[ROW][C]15[/C][C]39.1335[/C][C]1.24447592718113[/C][C]2.716[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=63652&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=63652&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
130.024750.4125630254882280.953
231.0410.1118838683635850.260999999999999
333.1991.119632975577272.457
433.823251.543201299247773.261
535.224750.7068136836064591.73000000000000
636.778251.129476095364572.72700000000000
735.56550.3869233343527730.931000000000004
835.7880.6687515732068731.42900000000000
937.8021.173257857420952.57
1034.344251.112126304277832.275
1134.890250.5095471028276011.12100000000000
1235.860251.150825030141422.637
1334.130750.4538152156990770.984000000000002
1435.403751.780636660485984.167
1539.13351.244475927181132.716






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha-2.46711194848567
beta0.096576894539679
S.D.0.0497723159357291
T-STAT1.94037373435443
p-value0.074337204251178

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & -2.46711194848567 \tabularnewline
beta & 0.096576894539679 \tabularnewline
S.D. & 0.0497723159357291 \tabularnewline
T-STAT & 1.94037373435443 \tabularnewline
p-value & 0.074337204251178 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=63652&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-2.46711194848567[/C][/ROW]
[ROW][C]beta[/C][C]0.096576894539679[/C][/ROW]
[ROW][C]S.D.[/C][C]0.0497723159357291[/C][/ROW]
[ROW][C]T-STAT[/C][C]1.94037373435443[/C][/ROW]
[ROW][C]p-value[/C][C]0.074337204251178[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=63652&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=63652&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)
alpha-2.46711194848567
beta0.096576894539679
S.D.0.0497723159357291
T-STAT1.94037373435443
p-value0.074337204251178






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha-22.1381547137850
beta6.15508092408769
S.D.2.38199212927824
T-STAT2.58400556762239
p-value0.0226837350964439
Lambda-5.15508092408769

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & -22.1381547137850 \tabularnewline
beta & 6.15508092408769 \tabularnewline
S.D. & 2.38199212927824 \tabularnewline
T-STAT & 2.58400556762239 \tabularnewline
p-value & 0.0226837350964439 \tabularnewline
Lambda & -5.15508092408769 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=63652&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-22.1381547137850[/C][/ROW]
[ROW][C]beta[/C][C]6.15508092408769[/C][/ROW]
[ROW][C]S.D.[/C][C]2.38199212927824[/C][/ROW]
[ROW][C]T-STAT[/C][C]2.58400556762239[/C][/ROW]
[ROW][C]p-value[/C][C]0.0226837350964439[/C][/ROW]
[ROW][C]Lambda[/C][C]-5.15508092408769[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=63652&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=63652&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)
alpha-22.1381547137850
beta6.15508092408769
S.D.2.38199212927824
T-STAT2.58400556762239
p-value0.0226837350964439
Lambda-5.15508092408769


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 4 ;
Parameters (R input):
par1 = 4 ;
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')