source('/home/pw/wessanet/cretab')



myrfcuid = 'r0481972'

x <- c(14,19,17,17,15,20,15,19,15,15,19,16,20,18,15,14,20,16,16,16,10,19,19,16,15,18,17,19,17,14,19,20,5,19,16,15,16,18,16,15,17,14,20,19,7,13,16,16,18,18,16,17,19,16,19,13,16,13,12,17,17,17,16,16,14,16,13,16,14,20,12,13,18,14,19,18,14,18,19,15,14,17,19,13,19,18,20,15,15,15,20,15,19,18,18,15,20,17,18,19,20,13,17,15,16,18,18,14,15,12,17,14,18,17,17,20,16,14,15,18,20,17,17,17,17,15,17,18,17,20,15,16,15,18,15,18,20,19,14,16,15,17,18,20,17,18,15,16,11,15,18,17,16,12,19,18,15,17,19,18,19,16,16,16,14)
#'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/14yv91482345427.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/2f5ey1482345427.tab") 

try(system("convert /home/pw/wessanet/rcomp/tmp/14yv91482345427.ps /home/pw/wessanet/rcomp/tmp/14yv91482345427.png",intern=TRUE))