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 03:41:15 -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/t12293377226rlsqistu1d33hi.htm/, Retrieved Wed, 15 May 2024 07:29:50 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=33656, Retrieved Wed, 15 May 2024 07:29:50 +0000
QR Codes:

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

Tree of Dependent Computations

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
F       [Standard Deviation-Mean Plot] [] [2008-12-15 10:41:15] [d41d8cd98f00b204e9800998ecf8427e] [Current]
Feedback Forum
2008-12-19 12:07:56 [006ad2c49b6a7c2ad6ab685cfc1dae56] [reply] 
De student feeft niet de juiste berekeningen gemaakt. Het was de bedoeling dat je ene tabel opnam in je document, zoals: http://www.freestatistics.org/blog/index.php?v=date/2008/Dec/13/t1229200568x9fasiirnug8ibn.htm

Uitleg bij deze tabel:
Kolom 1= tijdsindex
Kolom 2= Yt-waarden (gegevens)
Kolom 3= Ft= forecast= meest waarschijnlijke uitkomst, deze kolom moet je vergelijken met kolom 2
Kolom 4 en 5= lower en upper bound, met 95 % waarschijnlijkheid liggen de gegevens tussen deze grenzen
Kolom 6= p-waarde van een bepaalde toets (nulhypothese)= kans dat je je vergist bij het verwerpen van de nulhypothese (verschil tussen Yt en Ft)
Kolom 7= P= waarschijnlijkheid dat Ft groter is dan Yt-1 (stijging tov maand voordien)
Kolom 8= waarschijnlijkheid dat Ft groter is dan Yt-s (s= seizoenalitiet 12)= stijging tov dezelfde maand een jaar voordien
Kolom 9= waarschijnlijkheid dat Ft groter is dan de laatst gekende waarde

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
1961782
1787447
1953757
1828167
2078223
1777226
1538059
1678452
2262924
1616451
1829222
1763216
2138777
1565784
1781485
1962103
1775358
1837255
1917824
1878651
2124030
1948811
1947985
1719466
2139078
2086587
2020001
2405555
2152069
2791310
2397287
1885473
1978324
2165120
2052877
1726766
2267082
2501737
1916630
2177682
1859283
1718749
1819181
1463556
1979279
1723911
1528538
1635412
2255789
1698773
1635959
2054968
1794346
1938855
2112672
1446965
1610773
1576815
1509935
1769046

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'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=33656&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=33656&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=33656&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'George Udny Yule' @ 72.249.76.132






Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
11839577.16666667200789.105219125724865
21883127.41666667162624.566903739572993
32150037.25277491.6755914821064544
41882586.66666667308572.9572963351038181
51783741.33333333255691.743298034808824

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 1839577.16666667 & 200789.105219125 & 724865 \tabularnewline
2 & 1883127.41666667 & 162624.566903739 & 572993 \tabularnewline
3 & 2150037.25 & 277491.675591482 & 1064544 \tabularnewline
4 & 1882586.66666667 & 308572.957296335 & 1038181 \tabularnewline
5 & 1783741.33333333 & 255691.743298034 & 808824 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=33656&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]1839577.16666667[/C][C]200789.105219125[/C][C]724865[/C][/ROW]
[ROW][C]2[/C][C]1883127.41666667[/C][C]162624.566903739[/C][C]572993[/C][/ROW]
[ROW][C]3[/C][C]2150037.25[/C][C]277491.675591482[/C][C]1064544[/C][/ROW]
[ROW][C]4[/C][C]1882586.66666667[/C][C]308572.957296335[/C][C]1038181[/C][/ROW]
[ROW][C]5[/C][C]1783741.33333333[/C][C]255691.743298034[/C][C]808824[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=33656&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=33656&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
11839577.16666667200789.105219125724865
21883127.41666667162624.566903739572993
32150037.25277491.6755914821064544
41882586.66666667308572.9572963351038181
51783741.33333333255691.743298034808824






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha2694.23518199683
beta0.124928205078703
S.D.0.229402832429486
T-STAT0.544580046181875
p-value0.623913598087292

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 2694.23518199683 \tabularnewline
beta & 0.124928205078703 \tabularnewline
S.D. & 0.229402832429486 \tabularnewline
T-STAT & 0.544580046181875 \tabularnewline
p-value & 0.623913598087292 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=33656&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]2694.23518199683[/C][/ROW]
[ROW][C]beta[/C][C]0.124928205078703[/C][/ROW]
[ROW][C]S.D.[/C][C]0.229402832429486[/C][/ROW]
[ROW][C]T-STAT[/C][C]0.544580046181875[/C][/ROW]
[ROW][C]p-value[/C][C]0.623913598087292[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=33656&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=33656&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)
alpha2694.23518199683
beta0.124928205078703
S.D.0.229402832429486
T-STAT0.544580046181875
p-value0.623913598087292






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha-2.53459949191218
beta1.03057293912052
S.D.2.00797001637776
T-STAT0.51324119917866
p-value0.643188132298701
Lambda-0.0305729391205181

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & -2.53459949191218 \tabularnewline
beta & 1.03057293912052 \tabularnewline
S.D. & 2.00797001637776 \tabularnewline
T-STAT & 0.51324119917866 \tabularnewline
p-value & 0.643188132298701 \tabularnewline
Lambda & -0.0305729391205181 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=33656&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-2.53459949191218[/C][/ROW]
[ROW][C]beta[/C][C]1.03057293912052[/C][/ROW]
[ROW][C]S.D.[/C][C]2.00797001637776[/C][/ROW]
[ROW][C]T-STAT[/C][C]0.51324119917866[/C][/ROW]
[ROW][C]p-value[/C][C]0.643188132298701[/C][/ROW]
[ROW][C]Lambda[/C][C]-0.0305729391205181[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=33656&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=33656&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-2.53459949191218
beta1.03057293912052
S.D.2.00797001637776
T-STAT0.51324119917866
p-value0.643188132298701
Lambda-0.0305729391205181


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