*Unverified author*

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

Title produced by software: Standard Deviation-Mean Plot

Date of computation: Mon, 10 May 2010 16:16:39 +0000

Cite this page as follows:

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2010/May/10/t1273508224aex5a05i6k1miej.htm/, Retrieved Mon, 10 May 2010 18:17:06 +0200

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/2010/May/10/t1273508224aex5a05i6k1miej.htm/}, year = {2010}, } @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 = {2010}, note = {{ISBN} 3-900051-07-0}, url = {http://www.R-project.org}, }

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords:

KDGP2W83

Dataseries X:

» Textbox « » Textfile « » CSV «

15136 16733 20016 17708 18019 19227 22893 23739 21133 22591 26786 29740 15028 17977 20008 21354 19498 22125 25817 28779 20960 22254 27392 29945 16933 17892 20533 23569 22417 22084 26580 27454 24081 23451 28991 31386 16896 20045 23471 21747 25621 23859 25500 30998 24475 23145 29701 34365 17556 22077 25702 22214 26886 23191 27831 35406 23195 25110 30009 36242 18450 21845 26488 22394 28057 25451 24872 33424 24052 28449 33533 37351 19969 21701 26249 24493 24603 26485 30723 34569 26689 26157 32064 38870 21337 19419 23166 28286 24570 24001 33151 24878 26804 28967 33311 40226 20504 23060 23562 27562 23940 24584 34303 25517 23494 29095 32903 34379 16991 21109 23740 25552 21752 20294 29009 25500 24166 26960 31222 38641 14672 17543 25453 32683 22449 22316 27595 25451 25421 25288 32568 35110 16052 22146 21198 19543 22084 23816 29961 26773 26635 26972 30207 38687 16974 21697 24 etc...

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 1 seconds R Server 'Gwilym Jenkins' @ 72.249.127.135 Standard Deviation-Mean Plot Section Mean Standard Deviation Range 1 17398.25 2041.9674132235 4880 2 20969.5 2775.59236920698 5720 3 25062.5 3932.71496551682 8607 4 18591.75 2751.60684873403 6326 5 24054.75 4078.88746065231 9281 6 25137.75 4241.10460257702 8985 7 19731.75 2976.80873587807 6636 8 24633.75 2778.30312421089 5370 9 26977.25 3843.39124342032 7935 10 20539.75 2803.05611490982 6575 11 26494.5 3108.02214706824 7139 12 27921.5 5143.88397354891 11220 13 21887.25 3339.40188007773 8146 14 28328.5 5125.53167973821 12215 15 28639 5824.47382458994 13047 16 22294.25 3295.32865078230 8038 17 27951 3902.78558638997 8552 18 30846.25 5814.89577865904 13299 19 23103 2805.76644311912 6280 20 29095 4457.39924171035 9966 21 30945 5918.79114910018 12713 22 23052 3809.97226586581 8867 23 26650 4349.19862350142 9150 24 32327 5920.46129959482 13422 25 23672 2918.61656725694 7058 26 27086 4854.69429178261 10363 27 29967.75 4856.12845622244 10885 28 21848 3716.44319208568 8561 29 24138.75 3917.95077176832 8715 30 30247.25 6303.247595486 14475 31 22587.75 8128.80333033918 18011 32 24452.75 2546.30141905209 5279 33 29596.75 5007.54320167751 9822 34 19734.75 2680.46281762435 6094 35 25658.5 3460.51407549033 7877 36 30625.25 5610.55780797834 12052 37 21651.75 3301.63912574749 7205 38 25966.25 2947.26668570502 6326 39 29742.25 5638.8717769308 12084 40 22021.5 4305.48506752336 9339 41 27640.25 3166.47715661006 6499 42 30318.75 6185.18225950807 14474 43 21581 5499.66877790533 12671 44 26086 3038.03631753583 6304 Regression: S.E.(k) = alpha + beta * Mean(k) alpha -1419.06648402853 beta 0.219790711439122 S.D. 0.041964086545167 T-STAT 5.2375907480449 p-value 4.90668774041361e-06 Regression: ln S.E.(k) = alpha + beta * ln Mean(k) alpha -5.56541890162681 beta 1.36722762152273 S.D. 0.231824281189568 T-STAT 5.89768946767365 p-value 5.57384172493935e-07 Lambda -0.367227621522732

Charts produced by software:

http://www.freestatistics.org/blog/date/2010/May/10/t1273508224aex5a05i6k1miej/1dqg91273508197.png (open in new window) http://www.freestatistics.org/blog/date/2010/May/10/t1273508224aex5a05i6k1miej/1dqg91273508197.ps (open in new window) http://www.freestatistics.org/blog/date/2010/May/10/t1273508224aex5a05i6k1miej/2ozfu1273508197.png (open in new window) http://www.freestatistics.org/blog/date/2010/May/10/t1273508224aex5a05i6k1miej/2ozfu1273508197.ps (open in new window)

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

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