本实验初步探索6月数据
options(stringsAsFactors = F) rm(list=ls()) source("D:\\R\\packages\\Mreport\\scripts\\caculate.R", encoding = "utf-8") source("D:\\R\\packages\\Mreport\\scripts\\select.R", encoding = "utf-8")
library(Mreport) library(plyr) library(ggplot2) library(reshape2)
load_base() load_sample_base()
jd201806 <- read.csv("D:\\data\\sx_raw\\交调数据\\jd2018_06.csv") dim(jd201806)
当量计算,筛选属性,合并线路,合并样本库
jd201806 <- caculate_equivalent(jd201806) jd201806 <- select_atts(jd201806) jd201806 <- handle_mergeline(jd201806,station_line) jd201806 <- handle_mergesample(jd201806,sample_base) jd201806 <- merge(jd201806,roadlevel,by="index",all.x = T) dim(jd201806)
筛选交调站
jd201806s <- subset(jd201806,index %in% station_use) dim(jd201806s)
jd201805 <- read.csv("D:\\data\\sx_raw\\交调数据\\jd2018_05.csv") dim(jd201805)
当量计算,筛选属性,合并线路,合并样本库
jd201805 <- caculate_equivalent(jd201805) jd201805 <- select_atts(jd201805) jd201805 <- handle_mergeline(jd201805,station_line) jd201805 <- handle_mergesample(jd201805,sample_base) jd201805 <- merge(jd201805,roadlevel,by="index",all.x = T) dim(jd201805)
筛选交调站
jd201805s <- subset(jd201805,index %in% station_use) dim(jd201805s)
citygroup201806 <- caculate_level_carsmean(jd201806s,"citygroup2") citygroup201806
caculate_carsmean(jd201806s,"citygroup2") %>% gg_boxplot(xangle = 0,xlabname = "城市群",ylabname="6月机动车平均交通量")
citygroup201805 <- caculate_level_carsmean(jd201805s,"citygroup2") citygroup201805
同比
caculate_yearratio(citygroup201806,citygroup201805[-5,])
caculate_carsmean(jd201806s,"horizon10")
caculate_carsmean(jd201806s,"vertical10")
(portroad201806 <- caculate_frecarsmean(jd201806s,"portroad")) %>% tail(10)
(portroad201805 <- caculate_frecarsmean(jd201805s,"portroad")) %>% tail(10)
caculate_yearratio(portroad201806,portroad201805)
(bigcityio201805 <- caculate_frecarsmean(jd201805s,"bigcityio")) %>% tail(10)
(bigcityio201806 <- caculate_frecarsmean(jd201806s,"bigcityio")) %>% tail(10)
caculate_yearratio(bigcityio201806,bigcityio201805) %>% tail(10)
(airport201806 <- caculate_carsmean(jd201806s,"airport")) %>% tail(10)
(airport201805 <- caculate_carsmean(jd201805s,"airport")) %>% tail(10)
caculate_yearratio(airport201806,airport201805)
(scenery201806 <- caculate_passcarsmean(jd201806s,"scenery")) %>% tail(10)
(scenery201805 <- caculate_passcarsmean(jd201805s,"scenery")) %>% tail(10)
caculate_yearratio(scenery201806,scenery201805)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.