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

2015 Bike origin-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:10),]))
}

od(bike2015)
##    1  2  3  4  5  6  7 8 9 10
## 1 33  0  8  5  1  2  0 0 0  0
## 2  0 27  0  0  0  2  0 0 0  0
## 3  2  0 23  1  0  0  0 0 2  0
## 4  3  0  0 19  0  6  0 0 5  0
## 5  0  2  0  1 57  3  0 0 1  3
## 6  0  0  0 10  0 26 19 0 0  0
##        StartStation.Id EndStation.Id Freq
## 135977             191           191 3070
## 215050             307           307 1519
## 212094             303           303 1448
## 286732             407           407 1151
## 285993             406           406  898
## 174405             248           248  844