R version 3.6.2 (2019-12-12) -- "Dark and Stormy Night" Copyright (C) 2019 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-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. > source('/home/pw/wessanet/cretab') > > > > myrfcuid = '' > > x <- c(47.2,58.7,56.36,50.23,50.44,62.87,56.91,55.03,55.33,52.24,52.5,53.43,61.44,54.88,55.93,58.12,56.13,56.48,58.27,51.99,45.37,51.73,48.84,50.35,53.43,46.68,48.32) > #'GNU S' R Code compiled by R2WASP v. 1.2.327 (Tue, 16 Aug 2016 15:30:03 +0200) > #Author: root > #To cite this work: Wessa P., (2013), Tukey lambda PPCC Plot (v1.0.3) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_tukeylambda.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > # > 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="/home/pw/wessanet/rcomp/tmp/1q4rf1582648847.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 > > 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="/home/pw/wessanet/rcomp/tmp/27twc1582648847.tab") > > try(system("convert /home/pw/wessanet/rcomp/tmp/1q4rf1582648847.ps /home/pw/wessanet/rcomp/tmp/1q4rf1582648847.png",intern=TRUE)) convert-im6.q16: profile 'icc': 'RGB ': RGB color space not permitted on grayscale PNG `/home/pw/wessanet/rcomp/tmp/1q4rf1582648847.png' @ warning/png.c/MagickPNGWarningHandler/1654. character(0) > > proc.time() user system elapsed 0.310 0.075 0.390