Home on Website of ZhiYuan Wang

一篇新博客尝试

Wed, 05 Aug 2026 00:00:00 +0000

时隔五年，这是AI时代的一篇新博客

读《毛选（第一卷）》

Sun, 08 Aug 2021 00:00:00 +0000

《实践论》

马克思主义者认为人类的生产活动是最基本的实践活动，是决定其他一切活动的东西。

在阶级社会中，每一个人都在一定的阶级地位中生活，各种思想无不打上阶级的烙印。

ARCGIS常见操作汇总

Fri, 15 Jan 2021 00:00:00 +0000

1. 分乡镇出图

打开总图，选择要剪裁乡镇；
工具-分析-剪裁；
输入要素：需要剪裁的图层；剪裁要素：乡镇面；输出要素类：自己命名；XY容差，不填。

2. 矢量数据分面统计

分析工具-叠加分析-标识；
输入要素：公路网；标识要素：乡镇面；输出要素类：自己命名；XY容差，不填。
打开新标识图层属性表；
添加字段；长度-浮点型；
计算几何：长度，
汇总：

3. 按某一字段的属性值分成子元素集

分析-提取分析-按属性分割。

4. 将自然村点标记到乡镇面、行政村面

分析-叠加分析-标识。
输入要素：自然村点。
标识要素：乡镇面、行政村面。
输出要素：另存为。

5. 计算点到线的最近距离

分析-邻域分析-近邻分析。
输入要素：点。
临近要素：线（可为多图层的线）。

读《战神粟裕》

Sun, 10 Jan 2021 00:00:00 +0000

读书笔记

艰苦岁月不只给人难堪，还有作为战士必须经历的磨炼和考验，这样在进入下一个搏杀场时你才不致心虚胆怯。

粟裕技高一筹的地方，不是他没有失误，而是他能时时避免和纠正自己的失误。

写在新年的开始

Sat, 02 Jan 2021 00:00:00 +0000

为什么写在新年的开始不是1号而是2号呢？因为今天在翻1号日历的时候，突然感觉到了时间的流逝。之前看数字从1号到2号，没有什么感觉，人类的进化可能还没有对数字形成强烈的感受。而翻日历不同，每天都必须把前一天的日子翻到后面，这一天真的没有了，这种触感会虽然比撕下前一天的日历轻，但也足够让人震撼。

116号文-关于推动都市圈市域（郊）铁路加快发展意见的通知

Fri, 01 Jan 2021 00:00:00 +0000

国务院办公厅转发国家发展改革委等单位关于推动都市圈市域（郊）铁路加快发展意见的通知国办函〔2020〕116号

各省、自治区、直辖市人民政府，国务院各部委、各直属机构：

Geocomputation with R-2

Tue, 29 Dec 2020 00:00:00 +0000

载入相关包

library(sf)

## Linking to GEOS 3.8.0, GDAL 3.0.4, PROJ 6.3.1

library(raster)

## 载入需要的程辑包：sp

library(spData)
library(spDataLarge)

矢量数据：使用点，线和多边形；具有离散的，定义明确的边界，通常具有很高的精度。sf包处理矢量数据。
栅格数据：将表面分成恒定大小的方格，是底图的基础。raster包处理栅格数据。

2.2.1 简单特征（Simple features）简介

简单特征是一种开放标准，是一种表示大量几何图形的分层次的模型。

Geocomputation-with-R-3

Tue, 29 Dec 2020 00:00:00 +0000

Try add another post.

Geocomputation with R - Geographic data in R

Mon, 28 Dec 2020 00:00:00 +0000

2.2.1 简单特征（Simple features）简介

载入相关包

library(sf)

## Linking to GEOS 3.9.1, GDAL 3.2.1, PROJ 7.2.1

library(raster)

## 载入需要的程辑包：sp

library(spData)
library(spDataLarge)

矢量数据：使用点，线和多边形；具有离散的，定义明确的边界，通常具有很高的精度。sf包处理矢量数据。
栅格数据：将表面分成恒定大小的方格，是底图的基础。raster包处理栅格数据。

2.2.1 简单特征（Simple features）简介

简单特征是一种开放标准，是一种表示大量几何图形的分层次的模型。

城市轨道交通快线

Mon, 21 Dec 2020 00:00:00 +0000

1. 快线与市域铁路、城际铁路

国标《城市轨道交通线网规划标准》（GB/T50546-2018）中明确了市域城镇连绵地区或部分城市城区的快线承担以通勤出行为主的服务功能。快线主要服务于城市通勤圈范围。

又重返博客

Sun, 06 Dec 2020 00:00:00 +0000

上一次更新博客还是2020年的开头因为疫情困在家中记录蛋炒饭菜谱，一转眼再想着更新博客已经是2020年的末尾了。2020年是不平凡的一年，中间发生了很多变故，在这里就不多说了。总之，我慢慢找到了新的平衡，新的生活节奏。

通勤交通与职住关系

Sun, 06 Dec 2020 00:00:00 +0000

来源：《城市交通》2020年第5期

蛋炒饭

Mon, 10 Feb 2020 00:00:00 +0000

V1.0 黄瓜胡萝卜丝蛋炒饭

食材

黄瓜一个，擦丝；
胡萝卜一个，擦丝；
米饭，适量；
肉臊子或其他类型猪肉，适量；
鸡蛋2-3个，在小碗中打匀；
葱，切花；
火腿肠，切丁。
老抽/盐。

做法

热锅凉油，下葱花，炒香；
下鸡蛋，和葱花炒在一起，7-8成熟，盛出在小碗里；
热锅凉油，下肉臊子，炒香；
下米饭，压碎，翻炒；
下葱花鸡蛋，和米饭炒匀；
下黄瓜丝/胡萝卜丝/火腿肠丁，和米饭炒匀；
加少量老抽/盐，炒匀；
盛碗。

手撕包菜

Fri, 07 Feb 2020 00:00:00 +0000

食材

包菜（大的半颗，小的一颗）；
肉臊子适量；
干辣椒4个；
葱姜蒜；
盐；
料酒（50ml）；
生抽（15ml）；
蒸鱼豉油（10ml）；
蚝油（10ml）；

可选：

清炒芥菜

Fri, 07 Feb 2020 00:00:00 +0000

芥菜

芥菜俗称苦芥菜，被李时珍收入《本草纲目》中，能化痰平喘，消肿止痛。其种子磨粉就是芥末，榨出的油是芥子油。

芥菜味道苦，苦能降火，经常食用能解毒利尿。

六、交通分配与OD反推

Thu, 16 Jan 2020 00:00:00 +0000

交通分配是四阶段的最后一阶段，将各种出行方式的OD矩阵按一定的路径选择原则分配到交通网络中的各条道路上，求出各路段的流量及相关的交通指标。

三、出行生成预测

Wed, 15 Jan 2020 00:00:00 +0000

1. 出行的基本概念

1.1 出行的三个属性

每次出行有起讫点；
每次出行有一定目的；
每次出行使用一种或多种交通方式；

1.2 出行生成的表达

国内交通规划工作中对出行生成的定义：

五、方式划分预测

Wed, 15 Jan 2020 00:00:00 +0000

方式划分是四阶段的第三阶段，目的是将各小区间的出行分布量划分为各种交通方式的分布量。

在TransCAD中用Logit模型。

1. 基本原理与方法

1.1 基本概念

方式划分模型的基础：

四、出行分布预测

Wed, 15 Jan 2020 00:00:00 +0000

1. 基本原理与方法

出行生成预测：得到了未来年各小区出行的产生量与吸引量，需要进一步分析各校区之间的出行交换量。

出行分布预测：在已知分布矩阵各行、各列之和的条件下，求矩阵中每个元素具体数值的问题。

一、TransCAD——预备知识

Tue, 14 Jan 2020 00:00:00 +0000

1. 菜单项

File 文件
Edit 编辑
Map 地图
Dataview 数据视图
Selection 选择
Matrix 矩阵
Layout 布局
Tools 工具
Procedures 程序
Networks/Paths 网络/路径
Route Systems 路线系统
Planing 绘图
Transit 运输
Routing/Logistics 路线/物流
Statistics 统计数据
Window 窗口
Help 帮助

2. TransCAD的基本文件类型

文件类型	中文名	拓展名	归类	说明
Table	数据表	bin	数据	储存一般形式的表格数据
Matrix	矩阵	mtx	数据	储存矩阵形式的表格数据
Geographic File	地理文件	dbd	数据	储存点线面形式的空间数据
Route System	路线系统	rts	数据	储存公交线路和站点数据
Dataview	数据视图	dvw	视图	储存普通数据表的外观
Matrix View	矩阵视图	mvw	视图	储存矩阵数据表的外观
Map	地图	map	视图	储存地理文件、路线系等的外观
Chart/Figure	图表	fig	视图	储存基于数据表的统计图表
Layout	布局	lay	视图	将以上多种视图集合在一起
WorkSpace	工作空间	wrk	其他	同时储存多个窗口以便于管理

矩阵->矩阵视图
数据表->数据视图
数据表->图表
地理文件->地图
路线系->地图

3. 能够导入的文件类型

dbf文件，dBASE数据表
dxf文件，CAD图形文件
shp文件，ESRI地理文件

二、交通调查数据处理与分析

Tue, 14 Jan 2020 00:00:00 +0000

1. 交通调查的对象

社会经济基础资料
交通需求
交通设施
交通流量

2. 交通小区

不可能对每个交通源进行单独研究，要将分散的交通源按一定原则就近合并成若干个交通小区。交通小区在形式上表现为一个区域，但在逻辑上表现为一个点，这个点可称为小区的质心。

所里FTP设置

Thu, 09 Jan 2020 00:00:00 +0000

FXP 是File Exchange Protocol（文件交换协议）的英文简称，而中文简称为“文传协议”。用于Internet上的控制文件的双向传输。

打开flashfxp.exe文件。

铁路客运站点综合交通规划

Wed, 08 Jan 2020 00:00:00 +0000

11.4.1 铁路客运站点综合交通规划主要内容

铁路客运站点综合交通规划应从功能定位、场站布局、交通组织、用地布局等方面开展，主要内容包括：

功能定位及发展战略：根据铁路客运站点在城市总体规划和综合交通系统中的地位，确定铁路客运站点在城市发展中的定位，在综合交通系统中的定位以及所在片区的地位。根据城市社会经济发展目标与城市发展目标，制定铁路客运站点发展目标与发展策略。

交通指标汇总

Tue, 07 Jan 2020 00:00:00 +0000

0. 总

类别	建议指标	备选指标
交通需求	全市、关键区域或通道的：交通需求总量、交通出行结构	区域与对外交通、城市交通等各种交通方式的需求规模、结构、分布；每人或每车平均出行次数、出行距离。
交通供应	区域对外交通：可达性、运输能力	港口泊位等级、数量、航线数量、客货运吞吐能力等；公路等级、长度、密度、通行能力等；机场航线数量、客货运吞吐能力等；铁路长度、线路数量、客货运吞吐能力等；口岸客货运运输能力等；
	城市交通：可达性、交通承载能力	城市道路等级、长度、密度、通行能力等；公交线网数量、场站规模、车辆数、专用道长度，公交站点半径覆盖率等；轨道网络规模、密度、站点数量、车辆数，轨道站点半径覆盖率、周边覆盖人口岗位数量等；步行非机动车通道长度、密度，公共自行车租赁点数量等；停车设施供应规模、结构、供需缺口；
交通运行	基础指标：平均行程速度，服务水平/拥堵等级	交通流量、周转量、饱和度（负荷度）；平均延误；
	特征指标：各方式出行时耗，拥堵持续时长，拥堵里程比例	路段或路网行程时间比；不同拥堵等级的持续时间、比例；常发拥堵路段数量、长度、里程比例；道路行程时间可靠度；
	综合指标：交通运行指数	道路拥堵指数；拥堵延时指数；
