R, неправильный вывод из ggplot2 и facet_wrap при использовании пространственных данных

Я пытаюсь создать визуализацию Великобритании для разных периодов времени, где заполнение для каждого региона зависит от числовой переменной. Когда я разделяю кадр данных по периоду времени и создаю диаграммы независимо в ggplot2, они верны. Однако, когда я использую функцию facet_wrap, чтобы они отображались на одной диаграмме, результаты неверны (числа для каждого региона неверны).

(Обратите внимание, я старался максимально точно следовать этому руководству по пространственным данным в R: http://spatial.ly/2013/12/introduction-spatial-data-ggplot2/)

Команда, которую я пытаюсь использовать:

Map <- ggplot(data, aes(long, lat, group = group, fill = data$numericvariable)) 
+ geom_polygon() + coord_equal() + facet_wrap(~Date.ordered) + 
scale_fill_gradientn(limits = c(55,70),colours=c("white","red"))

Данные представляют собой пространственные данные, которые были преобразованы в фрейм данных с помощью команды fortify.

Вывод работает независимо, когда я использую следующие команды:

data.split<-split(data,data$Date.ordered)

Maps<-list()

for (i in 1:length(data.split){

    Maps[[i]]<- ggplot(data.split[[i]], aes(long, lat, group = group, 
    fill = data.split[[i]]$numericvariable)) + geom_polygon() + coord_equal() + 
    scale_fill_gradientn(limits = c(55,70),colours=c("white","red"))

    print(Maps[[i]])
}

Но это очень сложно красиво уместить в одну визуализацию, даже используя мультиплот.

Поскольку визуализации работают правильно, когда они выполняются независимо, я в настоящее время предполагаю, что это ошибка кодирования от моего имени, а не проблема с данными, но я могу опубликовать пример того, как выглядят данные, если это поможет.

Спасибо.