FreeStatistics.org

Free Statistics

of Irreproducible Research!

Description of Statistical Computation

Author's title

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

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

Tree of Dependent Computations

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 12:42:29] [5e28000efa8060aa7512f63d330b190a] [Current]
Feedback Forum

Post a new message

Dataset

Dataseries X:
Download CSV
Histogram
Boxplots 
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

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

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






Standard Deviation-Mean Plot
SectionMeanStandard DeviationRange
1867273.556561.0085983857120788
292115031977.632380566772448
3963160.25140672.478436912306706
41536769.5110576.972194335266042
51587361100396.242675378230999
61783177.75329359.737534685715627
71509584.75106194.971608437190594
81588769.2544947.8637191654100786
91790778.25345873.648704625752805
101523878.5101185.903117315218481
111548839.2566417.0261774042124334
121685325.75213131.561597330489163
13154214392395.025461331202147
141653498.567026.1870540164164050
151825436.25241336.211493945554525
16146587856532.4286346636129169
171586647.2552097.9915567769123548
181727308.5267651.096226038576585

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 867273.5 & 56561.0085983857 & 120788 \tabularnewline
2 & 921150 & 31977.6323805667 & 72448 \tabularnewline
3 & 963160.25 & 140672.478436912 & 306706 \tabularnewline
4 & 1536769.5 & 110576.972194335 & 266042 \tabularnewline
5 & 1587361 & 100396.242675378 & 230999 \tabularnewline
6 & 1783177.75 & 329359.737534685 & 715627 \tabularnewline
7 & 1509584.75 & 106194.971608437 & 190594 \tabularnewline
8 & 1588769.25 & 44947.8637191654 & 100786 \tabularnewline
9 & 1790778.25 & 345873.648704625 & 752805 \tabularnewline
10 & 1523878.5 & 101185.903117315 & 218481 \tabularnewline
11 & 1548839.25 & 66417.0261774042 & 124334 \tabularnewline
12 & 1685325.75 & 213131.561597330 & 489163 \tabularnewline
13 & 1542143 & 92395.025461331 & 202147 \tabularnewline
14 & 1653498.5 & 67026.1870540164 & 164050 \tabularnewline
15 & 1825436.25 & 241336.211493945 & 554525 \tabularnewline
16 & 1465878 & 56532.4286346636 & 129169 \tabularnewline
17 & 1586647.25 & 52097.9915567769 & 123548 \tabularnewline
18 & 1727308.5 & 267651.096226038 & 576585 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=41822&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]867273.5[/C][C]56561.0085983857[/C][C]120788[/C][/ROW]
[ROW][C]2[/C][C]921150[/C][C]31977.6323805667[/C][C]72448[/C][/ROW]
[ROW][C]3[/C][C]963160.25[/C][C]140672.478436912[/C][C]306706[/C][/ROW]
[ROW][C]4[/C][C]1536769.5[/C][C]110576.972194335[/C][C]266042[/C][/ROW]
[ROW][C]5[/C][C]1587361[/C][C]100396.242675378[/C][C]230999[/C][/ROW]
[ROW][C]6[/C][C]1783177.75[/C][C]329359.737534685[/C][C]715627[/C][/ROW]
[ROW][C]7[/C][C]1509584.75[/C][C]106194.971608437[/C][C]190594[/C][/ROW]
[ROW][C]8[/C][C]1588769.25[/C][C]44947.8637191654[/C][C]100786[/C][/ROW]
[ROW][C]9[/C][C]1790778.25[/C][C]345873.648704625[/C][C]752805[/C][/ROW]
[ROW][C]10[/C][C]1523878.5[/C][C]101185.903117315[/C][C]218481[/C][/ROW]
[ROW][C]11[/C][C]1548839.25[/C][C]66417.0261774042[/C][C]124334[/C][/ROW]
[ROW][C]12[/C][C]1685325.75[/C][C]213131.561597330[/C][C]489163[/C][/ROW]
[ROW][C]13[/C][C]1542143[/C][C]92395.025461331[/C][C]202147[/C][/ROW]
[ROW][C]14[/C][C]1653498.5[/C][C]67026.1870540164[/C][C]164050[/C][/ROW]
[ROW][C]15[/C][C]1825436.25[/C][C]241336.211493945[/C][C]554525[/C][/ROW]
[ROW][C]16[/C][C]1465878[/C][C]56532.4286346636[/C][C]129169[/C][/ROW]
[ROW][C]17[/C][C]1586647.25[/C][C]52097.9915567769[/C][C]123548[/C][/ROW]
[ROW][C]18[/C][C]1727308.5[/C][C]267651.096226038[/C][C]576585[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=41822&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=41822&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
1867273.556561.0085983857120788
292115031977.632380566772448
3963160.25140672.478436912306706
41536769.5110576.972194335266042
51587361100396.242675378230999
61783177.75329359.737534685715627
71509584.75106194.971608437190594
81588769.2544947.8637191654100786
91790778.25345873.648704625752805
101523878.5101185.903117315218481
111548839.2566417.0261774042124334
121685325.75213131.561597330489163
13154214392395.025461331202147
141653498.567026.1870540164164050
151825436.25241336.211493945554525
16146587856532.4286346636129169
171586647.2552097.9915567769123548
181727308.5267651.096226038576585






Regression: S.E.(k) = alpha + beta * Mean(k)
alpha-147649.984415763
beta0.187480635882917
S.D.0.07224204366866
T-STAT2.59517348018007
p-value0.0195313043404523

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & -147649.984415763 \tabularnewline
beta & 0.187480635882917 \tabularnewline
S.D. & 0.07224204366866 \tabularnewline
T-STAT & 2.59517348018007 \tabularnewline
p-value & 0.0195313043404523 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=41822&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-147649.984415763[/C][/ROW]
[ROW][C]beta[/C][C]0.187480635882917[/C][/ROW]
[ROW][C]S.D.[/C][C]0.07224204366866[/C][/ROW]
[ROW][C]T-STAT[/C][C]2.59517348018007[/C][/ROW]
[ROW][C]p-value[/C][C]0.0195313043404523[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=41822&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=41822&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)
alpha-147649.984415763
beta0.187480635882917
S.D.0.07224204366866
T-STAT2.59517348018007
p-value0.0195313043404523






Regression: ln S.E.(k) = alpha + beta * ln Mean(k)
alpha-11.5631194313630
beta1.62845573380642
S.D.0.670699998812246
T-STAT2.42799423988412
p-value0.0273497287049102
Lambda-0.628455733806419

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & -11.5631194313630 \tabularnewline
beta & 1.62845573380642 \tabularnewline
S.D. & 0.670699998812246 \tabularnewline
T-STAT & 2.42799423988412 \tabularnewline
p-value & 0.0273497287049102 \tabularnewline
Lambda & -0.628455733806419 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=41822&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-11.5631194313630[/C][/ROW]
[ROW][C]beta[/C][C]1.62845573380642[/C][/ROW]
[ROW][C]S.D.[/C][C]0.670699998812246[/C][/ROW]
[ROW][C]T-STAT[/C][C]2.42799423988412[/C][/ROW]
[ROW][C]p-value[/C][C]0.0273497287049102[/C][/ROW]
[ROW][C]Lambda[/C][C]-0.628455733806419[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=41822&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=41822&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-11.5631194313630
beta1.62845573380642
S.D.0.670699998812246
T-STAT2.42799423988412
p-value0.0273497287049102
Lambda-0.628455733806419


Figures (Output of Computation)

Input Parameters & R Code

Parameters (Session):
par1 = 4 ;
Parameters (R input):
par1 = 4 ;
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')