diff --git a/DESCRIPTION b/DESCRIPTION
index fb2c1ad6..35c4d35e 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -35,6 +35,6 @@ Imports:
       utils,
       raster,
       pbapply
-RoxygenNote: 7.0.2
+RoxygenNote: 7.1.0
 VignetteBuilder: knitr
 LinkingTo: Rcpp
diff --git a/NAMESPACE b/NAMESPACE
index 00fe7216..5e30404d 100644
--- a/NAMESPACE
+++ b/NAMESPACE
@@ -12,6 +12,7 @@ export(gps_as_sf)
 export(gtfs2gps)
 export(gtfs_shapes_as_sf)
 export(gtfs_stops_as_sf)
+export(merge_gtfs_feeds)
 export(read_gtfs)
 export(remove_invalid)
 export(test_gtfs_freq)
diff --git a/R/merge_gtfs_feeds.R b/R/merge_gtfs_feeds.R
new file mode 100644
index 00000000..3b693a07
--- /dev/null
+++ b/R/merge_gtfs_feeds.R
@@ -0,0 +1,56 @@
+#' @title Merge multiple GTFS feeds into a single one
+#' 
+#' @description Build a single GTFS by joinning together the elements of multiple GTFS feeds.
+#' @param gtfs_list A list of GTFS.zip files.
+#' @return A single list of data.tables, where each index represents the respective GTFS file name.
+#' @export
+#' @examples
+#' 
+#' # get a list of GTFS feeds
+#' spo <- system.file("extdata/saopaulo.zip", package = "gtfs2gps")
+#' poa <- system.file("extdata/poa.zip", package = "gtfs2gps")
+#' gtfs_list <- list(spo, poa)
+#' 
+#' new_gtfs <- merge_gtfs_feeds(gtfs_list)
+#' 
+
+merge_gtfs_feeds <- function(gtfs_list){
+  
+# read all fees separately
+  all_feeds <- lapply(gtfs_list, read_gtfs)
+
+  # separate 1st feed
+  new_gtfs <- all_feeds[[2]]
+
+  # function to extract elements in a series of lists
+  extract_list_element <- function(i, element){ all_feeds[[i]][[element]] }
+
+
+  ## piling up
+  
+  # 1/8 agency
+  new_gtfs$agency <- lapply(X=1:length(all_feeds), FUN = extract_list_element, 'agency') %>% rbindlist(fill=T)
+  
+  # 2/8 routes
+  new_gtfs$routes <- lapply(X=1:length(all_feeds), FUN = extract_list_element, 'routes') %>% rbindlist(fill=T)
+  
+  # 3/8 stops
+  new_gtfs$stops <- lapply(X=1:length(all_feeds), FUN = extract_list_element, 'stops') %>% rbindlist(fill=T)
+  
+  # 4/8 stop_times
+  new_gtfs$stop_times <- lapply(X=1:length(all_feeds), FUN = extract_list_element, 'stop_times') %>% rbindlist(fill=T)
+  
+  # 5/8 shapes
+  new_gtfs$shapes <- lapply(X=1:length(all_feeds), FUN = extract_list_element, 'shapes') %>% rbindlist(fill=T)
+  
+  # 6/8 trips
+  new_gtfs$trips <- lapply(X=1:length(all_feeds), FUN = extract_list_element, 'trips') %>% rbindlist(fill=T)
+  
+  # 7/8 calendar
+  new_gtfs$calendar <- lapply(X=1:length(all_feeds), FUN = extract_list_element, 'calendar') %>% rbindlist(fill=T)
+  
+  # 8/8 frequencies
+  new_gtfs$frequencies <- lapply(X=1:length(all_feeds), FUN = extract_list_element, 'frequencies') %>% rbindlist(fill=T)
+  
+  return(new_gtfs)
+}
diff --git a/man/merge_gtfs_feeds.Rd b/man/merge_gtfs_feeds.Rd
new file mode 100644
index 00000000..3be2db3c
--- /dev/null
+++ b/man/merge_gtfs_feeds.Rd
@@ -0,0 +1,27 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/merge_gtfs_feeds.R
+\name{merge_gtfs_feeds}
+\alias{merge_gtfs_feeds}
+\title{Merge multiple GTFS feeds into a single one}
+\usage{
+merge_gtfs_feeds(gtfs_list)
+}
+\arguments{
+\item{gtfs_list}{A list of GTFS.zip files.}
+}
+\value{
+A single list of data.tables, where each index represents the respective GTFS file name.
+}
+\description{
+Build a single GTFS by joinning together the elements of multiple GTFS feeds.
+}
+\examples{
+
+# get a list of GTFS feeds
+spo <- system.file("extdata/saopaulo.zip", package = "gtfs2gps")
+poa <- system.file("extdata/poa.zip", package = "gtfs2gps")
+gtfs_list <- list(spo, poa)
+
+new_gtfs <- merge_gtfs_feeds(gtfs_list)
+
+}
diff --git a/tests_rafa/merge_gtfs_feeds.R b/tests_rafa/merge_gtfs_feeds.R
new file mode 100644
index 00000000..3b693a07
--- /dev/null
+++ b/tests_rafa/merge_gtfs_feeds.R
@@ -0,0 +1,56 @@
+#' @title Merge multiple GTFS feeds into a single one
+#' 
+#' @description Build a single GTFS by joinning together the elements of multiple GTFS feeds.
+#' @param gtfs_list A list of GTFS.zip files.
+#' @return A single list of data.tables, where each index represents the respective GTFS file name.
+#' @export
+#' @examples
+#' 
+#' # get a list of GTFS feeds
+#' spo <- system.file("extdata/saopaulo.zip", package = "gtfs2gps")
+#' poa <- system.file("extdata/poa.zip", package = "gtfs2gps")
+#' gtfs_list <- list(spo, poa)
+#' 
+#' new_gtfs <- merge_gtfs_feeds(gtfs_list)
+#' 
+
+merge_gtfs_feeds <- function(gtfs_list){
+  
+# read all fees separately
+  all_feeds <- lapply(gtfs_list, read_gtfs)
+
+  # separate 1st feed
+  new_gtfs <- all_feeds[[2]]
+
+  # function to extract elements in a series of lists
+  extract_list_element <- function(i, element){ all_feeds[[i]][[element]] }
+
+
+  ## piling up
+  
+  # 1/8 agency
+  new_gtfs$agency <- lapply(X=1:length(all_feeds), FUN = extract_list_element, 'agency') %>% rbindlist(fill=T)
+  
+  # 2/8 routes
+  new_gtfs$routes <- lapply(X=1:length(all_feeds), FUN = extract_list_element, 'routes') %>% rbindlist(fill=T)
+  
+  # 3/8 stops
+  new_gtfs$stops <- lapply(X=1:length(all_feeds), FUN = extract_list_element, 'stops') %>% rbindlist(fill=T)
+  
+  # 4/8 stop_times
+  new_gtfs$stop_times <- lapply(X=1:length(all_feeds), FUN = extract_list_element, 'stop_times') %>% rbindlist(fill=T)
+  
+  # 5/8 shapes
+  new_gtfs$shapes <- lapply(X=1:length(all_feeds), FUN = extract_list_element, 'shapes') %>% rbindlist(fill=T)
+  
+  # 6/8 trips
+  new_gtfs$trips <- lapply(X=1:length(all_feeds), FUN = extract_list_element, 'trips') %>% rbindlist(fill=T)
+  
+  # 7/8 calendar
+  new_gtfs$calendar <- lapply(X=1:length(all_feeds), FUN = extract_list_element, 'calendar') %>% rbindlist(fill=T)
+  
+  # 8/8 frequencies
+  new_gtfs$frequencies <- lapply(X=1:length(all_feeds), FUN = extract_list_element, 'frequencies') %>% rbindlist(fill=T)
+  
+  return(new_gtfs)
+}