## Free Statistics

of Irreproducible Research!

Author's title
Author*The author of this computation has been verified*
R Software Modulerwasp_smp.wasp
Title produced by softwareStandard Deviation-Mean Plot
Date of computationThu, 12 Dec 2013 03:45:35 -0500
Cite this page as followsStatistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?v=date/2013/Dec/12/t13868379621hef552ave9kw17.htm/, Retrieved Sun, 05 Dec 2021 18:03:58 +0000
Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL https://freestatistics.org/blog/index.php?pk=232219, Retrieved Sun, 05 Dec 2021 18:03:58 +0000
QR Codes:

Original text written by user:
IsPrivate?No (this computation is public)
User-defined keywords
Estimated Impact81
Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)
-       [Standard Deviation-Mean Plot] [] [2013-12-12 08:45:35] [9e6a405f514733ea23d87e4507d39d29] [Current]
Feedback Forum

Post a new message
Dataseries X:
164
96
73
49
39
59
169
169
210
278
298
245
200
188
90
79
78
91
167
169
289
247
275
203
223
104
107
85
75
99
135
211
335
488
326
346
261
224
141
148
145
223
272
445
560
612
467
404
518
404
300
210
196
186
247
343
464
680
711
610
513
292
273
322
189
257
324
404
677
858
895
664
628
308
324
248
272

 Summary of computational transaction Raw Input view raw input (R code) Raw Output view raw output of R engine Computing time 9 seconds R Server 'Sir Ronald Aylmer Fisher' @ fisher.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 & 9 seconds \tabularnewline
R Server & 'Sir Ronald Aylmer Fisher' @ fisher.wessa.net \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=232219&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]9 seconds[/C][/ROW]
[ROW][C]R Server[/C][C]'Sir Ronald Aylmer Fisher' @ fisher.wessa.net[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=232219&T=0

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=232219&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 Input view raw input (R code) Raw Output view raw output of R engine Computing time 9 seconds R Server 'Sir Ronald Aylmer Fisher' @ fisher.wessa.net

 Standard Deviation-Mean Plot Section Mean Standard Deviation Range 1 154.083333333333 91.0938161091813 259 2 173 75.5079767008397 211 3 211.166666666667 134.386033726002 413 4 325.166666666667 166.10338146267 471 5 405.75 189.983312664321 525 6 472.333333333333 244.175173085871 706

\begin{tabular}{lllllllll}
\hline
Standard Deviation-Mean Plot \tabularnewline
Section & Mean & Standard Deviation & Range \tabularnewline
1 & 154.083333333333 & 91.0938161091813 & 259 \tabularnewline
2 & 173 & 75.5079767008397 & 211 \tabularnewline
3 & 211.166666666667 & 134.386033726002 & 413 \tabularnewline
4 & 325.166666666667 & 166.10338146267 & 471 \tabularnewline
5 & 405.75 & 189.983312664321 & 525 \tabularnewline
6 & 472.333333333333 & 244.175173085871 & 706 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=232219&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]154.083333333333[/C][C]91.0938161091813[/C][C]259[/C][/ROW]
[ROW][C]2[/C][C]173[/C][C]75.5079767008397[/C][C]211[/C][/ROW]
[ROW][C]3[/C][C]211.166666666667[/C][C]134.386033726002[/C][C]413[/C][/ROW]
[ROW][C]4[/C][C]325.166666666667[/C][C]166.10338146267[/C][C]471[/C][/ROW]
[ROW][C]5[/C][C]405.75[/C][C]189.983312664321[/C][C]525[/C][/ROW]
[ROW][C]6[/C][C]472.333333333333[/C][C]244.175173085871[/C][C]706[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=232219&T=1

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=232219&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 Section Mean Standard Deviation Range 1 154.083333333333 91.0938161091813 259 2 173 75.5079767008397 211 3 211.166666666667 134.386033726002 413 4 325.166666666667 166.10338146267 471 5 405.75 189.983312664321 525 6 472.333333333333 244.175173085871 706

 Regression: S.E.(k) = alpha + beta * Mean(k) alpha 14.5397976742765 beta 0.467419412979171 S.D. 0.0575820724704363 T-STAT 8.11744685325719 p-value 0.00125245910781976

\begin{tabular}{lllllllll}
\hline
Regression: S.E.(k) = alpha + beta * Mean(k) \tabularnewline
alpha & 14.5397976742765 \tabularnewline
beta & 0.467419412979171 \tabularnewline
S.D. & 0.0575820724704363 \tabularnewline
T-STAT & 8.11744685325719 \tabularnewline
p-value & 0.00125245910781976 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=232219&T=2

[TABLE]
[ROW][C]Regression: S.E.(k) = alpha + beta * Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]14.5397976742765[/C][/ROW]
[ROW][C]beta[/C][C]0.467419412979171[/C][/ROW]
[ROW][C]S.D.[/C][C]0.0575820724704363[/C][/ROW]
[ROW][C]T-STAT[/C][C]8.11744685325719[/C][/ROW]
[ROW][C]p-value[/C][C]0.00125245910781976[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=232219&T=2

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=232219&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 14.5397976742765 beta 0.467419412979171 S.D. 0.0575820724704363 T-STAT 8.11744685325719 p-value 0.00125245910781976

 Regression: ln S.E.(k) = alpha + beta * ln Mean(k) alpha -0.173310718003944 beta 0.914692457228613 S.D. 0.145035687678005 T-STAT 6.30667163284206 p-value 0.0032318403537262 Lambda 0.0853075427713865

\begin{tabular}{lllllllll}
\hline
Regression: ln S.E.(k) = alpha + beta * ln Mean(k) \tabularnewline
alpha & -0.173310718003944 \tabularnewline
beta & 0.914692457228613 \tabularnewline
S.D. & 0.145035687678005 \tabularnewline
T-STAT & 6.30667163284206 \tabularnewline
p-value & 0.0032318403537262 \tabularnewline
Lambda & 0.0853075427713865 \tabularnewline
\hline
\end{tabular}
%Source: https://freestatistics.org/blog/index.php?pk=232219&T=3

[TABLE]
[ROW][C]Regression: ln S.E.(k) = alpha + beta * ln Mean(k)[/C][/ROW]
[ROW][C]alpha[/C][C]-0.173310718003944[/C][/ROW]
[ROW][C]beta[/C][C]0.914692457228613[/C][/ROW]
[ROW][C]S.D.[/C][C]0.145035687678005[/C][/ROW]
[ROW][C]T-STAT[/C][C]6.30667163284206[/C][/ROW]
[ROW][C]p-value[/C][C]0.0032318403537262[/C][/ROW]
[ROW][C]Lambda[/C][C]0.0853075427713865[/C][/ROW]
[/TABLE]
Source: https://freestatistics.org/blog/index.php?pk=232219&T=3

Globally Unique Identifier (entire table): ba.freestatistics.org/blog/index.php?pk=232219&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 -0.173310718003944 beta 0.914692457228613 S.D. 0.145035687678005 T-STAT 6.30667163284206 p-value 0.0032318403537262 Lambda 0.0853075427713865

par1 <- as.numeric(par1)(n <- length(x))(np <- floor(n / par1))arr <- array(NA,dim=c(par1,np))j <- 0k <- 1for (i in 1:(np*par1)){j = j + 1arr[j,k] <- x[i]if (j == par1) {j = 0k=k+1}}arrarr.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.meanarr.sdarr.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')