Free Statistics

of Irreproducible Research!

Author's title

Standard Deviation Mean Plot - Verkoop vrouwenkleding VS 2000-2008 - Amelie...

Author*Unverified author*
R Software Modulerwasp_smp.wasp
Title produced by softwareStandard Deviation-Mean Plot
Date of computationFri, 05 Jun 2009 10:43:34 -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/05/t1244220421c6ftqzlhnmtbnjj.htm/, Retrieved Fri, 10 May 2024 02:43:52 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=41894, Retrieved Fri, 10 May 2024 02:43:52 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact115
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-05 16:43:34] [770bd392bad9be4bfacb32bbbd82a4eb] [Current]
Feedback Forum

Post a new message
Dataseries X:
461683
731993
522673
572709
382812
942567
802385
272643
992660
672651
122826
493878
801948
742156
612673
662804
572750
202510
652313
492663
292397
382618
872790
913865
221989
842162
142783
112683
862759
392482
602262
542540
982375
112533
622742
883970
562056
322091
832619
82723
872811
562534
112447
342593
422622
682787
602944
84298
942258
722391
662901
3007
452988
12759
942577
52626
882755
22968
433067
904437
812315
32428
53073
403145
243133
313018
32720
922852
192947
693108
33335
944749
962536
292541
953227
193417
73378
433219
962953
632986
333186
733212
463421
495024
752596
22640
123421
243440
623653
413304
852981
613232
493203
253344
803665
574904
132570
132785
393406
333452
583642
893230
272991
373159
133072
653091
373307




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

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=41894&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
1580956.666666667261114.437587442869834
2600123.916666667224194.516216472711355
3526773.333333333325313.54305384869842
4456877.083333333276112.611830904790088
5502894.5392686.57236448939570
6389735.25359878.333439139912321
7544091.666666667306759.019600471889575
8480865.25272004.374590866830341

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 580956.666666667 & 261114.437587442 & 869834 \tabularnewline
2 & 600123.916666667 & 224194.516216472 & 711355 \tabularnewline
3 & 526773.333333333 & 325313.54305384 & 869842 \tabularnewline
4 & 456877.083333333 & 276112.611830904 & 790088 \tabularnewline
5 & 502894.5 & 392686.57236448 & 939570 \tabularnewline
6 & 389735.25 & 359878.333439139 & 912321 \tabularnewline
7 & 544091.666666667 & 306759.019600471 & 889575 \tabularnewline
8 & 480865.25 & 272004.374590866 & 830341 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=41894&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]580956.666666667[/C][C]261114.437587442[/C][C]869834[/C][/ROW]
[ROW][C]2[/C][C]600123.916666667[/C][C]224194.516216472[/C][C]711355[/C][/ROW]
[ROW][C]3[/C][C]526773.333333333[/C][C]325313.54305384[/C][C]869842[/C][/ROW]
[ROW][C]4[/C][C]456877.083333333[/C][C]276112.611830904[/C][C]790088[/C][/ROW]
[ROW][C]5[/C][C]502894.5[/C][C]392686.57236448[/C][C]939570[/C][/ROW]
[ROW][C]6[/C][C]389735.25[/C][C]359878.333439139[/C][C]912321[/C][/ROW]
[ROW][C]7[/C][C]544091.666666667[/C][C]306759.019600471[/C][C]889575[/C][/ROW]
[ROW][C]8[/C][C]480865.25[/C][C]272004.374590866[/C][C]830341[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=41894&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=41894&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
1580956.666666667261114.437587442869834
2600123.916666667224194.516216472711355
3526773.333333333325313.54305384869842
4456877.083333333276112.611830904790088
5502894.5392686.57236448939570
6389735.25359878.333439139912321
7544091.666666667306759.019600471889575
8480865.25272004.374590866830341







Regression: S.E.(k) = alpha + beta * Mean(k)
alpha531316.110112213
beta-0.448878706127334
S.D.0.274562079562608
T-STAT-1.63488966445192
p-value0.153189565693403

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 531316.110112213 \tabularnewline
beta & -0.448878706127334 \tabularnewline
S.D. & 0.274562079562608 \tabularnewline
T-STAT & -1.63488966445192 \tabularnewline
p-value & 0.153189565693403 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=41894&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]531316.110112213[/C][/ROW]
[ROW][C]beta[/C][C]-0.448878706127334[/C][/ROW]
[ROW][C]S.D.[/C][C]0.274562079562608[/C][/ROW]
[ROW][C]T-STAT[/C][C]-1.63488966445192[/C][/ROW]
[ROW][C]p-value[/C][C]0.153189565693403[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=41894&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=41894&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)
alpha531316.110112213
beta-0.448878706127334
S.D.0.274562079562608
T-STAT-1.63488966445192
p-value0.153189565693403







Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha22.2110439641711
beta-0.73140277755222
S.D.0.442251694477019
T-STAT-1.65381565901547
p-value0.149249302551201
Lambda1.73140277755222

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & 22.2110439641711 \tabularnewline
beta & -0.73140277755222 \tabularnewline
S.D. & 0.442251694477019 \tabularnewline
T-STAT & -1.65381565901547 \tabularnewline
p-value & 0.149249302551201 \tabularnewline
Lambda & 1.73140277755222 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=41894&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]22.2110439641711[/C][/ROW]
[ROW][C]beta[/C][C]-0.73140277755222[/C][/ROW]
[ROW][C]S.D.[/C][C]0.442251694477019[/C][/ROW]
[ROW][C]T-STAT[/C][C]-1.65381565901547[/C][/ROW]
[ROW][C]p-value[/C][C]0.149249302551201[/C][/ROW]
[ROW][C]Lambda[/C][C]1.73140277755222[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=41894&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=41894&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)
alpha22.2110439641711
beta-0.73140277755222
S.D.0.442251694477019
T-STAT-1.65381565901547
p-value0.149249302551201
Lambda1.73140277755222



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