x <- c(6
,4
,6
,4
,2
,5
,1
,6
,3
,4
,3
,2
,4
,3
,5
,3
,3
,6
,1
,1
,1
,5
,2
,3
,3
,5
,5
,1
,2
,3
,5
,3
,5
,3
,2
,4
,3
,3
,6
,6
,3
,3
,4
,4
,2
,2
,3
,2
,5
,5
,6
,3
,3
,2
,5
,2
,4
,4
,1
,6
,2
,3
,4
,3
,4
,3
,6
,4
,5
,3
,4
,3
,3
,6
,6
,2
,2
,6
,5
,1
,2
,5
,3
,3
,4
,6
,4
,2
,4
,4
,2
,4
,3
,6
,4
,4
,3
,3
,3
,4
,2
,2
,3
,3
,4
,4
,2
,5
,5
,4
,4
,3
,4
,4
,5
,4
,5
,3
,2
,1
,5
,5
,3
,4
,4
,5
,5
,3
,4
,6
,4
,5
,4
,4
,5
,4
,3
,3
,6
,5
,2
,4
,3
,3
,5
,3
,3
,4
,5
,5
,5
,5
,3
,4
,4
,3
,4
,4
,3
,2
,4
,4
,6
,4
,2
,2
,2
,5
,1
,4
,5
,4
,1
,2
,4
,2
,6
,3
,2
,6
,2
,6
,1
,4
,1
,4
,2
,5
,5
,1
,1
,4
,4
,5
,5
,3
,1
,5
,2
,3
,6
,4
,4
,2
,3
,5
,5
,2
,2
,1
,4
,1
,4
,4
,5
,4
,5
,2
,3
,2
,1
,5
,2
,3
,4
,2
,5
,2
,6
,5
,4
,2
,4
,5
,4
,6
,6
,4
,4
,3
,6
,2
,2
,2
,3
,2
,4
,5
,4
,3
,4
,3
,2
,5
,3
,5
,3
,3
,4
,3
,3
,6
,3
,1
,5
,6
,4
,4
,5
,4
,4
,3
,1
,2
,3
,2
,5
,2
,2
,4
,2
,3
,4
,3
,1
,5
,1
,4
,5
,3
,4
,3
,3
,6
,2
,4
,5
,4
,4
,3
,4
,5
,3
,5
,2
,6
,6
,2
,3
,4
,4
,4
,4
,3
,7
,4
,4
,4
,4
,4
,3
,5
,3
,5
,5
,5
,5
,6
,6
,5
,5
,6
,7
,6
,7
,6
,6
,6
,6
,6
,7
,7
,7
,6
,5
,7
,5
,7
,6
,2
,2
,7
,4
,4
,5
,5
,4
,7
,7
,4
,7
,4
,6
,6
,6
,6
,6
,7
,6
,6
,4
,7
,6
,7
,6
,7
,7
,2
,6
,6
,5
,6
,6
,6
,6
,6
,4
,3
,4
,7
,6
,6
,7
,7
,7
,5
,6
,3
,6
,5
,5
,6
,6
,4
,4
,5
,6
,5
,5
,6
,7
,5
,7
,7
,2
,5
,4
,7
,6
,7
,5
,5
,4
,7
,5
,4
,6
,6
,5
,5
,3
,5
,6
,6
,7
,7
,5
,5
,6
,4
,7
,5
,6
,7
,7
,6
,6
,4
,7
,6
,5
,6
,5
,5
,6
,6
,4
,6
,7
,6
,7
,7
,6
,5
,5
,5
,7
,3
,6
,5
,7
,4
,6
,4
,6
,7
,4
,6
,7
,5
,6
,4
,4
,6
,4
,3
,5
,4
,5
,6
,5
,7
,7
,6
,7
,6
,7
,5
,3
,2
,7
,5
,5
,6
,7
,6
,5
,7
,5
,5
,6
,7
,7
,2
,6
,6
,7
,6
,6
,4
,6
,6
,7
,5
,7
,6
,5
,6
,7
,6
,6
,6
,3
,6
,4
,6
,4
,3
,7
,4
,7
,5
,5
,4
,6
,6
,7
,5
,4
,7
,4
,7
,6
,6
,6
,5
,5
,6
,5
,6
,6
,5
,5
,6
,7
,5
,7
,2
,5
,6
,7
,7
,5
,5
,1
,6
,5
,6
,7
,6
,6
,5
,6
,5
,6
,6
,6
,5
,5
,6
,6
,5
,6
,6
,1
,7
,7
,6
,4
,5
,6
,3
,7
,6
,6
,6
,5
,6
,6
,3
,7
,5
,6
,5
,5
,6
,5
,6
,5
,3
,7
,3
,7
,5
,6
,6
,4
,6
,7
,4
,6
,5
,6
,4
,5
,5
,5
,6
,7
,7
,6
,6
,5
,7
,7
,6
,4
,5)
par20 = ''
par19 = ''
par18 = ''
par17 = ''
par16 = ''
par15 = ''
par14 = ''
par13 = ''
par12 = ''
par11 = ''
par10 = ''
par9 = ''
par8 = ''
par7 = ''
par6 = ''
par5 = ''
par4 = ''
par3 = ''
par2 = ''
par1 = 'grey'
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/168lq1350395907.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()

#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/2omiz1350395907.tab") 

try(system("convert tmp/168lq1350395907.ps tmp/168lq1350395907.png",intern=TRUE))