交通运营	区域对外交通：客货运发送量，平均运送速度	各种方式客货运发送量的比例
	城市公共交通：全天或高峰客运量、交通分担率、平均运营车速、候车时间、准点率	公交平均换乘次数；公交、轨道交通车辆每平方米站立人数；公共交通平均步行到站距离；乘客满意度、投诉率；

1. 点

1.1 对外交通客运枢纽

1.1.0 换乘距离与时间限制

枢纽内主要换乘交通方式出入口之间旅客步行距离不宜超过200m。

城际铁路与综合交通

Tue, 07 Jan 2020 00:00:00 +0000

（出自《第七章城际铁路与城市群综合交通系统的协调和分工》）

城际铁路站点综合交通衔接规划

Tue, 07 Jan 2020 00:00:00 +0000

8.1 交通衔接规划主要内容

功能定位及发展战略

根据城际铁路站点在城市总体规划和综合交通系统中的地位，确定城际铁路站点在城市发展中的定位，在综合交通系统中的定位以及所在片区的地位。根据城市社会经济发展目标与城市发展目标，制定城际铁路站点发展目标与发展策略。

出差后联系商旅中心送机票

Mon, 23 Dec 2019 00:00:00 +0000

二院商旅部 86446868

读《亲密关系》

Tue, 12 Nov 2019 00:00:00 +0000

其他的人呢，纯粹是为恋爱而恋爱，一生中不断寻找能让我们热情不灭的特殊人物，而当热情熄灭时便转而寻找下一个可能目标。

亲密关系的过程，一开始是互相吸引，然后一步步发展出浪漫的火花。人一旦坠入情网，难免会有过高的期望，以致最后掉入幻灭的深渊。

读《巨人的陨落》

Tue, 12 Nov 2019 00:00:00 +0000

学校里的老师告诉男孩们，安全灯是矿主对雇员表示关切的方式之一。爸爸说：“就好像防止爆炸、避免停工或隧道受损对老板们没好处似的。”

他们看着那个可怕的洞穴，感到一种骇人的魔力，比利只觉得肚子里翻江倒海。黑暗似乎无限深远。他感受到了一种震颤，其中一半是侥幸，因为他不必进到里面，另一半是恐怖，因为总有一天他要下去的。

表头标准颜色

Tue, 29 Oct 2019 00:00:00 +0000

公司风格表格

表头颜色：224，240，235
表格边框颜色：153，204，255
表序号靠右，表名居中，小四黑体，Times New Roman
表内字体五号宋体，Times New Roman
图名：黑体五号
表名：黑体小四

找回路由器密码

Sun, 06 Oct 2019 00:00:00 +0000

今天突然发现我把屋里路由器的密码忘了，虽然笔记本、手机、kindle等设备都记着密码，但总觉得心里不踏实，还是找回密码的好。上网找到了解决方案：

读《局外人》

Wed, 02 Oct 2019 00:00:00 +0000

书评

《局外人》这本书是我毕业旅行时空闲时间边看的。毕业旅行后就到成都报到了，搬家和报到事情很杂很密，没时间梳理；再后来正式入职后笔记本电脑就搬到办公室用，所以此事竟一拖再拖。这次趁国庆放假要完结这件事，才好看后面的书。在整理书摘的过程中，我会一边分析主人公精神状态，一边回忆自己当时读过的感受，从我自己的观察、理解和思考来分析导致主人公最后悲剧的原因是什么。

读《怦然心动》

Sat, 28 Sep 2019 00:00:00 +0000

书评

之前在看过《怦然心动》电影后写过一篇博客看《怦然心动》。故事里女孩一家人的价值观帮助我更加清晰地认识了自己的价值观，看完后意犹未尽。于是在亚马逊上买下《怦然心动》的小说，想从小说中寻找更多关于女孩一家人价值观的解释，获取更多积极的能量。令我惊讶的是，小说的精彩程度丝毫不逊于改编的电影，甚至小说本身就是很好的作品，翻译的语言也简洁而真诚。故事分别从男孩和女孩的视角叙事，每人一章交替进行，最后汇合到一起，很有可读性。这仅仅是一篇中篇小说，而我勾画了185条摘抄，足以证明我对这部作品的喜爱和认可。

摘抄

（男孩视角）

朱莉没有马上被允许进入房间，因为妈妈注意到那双脏鞋并且要求她脱下来。等她脱下鞋子，妈妈又说她的脏袜子也不许穿进屋里。朱莉全然没觉得尴尬，一点儿也不。她只是拽下袜子，随手扔在我家门廊里。

PS里铁路线的画法

Sat, 24 Aug 2019 00:00:00 +0000

由于工作中经常要在底图上加铁路线，所以今天想攻克这个问题。

首先我尝试了用黑线块和白线块创建一个笔刷，用笔刷画铁路线。但实际效果和我想的不一样：只有在点笔刷时会出现铁路线，拖动笔刷时仍旧是黑色的，像是长了个铁路线外表的黑色刷子。

重返博客

Thu, 22 Aug 2019 00:00:00 +0000

已经有挺长一段时间没有更新博客和写日记了，原因是单位没有发电脑，就把笔记本电脑放在单位用了。因为既不可能在上班时间写博客写日记，也不想每天把电脑背来背去，所以就停更博客和日记将近两个月。而今天开始恢复写博客的原因也很简单：单位发电脑了。于是下班后我就把我的笔记本和显示器都带回来了。

姥爷说

Fri, 05 Jul 2019 00:00:00 +0000

姥爷是1929年出生的，今年已经九十岁了。俗话说“人活七十古来稀”，而从我懂事起姥爷就已经七十岁了，得益于坚持爬山运动和勤俭的生活习惯，一直身体硬朗、精神矍铄。可是岁月不饶人，姥爷近几年身体确实是越来越不好了，已有数次呼吸不畅被救护车送去急救，又都救回来了，但医生说两个肺中只有半个肺的功能了，每次发病得由姥姥立马插上氧气机输氧。所以每次放假回来必然回去看姥姥姥爷，因为见一次就少一次，谁知道这次不是最后一次呢？还记得爷爷去世之前的暑假，我坚持要去爷爷家几天，陪爷爷聊聊天。后面爷爷九月就去世了，还算是在去世前见了一面。所以这次从江南旅游回来去成都报到之前的短暂回家，也照例去看看姥姥和姥爷。

游江南（毕业旅行）

Fri, 05 Jul 2019 00:00:00 +0000

这不是一次单纯的旅行。
而是我给自己出的一道题目。
未来的路我得一人走，
能否依靠智慧做出判断，
依靠毅力坚持到底，
能否在完成目标的过程中，
欣赏美丽，体会快乐，同时保持体面。
我将用这道题试试自己。

读《边城》

Sun, 23 Jun 2019 00:00:00 +0000

书摘

他从不思索自己职务对于本人的意义，只是静静的很忠实的在那里活下去。

翠翠在风日里长养着，把皮肤变得黑黑的，触目为青山绿水，一对眸子清明如水晶，自然既长养她且教育她。

游北京世园会

Sat, 22 Jun 2019 00:00:00 +0000

从三月份知道北京要举办世园会，我就决定在离开北京之前一定要去一次，所以世园会和秦皇岛一样，都在我的毕业待完成清单里。今天是毕业典礼，而我昨晚毅然决定今天翘了毕业典礼，去游世园会，个中原因不再细说。Take full control of my life.

游山海关&北戴河

Fri, 21 Jun 2019 00:00:00 +0000

之前早打算在毕业这段时间去秦皇岛玩，但因为导师安排的编书任务和自己的拖延一直没抽时间去。眼看着距离离开北京越来越近了，列了个未完成事情清单，其中就包括去秦皇岛玩，总算把这件事提上了日程安排。

看《怦然心动》

Thu, 20 Jun 2019 00:00:00 +0000

“斯人若彩虹，遇上方知有。”最早知道这部电影是因为这句话。当时看到电影的封面是坐在树上的男孩女孩，以为故事讲述的是青梅竹马两小无猜的初恋。但我只猜对了一半，看完这部电影最大的感受不是懵懂初恋的美好，而是面对生活的态度。

读《世界文明史·中国&日本》

Mon, 17 Jun 2019 00:00:00 +0000

中国

中国和日本是《世界文明史》最后的一部分，看完了这部分，第一本《东方的遗产》就看完了。中间的印度部分由于过于晦涩，加上之前没有了解，看不太懂，所以直接跳过了。

离线安装chrome插件

Sat, 08 Jun 2019 00:00:00 +0000

上个月将笔记本重装了系统，但由于我的笔记本没法设置ipv6，所以无法登录谷歌商店下载、安装插件。昨天在学习OneNote时，接触到OneNote的chrome插件，突然有了灵感：能否离线安装chrome插件。

读《易经》（持续更新）

Thu, 16 May 2019 00:00:00 +0000

易之卦辞与爻辞含义十分灵活，从不同角度，不同层面，可以有不同理解，应仔细玩味，耐心探索。

经上

1. 乾卦

111111

乾，元亨利貞。

上乾下乾，也叫重天乾卦。天是阳气积聚而成，有纯阳之性。乾卦其象为天，其义为健，其德为元、亨、利、贞，六爻取龙为象，借潜、见、惕、跃、飞、亢，龙形之变，以明乾道变化之象。勉励人们效仿“天”之德，自强不息，顺天而行，相时而动。

读《动物农场》

Tue, 14 May 2019 00:00:00 +0000

例如他会说上帝赐给他尾巴以便驱赶苍蝇，然而他宁愿尾巴和苍蝇都不要。

我们还是实话实说吧：我们的生命是悲惨的，劳苦的，和短促的。我们生了下来，供给我们的食物仅仅够维持我们的躯体里始终有一口气，我们当中那些能活下来的，就被强迫干活，直到筋疲力尽；一旦我们的使用价值到了尽头，我们立马就会遭到骇人听闻的残酷杀戮。

读《多情剑客无情剑》

Mon, 29 Apr 2019 00:00:00 +0000

当他将“她”埋下去时，他自己的生命也就变得毫无意义。

但他一点也不着急，因为他知道这世上用金钱买不到的东西毕竟不多。

你看见过一匹狼在雪地上走路么？假如前面没有它的猎物，后面又没有追兵，它一定不肯走快的，因为它觉得光将力气用在走路上，未免太可惜了。

博客重新部署解决方案

Thu, 25 Apr 2019 00:00:00 +0000

前两天格式化电脑后，原来的Rstudio和各类设置都没有了，需要重新设置，本篇博客记录下收拾的顺序，以便下次参考。

下载并安装好R/Rstudio/Git。
设置Git
1. git config –global user.name “YOUR FULL NAME”
2. git config –global user.email “YOUR EMAIL ADDRESS”
生成SSH
1. 检查是否已有SSH：file.exists(“~/.ssh/id_rsa.pub”)
2. 使用Rstudio创建密钥，在Git面板，点击“Rreate RSA Key”
3. 点击“View public key”找到公钥
4. 将SSH公钥上传到Github
将博客等package克隆到本地
安装hugo
1. 尝试devtools::install_hugo()
2. 如果第一步连接不上，那么从github下载好hugo
3. 将hugo解压在D盘根目录下
4. 安装Hugo blogdown:::install_hugo_bin(“d:/hugo.exe”)
尝试修改并上传，理论上来说，设置好了SSH就不需要再在push时不需要输入用户名和密码了，证明连接成功

之后又遇到了Netlify部署网站出错的问题，具体表现在：

学习shiny包

Tue, 23 Apr 2019 00:00:00 +0000

1. Get Start

Shiny包有11个内置demo。

读《周恩来传》

