Free Statistics

of Irreproducible Research!

Author's title

Standard deviation -mean plot - werkloosheid in australie - Magali vn de Wi...

Author*Unverified author*
R Software Modulerwasp_smp.wasp
Title produced by softwareStandard Deviation-Mean Plot
Date of computationSat, 23 May 2009 06:45:21 -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/May/23/t1243082851t0t5renn49omeaa.htm/, Retrieved Sat, 04 May 2024 04:26:38 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=40323, Retrieved Sat, 04 May 2024 04:26:38 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact209
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-05-23 12:45:21] [946afef1f31655e4e01d6b8ef4c25a2b] [Current]
Feedback Forum

Post a new message
Dataseries X:
428800
424800
403400
398400
393500
380500
398300
387300
370400
372800
444600
449900
458100
424800
420600
400100
393000
387100
377500
400400
391400
363600
431000
441700
448500
415600
408000
416600
409300
387600
394500
407600
378500
359600
435700
433800
427700
413300
379500
379300
353700
378200
380600
394000
374000
375000
437600
443900
488800
463900
440000
453800
451600
453400
461400
509100
540600
555100
677400
694600
750100
733900
709300
720500
693200
687200
686800
720900
653100
624700
690000
717800
736500
699900
675600
635600
632500
594900
604000
620800
578400
571200
627400
657700
674100
672800
615300
609100
607600
566900
572700
589200
534800
543100
591100
624800
665300
642600
608700
594500
563800
596100
597600
633100
591000
584200
655800
670700
699700
712900
652000
635100
603100
610100
602000
597600
585400
567100
620600
646200
644800
645200
644800
593000
569100
518800
538700
554600
507900
488400
563300
592400
598100
546300
516100
518500
477400
483400
469400
501300
457400
446700
501900
550400
593700
548900
534200
550500
541800
569300
587400
627700
607000
629500
704600
767700
812200
824600
856300
812200
764100
801700
806000
867200
801600
817500
920900
959700
997700
949100
910900
920400
914200
926300
906400
926100
902500
895300
979900
1009700
1043800
979800
921600
923500
914500
891700
916000
931700
902400
893700
941500
980100
1006900
949200
883200
849900
839200
803900
797900
830800
753300
764100
807600
853700
886200
815700
743000
753600
724800
709600
721900




Summary of computational transaction
Raw Inputview raw input (R code)
Raw Outputview raw output of R engine
Computing time2 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24

