R version 2.13.0 (2011-04-13) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- c(293403 + ,277108 + ,264020 + ,260646 + ,246100 + ,244051 + ,241329 + ,234730 + ,234509 + ,233482 + ,233406 + ,228548 + ,223914 + ,223696 + ,223004 + ,213765 + ,210554 + ,202204 + ,199512 + ,195304 + ,191467 + ,191381 + ,191276 + ,190410 + ,188967 + ,188780 + ,185139 + ,185039 + ,184217 + ,181853 + ,181379 + ,181344 + ,179562 + ,178863 + ,178140 + ,176789 + ,176460 + ,175877 + ,175568 + ,174107 + ,173587 + ,173260 + ,172684 + ,167845 + ,167131 + ,167105 + ,166790 + ,164767 + ,162810 + ,162336 + ,161678 + ,158980 + ,157250 + ,156833 + ,155383 + ,154991 + ,154730 + ,151503 + ,146455 + ,143937 + ,142339 + ,142146 + ,142141 + ,142069 + ,141933 + ,139350 + ,139144 + ,137793 + ,136911 + ,136548 + ,135171 + ,134043 + ,131876 + ,131122 + ,130539 + ,130533 + ,130232 + ,129100 + ,128655 + ,128066 + ,127619 + ,127324 + ,126683 + ,126681 + ,125971 + ,125366 + ,122433 + ,121135 + ,119291 + ,118958 + ,118807 + ,118372 + ,116900 + ,116775 + ,115199 + ,114928 + ,114397 + ,113337 + ,111664 + ,108715 + ,107342 + ,107335 + ,106539 + ,105615 + ,105410 + ,105324 + ,103012 + ,102531 + ,101324 + ,100885 + ,100672 + ,99946 + ,99768 + ,99246 + ,98599 + ,98030 + ,94763 + ,93340 + ,93125 + ,91185 + ,90961 + ,90938 + ,89318 + ,88817 + ,84944 + ,84572 + ,84256 + ,80953 + ,78800 + ,78776 + ,75812 + ,75426 + ,74398 + ,74112 + ,73567 + ,69471 + ,68948 + ,67746 + ,67507 + ,65029 + ,64320 + ,61857 + ,61499 + ,50999 + ,46660 + ,43287 + ,38214 + ,35523 + ,32750 + ,31414 + ,24188 + ,22938 + ,21054 + ,17547 + ,14688 + ,7199 + ,969 + ,455 + ,203 + ,98 + ,0 + ,0 + ,0 + ,0) > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > gp <- function(lambda, p) + { + (p^lambda-(1-p)^lambda)/lambda + } > sortx <- sort(x) > c <- array(NA,dim=c(201)) > for (i in 1:201) + { + if (i != 101) c[i] <- cor(gp(ppoints(x), lambda=(i-101)/100),sortx) + } > postscript(file="/var/wessaorg/rcomp/tmp/146a71346340198.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot((-100:100)/100,c[1:201],xlab='lambda',ylab='correlation',main='PPCC Plot - Tukey lambda') > grid() > dev.off() null device 1 > > #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/wessaorg/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Tukey Lambda - Key Values',2,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Distribution (lambda)',1,TRUE) > a<-table.element(a,'Correlation',1,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Approx. Cauchy (lambda=-1)',header=TRUE) > a<-table.element(a,c[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Exact Logistic (lambda=0)',header=TRUE) > a<-table.element(a,(c[100]+c[102])/2) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Approx. Normal (lambda=0.14)',header=TRUE) > a<-table.element(a,c[115]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'U-shaped (lambda=0.5)',header=TRUE) > a<-table.element(a,c[151]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Exactly Uniform (lambda=1)',header=TRUE) > a<-table.element(a,c[201]) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/wessaorg/rcomp/tmp/2iacx1346340198.tab") > > try(system("convert tmp/146a71346340198.ps tmp/146a71346340198.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.094 0.184 1.265