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