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

2016 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(bike2016)
##     1   2   3   4   5   6  7  8  9 10
## 1 292   0   8  30   0  10  5  0  3  3
## 2   0 345   0   0  19   9  0  0  5  1
## 3  13   0 125  16   1   3  0  0 48 10
## 4   5   0  17 149   0  34  1  1  6 24
## 5   2  13   2   1 256  13  1  1  1  0
## 6  10   3   8  24  19 209 87 14  6  1
##        StartStation.Id EndStation.Id  Freq
## 148105             191           191 21192
## 599545             785           785 16309
## 232849             303           303 13051
## 236017             307           307 12590
## 191665             248           248  9012
## 602713             789           789  6874