R version 2.8.1 (2008-12-22) Copyright (C) 2008 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(4.154454406,4.146531122,4.20172477,4.155032229,4.168615322,4.167937313,4.143576832,4.12264177,4.194791758,4.193931191,4.179350682,4.16214604,4.179896333,4.173419384,4.208629438,4.149896207,4.145662471,4.153387979,4.143670433,4.093211492,4.189040908,4.208064948,4.179178229,4.164144582,4.16716959,4.161757195,4.232258962,4.181243161,4.18920949,4.23484617,4.202433822,4.176178109,4.260453018,4.251565224,4.255971639,4.23756887,4.222872477,4.224014811,4.276737459,4.249516316,4.223651667,4.256741793,4.223054413,4.224507165,4.290346301,4.284453294,4.299964669,4.279940573,4.29125782,4.28064665,4.35283829,4.275403503,4.296687124,4.30847904,4.259331025,4.280737674,4.318939392,4.334493962,4.349471799,4.308222819,4.326643036,4.312388949,4.364588514,4.303930064,4.326192438,4.335016361,4.318543299,4.31160554,4.344510166,4.383420373,4.39720978,4.33485569) > x <- c(3.971090213,3.970207359,4.006423253,3.990693961,3.989894564,3.996117476,3.974465636,3.943247125,4.008131562,4.022387126,3.996248915,3.990782692,3.983400738,3.988068203,4.024854948,3.982768577,3.972573081,3.978180517,3.986323777,3.926033597,3.993171605,4.025346835,3.997473759,4.009833179,3.972249136,3.983851719,4.050302404,3.999956568,4.003848122,4.065280871,4.022552084,3.998477303,4.063220736,4.05388479,4.067665882,4.069926969,4.025100961,4.031448862,4.086324231,4.067814511,4.040523226,4.070148736,4.034989216,4.021933166,4.085754244,4.077985291,4.092123869,4.091983332,4.085004999,4.082641778,4.156094631,4.102776615,4.141919874,4.132355762,4.104964846,4.09565744,4.126748142,4.144916527,4.148849343,4.118694508,4.127525925,4.132867788,4.187182178,4.136974096,4.166992231,4.172690482,4.145848757,4.141575183,4.157275396,4.195124423,4.201506367,4.146624089) > 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/rcomp/tmp/1svla1256664353.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/rcomp/tmp/2e1l91256664353.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/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/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/rcomp/tmp/3sik71256664354.tab") > > system("convert tmp/1svla1256664353.ps tmp/1svla1256664353.png") > system("convert tmp/2e1l91256664353.ps tmp/2e1l91256664353.png") > > > proc.time() user system elapsed 0.630 0.470 1.038