R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 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.005146,-0.016437,0.006005,0.008702,0.006065,0.016797,0.020605,0.048558,0.041351,0.028252,0.051707,0.037076,0.046126,0.038562,0.043101,0.040678,0.010128,0.028543,0.010615,-0.000725,0.020013,0.009351,0.004490,0.023993,0.009841,0.011276,0.004274,0.009074,0.014261,0.011233,0.013814,0.031786,0.020695,-0.006130,0.024601,0.025193,-0.021078,-0.016113,-0.014245,0.009134,-0.047135,-0.022520,-0.009361,-0.021052,-0.027148,-0.041973,-0.052227,-0.071104,-0.052942,-0.067569,-0.068452,-0.049375,-0.037657,-0.020642,-0.011183,-0.016302,-0.010016,-0.012496,-0.602865) > x <- c(0.461486645,-0.004084972,-0.022949858,0.314235793,0.719086015,0.519010603,0.214639302,0.724751373,0.947471981,0.942293662,1.06557672,0.290473814,0.466120842,0.844038098,0.614769703,0.262528017,-0.351824955,0.049966352,0.053616259,-0.321937027,-0.471965568,-0.57493886,-0.966618095,-0.889742211,-0.735032602,-0.678409402,-0.695641924,-0.471580389,-0.516139175,-0.668427732,-1.285444229,-1.261260411,-1.324871567,-1.625105921,-1.023530639,-0.655808986,-0.185705033,0.171733303,0.570691931,0.710504447,1.116546865,1.045570694,2.012784886,2.471584571,2.54198808,2.135458622,2.675297587,2.315305706,1.513051716,1.066878896,0.546485598,0.446071879,-0.873384,-1.141009819,-2.342678843,-3.455050932,-4.154431398,-3.353127131,-2.63931828) > ylab = 'Y' > xlab = 'X' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > library('Kendall') > k <- Kendall(x,y) > postscript(file="/var/www/html/rcomp/tmp/1wgdt1257933914.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > > plot(x,y,main='Scatterplot',xlab=xlab,ylab=ylab) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2jwau1257933914.ps",horizontal=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() null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/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/www/html/rcomp/tmp/3aovd1257933914.tab") > > system("convert tmp/1wgdt1257933914.ps tmp/1wgdt1257933914.png") > system("convert tmp/2jwau1257933914.ps tmp/2jwau1257933914.png") > > > proc.time() user system elapsed 0.545 0.327 0.654