Free Statistics

of Irreproducible Research!

Author's title

Standard deviation Mean Plot-omzetcijfers carrefour droge voeding- Angeliqu...

Author*Unverified author*
R Software Modulerwasp_smp.wasp
Title produced by softwareStandard Deviation-Mean Plot
Date of computationSat, 06 Jun 2009 04:55:43 -0600
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/Jun/06/t1244285795u85r04u0joy4s6u.htm/, Retrieved Sun, 28 Apr 2024 19:24:52 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=41955, Retrieved Sun, 28 Apr 2024 19:24:52 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact157
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Standard Deviation-Mean Plot] [Standard deviatio...] [2009-06-06 10:55:43] [5e28000efa8060aa7512f63d330b190a] [Current]
Feedback Forum

Post a new message
Dataseries X:
831581
808744
899237
929532
883165
908232
955613
937590
849396
978630
868513
1156102
1505713
1415151
1545021
1681193
1457973
1638575
1688972
1563924
1596359
1722061
1549332
2264959
1420268
1415099
1597279
1605693
1575400
1654752
1553966
1570959
1642414
1664774
1551560
2304365
1644081
1425600
1569344
1456489
1610786
1601519
1496600
1486452
1637939
1605759
1504221
1993384
1507620
1477037
1679184
1504731
1570141
1734191
1657498
1652164
1610941
1813765
1711573
2165466
1492778
1385488
1470589
1514657
1641395
1606185
1581162
1517847
1630080
1604623
1548973
2125558




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

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







Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
1917194.58333333390711.2292061095347358
21635769.41666667219057.449346661849808
31629710.75227031.354936283889266
41586014.5147908.767044541567784
51673692.58333333185073.569985490688429
61593277.91666667183302.354215630740070

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 917194.583333333 & 90711.2292061095 & 347358 \tabularnewline
2 & 1635769.41666667 & 219057.449346661 & 849808 \tabularnewline
3 & 1629710.75 & 227031.354936283 & 889266 \tabularnewline
4 & 1586014.5 & 147908.767044541 & 567784 \tabularnewline
5 & 1673692.58333333 & 185073.569985490 & 688429 \tabularnewline
6 & 1593277.91666667 & 183302.354215630 & 740070 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=41955&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]917194.583333333[/C][C]90711.2292061095[/C][C]347358[/C][/ROW]
[ROW][C]2[/C][C]1635769.41666667[/C][C]219057.449346661[/C][C]849808[/C][/ROW]
[ROW][C]3[/C][C]1629710.75[/C][C]227031.354936283[/C][C]889266[/C][/ROW]
[ROW][C]4[/C][C]1586014.5[/C][C]147908.767044541[/C][C]567784[/C][/ROW]
[ROW][C]5[/C][C]1673692.58333333[/C][C]185073.569985490[/C][C]688429[/C][/ROW]
[ROW][C]6[/C][C]1593277.91666667[/C][C]183302.354215630[/C][C]740070[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=41955&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=41955&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
1917194.58333333390711.2292061095347358
21635769.41666667219057.449346661849808
31629710.75227031.354936283889266
41586014.5147908.767044541567784
51673692.58333333185073.569985490688429
61593277.91666667183302.354215630740070







Regression: S.E.(k) = alpha + beta * Mean(k)
alpha-46383.9696761531
beta0.147348237940415
S.D.0.0456922659592005
T-STAT3.22479603160818
p-value0.0321338210901072

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & -46383.9696761531 \tabularnewline
beta & 0.147348237940415 \tabularnewline
S.D. & 0.0456922659592005 \tabularnewline
T-STAT & 3.22479603160818 \tabularnewline
p-value & 0.0321338210901072 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=41955&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-46383.9696761531[/C][/ROW]
[ROW][C]beta[/C][C]0.147348237940415[/C][/ROW]
[ROW][C]S.D.[/C][C]0.0456922659592005[/C][/ROW]
[ROW][C]T-STAT[/C][C]3.22479603160818[/C][/ROW]
[ROW][C]p-value[/C][C]0.0321338210901072[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=41955&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=41955&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-46383.9696761531
beta0.147348237940415
S.D.0.0456922659592005
T-STAT3.22479603160818
p-value0.0321338210901072







Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha-6.66070452219463
beta1.31599682510676
S.D.0.302820679962041
T-STAT4.34579575368409
p-value0.0121972739782109
Lambda-0.315996825106764

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & -6.66070452219463 \tabularnewline
beta & 1.31599682510676 \tabularnewline
S.D. & 0.302820679962041 \tabularnewline
T-STAT & 4.34579575368409 \tabularnewline
p-value & 0.0121972739782109 \tabularnewline
Lambda & -0.315996825106764 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=41955&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-6.66070452219463[/C][/ROW]
[ROW][C]beta[/C][C]1.31599682510676[/C][/ROW]
[ROW][C]S.D.[/C][C]0.302820679962041[/C][/ROW]
[ROW][C]T-STAT[/C][C]4.34579575368409[/C][/ROW]
[ROW][C]p-value[/C][C]0.0121972739782109[/C][/ROW]
[ROW][C]Lambda[/C][C]-0.315996825106764[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=41955&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=41955&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-6.66070452219463
beta1.31599682510676
S.D.0.302820679962041
T-STAT4.34579575368409
p-value0.0121972739782109
Lambda-0.315996825106764



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