Mon, 22 Apr 2019 00:00:00 +0000

摘抄

文雅、非常富有忍耐性，极为聪慧、机敏。

周的家庭长辈们是一些有修养的、遭受了艰难时世磨炼的绅士。

养母是淮安陈氏家族中最了不起的妇女，她有自信心，具有高度的智慧和社会良知。周深受这位妇女的影响，以至于在其后来的生活中，他完全把她当成了自己的“母亲”，正如他曾用此来称呼他的生母。

读《月亮与六便士》

Tue, 09 Apr 2019 00:00:00 +0000

摘抄

《圣经》反复告诫凡人不能彼此论断是有道理的，因为人和人的相互了解往往肤浅、局限而片面。

月亮象征着崇高的理想追求和美妙的精神境界。

六便士这种小面额的硬币代表着世俗的鸡虫得失与蝇头小利。

给每一行添加一个引用符

Tue, 02 Apr 2019 00:00:00 +0000

在用Kindle + Kindle Mate整理读书笔记，并转化为博客时候，经常需要给每一行一个符号“>”，来在markdown中表示引用。而当笔记较多时，一个一个添加引用符太过繁琐，所以想着写一个程序解决，如下：

读《战后日本经济史》

Tue, 02 Apr 2019 00:00:00 +0000

书摘

和我们一同躲在防空洞里的人大部分都窒息而死。

以后如果要躲在防空洞里，一定要躲在靠洞口的位置，氧气比里面充足。

“国家总动员体制”，即“1940”年体制：将全国所有资源都用来为战争服务。其理念是对产业实行国家统制。他们认为企业必须为公共利益做奉献，而不得追求私利。此外，也不允许不劳而获的特权阶级的存在。税收不经过地方政府直接交给国家，再由国家拨款给地方的模式初步形成。

读《追风筝的人》

Wed, 20 Mar 2019 00:00:00 +0000

书评

很早就听说过这本书了，因为它常年蝉联豆瓣读书排行榜前列，前段时间终于将其一口气读完。读的过程中很有触动，曾一度读到泪目，在候车厅盯着天花板尽量不让眼泪流出来。虽然读的是翻译版，也能感受到作者文笔很棒。不论是阿富汗曾经的富庶美好、还是现在的惨败零落，对景色的描写都很有画面感；友情、亲情、爱情，对人与人之间的羁绊渲染和刻画都很有代入感；情节的推动犹如电视剧，甚至有些语言描写颇有镜头感。虽然有以上种种优点，算是一本好书；但这本书仍有些过誉的可能，还算不上是本“伟大的书”。

《航拍中国》中的旅游景点

Wed, 06 Mar 2019 00:00:00 +0000

四川

九曲黄河第一湾；
黄龙、九寨沟；
都江堰；
龙泉山；
宽窄巷子；
人民公园；
汶川-映秀；
剑门关；
阆中古城；
渠县；
安岳石刻；
蜀南竹海；
乐山大佛；
峨眉山；
青城山；
上里古镇；
雅西高速；
大凉山；
德格印经院；
四姑娘山；
卧龙大熊猫基地；
川藏公路；
稻城亚丁；
海螺沟；

上海

外滩；
南京路；

R包中的相对路径

Wed, 06 Mar 2019 00:00:00 +0000

今天无意中解决了一个困扰我多年的问题——R中的路径究竟是怎么回事。我发现之前的我在引用路径时有两个问题：

使用了反斜杠\；这样为了在R中用路径，就必须打N个反斜杠。
使用了绝对路径；

所以正确的做法是：

利用插值法填补缺失值

Tue, 05 Mar 2019 00:00:00 +0000

今天解决了一个很重要的问题——对交通流序列插补缺失值。去年此时，我是用手工方法插补的，非常笨拙。这次要实现对所有交调站的自动化插补。

该问题可以用“分而治之”的思路切分，先用“桐木岭站”实验，成功后用lapply推广到所有站。

分组求和

Mon, 04 Mar 2019 00:00:00 +0000

今天下午在做数据预处理时，碰到了去年此时难到我的一个问题——将不同车道号或上下行方向的交通流量相加。我当时的解法很笨很粗暴，这次寻求更加简单的解决方案。

将字符串分拆成数据框

Mon, 04 Mar 2019 00:00:00 +0000

今天要解决一个长期困扰我的问题：使用strsplit或stringr::str_split，都是将元素分拆为列表的形式，虽然这样的通用性比较好，但是在大多数情形下，分拆元素的个数是固定的，更希望以数据框的形式返回，这样更具有通用性。

读《世界文明史·近东》

Fri, 01 Mar 2019 00:00:00 +0000

从这里开始，这本书就进入了正题——世界各个时期，各地的文明。文字也变得具体、生动、有趣起来。

苏美尔

一开始介绍的当然是人类历史上最古老的文明——苏美尔文明。它起源于两河流域（幼发拉底河和底格里斯河之间，也叫美索不达米亚平原）。历史从公元前40世纪-公元前24世纪，后被阿卡德（两河流域的另一文明）所灭，其文明被苏美尔和阿卡德的后裔——巴比伦人所继承。

读《咏远有李》

Fri, 15 Feb 2019 00:00:00 +0000

书摘

关于爱情

包容哈文的直接与脾气：

我给她起了个绰号，叫“劈头士”——劈头盖脸谩骂的人士。

越是直白，越是在乎。

怎样处理矛盾：

矛盾不过夜，过夜就是仇。”

读《世界文明史·文明的条件》

Tue, 12 Feb 2019 00:00:00 +0000

文明的条件是《世界文明史·东方的智慧》的第一部分，主要讲了文明的定义和内涵，以及其诞生的经济、政治、伦理条件。这一部分说理性的文字较多，读起来较为枯燥，描述案例也很简略，看起来较为费劲。不过能理解今天人类社会的种种现象，背后都有其起源，都有其原因，也饶有趣味。

读《中国历代政治得失》

Fri, 08 Feb 2019 00:00:00 +0000

本书的作者是史学大家钱穆先生。本书从历史学的观点，细数中国五个重要朝代：汉、唐、宋、明、清，的政治组织形态，特别是从“政治制度”的层面，陈述事实，分析利弊，论述观点。读完这本书收获很大，第一是理解了政府是怎么组织的，第二是纠正了过去对中国历史政治的很多偏见，第三是对很多政治的概念理解更加深刻。但本书仍有其自身的历史局限性，对于部分观点，我并非全部赞同。

学习《R for Data Science》(8)——Stringr

Wed, 30 Jan 2019 00:00:00 +0000

library(tidyverse)
library(stringr)

stringr没有被包括在tidyverse中

basic

在字符串中引用符号，需要加\

特殊字符：

\n newline
\t tab

String length 字符长度

str_length(c("a", "R for data science", NA))

## [1]  1 18 NA

Combining strings 字符合并

str_c("x", "y")

## [1] "xy"

str_c("x", "y", "z")

## [1] "xyz"

str_c("x", "y", sep = ", ")

## [1] "x, y"

将所要合并的向量作为第一个参数传入，需要用到collapse参数

读《21世纪资本论》（精选本）

Wed, 30 Jan 2019 00:00:00 +0000

摘抄

托马斯·马尔萨斯在1798年发表的《人口原理评述》一文中明确指出：人口过剩是影响财富分配的首要因素。

农业时代的古老观点，但现在读到仍然脊背发凉。精英阶层（地主阶层）为了解决人口过剩，会不会主动发起战争？“天地不仁以万物为刍狗，圣人不仁以百姓为刍狗。”

学习《R for Data Science》(7)——Relation Data

Tue, 29 Jan 2019 00:00:00 +0000

涉及到多个表的数据分析任务。

关系型数据的三个动词：

mutating joins，通过一个表中相匹配的观测向另一个表中添加变量；
filtering joins，通过一个表中的观测筛选另一个表中的观测；
set operations，将观测视为集合的元素。

library(tidyverse)
library(nycflights13)

nycflights13数据集

flights查询航班信息

学习《R for Data Science》(6)——Tidy

Mon, 28 Jan 2019 00:00:00 +0000

本章需要用到tidyr包

library(plyr)
library(tidyverse)
library(ggplot2)
library(pryr)

tidy data

同一组数据有多种表达方式

tidy data的定义：

Each variable must have its own column.
Each observation must have its own row.
Each value must have its own cell.

为什么要tidy data？

读《围城》

Mon, 28 Jan 2019 00:00:00 +0000

摘抄

我们对采摘不到的葡萄，不但想像它酸，也很可能想像它是分外地甜。

忠厚老实人的恶毒，像饭里的砂砾或者出骨鱼片里未净的刺，会给人一种不期待的伤痛。

张小龙在微信公开课上的演讲摘抄

Sun, 20 Jan 2019 00:00:00 +0000

讲到设计微信时的原动力：

原动力应该是你内心深处的一种很深的认知或者希望。它很强大，以至于说它可以坚持很久，并且你克服很多困难都要去做到它。

微信的原动力首先是一个好的工具，一个好的朋友。

读《中国哲学简史》

Sat, 19 Jan 2019 00:00:00 +0000

由于本书前半部分（先秦哲学）是在纸质书上看的，后半部分是在kindle上看的，所以摘抄只有后半部分的。前半部分已看过多边，印象深刻，所留到以后再重读。

学习链接

Tue, 01 Jan 2019 00:00:00 +0000

