source('/home/pw/wessanet/cretab')



myrfcuid = 'r0702855'

y <- c(13,0,6,0,0,1,3,6,0,0,5,1,1,0,2,6,1,0,0,0,7,4,1,4,0,4,1,1,0,5,0,3,2,0,0,3,0,0,0,1,0,0,4,0,0,0,0,0,0,4,0,0,0,0,5,20,4,1,0,1,10,3,0,1,0,4,4,16,5,0,1,1,0,0,4,0,0,1,0,2,0,0,0,0,6,16,1,1,3,1,1,35,0,2,0,0,4,1,0,3,8,4,0,0,0,1,0,1,1,0,9,1,1,2,0,0,0,1,3,1,0,8,0,0,5,0,0,1,1,0,2,1,4,1,0,0,0,1,1,2,0,0,3,1,1,1,3,0,1,0,1,1,0,2,2,0,22,4,0,25,0,6,7,10,2,1,0,0,0,0,0,0,0,2,3,7,0,0,6,0,0,0,0,0,0,3,8,0,6,3,3,0,6,1,0,2,0,25,14,0,0,0,1,1,0,0,1,2,10,0,0,0,12,0,0,0,0,0,3,29,0,1,1,5,9,6,1,5,2,5,0,1,1,0,1,3,9,9,3,5,5,0,0,0,0,2,2,0,4,0,0,2,2,4,1,1,10,5,0,0,5,1,4,1,0,1,6,9,4,3,0,0,108,1,2,0,2,7,13,2,0,3,2,1,0,8,0,7,1,0,4,24,4,9,0,2,7,1,0,17,0,0,5,0,0,0,1,2,8,0,0,2,0,1,0,0,0,6,0,0,0,2,0,0,2,4,5,0,4,4,0,0,1,0,2,10,0,1,0,0,0,0,5,1,9,0,0,0,0,0,32,4,0,7,0,12,0,0,4,0,2,3,3,1,6,7,1,0,0,0,2,0,24,6,12,6,10,4,3,4,1,8,0,5,5,6,5,3,6,0,11,3,0,8,0,0,0,1,2,0,3,2,7,5,0,4,0,3,0,0,1,3,0,3,0,1,2,1,1,3,0,0,1,0,0,0,2,3,0,1,1,8,0,4,0,2,3,0,1,0,0,0,6,1,2,3,3,0,0,0,19,0,2,6,0,4,1,0,0,4,1,0,5,0,0,0,0,4,3,0,1,0,0,10,3,3,6,4,1,3,3,0,5,0,0,8,9,10,0,1,0,10,0,0,1,0,3,0,0,0,13,7,19,0,7,25,0,0,5,3,30,0,0,3,0,4,16,0,7,0,0,0,0,6,0,0,0,0,0,0,0,0,2,1,0,0,0,1,15,6,2,0,9,5,0,2,1,1,4,7,6,1,0,6,0,0,0,0,0,9,0,2,1,29,3,0,2,4,6,0,3,3,2,6,0,0,5,0,0,0,4,8,14,0,0,9,0,3,0,0,2,3,7,3,3,2,0,8,6,3,14,6,12,6,4,2,1,0,5,6,0,1,6,0,0,2,3,0,5,1,0,1,1,1,3,0,0,0,3,0,0,2,7,0,0,6,1,3,9,0,2,2,0,9,4,0,8,0,0,1,6,0,9,5,1,38,4,1,0,0,0,1,22,1,3,2,3,2,2,3,0,0,0,0,0,0,0,0,8,0,0,0,2,6,5,0,0,1,0,12,27,1,3,0,0,12,0,0,2,0,9,0,7,1,0,8,2,0,7,0,0,1,2,0,1,3,1,0,0,4,0,0,14,6,0,0,0,0,1,0,15,3,10,0,3,27,1,1,2,5,0,1,0,2,6,3,0,0,1,0,0,0,1,7,0,2,3,3,3,1,0,0,14,7,0,0,4,11,0,5,0,1,0,3,4,0,0,0,6,1,6,2,0,2,0,2,4,11,6,0,0,5,4,1,0,1,0,11,1,0,4,13,0,5,0,0,3,0,5,1,0,0,0,0,4,0,1,0,0,0,0,0,1,0,7,3,4,1,0,4,39,0,0,4,0,2,0,1,0,6,4,0,4,0,2,11,3,2,2,0,3,5,0,11,0,5,0,0,0,0,0,3,0,1,1,6,5,21,0,1,0,1,0,0,0,0,33,4,0,0,0,3,11,0,0,1,0,1,5,0,0,2,1,0,0,5,3,0,1,0,0,2,3,3,0,1,0,0,4,1,14,0,1,0,2,1,5,3,4,2,0,4,7,4,0,0,3,1,0,0,4,6,1,0,0,12,0,8,8,0,0,0,0,4,0,0,0,0,0,5,6,13,5,0,18,7,1,3,1,0,0,0,0,4,2,4,0,0,0,0,0,1,9,2,0,3,0,0,0,1,3,0,0,1,2,0,5,1,0,0,0,2,1,1,3,7,0,0,1,0,0,0,6,2,5,0,5,12,0,2,2,3,1,0,3,6,1,0,0,0,0,0,0,6,1,4,7,1,0,0,1,6,1,0,0,1,5,1,21,28,5,3,3,3,5,5,0,0,0,0,0,0,2,6,0,10,5,1,5,1,1,3,2,0,0,0,3,2,0,0,3,0,0,4,0,2,0,4,0,8,7,5,0,1,2,0,5,0,4,3,2,0,0,10,5,0,17,2,3,11,0,0,7,0,18,4,4,0,0,0,0,7,1,5,0,5,11,2,9,0,0,1,0,15,3,15,8,1,0,3,3,14,4,2,0,0,1,0,1,0,0,5,1,5,4,1,0,0,4,3,8,0,0,0,0,0,8,24,0,0,0,15,5,4,5,0,1,0,0,8,5,3,0,0,0,7,1,0,5,0,10,0,28,3,4,0,0,0,0,2,0,1,3,6,0,0,0,8,6,0,0,3,1,1,0,10,41,3,1,0,0,2,0,0,1,1,3,7,4,1,0,4,7,10,4,1,2,2,0,0,6,0,4,11,0,0,1,0,1,0,11,0,0,0,0,0,0,1,1,1,2,12,1,2,1,2,0,0,0,1,2,7,8,1,7,0,6,0,0,0,16,7,0,18,0,1,15,3,0,0,0,6,3,0,1,0,6,2,2,0,3,0,6,1,0,2,2,0,0,0,3,0,10,2,0,4,7,4,2,1,5,0,0,0,0,0,2,1,5,1,17,1,0,0,0,0,1,0,0,0,0,3,0,6,6,4,0,2,7,5,3,5,0,12,7,0,0,0,0,4,4,6,1,0,7,0,2,0,0,1,0,1,4,1,2,0,5,6,3,0,2,0,5,0,9,0,0,4,0,3,11,2,2,1,0,22,0,0,0,0,1,0,6,1,0,5,0,0,0,8,0,1,2,7,3,31,7,0,7,0,0,4,2,44,0,6,0,4,0,0,0,0,1,5,0,5,6,1,0,3,0,0,0,0,5,4,0,8,0,0,0,2,0,0,11,2,3,0,2,0,11,3,1,5,1,4,7,0,1,5,7,11,4,2,0,4,0,7,0,4,3,0,1,0,2,5,0,0,8,5,1,0,9,7,0,2,0,16,0,0,8,3,1,0,1,16,1,3,6,5,1,1,17,26,0,0,4,7,0,0,5,0,0,0,1,14,0,0,0,2,0,5,2,1,4,19,2,3,0,6,12,2,0,1,0,3,0,0,1,6,0,10,0,0,1,2,6,2,0,11,18,0,0,1,22,0,0,1,0,0,12,43,2,0,4,0,0,0,1,1,7,5,7,4,1,10,11,31,4,1,0,0,1,16,12,3,0,1,0,1,0,0,3,5,4,1,0,5,0,0,0,8,13,2,6,5,0,0,0,4,0,1,0,2,8,3,0,3,0,1,0,0,2,2,8,1,6,3,4,2,0,0,6,4,0,0,29,0,1,1,1,0,0,0,4,0,0,0,4,8,4,16,3,0,1,2,0,0,0,0,6,2,0,3,2,0,1,0,6,0,0,0,1,2,0,0,0,7,1,2,2,1,0,2,0,0,7,12,1,0,3,2,8,9,7,27,0,0,4,0,1,1,1,0,0,0,6,5,2,3,0,0,0,0,7,5,0,0,4,0,0,3,5,0,0,3,0,2,9,6,1,9,7,5,8,0,0,1,1,5,0,5,0,0,3,1,1,0,0,0,3,0,9,0,0,2,0,0,4,0,4,0,6,0,3,0,4,0,1,0,3,1,1,0,9,4,4,3,3,1,1,0,7,0,0,3,2,4,3,2,0,0,1,6,3,1,5,0,0,0,0,10,0,11,0,7,0,6,4,4,0,2,1,0,2,2,4,0,2,10,0,1,6,1,3,0,1,20,3,3,8,1,2,0,1,0,0,1,16,0,3,19,0,3,16,1,1,8,14,40,23,0,0,0,0,2,3,7,5,0,0,2,0,0,0,4,0,0,11,3,3,1,1,1,7,9,4,8,4,1,0,2,0,0,1,1,3,6,1,0,2,0,0,2,1,6,3,0,3,3,0,1,5,0,4,1,5,14,3,0,1,0,0,0,0,0,3,0,8,0,5,13,2,0,0,0,5,0,0,0,0,2,3,0,0,0,2,2,2,4,2,4,0,0,4,1,4,0,0,0,0,0,0,3,0,0,4,2,3,0,2,3,0,2,1,0,8,2,0,0,1,9,4,2,2,0,0,4,0,0,1,0,0,2,6,0,9,5,1,0,5,0,2,0,5,7,5,4,0,1,1,0,0,2,5,5,11,25,1,4,2,1,3,2,8,0,3,8,0,0,0,0,0,0,0,3,0,0,2,13,0,0,5,0,6,0,3,2,1,1,0,0,0,0,0,0,5,23,0,3,5,0,0,0,0,12,1,3,2,0,0,1,0,2,9,0,17,0,0,4,2,0,5,0,2,3,0,1,2,2,6,0,0,0,0,0,14,1,0,2,19,1,4,0,2,2,0,2,3,2,0,0,0,2,0,2,3,2,0,0,0,1,6,0,0,0,1,2,5,0,1,0,0,3,0,0,4,4,14,0,2,8,4,0,8,0,0,0,1,0,4,0,6,0,0,0,0,0,10,9,0,1,40,5,0,5,3,0,1,0,7,11,1,2,0,3,4,0,13,0,0,0,0,0,6,0,2,4,1,0,4,1,1,0,1,0,0,0,4,4,3,1,1,4,2,3,0,0,0,1,9,0,3,5,0,0,1,1,2,2,16,2,3,0,16,5,1,5,24,1,0,0,0,1,0,2,4,2,0,6,1,0,0,3,2,0,0,0,0,1,1,1,2,0,0,4,4,0,18,3,0,0,0,6,4,0,0,0,0,6,0,0,1,1,0,14,7,0,0,0,7,0,0,4,1,0,1,1,0,0,5,0,0,6,0,0,0,3,0,0,4,0,2,0,0,0,0,2,0,0,0,4,10,0,0,10,2,0,0,3,0,0,0,0,5,5,5,0,0,0,0,4,0,1,4,10,2,18,3,0,3,8,0,6,0,0,26,4,0,0,1,1,0,4,3,4,1,0,0,0,0,1,0,6,16,0,0,9,5,0,0,1,6,38,0,1,0,0,1,10,8,1,11,12,6,4,0,0,11,0,0,5,3,19,0,0,5,1,0,0,50,6,0,0,1,1,0,6,14,4,1,7,0,5,0,1,0,0,0,0,0,0,0,0,1,1,0,0,0,11,8,4,5,4,0,6,1,4,0,0,0,0,0,0,1,1,1,0,0,0,7,0,0,0,1,1,2,0,0,0,5,1,3,0,0,2,0,0,1,1,2,1,4,0,8,0,3,0,0,0,0,0,0,0,2,0,0,0,0,0,19,0,3,0,2,1,2,9,0,0,3,0,0,2,0,1,1,0,21,0,5,0,0,8,2,0,0,12,3,0,5,5,5,1,9,6,6,0,0,4,17,0,0,2,1,0,0,4,3,0,0,0,7,2,0,7,0,0,1,1,0,0,1,0,8,0,8,13,4,7,1,9,0,0,2,1,5,2,0,2,2,0,2,4,5,0,5,0,3,0,1,3,11,1,3,5,14,0,1,1,0,0,5,0,1,59,2,2,10,8,5,7,0,0,0,2,1,0,4,2,0,0,0,1,0,3,11,4,2,3,0,1,3,0,7,0,0,0,0,16,0)
x <- c(24,0,10,4,45,2,4,9,3,10,7,40,5,4,8,6,1,2,0,3,44,5,0,14,37,3,4,0,8,11,6,4,26,0,37,5,0,4,39,8,0,3,5,5,45,2,0,15,5,4,7,0,7,0,19,1,3,4,57,4,47,8,42,0,52,5,5,29,5,10,10,7,79,38,5,0,43,3,60,7,79,0,8,3,4,9,4,9,2,0,2,26,40,35,79,0,5,0,10,10,2,4,0,44,41,43,0,6,5,55,11,3,0,5,76,0,0,2,20,0,15,5,18,20,6,23,0,0,14,12,3,0,6,17,47,32,0,0,0,3,0,4,2,0,6,1,8,0,5,5,3,9,45,8,5,4,47,10,32,46,5,5,9,34,0,3,38,0,2,30,3,70,0,5,4,22,54,0,3,3,50,0,50,0,15,9,5,40,5,7,5,0,6,24,0,3,45,45,6,5,44,4,50,30,30,0,10,3,7,50,0,7,7,0,79,0,55,45,2,14,17,3,0,7,30,10,11,7,5,5,4,0,0,15,0,4,6,7,3,3,6,4,3,0,19,7,6,28,4,5,10,3,8,9,6,13,10,9,60,0,30,12,3,4,0,30,7,2,9,0,0,15,65,20,8,22,4,6,24,10,3,9,9,0,4,10,0,8,3,0,8,48,11,9,0,27,10,9,30,2,3,6,9,0,25,0,1,6,5,70,0,5,36,6,79,0,5,4,0,8,26,5,2,37,5,7,4,21,4,3,0,0,4,44,5,4,0,7,3,80,0,5,4,7,26,0,3,1,0,4,48,8,20,25,24,7,17,0,5,0,3,6,2,22,10,0,5,0,2,9,10,20,49,10,4,10,51,9,7,14,3,3,2,2,2,23,2,2,20,0,15,5,7,6,0,2,0,12,28,20,9,2,5,4,4,9,0,7,5,5,2,4,0,8,3,6,5,0,1,14,8,23,7,0,18,108,5,3,2,2,15,4,55,4,0,5,7,2,46,12,0,25,4,3,3,35,50,121,55,0,50,48,4,3,8,10,10,9,0,4,7,10,8,0,45,58,0,8,4,0,7,68,0,3,8,9,5,5,5,6,3,14,3,25,0,10,4,7,17,0,49,43,6,6,8,0,10,3,3,7,6,10,12,15,8,38,4,8,17,5,43,59,0,4,0,7,6,2,4,4,40,0,0,8,0,3,3,13,4,0,0,5,4,8,45,1,5,7,42,6,6,0,24,4,41,4,45,5,10,10,4,3,4,4,0,3,3,60,5,3,0,3,42,47,50,2,7,6,28,0,4,7,4,12,0,43,4,40,12,0,25,4,0,4,0,4,0,8,42,6,7,5,7,25,5,8,2,8,19,7,51,22,9,3,5,8,30,0,8,8,2,2,7,40,58,2,5,3,7,8,3,8,4,6,10,0,2,0,10,40,55,2,6,148,0,4,4,12,8,0,10,6,0,28,9,5,10,0,13,12,8,24,24,15,6,48,7,7,0,0,2,2,22,30,10,10,5,9,6,2,0,0,18,3,0,4,59,0,23,5,4,16,4,14,5,120,0,9,5,9,31,5,7,120,51,23,12,12,12,12,6,15,24,28,0,6,4,0,10,0,6,4,6,4,0,7,5,0,0,4,0,1,22,4,4,4,49,0,10,3,8,5,8,30,8,28,3,6,5,20,59,10,40,8,5,34,3,3,24,0,9,0,7,6,5,4,4,21,5,0,0,0,10,6,60,0,5,3,47,7,4,1,60,19,4,8,7,0,4,0,5,9,1,9,2,45,9,5,7,0,41,10,5,0,8,7,8,20,2,8,5,22,65,10,8,0,5,7,5,4,9,24,0,4,4,0,4,2,2,10,0,3,6,6,15,8,6,5,5,7,44,10,15,8,8,3,35,3,0,32,0,25,5,4,3,10,8,6,6,0,3,4,3,5,0,7,15,75,0,56,30,7,43,4,25,1,4,35,39,18,0,12,0,7,4,0,46,3,5,80,0,8,31,5,0,32,7,8,6,1,6,2,47,0,20,3,4,7,3,3,0,7,3,4,35,24,34,30,5,0,23,4,6,49,0,9,6,8,5,9,1,5,9,7,3,0,2,8,3,108,7,10,2,8,30,9,20,23,9,3,40,9,0,8,1,60,2,1,4,24,4,13,5,100,5,2,46,5,37,5,3,33,25,30,7,4,3,7,35,0,2,3,25,9,0,15,22,17,0,5,59,0,38,7,13,8,7,10,0,3,0,4,5,20,7,0,6,7,1,40,69,0,9,6,7,5,8,20,54,6,3,8,5,0,11,7,4,125,3,16,10,25,7,8,4,8,6,7,2,2,7,10,5,1,28,35,3,10,48,15,6,6,9,9,9,2,4,0,60,0,2,6,1,7,1,25,4,0,8,0,0,4,4,5,55,0,2,10,60,0,10,0,0,8,0,6,5,9,49,10,5,3,8,3,3,0,48,0,8,7,7,0,3,4,5,0,45,5,60,3,3,0,6,5,13,15,7,0,0,5,22,18,3,8,0,7,3,5,9,5,5,8,0,0,10,16,20,29,0,9,6,36,4,0,0,5,5,5,5,0,5,23,5,2,8,4,7,5,10,5,7,0,0,0,5,50,25,47,25,0,52,4,7,9,12,0,40,1,0,0,5,5,5,3,0,7,0,1,9,0,7,24,36,8,5,69,0,70,0,4,47,13,0,9,0,11,19,7,6,0,2,3,5,5,0,37,60,2,0,39,7,10,0,2,5,0,5,13,4,5,35,2,0,4,10,4,0,0,45,45,7,10,9,25,9,8,17,0,14,1,24,6,14,0,5,2,3,10,0,20,8,23,10,10,6,0,50,52,23,4,4,25,10,5,4,1,8,0,70,0,16,5,17,22,62,25,3,4,37,5,1,3,10,6,5,0,5,5,5,42,3,10,6,25,0,5,0,0,5,29,4,0,4,10,47,5,7,6,7,11,10,79,0,6,0,12,7,7,8,6,25,8,3,0,2,4,7,70,0,78,0,8,1,5,10,7,35,5,10,5,10,3,0,42,5,0,5,40,0,9,7,10,6,0,7,0,8,3,10,0,10,7,2,2,6,0,0,8,4,50,32,13,7,15,4,6,0,4,0,8,18,10,10,2,2,33,60,0,79,0,5,45,17,2,0,4,0,5,3,4,0,8,3,7,6,39,7,49,7,60,40,3,4,46,7,9,25,10,0,0,0,56,10,5,42,5,6,0,90,8,0,38,0,3,5,7,0,8,6,0,5,7,0,50,21,25,1,30,5,0,17,5,25,2,4,5,5,110,4,9,17,19,6,29,61,2,55,10,10,4,3,0,3,4,4,10,20,5,5,8,3,20,10,4,50,9,80,42,1,3,2,5,3,0,40,5,18,8,9,9,3,0,0,0,0,0,7,6,0,0,10,0,5,1,1,23,8,8,3,2,4,2,6,0,8,49,6,4,39,6,50,12,12,2,3,8,57,31,0,8,0,7,10,12,5,4,9,4,1,9,33,76,0,4,10,0,4,23,2,4,9,46,13,28,5,0,5,5,5,20,9,5,10,2,3,9,3,39,6,0,0,4,1,2,21,9,15,4,10,10,2,2,10,4,10,6,35,41,15,0,18,19,20,3,5,7,46,0,5,10,3,77,0,7,2,10,0,6,0,10,0,6,7,3,9,10,8,7,4,0,44,37,10,5,4,0,43,78,0,0,12,28,6,0,3,5,7,2,2,5,5,21,4,0,12,0,10,5,8,0,2,8,5,3,7,0,2,0,5,0,24,44,4,7,42,75,0,10,5,7,10,14,0,2,10,0,7,10,47,2,8,6,7,46,9,10,75,0,10,3,14,8,5,3,3,45,7,10,6,3,0,8,75,0,8,7,0,43,0,14,4,4,5,0,3,6,4,4,30,4,4,7,7,5,4,0,3,2,8,3,20,5,8,8,7,0,4,5,45,59,24,0,5,63,0,0,2,0,9,0,6,150,10,2,5,10,8,0,48,50,10,0,10,0,8,0,9,6,0,5,4,5,8,0,0,8,6,2,2,0,35,0,4,6,4,10,5,6,0,7,3,5,12,5,0,10,39,1,9,1,0,4,30,9,30,4,4,50,5,1,6,4,4,4,2,4,12,25,24,4,8,13,10,0,3,0,5,30,2,4,7,28,0,25,0,13,5,10,4,49,48,25,0,50,0,7,12,26,8,66,0,4,3,7,78,6,9,0,7,5,4,8,3,10,9,54,9,7,5,6,0,8,10,10,5,8,8,25,0,5,10,7,7,7,8,8,2,15,9,7,2,4,4,0,4,5,47,8,4,0,0,4,2,0,2,8,3,3,6,38,8,25,7,2,0,2,5,7,0,18,5,2,10,48,7,4,3,5,6,1,4,4,0,0,2,6,7,6,0,5,45,77,0,6,50,3,10,5,6,0,4,5,5,38,5,3,5,7,5,29,10,1,4,8,15,11,4,13,0,20,10,35,66,5,6,0,14,29,2,4,5,27,11,0,10,7,10,7,9,5,0,5,55,4,5,6,9,52,3,8,9,0,6,19,10,17,16,13,40,8,0,4,0,8,3,5,37,56,21,0,0,1,8,42,7,5,0,6,8,1,1,4,23,71,0,6,40,28,5,5,6,0,0,1,4,9,8,4,0,79,0,22,5,4,6,16,5,0,125,7,48,0,3,14,0,8,3,4,4,5,2,110,40,0,5,1,1,28,45,4,31,10,6,2,5,9,0,46,3,0,0,4,2,4,0,57,4,4,35,60,0,2,5,0,40,3,3,7,8,60,25,0,0,4,30,0,5,18,6,0,12,5,4,7,19,67,69,0,2,2,3,50,9,0,4,58,0,2,18,9,5,2,53,5,0,5,4,0,4,0,11,0,49,13,0,10,6,0,32,0,10,2,13,1,9,45,3,6,6,5,2,5,9,0,4,10,5,10,3,1,0,6,0,5,3,6,5,0,7,3,10,40,5,25,8,7,8,8,7,4,5,7,9,0,47,5,20,3,34,5,4,0,48,0,5,5,17,4,0,4,5,3,0,5,7,36,6,0,3,5,5,9,5,0,53,4,12,5,33,10,75,75,75,0,4,4,1,5,3,3,0,4,4,4,0,4,44,2,48,28,14,0,2,10,0,37,5,2,1,2,7,0,0,7,0,5,0,6,0,30,5,10,0,5,5,7,0,3,0,30,0,3,7,0,10,5,4,0,4,5,0,0,16,0,6,7,4,0,3,4,0,3,0,3,6,22,8,5,3,0,6,4,4,7,7,5,41,4,0,8,3,3,10,10,5,1,4,50,0,11,43,20,3,3,4,2,5,23,6,5,43,5,7,48,4,3,75,40,10,8,34,0,7,23,9,5,0,0,14,9,7,5,3,52,9,0,43,7,0,6,41,9,1,135,10,0,45,25,5,4,4,10,0,5,0,5,20,0,60,0,30,0,50,0,42,42,22,12,0,9,7,20,3,9,0,4,5,10,79,30,60,0,19,0,40,8,0,7,7,15,2,10,0,50,1,1,50,2,0,40,4,8,15,0,0,10,75,6,5,2,2,7,9,0,4,3,21,0,35,25,0,60,52,5,5,4,0,125,20,0,4,0,5,45,5,3,0,18,40,32,7,0,66,8,0,24,8,26,27,28,11,5,9,9,5,20,0,7,6,4,10,8,8,0,26,4,4,0,0,7,40,5,0,5,41,54,0,3,3,48,0,34,10,3,0,10,0,8,5,3,0,5,8,17,9,0,8,21,8,9,21,13,0,8,5,0,3,14,0,7,8,40,8,6,10,25,1,44,2,60,2,3,21,7,4,0,20,0,60,4,70,0,6,0,10,48,2,5,14,2,3,8,3,0,7,5,5,0,8,10,10,2,0,2,8,5,9,7,2,10,0,6,4,10,8,4,4,4,18,24,39)
ylab = 'cars derailed'
xlab = 'speed'
#'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/1w5zt1513423549.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/2d2mv1513423549.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/3nwea1513423549.tab") 

try(system("convert /home/pw/wessanet/rcomp/tmp/1w5zt1513423549.ps /home/pw/wessanet/rcomp/tmp/1w5zt1513423549.png",intern=TRUE))
try(system("convert /home/pw/wessanet/rcomp/tmp/2d2mv1513423549.ps /home/pw/wessanet/rcomp/tmp/2d2mv1513423549.png",intern=TRUE))
