In ahorawzy/Mreport: Monthly Report of traffic condition

knitr::opts_chunk$set(echo = F,message = F)

本实验探索分日数据的分析框架

环境准备

options(stringsAsFactors = F)
rm(list = ls())
source("D:\\R\\packages\\Mreport\\scripts\\source_toglobal.R", encoding = "utf-8")

library(Mreport)
library(plyr)
library(ggplot2)
library(reshape2)
library(knitr)

load_base()
load_sample_base()

数据准备

2018年端午

jddwnew <- read.csv("D:\\data\\sx_raw\\交调数据\\jd2018_06_DuanWu.csv")
dim(jddwnew)

拆分日期

x <- strsplit(jddwnew$观测日期,split = "-")
y <- sapply(x,'[',3)
jddwnew$日 <- y

jddwnew <- split_day(jddwnew)
jddwnew <- caculate_equivalent(jddwnew)
jddwnew <- select_atts_forday(jddwnew)
jddwnew <- handle_mergeline(jddwnew,station_line)
jddwnew <- handle_mergesample(jddwnew,sample_base)
jddwnew <- merge(jddwnew,roadlevel,by="index",all.x = T)
jddwnew <- subset(jddwnew,index %in% station_use)

2017年端午

jddwprevious <- read.csv("D:\\data\\sx_raw\\交调数据\\jd2017_06_DuanWu.csv")
dim(jddwprevious)

jddwprevious <- split_day(jddwprevious)
jddwprevious <- caculate_equivalent(jddwprevious)
jddwprevious <- select_atts_forday(jddwprevious)
jddwprevious <- handle_mergeline(jddwprevious,station_line)
jddwprevious <- handle_mergesample(jddwprevious,sample_base)
jddwprevious <- merge(jddwprevious,roadlevel,by="index",all.x = T)
jddwprevious <- subset(jddwprevious,index %in% station_use)

1. 总体分析

1.1 分道路等级

(dwallnewlevel <- ddply(jddwnew,"level",summarise,Wmean = weighted.mean(cars,w=mileage))) %>% kable()

(dwallpreviouslevel <- ddply(jddwprevious,"level",summarise,Wmean = weighted.mean(cars,w=mileage))) %>% kable()

caculate_increaseratio(dwallnewlevel,dwallpreviouslevel) %>% kable()

1.2 分日期

(dwallnewday <- ddply(jddwnew,"day",summarise,Wmean = weighted.mean(cars,w=mileage))) %>% kable()

2. 分省分析

dwprovincenew <- caculate_carsmean(jddwnew,"province")

dwprovinceprevious <- caculate_carsmean(jddwprevious,"province")
x1 <- caculate_increaseratio(dwprovincenew,dwprovinceprevious)

merge(dwprovincenew,x1,by="province") %>% kable()

3. 旅游景区

dwscenerynew <- caculate_carsmean(jddwnew,"scenery")

同比

dwsceneryprevious <- caculate_carsmean(jddwprevious,"scenery")
x1 <- caculate_increaseratio(dwscenerynew,dwsceneryprevious)

merge(dwscenerynew,x1,by="scenery") %>% kable()

4. 大城市出入口

dwbigcityionew <- caculate_carsmean(jddwnew,"bigcityio")

dwbigcityioprevious <- caculate_carsmean(jddwprevious,"bigcityio")
x1 <- caculate_increaseratio(dwbigcityionew,dwbigcityioprevious)

merge(dwbigcityionew,x1,by="bigcityio") %>% kable()

5. 机场

dwairportnew <- caculate_carsmean(jddwnew,"airport")

dwairportprevious <- caculate_carsmean(jddwprevious,"airport")
x1 <- caculate_increaseratio(dwairportnew,dwairportprevious)

merge(dwairportnew,x1,by="airport") %>% kable()

ahorawzy/Mreport documentation built on May 3, 2019, 3:40 p.m.