R version 3.4.0 (2017-04-21) -- "You Stupid Darkness" Copyright (C) 2017 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. > y <- c(0.007142857,0.008695652,0.008450704,0.008450704,0.008219178,0.008275862,0.008053691,0.010457516,0.009302326,0.00952381,0.009876543,0.01863354,0.01875,0.019379845,0.021126761,0.021276596,0.021337127,0.028409091,0.028368794,0.028368794,0.028368794,0.02739726,0.028169014,0.032531825,0.032624113,0.035360679,0.035211268,0.035310734,0.03943662,0.036781609,0.03902439,0.041131105,0.035874439,0.025062657,0.015026296,0.010883169,0.015076459,0.017193691,0.014276847,0.012734262,0.006318083,0.00521096,0.008375,0.007494407,0.008815789,0.010876623,0.010244648,0.008569967,0.00688592,0.008167744,0.008354115,0.009839461,0.012034546,0.012762763,0.010848756,0.011089367,0.010981912,0.011517615,0.01654764,0.017838587,0.017743325,0.018888685,0.018625678,0.015025529,0.01234272,0.011822773,0.010038986,0.008101306,0.006156605,0.006726071,0.006022989,0.005104735,0.004735467,0.004356971,0.006019095,0.006045445,0.006839623,0.005410448,0.005809295) > x <- c(0.25,0.3,0.3,0.3,0.3,0.3,0.3,0.4,0.4,0.4,0.4,0.75,0.75,0.75,0.75,0.75,0.75,1,1,1,1,1,1,1.15,1.15,1.25,1.25,1.25,1.4,1.6,1.6,1.6,1.6,1.6,1.6,2,2.1,2.3,2.3,2.65,2.9,3.1,3.35,3.35,3.35,3.35,3.35,3.35,3.35,3.35,3.35,3.8,4.25,4.25,4.25,4.25,4.25,4.25,4.75,5.15,5.15,5.15,5.15,5.15,5.15,5.15,5.15,5.15,5.15,5.85,6.55,7.25,7.25,7.25,7.25,7.25,7.25,7.25,7.25) > ylab = 'Minium Wage Buying Power' > xlab = 'Nominal Minimum Wage' > #'GNU S' R Code compiled by R2WASP v. 1.2.291 () > #Author: root > #To cite this work: Wessa, (2012), Kendall tau Rank Correlation (v1.0.11) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_kendall.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > # > library('Kendall') > k <- Kendall(x,y) > postscript(file="/var/wessaorg/rcomp/tmp/1htpf1499868251.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=8.3333333333333) > > plot(x,y,main='Scatterplot',xlab=xlab,ylab=ylab) > grid() > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/2v4fg1499868251.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=8.3333333333333) > > plot(rank(x),rank(y),main='Scatterplot of Ranks',xlab=xlab,ylab=ylab) > grid() > dev.off() null device 1 > > #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,'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="/var/wessaorg/rcomp/tmp/3omv51499868251.tab") > > try(system("convert tmp/1htpf1499868251.ps tmp/1htpf1499868251.png",intern=TRUE)) character(0) > try(system("convert tmp/2v4fg1499868251.ps tmp/2v4fg1499868251.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.722 0.185 1.929