library(readr)
bike2013<-read.csv("bike2013.csv",sep = ",")

2013 Bike orgin-destination matrix and relative frequencies

od<- function(data){
  odmatrix<-as.data.frame.matrix(table(data$StartStation.Id,data$EndStation.Id))
  frq<-as.data.frame(table(data$StartStation.Id,data$EndStation.Id))
  frq<-frq[with(frq, order(Freq,decreasing = T)),] 
  colnames(frq)<- c("StartStation.Id","EndStation.Id","Freq")
  print(head(odmatrix[c(1:10),c(1:10)]))
  print(head(frq[c(1:20),]))
}
od(bike2013)
##     0   1   2   3   4   5   6  7 8  9
## 1 123 343   0  45  32   0   5  0 1  5
## 2 162   0 295   0   0  13   2  2 1  2
## 3 178  18   0 112  13   2   6  1 0 37
## 4 213   7   0   4 148   1  40  5 1 20
## 5 208   1  12   0   3 293   6  4 0  4
## 6 243  11   2   7  39   9 235 53 8  1
##        StartStation.Id EndStation.Id  Freq
## 138348             191           191 18385
## 221796             307           307 10561
## 218868             303           303  9313
## 293531             407           407  6638
## 180072             248           248  5906
## 292799             406           406  5521