[谢益辉的博客](https://yihui.name/) [Rstudio Webinars](https://resources.rstudio.com/webinars) [bookdown书单-精华版](https://bookdown.org/) [bookdown书单-完全版](https://bookdown.org/home/archive/) [Geocomputation with R](https://geocompr.robinlovelace.net/) [blogdown](https://bookdown.org/yihui/blogdown/) [bookdown](https://bookdown.org/yihui/bookdown/) [rmarkdown](https://bookdown.org/yihui/rmarkdown/) [R for Data Science](http://r4ds.had.co.nz/) [Advanced R](https://adv-r.hadley.nz/) [统计之都](https://cosx.org/) [w3cschool](https://www.w3cschool.cn/)

学习《R for Data Science》(4)——Tibbles

Sun, 30 Dec 2018 00:00:00 +0000

学习更多的tibble知识，vignette(“tibble”)

（除了打印好看点，好像没发现tibble多有用？）

library(tidyverse)

## Warning: 程辑包'tidyverse'是用R版本3.5.1 来建造的

## -- Attaching packages -------------------------------------------------------- tidyverse 1.2.1 --

## √ ggplot2 2.2.1     √ purrr   0.2.5
## √ tibble  1.4.2     √ dplyr   0.7.6
## √ tidyr   0.8.1     √ stringr 1.3.1
## √ readr   1.1.1     √ forcats 0.3.0

## Warning: 程辑包'tidyr'是用R版本3.5.1 来建造的

## Warning: 程辑包'readr'是用R版本3.5.1 来建造的

## Warning: 程辑包'forcats'是用R版本3.5.1 来建造的

## -- Conflicts ----------------------------------------------------------- tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag()    masks stats::lag()

创建Tibble

tidyverse包里所有的函数，创建的都是tibble；

学习《R for Data Science》(5)——Data Import

Sun, 30 Dec 2018 00:00:00 +0000

概述
解析函数
写出文件

本文需要用到readr包

library(tidyverse)
library(hms)

概述

read_csv()读取逗号分隔；read_csv2()读取分号分隔；
read_tsv()读取tab分隔；
read_delim()读取任意分隔；
read_fwf()读取固定宽度文件；
read_log()读取Apache的log文件；

readr::read_csv相比base::read.csv的好处：

学习《R for Data Science》(2)——Transformation

Tue, 25 Dec 2018 00:00:00 +0000

使用dplyr包进行数据变形

library(tidyverse)

## Warning: 程辑包'tidyverse'是用R版本3.5.1 来建造的

## -- Attaching packages -------------------------------------------------------- tidyverse 1.2.1 --

## √ ggplot2 2.2.1     √ purrr   0.2.5
## √ tibble  1.4.2     √ dplyr   0.7.6
## √ tidyr   0.8.1     √ stringr 1.3.1
## √ readr   1.1.1     √ forcats 0.3.0

## Warning: 程辑包'tidyr'是用R版本3.5.1 来建造的

## Warning: 程辑包'readr'是用R版本3.5.1 来建造的

## Warning: 程辑包'forcats'是用R版本3.5.1 来建造的

## -- Conflicts ----------------------------------------------------------- tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag()    masks stats::lag()

library(nycflights13)

## Warning: 程辑包'nycflights13'是用R版本3.5.1 来建造的

tibble是一种基于data.frame的数据类型，为了更好地适应于tidyverse包。

学习《R for Data Science》(3)——EDA

Tue, 25 Dec 2018 00:00:00 +0000

问题
Variation
缺失值
Covariation
模式和模型
Project式工作流

系统地综合使用可视化和数据变形来探索数据。EDA(Exploratory Data Analysis)

思维模式(2)——系统式思维

Tue, 11 Dec 2018 00:00:00 +0000

本文主要来源于《系统之美》一书。

《系统之美》提供了一套“系统理论”，指导我们在面临混乱不堪、纷繁复杂且迅速变化的局面时，重新找回对整个系统的直觉，理解系统各个构成部分的能力，看清系统跟各个构成部件的关系，分析未来可能的行为趋势，创造性地得出“破局之道”。

思维模式(1)——供应链与供应商

Mon, 10 Dec 2018 00:00:00 +0000

本思维模式主要来源于《采购与供应链管理——一个实践者的角度》这本书。

供应链视角

供应链管理是对贯穿其中的产品流、信息流和资金流的集成。三者结合紧密。必须从整体来看待三条流。供应链的改善，最终都要落实在三条流的改善上。
通过外包战略改变系统边界、系统要素间结构和系统与环境的关系。内部职能外包给供应商，原来由运营部门负责的任务转移到采购部门。专业化程度更强，核心竞争力更强，但需要更强的供应（采购）管理能力。从集成式设计到模块化设计，所以外包也叫模块化。建立长期协作关系成本高则纵向集成，内部协调成本高则纵向剥离。剥离重资产运作，还可提高母公司的投资回报率，以提高股价。投资回报率=所得/付出。
物流管理：库存->仓储->物流网络->运输->3PL->订单履行->供给与需求规划。
运营管理：把原材料、人力、技术、资金、设备等转化为产品、服务的增值过程。需要有成套的系统、流程、培训和监控。运营管理的价值体现在对这些系统的设计、运营和改进。产品流还包括增值的生产过程，即运营的部分。
将销售计划，转变为生产总体计划，转变为物料供应计划、生产进度规划和生产的实施和控制。
包裹本身丢掉的概率很低，丢掉的是单据，或者单据和包裹分离，是个信息流问题。
如何确保数据的准确性，并从中提炼出合适的信息，是管理者的一项重要任务。数据、信息才是对供应链的综合掌握。职位越高，工作对象越是信息流。
提高信息透明度是改善系统运作的杠杆点之一。信息流的优化相对容易。信息流集成不仅是信息技术问题，还有制度问题和商业问题。关系：愿不愿意集成（合作伙伴层），是根本；连接：能不能够集成（订单层），是方法。
集成化供应链（ISC），供应链所有成员基于共同目标而组成的“虚拟组织”，组织内的成员通过三流的合作，在战略制定、业务流程、信息共享、联合预测、产品开发上全面合作，做到供应链计划协调、业务顺畅。其实质是通过虚拟化，对系统边界的模糊化和扩大化并对系统要素重新优化。
库存周转率：衡量材料在价值链中流通快慢。库存周转率是运营效率的体现。库存是企业运营水平的标尺。库存是供应链系统中的存量，容易监控和调节整个系统的指标。
深入基层，一线工人参与，群策群力。
不同的利益诉求来源于系统中各要素的有限理性。给不同部门设立统一的指标，给员工更具体的绩效考核目标，并适当组织调整。
复杂度是供应链管理的大敌。重要的是在控制成本的同时，对复杂度的控制与管理。设计独特、缺乏产品标准、缺乏对产品线的协调管理，是增加产品复杂度的根源之一。部门之间的壁垒和目标冲突，是复杂度泛滥的催化剂。谁能降低复杂度，谁就能更好地控制成本，谁就能领先对手。（但长尾理论，引用互联网管理方式，以较低的成本管理了高复杂度的产品线，并利用长尾效应，产生长尾盈利。）复杂度大增成了企业的通病。营收一直在长，利润在走下坡路，有了市场，没了利润。
非标件盛行，一个原因是缺乏反馈机制，让设计、营销、客户尽早知道非标设计的成本。需求管理需要跨职能协作。
真正的自由都是戴着某种镣铐跳舞。对成熟的流程、框架、模板的需要严格拷贝，从而减少复杂度。
纳入过多供应商，造成管理资源摊薄，管理深度不够；同时采购份额下降，规模效益下降，采购成本高。是供应链管理能力不足的体现。
复杂度的控制。产品的复杂度决定了组织、流程和供应链的复杂度。产品复杂度，精简产品线并产品标准化；流程复杂度，精益生产；组织复杂度，组织重构，机构精简；需求复杂度，整合需求；供应复杂度，整合供应商。
精益生产。拉动式准时化生产，强调物流平衡，由看板传递各层信息，对每一道工序全面质量管理，团队工作法（不安行政组织划分，按业务关系划分），并行工程（在产品设计开发期间，将概念设计、结构设计、工艺设计统一，项目小组定期或随时反馈，利用适当信息工具），追求零库存，追求快速应对市场变化，内外环境和谐统一，人本主义（尊重员工，重视培训，共同协作）。用最少的投入实现最大产出的目的。
复杂问题的简单解决方案——利用系统中的杠杆点。并且利用系统的自适应能力，系统有一定纠偏功能，时间长了总会达到某种平衡。
常年如一日，把一件件小事做好、做到位的结果。文火熬汤，慢药调理，大处着眼，小处着手。把小事做好，小事做扎实。治大国如烹小鲜。给飞驰的车换轮子，采取演进的方法，看上去慢，却能系统地解决问题。
大企业看上去强大，但在根本问题上，往往没有犯错、纠错的机会。复杂度控制是全民战争，预防为主。将复杂度与设计人员、销售人员的绩效挂钩。股东驱使公司多元化发展，使得公司落入大公司陷阱，从伟大到庞大。纵向一体化需要考虑：产业链上的利润分布，竞争优势和核心竞争力。动态多元化：投资公司。一个公司对外来资源的管理能力越低，越依赖竖向集成。巨无霸做事，速度慢，决策保守，运营水平低。后工业时代的纵向并购：获得知识产权和管理能力。国企也在转变思路，从纵向并购变为吃苦，利用资本增值，结果是投资现金流入可观。啥都做，啥都做不好，企业不堪自重而垮掉。
在专注和多元化上找到一个平衡点，不仅需要对未来的判断，还需要有能破能立的勇气。
生产型企业的三个设计：产品设计、工艺设计和供应链设计。产品设计和工艺设计交叉并行，相互应征。供应链设计围绕产品和工艺，设计和优化三流，以降低供应链成本，提高供应链效率。
推拉模型。按预测生产为推，按订单生产为拉。推的优点是缩短交货周期，充分利用产能，缺点是库存风险；拉的优点是控制库存积压风险，缺点是供应风险大增。推到一定地步，就得变成拉，即由客户的订单来驱动，以控制库存积压风险。定制化程度越高，预测准确度越低，推拉结合点离消费点越远。对产品实时性要求越高，推拉结合点离消费点越近。产品共同的部分好预测，用推式生产，以取得规模效益、降低成本；差异化部分用拉式生产，以降低需求变动带来库存风险。推拉结合点选择适当，会有效地平衡供应链的响应速度、成本和服务水平。供应链推拉点还要和产品战略相匹配。总体来说，推拉结合点在预测的边界上。
大公司的很多精力都花在整合供应商、集中采购额上。整合多方需求做预测，预测精度高，安全库存低。
一级供应商的选择不一定是工序，因为工序只考虑了产品流，没有考虑资金流和信息流。苹果开发和管理下级供应商，增强了自己对整个供应链的管控，同时也限制了富士康的议价能力。技术转由一级供应商管理，自己重营销而轻制造。一级供应商主要干力气活，靠劳苦谋口饭吃。
牛鞭效应的四大成因： • 多重需求预测 • 批量生产与订购 • 价格浮动和促销 • 理性预期三种应对措施： • 跨行 • 外包 • 信息共享
CPFR协同计划、预测和补货，在同一个预测下运作。销售与运营计划（S&OP），让公司在一个数字上运作。
从库存角度着眼供应链。供应链上的各种问题，都会在库存上体现出来。库存是此系统中最重要的存量。库存取决于企业和供应链的整体运营水平，是结果，而不是根源。库存需要系统地解决。
缩短周转周期，降低周转库存。周转周期包括了采购前置期，生产调试及安装期，信息周转周期。可以通过增加过程库存的手段，增加系统的平稳性，降低安全库存。精益生产的前提是稳定的生产计划。需求信息沿着供应链传导。信息周期增大，则采购前置期增大，则库存增多，供应链成本增加。流程改进小组，从跨职能层面审视信息流，进行公司层面的流程优化。信息流驱动产品流和资金流。成本是现象，根源是时间。
控制不确定因素，降低安全库存。由于需求和供给的不确定性，导致了安全库存，以平衡短缺风险和库存风险。信息不对称，有关系的问题也有连接的问题。成本沿着供应链传导提升了供应链的成本。库存意味着成本，供应商的成本总是会以各种方式转嫁给客户。信息不对称，执行不力，需求变动导致的不确定性。将噪声与模式分解，由系统自身纠正噪声，造成通过改进系统优化模式。
管理需求，管好供应链库存。需求管理不到位，导致供应链的波动。供应链管理的根本思路是突破部门壁垒，从局部优化上升到全局优化。
供应链降本三个台阶。第一个台阶是通过谈判来降低价格，谈判降价做的是利润的转移，太过激烈，则不利于整体供应链的可持续发展。第二个台阶是流程改进，精益生产，从而降低生产成本，电子商务，从而降低交易成本。第三个台阶是价值工程和价值分析。低附加值的事务由系统去做，高附加值的事务，如管理需求，理顺关系等，由采购人员做。采购必须发挥领袖作用，把公司内部的力量整合起来，与供应商一起推动流程改进，设计优化来降本。在集成产品设计与工艺设计上，采购扮演着重要角色，对外要管理，供应商，对内要管理内部需求。政府的繁文缛节增加了系统性成本。成本最终都属于整个供应链。供应商的库存增加，成本上升，最后还得转嫁给采购方。

供应商视角

语重心长的话

给飞驰的车换轮子。难就难在从一个烂摊子着手，在现有的框架上修修补补。采取演进的方式，看起来慢，却能更系统地解决问题。
系统工程需要漫长的微调和尝试。
多来点实在的修修补补，把一件件小事做好，把一个个小环节理清，把一个个细节做好、做扎实。
啥都做，啥都做不好；最后不堪自重而垮掉。
每条狗都有走运的一天。
决心永远大于技巧。
如何在专注和多元化上找到一个平衡点，不仅要基于对未来的判断，还要有能破能立的勇气。
复杂的问题大都是系统，而系统有一定的纠偏功能，时间长了总会达到某种平衡。系统地解决问题：治大国如烹小鲜。抓住复杂背后的简单原则——系统中的杠杆点。
注意各子系统之间的信息交流。提高信息透明度是改善系统运作的杠杆点之一。信息流不仅仅是连接问题，还有关系问题。
真正的自由都是戴着镣铐跳舞。
无谓的多元化不是创新。
苟延残喘的结果，无非是让让公司价值更低。
没法度量的没法管理。
与其埋头在沙子里，奢望奇迹发生，不如勇于承担责任，积极主动，将其转变为竞争优势。
供应链的能力由其最薄弱的环节决定。
大幅衰退以后是大幅上扬。
谁干得最好、最经济，就该由谁来干。
供应链各种问题在公司内部也存在。

读《在困境中更要发奋求进》(华罗庚演讲)

Mon, 03 Dec 2018 00:00:00 +0000

真正肯钻研的人在什么场合都可以钻研。

学习要自己努力，努力就可以很快上去。

白天紧张地搞优选法…我的理论研究室晚上进行的。

于其讲我自己所长的，不如将自己所短的。弄斧到班门，下棋找高手。

看《无敌破环王2·大闹互联网》

Sat, 01 Dec 2018 00:00:00 +0000

应豪杰邀，今天下午和豪杰、兴东一起去看《无敌破环王2·大闹互联网》（迪士尼），据说是年度动画片。上一年度是《寻梦环游记》（皮克斯），再上一年度是《你的名字》（新海诚），再上一年度是《疯狂动物城》（迪士尼）。

捕捉幸运之风

Sun, 25 Nov 2018 00:00:00 +0000

1

TED上一个关于“捕捉幸运之风”中提到的三点：

走出舒适区，承担一定的风险；
常怀感恩之心；
尊重各种各样的想法。

2

另外一篇讲运气的专栏文章

提升好运气的催化剂：

25岁做什么可以让五年后受益匪浅

Fri, 23 Nov 2018 00:00:00 +0000

知乎上关于25岁做什么可以在五年后受益匪浅的讨论：

学霸猫

作者：学霸猫链接：https://www.zhihu.com/question/52178718/answer/131556583 来源：知乎著作权归作者所有。商业转载请联系作者获得授权，非商业转载请注明出处。

当代中国与社会的书单

Thu, 22 Nov 2018 00:00:00 +0000

景跃进等《当代中国政府与政治》
李侃如《治理中国》
吴毅《小镇喧嚣》
项飙《跨越边界的社区：北京“浙江村”的生活史》
谭同学《双面人：转型农村中的人生、欲望与社会心态》
吴晓林《房权政治：中国城市社区的业主维权》
乔·史塔威尔《亚洲教父：香港、东南亚的金钱与权力》
卡尔·沃尔特《红色资本》
萧东连《荜路维艰》
连阔如《江湖丛谈》

交通流问题文献综述(2)

Wed, 21 Nov 2018 00:00:00 +0000

参考文献引用

Ahmed, M. S. and Cook, A. R. (1979) Analysis of freeway traffic time-series data by using Box–Jenkins techniques, Transportation Research Board, 722, pp. 1–9.
Vlahogianni E , Golias J , Karlaftis M . Short-term traffic forecasting: Overview of objectives and methods[J]. Transport Reviews, 2004, 24(5):533-557.
Vlahogianni E I , Karlaftis M G , Golias J C . Short-term traffic forecasting: Where we are and where we’re going[J]. Transportation Research Part C Emerging Technologies, 2014, 43:3-19.
Van Lint H , Van Hinsbergen C . Short-Term Traffic and Travel Time Prediction Models[J]. Transportation Research E-Circular, 2012.
Smith, B. L. and Demetsky, M. J. (1996) Multiple-interval freeway traffic flow forecasting, Transportation Research Record, 1554, pp. 136–141.
Habtemichael F G , Cetin M . Short-term traffic flow rate forecasting based on identifying similar traffic patterns[J]. Transportation Research Part C, 2016, 66:61-78.
Aghabozorgi S , Shirkhorshidi A S , Wah T Y . Time-series clustering - A decade review[M]. Elsevier Science Ltd. 2015.
Xia J , Chen M . A Nested Clustering Technique for Freeway Operating Condition Classification[J]. Computer-Aided Civil and Infrastructure Engineering, 2010, 22(6):430-437.
Xia J , Huang W , Guo J . A clustering approach to online freeway traffic state identification using ITS data[J]. KSCE Journal of Civil Engineering, 2012, 16(3):426-432.
Cheslow M , Hatcher S G , Patel V M . AN INITIAL EVALUATION OF ALTERNATIVE INTELLIGENT VEHICLE HIGHWAY SYSTEMS ARCHITECTURES[J]. System Architecture, 1992.

待补充论文

Vlahogianni 2004

关于短时交通流预测和ITS关系：

Cheslow et al. (1992) underlined the relation of dynamic traffic control to the ability to make and continuously update traffic predictions.

leaflet for R(4)——leafletCN色块图

Tue, 20 Nov 2018 00:00:00 +0000

(本博客中的数据已替换为随机数)

之前发现了一个leafletCN中的一个bug，就是geojsonMap函数的颜色映射有问题。

## Using Wmean as value column: use value.var to override.

	province	diff
1	安徽省	0.009
2	北京市	-0.007
3	福建省	-0.021
4	甘肃省	0.012
5	广东省	0.008
6	广西壮族自治区	0.078
7	贵州省	0.035
8	海南省	0.003
9	河北省	0.004
10	河南省	-0.018
11	黑龙江省	-0.006
12	湖北省	0.006
13	湖南省	0.000
14	吉林省	-0.001
15	江苏省	0.002
16	江西省	0.009
17	辽宁省	0.008
18	内蒙古自治区	-0.051
19	宁夏回族自治区	0.029
20	青海省	0.030
21	山东省	-0.012
22	山西省	-0.013
23	陕西省	0.009
25	四川省	0.006
26	天津市	-0.027
27	西藏自治区	-0.082
28	新疆维吾尔自治区	0.004
29	云南省	0.010
30	浙江省	0.016
31	重庆市	0.026

如果用默认的“Blues”，“Reds”这类颜色，是没有问题的

线性代数的本质（更新中）

Tue, 20 Nov 2018 00:00:00 +0000

在知乎看矩阵乘法的意义时，无意间发现B站上讲线性代数本质的系列视频，看了几集后觉得很好。重新看一遍，把重要观点记录下来。

0. 序言

一般课堂上讲线性代数总是从如何运算的角度讲，但学生往往不理解其中意义。几何水平的理解能够让你判断出：

文件的批量处理(2)——ldply

Mon, 19 Nov 2018 00:00:00 +0000

(本博客中使用的数据已全部替换为随机数)

在上一次的工作中，我将文件的批量处理执行地更加彻底，基本上驾轻就熟了。并且发现了plyr里很有用的工具：ldply来配合我的文件批处理。

川普的几句话

Sun, 18 Nov 2018 00:00:00 +0000

今天无意看到了个川普关于“成功学”的一个视频，虽然有些鸡汤（或者说是鸡血），但想想还是挺有启发的。摘录如下。

Thinking Big. 要开始思考任何事时，要大胆设想。大胆设想始于细节，成于细节。
Momentum 行动起来。
Stay focus 保持专注。失去专注的同时也就失去了行动。
Look at the solution, not the problem. 不管做什么，都会遇到麻烦。无论事情有多糟糕，持续专注于在解决问题上。会发现问题往往都能化险为夷，最终被解决。
see opportunity for what it is. 探寻机会的本质？每天问问自己，今天能学到什么以前所不知道的？把今天看作是一个机会。be open to ideas.
Know everything as you can about what you’re doing. 清楚掌握一切你正在做的事。这样能打好基础。
Be lucky. 越勤劳，越幸运。保持热情 Be passionate。热爱所从事的事业。
See yourself as a victor.这样能够消除消极的态度，积极地面对问题。将问题视作挑战。因为无论做什么，都会遇到挑战。要相信自己有能力去解决问题。
Be samrt. 天生我材必有用。善于利用手上的资源，而非被缺乏的资源所拘束。
Never give up. 最大的输家都是那些放弃的人。

Advanced R(1)——basic

Thu, 08 Nov 2018 00:00:00 +0000

library(magrittr)

第2章数据结构

2.1 向量

两种向量：原子向量，列表，三个共同属性：

typeof()
length()
attributes() 附加任意元数据

a <- 1:3
typeof(a)

## [1] "integer"

length(a)

## [1] 3

attributes(a)

## NULL

l <- list(a)
typeof(l)

## [1] "list"

length(l)

## [1] 1

attributes(l)

## NULL

2.1.1 原子向量

四种常见类型：

Advanced R(2)——debug

Thu, 08 Nov 2018 00:00:00 +0000

调试是一门用来解决代码中意想不到的错误的艺术和科学。

致命错误由stop()引起
警告由warning()产生
消息由message()产生

9.1 调试技巧

确认那些你认为对的的确是对的，知道找到一个错的（本来你认为是对的）。

Advanced R(3)——OO

Thu, 08 Nov 2018 00:00:00 +0000

面向对象

A class defines the behaviour of objects by describing their attributes and their relationship to other classes.

The class is also used when selecting methods, functions that behave differently depending on the class of their input.

基于消息传递的OO系统：With message-passing, messages (methods) are sent to objects and the object determines which function to call.
基于泛型函数的OO系统：a special type of function called a generic function decides which method to call.

library(pryr)

Base types

functions
environments
names
calls
promises

f <- function() {}
typeof(f)

## [1] "closure"

is.function(f)

## [1] TRUE

typeof(sum)

## [1] "builtin"

is.primitive(sum)

## [1] TRUE

S3对象

S对象统是最普遍使用的对象系统。是非正式的并且简洁的。

Advanced R(4)——FP

Thu, 08 Nov 2018 00:00:00 +0000

library(magrittr)

rm(list=ls())

第10章：函数式编程

函数式编程：创建和操作函数的工具。

适用于向量的所有操作也都适用于函数：

将函数赋值给变量；
将函数存储在列表；
函数作为参数传递给其他函数；
函数内创建函数；
把函数作为函数的结果返回；

10.1 动机

DRY: Do Not Repeat Yourself

思维模式(0)——序

Thu, 08 Nov 2018 00:00:00 +0000

一直以来我都有一种隐隐的思考，但直到今天我才能较为清楚地表达出来。那就是我发现，一个人通过书本学习或实践经历，将理论或方法抽象并提升到更高的维度，从而构建起将来思考问题的某一独特角度，或解决问题的某一独特方法，即思维模式；而各种思维模式通过独特的连接方式构成属于个人的独特思维体系，对每个人来说都是一件非常重要的事情。

画因子型日期的点画线时存在的问题

Thu, 08 Nov 2018 00:00:00 +0000

今天在画因子型日期的点画线时遇到了两个问题。

数据如下：

head(xzyk_mds[c("md","cars")],10)

##       md  cars
## 1   10-1 52356
## 2   10-2 53108
## 3   10-3 53738
## 4   10-4 50594
## 5   10-5 50012
## 6   10-6 49111
## 7   10-7 49118
## 8   10-8 49744
## 9   10-9 49855
## 10 10-10 48507

其中md为有序因子

《西游记》中最难的三难

Wed, 07 Nov 2018 00:00:00 +0000

今天看了“键客行”的视频，点评西游记中最艰难的三难。首先他说难不仅仅是肉体上的，更多是内心跟自己的斗争，或是一些无形的斗争。因为有形的斗争都容易警觉，想想办法总能克服；但能战胜自己的内心就难了；能识破软刀子，在无形的斗争中获胜就更难了。

【转载】中国进口博览会纪实

Wed, 07 Nov 2018 00:00:00 +0000

12个主宾国：

印度尼西亚、越南、巴基斯坦、南非、埃及、俄罗斯、英国、匈牙利、德国、加拿大、巴西、墨西哥

展馆展区

国家贸易投资综合展区：展示贸易投资领域有关情况，包括货物贸易、服务贸易、产业状况、投资旅游，以及各国有特色的产品。只展示不成交。
消费电子及家电展区：移动设备，智能家居，智能家电，虚拟现实与增强现实，电子游戏，健康运动产品，音频产品，视频与高清设备，生活科技，显示技术，在线与家庭娱乐，产品与系统解决方案等。
服装服饰及日用消费品展区：服装、纺织品、丝绸产品、餐厨用品、家居用品、礼品、家居装饰品、节日用品、珠宝首饰、家具、婴童用品、玩具、文化用品、美容美发护理产品、运动及休闲产品、箱包、鞋、钟表、陶瓷和玻璃制品等。汽车展区：智能驾驶汽车与技术，互联网汽车与技术，新能源汽车与技术，品牌汽车等。
智能及高端装备展区：人工智能，工业自动化与机器人，数字化工厂，及物联网，材料加工及成型装备，工业零部件，信息通信技术装备，节能环保装备，新能源电力电工装备，航空航天技术装备，动力传动与控制技术装备，3D 打印等。食品及农产品展区：乳制品、肉制品、水产品、蔬果、茶和咖啡、饮料及酒类、甜食及休闲食品、调味品、罐头及方便食品等。
医疗器械及医药保健展区：医疗器械，药品，健康及保健品，传统医学产品，医美产品，养老与康复，制药机械与设备等。
服务贸易展区：新兴技术，服务外包，创意设计，文化教育，旅游服务，物流服务，综合服务等。

纪实

【转自新华社《人间正道浩荡前行——习近平主席出席首届中国国际进口博览会纪实》】

ggplot2一页多图的解决方案

Tue, 06 Nov 2018 00:00:00 +0000

在网上找到一个函数，能够很好地解决ggplot2一页多图，函数如下：

# Multiple plot function
#
# ggplot objects can be passed in ..., or to plotlist (as a list of ggplot objects)
# - cols:   Number of columns in layout
# - layout: A matrix specifying the layout. If present, 'cols' is ignored.
#
# If the layout is something like matrix(c(1,2,3,3), nrow=2, byrow=TRUE),
# then plot 1 will go in the upper left, 2 will go in the upper right, and
# 3 will go all the way across the bottom.
#
multiplot <- function(..., plotlist=NULL, file, cols=1, layout=NULL) {
  library(grid)
 
  # Make a list from the ... arguments and plotlist
  plots <- c(list(...), plotlist)
 
  numPlots = length(plots)
 
  # If layout is NULL, then use 'cols' to determine layout
  if (is.null(layout)) {
    # Make the panel
    # ncol: Number of columns of plots
    # nrow: Number of rows needed, calculated from # of cols
    layout <- matrix(seq(1, cols * ceiling(numPlots/cols)),
                    ncol = cols, nrow = ceiling(numPlots/cols))
  }
 
 if (numPlots==1) {
    print(plots[[1]])
 
  } else {
    # Set up the page
    grid.newpage()
    pushViewport(viewport(layout = grid.layout(nrow(layout), ncol(layout))))
 
    # Make each plot, in the correct location
    for (i in 1:numPlots) {
      # Get the i,j matrix positions of the regions that contain this subplot
      matchidx <- as.data.frame(which(layout == i, arr.ind = TRUE))
 
      print(plots[[i]], vp = viewport(layout.pos.row = matchidx$row,
                                      layout.pos.col = matchidx$col))
    }
  }
}

我已将其收藏到usefulr里。

《天龙八部》中的主角们

Mon, 05 Nov 2018 00:00:00 +0000

最近在看一个Up主叫“键行客”推的关于金庸小说介绍的系列视频，每段视频先用非常浅显的语言叙述故事梗概，后半部分有其对作品和人物的相关理解和思考。我认为思考和点评挺有水平的，值得反复品味，所以记下来。

缅怀金庸先生

Wed, 31 Oct 2018 00:00:00 +0000

昨晚看新闻说金庸先生去世了，享年94岁。金庸的“飞雪连天射白鹿，笑书神侠倚碧鸳”14部小说，对中国的70后、80后、90后影响太深，对我的影响也太深。我从完全记事起看过的第一部电视剧就是奶奶给我买的83版的《射雕英雄传》，那是我看过最好看的电视剧。后来《神雕侠侣》、《天龙八部》、《倚天屠龙记》等电视剧也是看过不止三四遍。

面试

Wed, 31 Oct 2018 00:00:00 +0000

无领导小组

注意要点：

询问结果呈现形式的要求；
服务组员；
征求大家的讨论方式，是直接讨论还是先个人思考（5 min）后再讨论；
选出时间控制者；
每个人的结论呈现方式；
用纸笔记录每个人所陈述的意见和理由；
按顺序陈述，根据之前陈述人的观点调整自己观点中的不足；
分析题目中的判别标准；
先求同，再存异；
细化他人意见；
跳出原有思维框架，进行补充；
最后陈述者先向小组成员陈述一遍，查漏补缺，取得大家的共同同意；
共同鼓励；
最后陈述时要有对讨论过程的总结；
哪位组员贡献最大？集体智慧，从不同角度给出建议，并最终统一形成答案。
认为讨论环节哪些是做的好的，哪些是无效的？
结束后有头有尾，收拾好桌面细节。
在自己取得优势后对其他人的关怀。
做有把握的事情、做无争议的事情、出现争议时给出理由、陈述清楚即可，不要纠缠讨论。
在组内有人获得领导力，得到多人支持时，如果不是只有一人可以胜出的测试，此人又是一个懂“规则”的人，尽量更换跟随策略，做好一个“可合作”者，同时将自己不擅长展现的内容都尽量让此人来展现。

具体企业面试

1. 移动研究院

一面：

申论

Mon, 29 Oct 2018 00:00:00 +0000

1. 总论

考察对给定材料的分析、概括、提炼、加工。

在把握材料主旨和精神的基础上，形成并提出自己的观点、思路或解决方案，准确流畅地用文字形式表达出来。

港珠澳大桥的脆弱性假想

Sun, 28 Oct 2018 00:00:00 +0000

背景

今天重庆发生了一起交通事故：重庆万州一辆公交车在桥上撞上一辆逆行的小客车，掉入了江中，目前遇难人数不明。我在脑海里模拟这件事故发生的始末，发现除了逆行的小客车是可以避免的以外，其他可控的因素非常少，而仅仅是一辆逆行的小客车就能立马结束人数未知的一车人的生命，几乎没有能够生还的可能。由此可以提出疑问：桥，这种古老的、常见的又在近年来大规模建设的交通基础设施，在整个交通系统中，是否具有显著的脆弱性？

两首神似的歌

Sat, 27 Oct 2018 00:00:00 +0000

今天周六，一个人来到单位刷行测题，今天的主要任务是材料分析题。边刷题边听网易云音乐的日推，其中有这么一首歌。

Tennessee Waltz。

我听怎么感觉这么耳熟呢？想起来后弦有首类似的，一听果然很像，节奏和旋律非常神似。

看《American Geinus》（未完）

Fri, 26 Oct 2018 00:00:00 +0000

今晚终于把这部8集的纪录片看完第一遍了。这部纪录片是我看过为数不多觉得很精彩的，能推动我完整看完一遍的。这么喜欢的之前只看过《河西走廊》和《鸟瞰中国》。这部纪录片好看的原因大体有三点：一是美剧版的镜头语言和叙事手法，类似于《河西走廊》，感觉就像在历史场景中身临其境一般，非常有代入感，画面也很精致典雅；二是每集讲的都不是一个天才，二是两个天才的竞争和“对决”，这样的竞争氛围产生了刺激感，激发了观众的好奇心；三是和我个人的兴趣相吻合，我从内心里也是崇拜这样的天的，也是想通过发挥创造力达到一些非凡的成就。

blogdown中的公式支持

Thu, 25 Oct 2018 00:00:00 +0000

最近在总结行测题时，发现我的xmin主题不支持公式了，加了反双引号也没用。用knit渲染出来是带公式的，但博客里怎么都显示不了公式。所以判断不是rmarkdown的问题，而是blogdown的问题。虽然是个小问题，但如鲠在喉，怎么都解决不了，快把我的强迫症逼出来了。今天实在忍无可忍，下手去解决。

行测(3)——材料分析题

Thu, 25 Oct 2018 00:00:00 +0000

注意

注意审题，严格对照清楚问的是谁，别计算错了对象；

例题

1. 已知今年水平，问去年（前年）水平

去年=今年/(1+增长率)

2015年蔬菜产量15.79万吨，下降3.4%；水果产量7.84万吨，增长7.4%；2014年蔬菜比水果约高多少万吨？

行测(4)——专业知识

Thu, 25 Oct 2018 00:00:00 +0000

1. 金融学

1.1 货币

货币是固定充当一般等价物的特殊商品。

五个职能：价值尺度、流通手段、支付手段（清偿债务）、贮藏手段、世界货币。

充当世界货币的条件：

学习《R for Data Science》(1)——visualization

Wed, 24 Oct 2018 00:00:00 +0000

3. 数据可视化

library(tidyverse)

这本书是哈德利的入门之作。虽然对于我来说有点过于简单，不过这本书是以系统的观点阐述数据分析中的各个问题，并且都是用tidyverse工具箱中的包解决的。利用闲暇时间读一遍能够系统学习一边tidyverse，也不算浪费时间。不过这本书过的速度可以稍微快一点，毕竟很多都是太基础的内容。

小米手环3

Wed, 24 Oct 2018 00:00:00 +0000

【这篇博客是没有收广告费的“硬文”，只因为非常喜欢，想记录下来】

小米手环3发布已经快半年了，之前看评测说屏幕大了，多了很多功能，也没太在意。最近觉得手环2不能看信息比较麻烦，这一年各种智联招聘的短信又很多，所以想买个能看信息内容的手环，就是小米手环3了。NFC版放开购买已经是八九月份了，所以我还算是较早买的。周五订了小米手环3，周六就到了，晚上开盒，看到水滴般的屏幕，戴上以后腕带很柔软，手环的存在感比2代更低，非常喜欢。

看《超级工程·港珠澳大桥》

Tue, 23 Oct 2018 00:00:00 +0000

今天是港珠澳大桥正式通车的日子，但之前对这座桥知之甚少。于是和豪杰看了《超级工程·港珠澳大桥》这部纪录片。

内容

纪录片里主要讲了大桥在规划、设计、施工建设方面遇到的一些工程难题，并介绍是如何解决的，以下几个印象比较深刻。

行测(1)——数量关系题

Tue, 23 Oct 2018 00:00:00 +0000

数字推理题

1. 多级数列

对相邻两项进行四则运算（做差，做商，做和），进行一次或两次运算。

若呈平稳递增趋势，可尝试使用倒三角法寻找规律。

有时做差和做商交替使用。

行测(2)——推理题

Tue, 23 Oct 2018 00:00:00 +0000

1. 图形推理题

注意：

九宫格问题一般横着找规律；
有些存在交替规律或递增递减规律；
九宫格问题：第一与第二存在相同规律，第二与第三存在相同规律

规律总结：

1.1 图形元素数量规律

规律：相同，递增，递减

看《萤火之森》

Sun, 21 Oct 2018 00:00:00 +0000

《萤火之森》是一部2011年的日本动漫，导演是大森贵弘。

故事

故事情节非常简单，讲一个小女孩萤在6岁那年暑假在乡下爷爷家玩耍时，误入山上的森林迷路了。这时候遇到一个戴着面具的小哥哥银。萤觉得能获救了，扑向银，结果他躲开了让她扑了个空。原来银是森林里的“妖怪”，不能让人碰，碰了就会消失，也就意味着灭亡。两个人用树枝牵着，走出了森林。第二天，萤又来了，带了吃的表达感谢。两个人愉快地玩了一个暑假，玩的时候也会碰到其他的“妖怪”，警告萤不要碰银。暑假结束的时候，萤要回城市了，银邀请她明年暑假继续来玩。就这样，玩过了很多个暑假。小姑娘慢慢长成了高中生，对银的感情也不仅仅是童年时的玩伴，听说冬天会冷，给银送了围巾，也会在冬天想念银。

听《言叶之庭》片尾曲《Rain》

Fri, 19 Oct 2018 00:00:00 +0000

《言叶之庭》是好几年前看过的动漫了。昨天在网易云音乐的视频推荐里看到了这首歌的MV，感觉色彩很好看，歌手的声音很温暖温柔，所以收藏到博客里。

有时候感觉北京天空的云和夕阳也很有新海诚漫画中的感觉。

【转载】阿里云的这群疯子

Wed, 17 Oct 2018 00:00:00 +0000

链接《阿里云的这群疯子》

这篇文章很长，不过我是一口气读完的。也许是故事本身很坎坷，也许是作者描写得过于精彩。读完后，我为阿里云创造出这样的奇迹感到惊讶，也为自己的渺小和之前浪费的时间感到惭愧。就如同雷军看完《硅谷之火》后，下决心创造出一家世界级的科技公司。我觉得24岁已过的我，也是时候去解放自己给自己的束缚，放手去做一些事了。

交通流问题文献综述(1)

Wed, 17 Oct 2018 00:00:00 +0000

本篇文章服务于我的硕士毕业论文以及与之相关的准备投寄的小论文，也许可能也可以为未来在这方面更多的研究做积累。

英文文献

0. Analysis of freeway traffic time-series data by using Box-Jenkins techniques.

开篇之作

log函数在图形映射中的妙用

Tue, 16 Oct 2018 00:00:00 +0000

library(leaflet)
library(leafletCN)

今天在工作中遇到这样一个问题：将某数值映射到leaflet上点的颜色，比如数值越大红色越深，数值越小红色越浅。

我第一次是这样做的

reg <- leafletCN::leafletGeo("China")
pal <- colorNumeric("Reds",y$cars)
leaflet(y) %>% 
  addProviderTiles("OpenStreetMap.BlackAndWhite") %>% 
  addPolygons(data = reg, stroke = TRUE, smoothFactor = 1, fillOpacity = 0, weight = 1) %>% 
  addCircles(lng = ~lng,lat = ~lat,radius = ~cars/5,color = ~pal(y$cars))

## Warning in validateCoords(lng, lat, funcName): Data contains 61 rows with
## either missing or invalid lat/lon values and will be ignored

操作没啥大问题，但是由于数值差距太大，使得大数值的被映射为正常的红色，而小数值的是浅肉色，接近于底图的颜色，很不好分辨。

28字箴言

Sun, 14 Oct 2018 00:00:00 +0000

邓公针对1990年前的国际局势提出的对外总方针。

冷静观察，
稳住阵脚，
沉着应付，
韬光养晦，
善于藏拙，
决不当头，
有所作为。

交通流问题相关表达总结(1)

Sat, 13 Oct 2018 00:00:00 +0000

1. 套话

结论性 This research provides strong evidence suggesting that
得到良好结果 provides promising results
良好的实用性 it can be easily incorporated with real-time traffic control for proactive freeway traffic management.
最原始的目标 The prime objective of traffic management strategies is to
局限于 has been limited to
随着科技进步 with the advancement in technology and the wide deployment of intelligent transportation systems
该领域已有大量的研究 , considerable amount of research has been focused on the topic
提供方法解决问题 In this paper, we present a methodology for short-term traffic forecast based on learning similar traffic patterns as a reference for providing the predictions on future traffic.

2. 短时交通流

2.1 预测

短时交通流预测 short-term traffic flow rate forecasting
短时交通预测 short-term traffic forecasting
及时可靠准确的预测 timely, reliably, and accurately
交通管理和控制应用 traffic management and control applications
动态选择交通管理策略 dynamically apply alternative strategies proactively in response to predicted traffic conditions.
短时交通流预测的必要性 Short-term traffic forecasting models, therefore, are an integral element of the toolset needed for real-time traffic control and management. Moreover, such tools are important in providing travelers with reliable travel time information, optimizing traffic signals, and deployment of emergency management systems.
提前预测交通量 predicting the expected volume of traffic ahead of time
短时交通流预测的原理 The availability of a vast amount of spatial and temporal traffic data coupled with advancements in statistics and data analysis techniques have created an opportunity to perform short-term traffic forecast with a reasonable prediction accuracy and short processing time.
短时交通流预测的定义 Short-term traffic forecast aims at predicting the evolution of traffic over time horizons ranging from few seconds to few hours
短时交通流预测的应用领域 This interest is the direct result of an increasing need for developing user friendly applications which can both provide accurate information to drivers and be used for signal optimization.
交通管理部门 transportation agencies
交通管理与控制 traffic management and control application
时间粒度 data resolution (e.g. 5mins data)
交通拥堵 traffic congestion
在线应用 online implementation
高速公路 freeway

2.2 相似交通流识别

识别相似交通模式 identifying similar traffic pattern
交通流模式识别 traffic flow pattern analysis
Identifying spatial and temporal flow patterns has been an important consideration in short-term traffic forecasting research.
在线交通流状态识别 Online traffic state identification

聚类：

读《聊斋志异·画皮》

Fri, 12 Oct 2018 00:00:00 +0000

《画皮》是《聊斋》中的经典名篇，鬼画人皮的形象也多次在荧幕上呈现，被大众所熟知。但《画皮》的原版故事我一直没有读过，今天正好读到。

故事原文见链接《画皮》原文，白话文见链接《画皮》白话文。

文件的批量处理

Thu, 11 Oct 2018 00:00:00 +0000

今天我要面对之前一直想正面解决的一个问题：文件的批量操作。但之前对列表这种数据结构认识不深刻，不敢用也不会用，所以一直没能真正解决。

通过《Advanced R》里关于列表和泛函的操作，我现在已经有信心通过列表来进行文件的批量操作了。这篇文章仅展示批量操作的代码，并不运行。

合并列表

Wed, 10 Oct 2018 00:00:00 +0000

今天在工作中遇到一个问题：当我用lapply处理完一个列表后，会得到结果列表。结果列表的某一列都是相同的索引，我需要把列表中每一个元素都根据索引合并到一个DataFrame中，该怎么操作？

看《影》

Thu, 04 Oct 2018 00:00:00 +0000

今天和豪杰去看了张艺谋导演的电影《影》。看之前就有评价说，这是张艺谋电影的新高度。买的中影电影院（小西天店）的巨幕厅，还是有些小期待的。

故事

故事是由《三国·荆州》的故事架空后改编的。故事讲一国的主将都督（子虞）一年前被敌国主将杨苍砍伤，从此容貌大变，于是启用了一直被训练了替身“影子”境州（取此名是子虞想让他牢记自己的使命——破境州）。由于境州和子虞长得很像，而且多年来训练有素，一直没有被识破。但战败于杨苍，丢掉境州的失败，一直让子虞有着复仇的计划。可惜大王“沛公”沉迷于酒色，一心求和，不想拿下境州。于是（子虞让）境州去给杨苍父子下战书，之后回来后辞官。在朝堂上，大王借给境州敷药，发现境州身上的是新伤。

读《心流》

Thu, 04 Oct 2018 00:00:00 +0000

0.序言
- 0.1 自造内心秩序之途
第1章心流，快乐的源泉
第2章控制意识，改善体验的品质
第3章享乐与乐趣
第4章如何在日常生活中寻找心流
第5章感官之乐
第6章思维之乐
第7章工作之乐
第8章人际之乐

本书的一大目标就是发掘各种方法，以便把理性的细节转变成具有个人意义的游戏，导向最优体验。

谢益辉精彩博客（持续更新）

Thu, 04 Oct 2018 00:00:00 +0000

韩愈的《师说》中说：“古之学者必有师。师者，所以传道受业解惑也。”从小学到中学到大学，拥有老师似乎是理所当然的事。但一个人一旦脱离了校园环境，或是在校园环境中无法找到能够指引自己的老师（比如硕士阶段的我），便需要向外寻求能够在专业、技能，甚至是观念上能够指引自己的老师。好在互联网时代，人与人、人与社区之间的连接不再被距离或连接方式或连接成本所约束。博客、webinar（网络研讨会）就是非常好的找到老师、跟老师学习的途径。就如同古代佛学大师们开坛讲经一样，现在拥有一技之长或一家之言的人们开博客，不仅讲技术，也讲自己的观念，让思想能够广泛传播出去。

看《海上钢琴师》

Wed, 03 Oct 2018 00:00:00 +0000

背景

《海上钢琴师》这部电影我很早之前就已经听人推荐，知道是一部豆瓣评分很高（9.1分）的电影。导演是意大利人吉赛贝·托纳多雷，他还导演了《天堂电影院》和《西西里的美丽传说》。近期在网易云上听了一首《海上钢琴师》的背景音乐（The Crisis），觉得很喜欢，于是去看电影。

读《吕不韦列传》

Tue, 02 Oct 2018 00:00:00 +0000

往来贩贱卖贵，家累千金。

吕不韦是个低买高卖的渠道商，善于发现处于价值洼地的资产。

“此奇货可居。”

吕不韦知道子楚时，发现了子楚的价值，认为“此奇货可居也”。之后又发现子楚生活条件不佳，抓住了子楚内心的需求，用千金作为投入资本，获得了子楚的信任，子楚用自己皇子的身份入股，吕不韦达成了和子楚的利益联盟。

读《屈原贾生列传》

Mon, 01 Oct 2018 00:00:00 +0000

博闻强识，明于治乱，娴于辞令。

这是描写屈原的个人能力：记忆力好，管理能力强，系统性思维能力强，文笔好、讲话水平高。

王使屈平为令，众莫不知，每一令出，平伐其功，以为“非我莫能为”也。

stringr学习总结

Sun, 30 Sep 2018 00:00:00 +0000

library(stringr)

1. stringr包函数一览

1.1 字符串拼接函数

str_c: 字符串拼接。
str_join: 字符串拼接，同str_c。
str_trim: 去掉字符串的空格和TAB()
str_pad: 补充字符串的长度
str_dup: 复制字符串
str_wrap: 控制字符串输出格式
str_sub: 截取字符串
str_sub<- 截取字符串，并赋值，同str_sub

1.2 字符串计算函数

str_count: 字符串计数
str_length: 字符串长度
str_sort: 字符串值排序
str_order: 字符串索引排序，规则同str_sort

1.3 字符串匹配函数

str_split: 字符串分割
str_split_fixed: 字符串分割，同str_split
str_subset: 返回匹配的字符串
word: 从文本中提取单词
str_detect: 检查匹配字符串的字符
str_match: 从字符串中提取匹配组。
str_match_all: 从字符串中提取匹配组，同str_match
str_replace: 字符串替换
str_replace_all: 字符串替换，同str_replace
str_replace_na:把NA替换为NA字符串
str_locate: 找到匹配的字符串的位置。
str_locate_all: 找到匹配的字符串的位置,同str_locate
str_extract: 从字符串中提取匹配字符
str_extract_all: 从字符串中提取匹配字符，同str_extract

1.4 字符串变换函数

str_conv: 字符编码转换
str_to_upper: 字符串转成大写
str_to_lower: 字符串转成小写,规则同str_to_upper
str_to_title: 字符串转成首字母大写,规则同str_to_upper

1.5 参数控制函数，仅用于构造功能的参数，不能独立使用。

boundary: 定义使用边界
coll: 定义字符串标准排序规则。
fixed: 定义用于匹配的字符，包括正则表达式中的转义符
regex: 定义正则表达式

2. 详细函数

2.1 字符串拼接函数

2.1.1 str_c

str_c(…, sep = “”, collapse = NULL)

leaflet for R(3)——leaflet.extras

Sat, 29 Sep 2018 00:00:00 +0000

11 leaflet.extras

11.1 addHeatmap

library(leaflet)
library(leafletCN)
library(leaflet.extras)
library(magrittr)

leaflet(quakes) %>% addProviderTiles(providers$Esri.WorldImagery) %>%
  setView( 178, -20, 5 ) %>%
    addHeatmap(lng = ~long, lat = ~lat, intensity = ~mag,
               blur = 20, max = 0.05, radius = 15)

library(sp)
jsURL <- 'https://rawgit.com/Norkart/Leaflet-MiniMap/master/example/local_pubs_restaurant_norway.js'
v8 <- V8::v8()
v8$source(jsURL)
geoJson <- geojsonio::as.json(v8$get('pubsGeoJSON'))
spdf <- geojsonio::geojson_sp(geoJson)

leaflet(spdf) %>%
  addProviderTiles(providers$Thunderforest.TransportDark) %>%
  addHeatmap(blur = 20, max = 0.05, radius = 15)
  
leaflet(spdf) %>%
  addProviderTiles(providers$Thunderforest.Transport) %>%
  addHeatmap(blur = 20, max = 0.05, radius = 15, gradient = 'Greys')
  
london.crimes.files <- Sys.glob(
  paste0(system.file('examples/data/London-Crimes', package='leaflet.extras'),
         '/*/*-city-of-london-street.csv.zip'))
london.crimes <- suppressMessages(
  purrr::map(
    london.crimes.files,
    ~readr::read_csv(.) %>%
      dplyr::select(Latitude, Longitude) %>%
      dplyr::filter(!is.na(Latitude))) %>%
  magrittr::set_names(basename(Sys.glob(
    paste0(system.file('examples/data/London-Crimes', package='leaflet.extras'),
           '/2016*')))))

leaf <- leaflet() %>%
  addProviderTiles(providers$CartoDB.Positron)

purrr::walk(
  names(london.crimes),
  function(month) {
    leaf <<- leaf %>%
      addHeatmap(
        data = london.crimes[[month]],
        layerId = month, group = month,
        lng=~Longitude, lat=~Latitude,
        blur = 20, max = 0.05, radius = 15)
  })

leaf %>%
  setView(-0.094106, 51.515, 14) %>%
  addLayersControl(
    baseGroups = names(london.crimes),
    options = layersControlOptions(collapsed = FALSE)
  )

11.2 addWebGLHeatmap

案例一

Just by number of quakes

leaflet for R(2)—— learletCN

Fri, 28 Sep 2018 00:00:00 +0000

library(leaflet)
library(leafletCN)

10. LeafletCN

10.1 函数简介

regionNames：找地名函数，极其方便
demomap：有地名就给你输出一个交互、动态的区域图，真方便
geojsonMap：作用是分区块标色

leaflet for R(1)

Thu, 27 Sep 2018 00:00:00 +0000

library(leaflet)
library(leafletCN)

1. 简介

1.1 leaflet的功能

交互地图浏览（缩放、平移）
使用多种底图进行任意组合
加载地图瓦片（WMTS）
点要素定位标记
多边形要素标记
线要素标记
弹出窗口
解析加载GeoJson
从R或者RSutido创建地图窗口
可以把地图嵌入 knitr/R包所生成的Markdown文档中，或者是Shiny制作的APP中。
可以直接获取通过SP包生成就（加载）的空间对象以及包含经纬度的数据框进行展示。
可以设定地图范围以及封装自定义的鼠标事件。

1.2 leaflet的使用步骤

加载leaflet包
通过leaflet包创建地图控件。
通过图层操作的方法（如addTiles、 addMarkers、 addPolygons）来处理图层数据，并且修改地图插件的各种参数，来把图层显示在地图控件上。
可以重复第三步，可以增加更多的图层数据。
把地图部件显示出来。完成绘图。

m <- leaflet()
at <- addTiles(m)
addMarkers(at,lng = 116.391,lat=39.912,popup = "这里是北京")

管道式写法

Rstudio的Git崩溃解决方案

Thu, 27 Sep 2018 00:00:00 +0000

Rsutdio是我最喜爱的软件之一。它已经不仅仅是一个IDE了，而且能够在里面做实验、写报告、开发报告模板、写博客、做幻灯片（目前还不会）。Rstudio的Git也相当好用，简化了在Git命令行里面的操作。

读《人性的弱点》

Wed, 26 Sep 2018 00:00:00 +0000

最近闲暇时候在读《人性的弱点》一书，作者卡耐基。

这本书在现在人的观点是“成果学”、“鸡汤学”、“读了很多大道理仍然过不好一生系列”。其实我一直对现在社会主流批判“成功学”的书的风气保留意见。确实有很多粗制滥造内容浅薄的“成功学”书，但对于一些内容正统、说理透彻、案例丰富的书，为什么要批评这些书呢？不能按书上说的做，这是个人的注意力容易分散、意志力薄弱的原因，并不是书的问题。如果说现在的成功学是“鸡汤”，那么指导了两千年中华文化方向的《论语》算不算鸡汤？任何宗教的教义算不算鸡汤？“纸上得来终觉浅，绝知此事要躬行”，书能够直接帮助人的确实很有限。每个人不经历扒皮般的痛苦，苦苦想寻找一根稻草，又怎能在心平气和、气定神闲下靠读一本书就大彻大悟，从此不一样？相反，我认为人脑不是那么机械的。通过反复阅读正能量的、有启发的句子，会在脑海里留下一些光斑，在某些场景下、在经历某些事件后被点亮，这便是我们能够一点点改变自我的机会。

听《给电影人的情书》

Tue, 25 Sep 2018 00:00:00 +0000

给电影人的情书（网易云音乐）

多少人爱你遗留银幕的风采
多少人爱你遗世独立的姿态
你永远的童真，赤子的期待，孤芳自赏的无奈
谁明白你悉心隐藏的悲哀
谁了解你褪色脸上的缅怀
你天衣无缝的潇洒，心底的害怕，慢慢渗出了苍白
你苦苦地追求永恒
生活却颠簸无常遗憾
你傻傻地追求完美
却一直给误会给伤害给放弃给责备
何悲何爱，何必去愁与苦，何必笑骂恨与爱
人间不过是你寄身之处
银河里才是你灵魂地徜徉地
人间不过是你无形的梦，偶然留下的梦，尘世梦
以身外身做银亮色的梦
以身外身做梦中梦

《新京报》40书单

Mon, 24 Sep 2018 00:00:00 +0000

《新京报》举办了“大民大国40本书”（1978-2018）影响中国的40本书投票活动，希望以40本书来概括40年的精神潮流，呈现国家与私人生活的变迁。以下是40本书的书单。

分而治之的策略(3)——再思考

Mon, 24 Sep 2018 00:00:00 +0000

在之前的两篇文章中，我首先用手动版的分而治之解决了工作中遇到的问题，之后根据Hadley的论文学习了plyr的应用方法和思想。在本文中，我将探索使用更多方法来解决第一篇中的问题。

分而治之的策略(2)——plyr包

Sun, 23 Sep 2018 00:00:00 +0000

plyr包是Hadley Wickham在早期创作的一个R包，创作时间大约是2011年左右。这个包的主题思想是“The Split-Apply-Combine Strategy”，我将它翻译为“分而治之的策略”。它是指将一个大问题切分为若干个相似的、并列的小问题。如果能找到小问题的解决办法，那么将其平行的运用在其他小问题上，就解决了整个大问题。Wickham为阐述plyr的使用方法，写了一篇论文，投在《Journal of Statistical Software》开源期刊上，论文名为《The Split-Apply-Combine Strategy for Data Analysis》。本文便是这篇论文的学习笔记。

分而治之的策略(1)——手动版

Sat, 22 Sep 2018 00:00:00 +0000

昨天在工作中遇到了这样一个问题：我希望按年月，按省份，按公路等级计算断面交通量，进而计算总行驶量。

PS：以下交通量计算数据均用随机数代替。

使用我之前写的caculate_frecarsmean函数，可以调用plyr::ddply函数，这里便使用了plyr的思想：分而治之。关于plyr的思想，具体会在之后的《分而治之的策略》系列中详细写。考虑到日期必须是列的形式，公路等级必须是行的形式（因为各省公路里程是行的形式，如下表）。

矩阵乘向量

Fri, 21 Sep 2018 00:00:00 +0000

今天在工作中遇到一个问题：有一个矩阵，有一个向量，列数是相同的，需要矩阵的每行乘以向量。

例如：

A <- matrix(1:20,nrow=5)
A

##      [,1] [,2] [,3] [,4]
## [1,]    1    6   11   16
## [2,]    2    7   12   17
## [3,]    3    8   13   18
## [4,]    4    9   14   19
## [5,]    5   10   15   20

b <- c(3,4,5,6)
b

## [1] 3 4 5 6

我一开始是这么做的：

因子的顺序

Thu, 20 Sep 2018 00:00:00 +0000

因子的自然排序

今天在工作中，要处理一个之前已经遇到，但一直没有正面解决的问题：当所要分析的时间是（日度）跨月份或（月度）跨年份，该怎样操作。

关于我

Wed, 19 Sep 2018 00:00:00 +0000

这是王致远的个人主页，既用作技术博客也用作文字博客。用作技术博客主要是记录我学习R语言（或未来学习其他技术）的学习里程、心得和解决方法。文字博客是用来记录我的所读、所观、所思、所感。

创建R包

Wed, 19 Sep 2018 00:00:00 +0000

第0章：工具汇总

0.1 devtools工具汇总

devtools::has_devel() 检查开发环境是否完整；
devtools::create(“path//to//package//pkgname”) 创建包；
devtools::install() 从源码包到安装包；
devtools::build() 从源码包到压缩包，同时生成使用指南；
devtools::burld(binary=TRUE) 制作一个二进制包；
devtools::install_github() 从github上下载源码包，编译，安装，生成安装包；
devtools::load_all() 从源码包到内存中的附加包；
devtools::use_build_ignore(“notes”) 忽略特定的文件或目录；
devtools::use_package(“dplyr”) 添加dplyr到DESCRIPTION中的Imports字段，使用dplyr::fun()引用；
devtools::use_package(“dplyr”,“Suggests”) 添加dplyr到Suggests字段，使用requirenamespace(“dplyr”,quietly=TRUE)来测试包是否安装，使用dplyr::fun()引用；
devtools::document() 使用roxygen2创建对象文档；
devtools::use_vignette(“my-vignette”) 创建使用指南；
devtools::build_vignettes() 编译所有使用指南；
devtools::use_testthat() 创建testthat测试相关文件夹和文件;
devtools::test() 测试包；
devtools::use_data() 使用导出数据；
devtools::use_data_raw() 创建原始数据文件夹，用于存放原始数据和预处理数据的代码；

0.2 formatR工具

formatR::tidy_dir(“R”)

第1章：简介

1.1 为什么开发R包

共享R代码
组织自己的R代码，培养系统性思维

1.2 理念

任何可以自动化的都应该自动化

博客！出发！

Tue, 18 Sep 2018 00:00:00 +0000

今天开始在这里写博客的契机很偶然，是看到微信公众号“R语言中文网”上关于blogdown的一篇简单介绍。中午用10分钟看完后，觉得用blogdown搭建一个博客非常简单，于是乎，晚上回到宿舍尝试。