\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 & 2 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ 193.190.124.24 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=40323&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]2 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Ronald Aylmer Fisher' @ 193.190.124.24[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=40323&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=40323&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 time2 seconds
R Server'Sir Ronald Aylmer Fisher' @ 193.190.124.24







Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
1404391.66666666726814.190894623979500
2407441.66666666727914.006205205394500
3407941.66666666725258.534411284888900
4394733.33333333328814.968753105190200
5515808.33333333387529.4387707082254600
6698958.33333333334720.4539751776125400
7636208.33333333349331.3556626529165300
860012544009.8775070077139300
961695035170.9461707452106900
1062765044132.0230053581145800
1157175054136.2674604138156800
1250557543609.4052616259151400
13605191.66666666770289.5499582501233500
1483700055528.6821872674195600
15936541.66666666738764.427421714114400
16936691.66666666744347.4803698194152100
1784497573313.487845007253600

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 404391.666666667 & 26814.1908946239 & 79500 \tabularnewline
2 & 407441.666666667 & 27914.0062052053 & 94500 \tabularnewline
3 & 407941.666666667 & 25258.5344112848 & 88900 \tabularnewline
4 & 394733.333333333 & 28814.9687531051 & 90200 \tabularnewline
5 & 515808.333333333 & 87529.4387707082 & 254600 \tabularnewline
6 & 698958.333333333 & 34720.4539751776 & 125400 \tabularnewline
7 & 636208.333333333 & 49331.3556626529 & 165300 \tabularnewline
8 & 600125 & 44009.8775070077 & 139300 \tabularnewline
9 & 616950 & 35170.9461707452 & 106900 \tabularnewline
10 & 627650 & 44132.0230053581 & 145800 \tabularnewline
11 & 571750 & 54136.2674604138 & 156800 \tabularnewline
12 & 505575 & 43609.4052616259 & 151400 \tabularnewline
13 & 605191.666666667 & 70289.5499582501 & 233500 \tabularnewline
14 & 837000 & 55528.6821872674 & 195600 \tabularnewline
15 & 936541.666666667 & 38764.427421714 & 114400 \tabularnewline
16 & 936691.666666667 & 44347.4803698194 & 152100 \tabularnewline
17 & 844975 & 73313.487845007 & 253600 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=40323&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]404391.666666667[/C][C]26814.1908946239[/C][C]79500[/C][/ROW]
[ROW][C]2[/C][C]407441.666666667[/C][C]27914.0062052053[/C][C]94500[/C][/ROW]
[ROW][C]3[/C][C]407941.666666667[/C][C]25258.5344112848[/C][C]88900[/C][/ROW]
[ROW][C]4[/C][C]394733.333333333[/C][C]28814.9687531051[/C][C]90200[/C][/ROW]
[ROW][C]5[/C][C]515808.333333333[/C][C]87529.4387707082[/C][C]254600[/C][/ROW]
[ROW][C]6[/C][C]698958.333333333[/C][C]34720.4539751776[/C][C]125400[/C][/ROW]
[ROW][C]7[/C][C]636208.333333333[/C][C]49331.3556626529[/C][C]165300[/C][/ROW]
[ROW][C]8[/C][C]600125[/C][C]44009.8775070077[/C][C]139300[/C][/ROW]
[ROW][C]9[/C][C]616950[/C][C]35170.9461707452[/C][C]106900[/C][/ROW]
[ROW][C]10[/C][C]627650[/C][C]44132.0230053581[/C][C]145800[/C][/ROW]
[ROW][C]11[/C][C]571750[/C][C]54136.2674604138[/C][C]156800[/C][/ROW]
[ROW][C]12[/C][C]505575[/C][C]43609.4052616259[/C][C]151400[/C][/ROW]
[ROW][C]13[/C][C]605191.666666667[/C][C]70289.5499582501[/C][C]233500[/C][/ROW]
[ROW][C]14[/C][C]837000[/C][C]55528.6821872674[/C][C]195600[/C][/ROW]
[ROW][C]15[/C][C]936541.666666667[/C][C]38764.427421714[/C][C]114400[/C][/ROW]
[ROW][C]16[/C][C]936691.666666667[/C][C]44347.4803698194[/C][C]152100[/C][/ROW]
[ROW][C]17[/C][C]844975[/C][C]73313.487845007[/C][C]253600[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=40323&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=40323&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
1404391.66666666726814.190894623979500
2407441.66666666727914.006205205394500
3407941.66666666725258.534411284888900
4394733.33333333328814.968753105190200
5515808.33333333387529.4387707082254600
6698958.33333333334720.4539751776125400
7636208.33333333349331.3556626529165300
860012544009.8775070077139300
961695035170.9461707452106900
1062765044132.0230053581145800
1157175054136.2674604138156800
1250557543609.4052616259151400
13605191.66666666770289.5499582501233500
1483700055528.6821872674195600
15936541.66666666738764.427421714114400
16936691.66666666744347.4803698194152100
1784497573313.487845007253600







Regression: S.E.(k) = alpha + beta * Mean(k)
alpha26894.9288693815
beta0.0309512105133214
S.D.0.0239476930894724
T-STAT1.29245060881993
p-value0.215756932034609

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 26894.9288693815 \tabularnewline
beta & 0.0309512105133214 \tabularnewline
S.D. & 0.0239476930894724 \tabularnewline
T-STAT & 1.29245060881993 \tabularnewline
p-value & 0.215756932034609 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=40323&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]26894.9288693815[/C][/ROW]
[ROW][C]beta[/C][C]0.0309512105133214[/C][/ROW]
[ROW][C]S.D.[/C][C]0.0239476930894724[/C][/ROW]
[ROW][C]T-STAT[/C][C]1.29245060881993[/C][/ROW]
[ROW][C]p-value[/C][C]0.215756932034609[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=40323&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=40323&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)
alpha26894.9288693815
beta0.0309512105133214
S.D.0.0239476930894724
T-STAT1.29245060881993
p-value0.215756932034609







Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha2.55766394557760
beta0.61035401762357
S.D.0.28429109163266
T-STAT2.14693332147116
p-value0.0485561935380429
Lambda0.38964598237643

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & 2.55766394557760 \tabularnewline
beta & 0.61035401762357 \tabularnewline
S.D. & 0.28429109163266 \tabularnewline
T-STAT & 2.14693332147116 \tabularnewline
p-value & 0.0485561935380429 \tabularnewline
Lambda & 0.38964598237643 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=40323&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]2.55766394557760[/C][/ROW]
[ROW][C]beta[/C][C]0.61035401762357[/C][/ROW]
[ROW][C]S.D.[/C][C]0.28429109163266[/C][/ROW]
[ROW][C]T-STAT[/C][C]2.14693332147116[/C][/ROW]
[ROW][C]p-value[/C][C]0.0485561935380429[/C][/ROW]
[ROW][C]Lambda[/C][C]0.38964598237643[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=40323&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=40323&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)
alpha2.55766394557760
beta0.61035401762357
S.D.0.28429109163266
T-STAT2.14693332147116
p-value0.0485561935380429
Lambda0.38964598237643



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