Taak 8 - Step 1 (5)

*The author of this computation has been verified*

R Software Module: rwasp_smp.wasp (opens new window with default values)

Title produced by software: Standard Deviation-Mean Plot

Date of computation: Fri, 05 Dec 2008 03:40:23 -0700

Cite this page as follows:

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2008/Dec/05/t1228473660jsxk24w6x2g9lxd.htm/, Retrieved Fri, 05 Dec 2008 10:41:00 +0000

BibTeX entries for LaTeX users:

@Manual{KEY, author = {{YOUR NAME}}, publisher = {Office for Research Development and Education}, title = {Statistical Computations at FreeStatistics.org, URL http://www.freestatistics.org/blog/date/2008/Dec/05/t1228473660jsxk24w6x2g9lxd.htm/}, year = {2008}, } @Manual{R, title = {R: A Language and Environment for Statistical Computing}, author = {{R Development Core Team}}, organization = {R Foundation for Statistical Computing}, address = {Vienna, Austria}, year = {2008}, note = {{ISBN} 3-900051-07-0}, url = {http://www.R-project.org}, }

Family? (F = Feedback message, R = changed R code, M = changed R Module, P = changed Parameters, D = changed Data)

Feedback Forum:

2008-11-27 13:41:43 [a2386b643d711541400692649981f2dc]  [reply]  test Post a new message

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords:

Dataseries X:

» Textbox « » Textfile « » CSV «

267.850 255.468 246.717 238.276 229.336 231.047 264.440 261.308 257.844 253.687 248.286 248.641 243.753 233.495 222.356 211.333 202.478 202.746 238.802 235.385 224.820 219.386 213.106 213.166 208.201 197.775 189.191 178.543 171.809 172.365 204.140 205.467 193.079 190.224 185.553 185.184 183.059 175.496 168.120 160.374 155.290 156.152 189.784 192.250 184.896 184.835 180.172 181.875 182.412 180.627 174.303 169.431 163.902 166.114 198.414 205.626 199.333 199.588 196.569 200.880 201.579 195.483 190.617 187.576 183.968 186.998 216.617 224.692 222.476 223.247 225.618 232.758 235.868 232.863 227.564 226.822 223.864 227.155 260.300 273.944 270.779 268.104 268.703 273.413 275.597 270.111 262.272 255.823 250.753 250.512 277.888 289.694 281.310 275.425 271.287 274.059 274.113 267.546 257.622 250.612 243.829 243.180 275.362 287.027 279.175 282.416 275.424 277.862 284.998 272.182 258.613 253.046 243.315 234.312 265.912 279.384 262.547 256.102 251.133 249.598 251.592 244.976 237.527 232.790 224.726 223.918 252.637 263.736 251.143 239.530 232.401 233.749 232.055 224.473 215.866 207.808 199.440 193.330 222.787 241.434 221.263 207.448 200.241 205.009 206.230 198.253 194.660 185.847 180.314 176.282 203.541 222.042 197.519 185.142 176.355 180.448 180.143 173.666 165.687 162.719 157.079 153.730 182.698 200.765 176.512 166.618 158.644 159.585 163.095 159.044 155.511 153.745 150.569 150.605 179.612 194.690 189.917 184.128 175.335 179.566 181.140 177.876 175.041 169.292 166.070 166.972 206.348 215.706 202.108 195.411 193.111 195.198 198.770 194.163 190.420 189.733 186.029 191.531 232.571 243.477 227.247 217.859 208.679 213.188 216.234 213.587 209.465 204.045 200.237 203.666 241.476 260.307 243.324 244.460 233.575 237.217 235.243 230.354 227.184 221.678 217.142 219.452 256.446 265.845 248.624 241.114 229.245 231.805 219.277 219.313 212.610 214.771 211.142 211.457 240.048 240.636 230.580 208.795 197.922 194.596 194.581 185.686 178.106 172.608 167.302 168.053 202.300 202.388 182.516 173.476 166.444 171.297

Output produced by software:

Summary of computational transaction Raw Input view raw input (R code) Raw Output view raw output of R engine Computing time 2 seconds R Server 'Herman Ole Andreas Wold' @ 193.190.124.10:1001 Standard Deviation-Mean Plot Section Mean Standard Deviation Range 1 250.241666666667 12.4263631286898 38.514 2 221.7355 13.8288852932746 41.275 3 190.127583333333 12.2734839942072 36.392 4 176.02525 12.9417145563905 36.96 5 186.43325 15.2834570168651 41.724 6 207.63575 18.2221032421667 48.79 7 249.114916666667 21.4613089385541 50.08 8 269.560916666667 12.3049008008881 39.182 9 267.847333333333 15.2223297185617 43.847 10 259.261833333333 14.6923140316990 50.686 11 240.727083333333 12.2367385840445 39.818 12 214.262833333333 14.4962906153975 48.104 13 192.219416666667 14.0161952200649 45.76 14 169.8205 13.5452537511048 47.035 15 169.651416666667 15.9897275743371 44.121 16 187.02275 16.5564737249375 49.636 17 207.805583333333 19.1412672491012 57.448 18 225.63275 20.0144900787813 60.07 19 235.344333333333 15.0668297644699 48.703 20 216.76225 14.4808854458370 46.04 21 180.396416666667 13.1845709997318 35.944 Regression: S.E.(k) = alpha + beta * Mean(k) alpha 14.0434030040856 beta 0.0049747597655878 S.D. 0.0186340672745821 T-STAT 0.266971224922733 p-value 0.792367481159726 Regression: ln S.E.(k) = alpha + beta * ln Mean(k) alpha 2.37185874740021 beta 0.0616260555142846 S.D. 0.250286235080674 T-STAT 0.246222312203549 p-value 0.808150577873492 Lambda 0.938373944485715

Charts produced by software:

http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228473660jsxk24w6x2g9lxd/17uyp1228473617.png (open in new window) http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228473660jsxk24w6x2g9lxd/17uyp1228473617.ps (open in new window) http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228473660jsxk24w6x2g9lxd/2gfcd1228473617.png (open in new window) http://127.0.0.1/wessadotnet/public_html/freestatisticsdotorg/blog/date/2008/Dec/05/t1228473660jsxk24w6x2g9lxd/2gfcd1228473617.ps (open in new window)

Parameters (Session):

par1 = 36 ; par2 = 1.2 ; par3 = 1 ; par4 = 1 ; par5 = 12 ;

Parameters (R input):

par1 = 12 ; par2 = 1.2 ; par3 = 1 ; par4 = 1 ; par5 = 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')

Copyright

This work is licensed under a Creative Commons Attribution-Noncommercial-Share Alike 3.0 License.

Software written by Ed van Stee & Patrick Wessa

Disclaimer

Information provided on this web site is provided "AS IS" without warranty of any kind, either express or implied, including, without limitation, warranties of merchantability, fitness for a particular purpose, and noninfringement. We use reasonable efforts to include accurate and timely information and periodically update the information, and software without notice. However, we make no warranties or representations as to the accuracy or completeness of such information (or software), and we assume no liability or responsibility for errors or omissions in the content of this web site, or any software bugs in online applications. Your use of this web site is AT YOUR OWN RISK. Under no circumstances and under no legal theory shall we be liable to you or any other person for any direct, indirect, special, incidental, exemplary, or consequential damages arising from your access to, or use of, this web site.

Privacy Policy

We may request personal information to be submitted to our servers in order to be able to:

• personalize online software applications according to your needs
• enforce strict security rules with respect to the data that you upload (e.g. statistical data)
• manage user sessions of online applications
• alert you about important changes or upgrades in resources or applications

We NEVER allow other companies to directly offer registered users information about their products and services. Banner references and hyperlinks of third parties NEVER contain any personal data of the visitor.

We do NOT sell, nor transmit by any means, personal information, nor statistical data series uploaded by you to third parties.

We carefully protect your data from loss, misuse, alteration, and destruction. However, at any time, and under any circumstance you are solely responsible for managing your passwords, and keeping them secret.

We store a unique ANONYMOUS USER ID in the form of a small 'Cookie' on your computer. This allows us to track your progress when using this website which is necessary to create state-dependent features. The cookie is used for NO OTHER PURPOSE. At any time you may opt to disallow cookies from this website - this will not affect other features of this website.

We examine cookies that are used by third-parties (banner and online ads) very closely: abuse from third-parties automatically results in termination of the advertising contract without refund. We have very good reason to believe that the cookies that are produced by third parties (banner ads) do NOT cause any privacy or security risk.

FreeStatistics.org is safe. There is no need to download any software to use the applications and services contained in this website. Hence, your system's security is not compromised by their use, and your personal data - other than data you submit in the account application form, and the user-agent information that is transmitted by your browser - is never transmitted to our servers.

As a general rule, we do not log on-line behavior of individuals (other than normal logging of webserver 'hits'). However, in cases of abuse, hacking, unauthorized access, Denial of Service attacks, illegal copying, hotlinking, non-compliance with international webstandards (such as robots.txt), or any other harmful behavior, our system engineers are empowered to log, track, identify, publish, and ban misbehaving individuals - even if this leads to ban entire blocks of IP addresses, or disclosing user's identity.

FreeStatistics.org is powered by