Free Statistics

of Irreproducible Research!

Author's title

Author*Unverified author*
R Software Modulerwasp_smp.wasp
Title produced by softwareStandard Deviation-Mean Plot
Date of computationMon, 18 Aug 2014 07:30:50 +0100
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2014/Aug/18/t14083436112ut7k6uffm3v272.htm/, Retrieved Thu, 31 Oct 2024 23:31:47 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=235675, Retrieved Thu, 31 Oct 2024 23:31:47 +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)
-     [Notched Boxplots] [] [2013-02-26 20:31:34] [f974b105a61ab974a820d469d59cfaf7]
- RMPD  [Harrell-Davis Quantiles] [] [2014-08-18 05:22:13] [f974b105a61ab974a820d469d59cfaf7]
- RMP     [(Partial) Autocorrelation Function] [] [2014-08-18 06:01:25] [f974b105a61ab974a820d469d59cfaf7]
- RMP         [Standard Deviation-Mean Plot] [] [2014-08-18 06:30:50] [8f84a338303fe8d74ac0d8ad91c8b331] [Current]
Feedback Forum

Post a new message
Dataseries X:
95870
95523
95208
94541
101097
100781
95870
92581
92928
92928
93244
93910
93559
96190
97172
96190
99799
97488
92261
90964
90964
91599
89004
90964
89319
90964
93559
94541
96852
95870
89981
87670
86692
87670
86057
86692
84728
88021
89635
89981
96190
96190
88021
86057
86057
87039
82764
80799
78524
79186
82133
79821
86057
87039
80799
78524
77221
78524
74910
73613
68390
69688
70004
70355
76559
75893
68390
65097
63799
65444
59208
54964
47115
47777
47777
47115
52684
53004
46448
45151
42524
46133
39577
35653
28146
29764
27800
28462
33373
34355
31093
30742
30742
35017
27484
22573
14058
20929
19946
20293
28146
27164
23555
25204
25204
31093
24222
20293
14058
22258
21595
21911
28782
28146
25835
26186
27800
31409
25835
21275




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' @ yule.wessa.net

\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' @ yule.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=235675&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' @ yule.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=235675&T=0

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







Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
195373.41666666672855.153696176078516
293846.16666666673394.6436372972110795
390488.91666666673831.8009878653710795
487956.83333333334662.9917012643915391
579695.91666666673953.9202056217913426
667315.91666666676197.5169581833121595
745913.16666666674896.9416027359417351
829962.58333333333437.5131195341912444
923342.254549.3136370426117035
1024590.83333333334633.5305375517117351

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 95373.4166666667 & 2855.15369617607 & 8516 \tabularnewline
2 & 93846.1666666667 & 3394.64363729721 & 10795 \tabularnewline
3 & 90488.9166666667 & 3831.80098786537 & 10795 \tabularnewline
4 & 87956.8333333333 & 4662.99170126439 & 15391 \tabularnewline
5 & 79695.9166666667 & 3953.92020562179 & 13426 \tabularnewline
6 & 67315.9166666667 & 6197.51695818331 & 21595 \tabularnewline
7 & 45913.1666666667 & 4896.94160273594 & 17351 \tabularnewline
8 & 29962.5833333333 & 3437.51311953419 & 12444 \tabularnewline
9 & 23342.25 & 4549.31363704261 & 17035 \tabularnewline
10 & 24590.8333333333 & 4633.53053755171 & 17351 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=235675&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]95373.4166666667[/C][C]2855.15369617607[/C][C]8516[/C][/ROW]
[ROW][C]2[/C][C]93846.1666666667[/C][C]3394.64363729721[/C][C]10795[/C][/ROW]
[ROW][C]3[/C][C]90488.9166666667[/C][C]3831.80098786537[/C][C]10795[/C][/ROW]
[ROW][C]4[/C][C]87956.8333333333[/C][C]4662.99170126439[/C][C]15391[/C][/ROW]
[ROW][C]5[/C][C]79695.9166666667[/C][C]3953.92020562179[/C][C]13426[/C][/ROW]
[ROW][C]6[/C][C]67315.9166666667[/C][C]6197.51695818331[/C][C]21595[/C][/ROW]
[ROW][C]7[/C][C]45913.1666666667[/C][C]4896.94160273594[/C][C]17351[/C][/ROW]
[ROW][C]8[/C][C]29962.5833333333[/C][C]3437.51311953419[/C][C]12444[/C][/ROW]
[ROW][C]9[/C][C]23342.25[/C][C]4549.31363704261[/C][C]17035[/C][/ROW]
[ROW][C]10[/C][C]24590.8333333333[/C][C]4633.53053755171[/C][C]17351[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=235675&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=235675&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
195373.41666666672855.153696176078516
293846.16666666673394.6436372972110795
390488.91666666673831.8009878653710795
487956.83333333334662.9917012643915391
579695.91666666673953.9202056217913426
667315.91666666676197.5169581833121595
745913.16666666674896.9416027359417351
829962.58333333333437.5131195341912444
923342.254549.3136370426117035
1024590.83333333334633.5305375517117351







Regression: S.E.(k) = alpha + beta * Mean(k)
alpha4872.01800622897
beta-0.00987782657570734
S.D.0.0107063983005607
T-STAT-0.922609667453716
p-value0.383192441117397

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 4872.01800622897 \tabularnewline
beta & -0.00987782657570734 \tabularnewline
S.D. & 0.0107063983005607 \tabularnewline
T-STAT & -0.922609667453716 \tabularnewline
p-value & 0.383192441117397 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=235675&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]4872.01800622897[/C][/ROW]
[ROW][C]beta[/C][C]-0.00987782657570734[/C][/ROW]
[ROW][C]S.D.[/C][C]0.0107063983005607[/C][/ROW]
[ROW][C]T-STAT[/C][C]-0.922609667453716[/C][/ROW]
[ROW][C]p-value[/C][C]0.383192441117397[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=235675&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=235675&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)
alpha4872.01800622897
beta-0.00987782657570734
S.D.0.0107063983005607
T-STAT-0.922609667453716
p-value0.383192441117397







Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha9.48268417013394
beta-0.105413692737901
S.D.0.130872787069247
T-STAT-0.805466859066163
p-value0.443829715932083
Lambda1.1054136927379

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & 9.48268417013394 \tabularnewline
beta & -0.105413692737901 \tabularnewline
S.D. & 0.130872787069247 \tabularnewline
T-STAT & -0.805466859066163 \tabularnewline
p-value & 0.443829715932083 \tabularnewline
Lambda & 1.1054136927379 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=235675&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]9.48268417013394[/C][/ROW]
[ROW][C]beta[/C][C]-0.105413692737901[/C][/ROW]
[ROW][C]S.D.[/C][C]0.130872787069247[/C][/ROW]
[ROW][C]T-STAT[/C][C]-0.805466859066163[/C][/ROW]
[ROW][C]p-value[/C][C]0.443829715932083[/C][/ROW]
[ROW][C]Lambda[/C][C]1.1054136927379[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=235675&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=235675&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)
alpha9.48268417013394
beta-0.105413692737901
S.D.0.130872787069247
T-STAT-0.805466859066163
p-value0.443829715932083
Lambda1.1054136927379



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