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()
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)
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)
(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()
(dwallnewday <- ddply(jddwnew,"day",summarise,Wmean = weighted.mean(cars,w=mileage))) %>% kable()
dwprovincenew <- caculate_carsmean(jddwnew,"province")
dwprovinceprevious <- caculate_carsmean(jddwprevious,"province") x1 <- caculate_increaseratio(dwprovincenew,dwprovinceprevious)
merge(dwprovincenew,x1,by="province") %>% kable()
dwscenerynew <- caculate_carsmean(jddwnew,"scenery")
同比
dwsceneryprevious <- caculate_carsmean(jddwprevious,"scenery") x1 <- caculate_increaseratio(dwscenerynew,dwsceneryprevious)
merge(dwscenerynew,x1,by="scenery") %>% kable()
dwbigcityionew <- caculate_carsmean(jddwnew,"bigcityio")
dwbigcityioprevious <- caculate_carsmean(jddwprevious,"bigcityio") x1 <- caculate_increaseratio(dwbigcityionew,dwbigcityioprevious)
merge(dwbigcityionew,x1,by="bigcityio") %>% kable()
dwairportnew <- caculate_carsmean(jddwnew,"airport")
dwairportprevious <- caculate_carsmean(jddwprevious,"airport") x1 <- caculate_increaseratio(dwairportnew,dwairportprevious)
merge(dwairportnew,x1,by="airport") %>% kable()
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.