source('/home/pw/wessanet/cretab') myrfcuid = '' y <- c(29.9,47,15.3,5.4,8,14.9,17,11.9,3.7,5.2) x <- c(5.3,8.48,5.36,7.44,9.83,8.47,5.8,7.15,8.05,8.2) ylab = 'y' xlab = 'x' #'GNU S' R Code compiled by R2WASP v. 1.2.327 (Sat, 15 Jul 2017 00:25:15 +0200) #Author: root #To cite this work: Wessa, (2017), Kendall tau Rank Correlation (v1.0.13) in Free Statistics Software (v$_version), Office for Research Development and Education, URL https://www.wessa.net/rwasp_kendall.wasp/ #Source of accompanying publication: Office for Research, Development, and Education # library('Kendall') k <- Kendall(x,y) postscript(file="/home/pw/wessanet/rcomp/tmp/11f4b1562668237.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) plot(x,y,main='Scatterplot',xlab=xlab,ylab=ylab) grid() dev.off() postscript(file="/home/pw/wessanet/rcomp/tmp/2ec4w1562668237.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) plot(rank(x),rank(y),main='Scatterplot of Ranks',xlab=xlab,ylab=ylab) grid() dev.off() a<-table.start() a<-table.row.start(a) a<-table.element(a,'Kendall tau Rank Correlation',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Kendall tau',header=TRUE) a<-table.element(a,k$tau) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'2-sided p-value',header=TRUE) a<-table.element(a,k$sl) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Score',header=TRUE) a<-table.element(a,k$S) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Var(Score)',header=TRUE) a<-table.element(a,k$varS) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Denominator',header=TRUE) a<-table.element(a,k$D) a<-table.row.end(a) a<-table.end(a) table.save(a,file="/home/pw/wessanet/rcomp/tmp/3vkwi1562668237.tab") try(system("convert /home/pw/wessanet/rcomp/tmp/11f4b1562668237.ps /home/pw/wessanet/rcomp/tmp/11f4b1562668237.png",intern=TRUE)) try(system("convert /home/pw/wessanet/rcomp/tmp/2ec4w1562668237.ps /home/pw/wessanet/rcomp/tmp/2ec4w1562668237.png",intern=TRUE))