该模型基于坐标变换。基本假设是:不同地区的整个疫情进程整体形态相似,只是曲线的高低长短形态有所不同。所以我们通过一个已经经历整个疫情进程地区的曲线,进行坐标映射和尺度放缩得到其它地区的曲线。
形式化表示即为,设已知完整进程曲线f(x)
(下称基准函数),引入坐标变换函数g
和尺度变换函数h
,得新曲线h(f(g(x')))
,x'
为新地区的自变量(我们的例子中为天数)。
当前基准函数已知,需要求g
与h
的参数。这里我们使用反向传播算法[1],希望基准函数以最小的变换代价贴合所求地区的已知曲线,这个变换通过g
与h
完成,达成最小代价即代表着二者取到最优参数。这样,在所得曲线贴合当前地区已知曲线的同时,通过基准函数提供的先验形态,该地区样本外的曲线值也可以得到估计。
因为要使用反向传播算法,所以必须得到f
、g
与h
三者可微的表示。f(x)
可以通过任意函数拟合的方式求得,三次样条插值[2]是我们推荐的手段。g
与h
因为需要在反向传播的过程中调整,因此这里我们使用两个简单的神经网络模型来表示。为了方便编程,基准函数也就直接使用神经网络拟合得到。如果你使用其它插值方法求得基准函数的解析式,那么可以使用keras的自定义层来替换我们例子中表示基准函数的网络。
我们的例子模型示意图如下:
需要注意的是g
与h
的函数原型。h
表示对基准函数的尺度变换,因此需要控制其表达能力,我们使用了简单的单线性层。如果将h
设计的过于复杂,反向传播过程将直接通过h
函数拟合当前地区的形态,而不是试图对基准函数进行变换(因为这相比直接调整h
代价更大)。g
的原型取决于你所期待的坐标变换形式,在我们的例子是一个简单的时间-人数
的一维函数,那么从x轴来说,我只希望它进行左右平移和放缩,即一维仿射变换。所以在设计中只尝试使用数个带偏置的线性层组合(在后来因为效果不理想,克制地更换了一些relu)。如果觉得你所要达到的效果有必要使用更复杂的映射函数,那么可以酌情进行修改。不过值得庆幸的是,即使g
相对较深也不会像h
那样对结果直接造成严重影响,因为g
只决定了f
的自变量,最终的曲线值还是由f
计算得到。
我们的例子只针对了一维函数,那么在高维情况下这个方法还能否适用呢?如添加多个与传染病相关的时序变量时?首先,坐标和尺度变换在高维情况下当然适用,但需要注意的是,二者都是针对高维图形整体。如果你的指标仅关注这个高维函数的某个截面(如添加时序变量时),那么该模型将无法直接套用。
举个例子,假设你希望研究x
、y
坐标与函数值的关系(即感染人数随区域的变化,假设二者有函数关系),并希望对样本外的x
、y
区间的函数值进行预测,那么你将会关注的是基准函数f(x,y)
所描述的三维图形整体。在这种情况下,坐标变换和尺度变换依然能有很好的表现,模型也可以使用。但假如你要引入医疗资源(随时间变化)作为外生变量,希望使用二者对未来的感染人数进行预测,这将无法得到好的结果。
这是为什么呢?
30×30
网格内的x
、y
值已知,那么我们将通过900个数据来了解整个区间的分布。y=g(x)
(表示医疗资源随时间变化)的曲面截基准函数f(x,y)
得到的截面高,而不是基准函数整体虽然当我们知道f(x,y)
整体的分布时,当然可以知道其任意一个截面的信息。但我们所拥有的数据全部都在一个截面上(30个)而不是均匀分布在整个二维区间内(900个),自然就没有足够的信息对整个基准函数进行建模。既然无法保证基准函数的准确性,坐标变换与尺度变换也就失去了意义。
如果需要引入时序变量,怎样对这个模型进行改进呢?
首先,我们认为,其它时序变量带来的影响即为当它们取不同值时曲线会展现出不同的形态,通过我们对疫情数据的分析,曲线具有几种典型形态,大部分地区的疫情进程都可以落在这几种形态之内。
因此我们考虑的方案是,分离时间变量与其它时序变量。首先通过已知的时序变量对数据进行分类,判断其对应哪种形态,之后再使用对应形态的模型采取上面所说的方案进行训练和预测。这样可以将多维回归问题转化成小规模的分类和简单回归,有效利用已知信息,提高训练的稳定性。
一个也不能少
check.py
:画出指定城市所有图像to_submission.py
:整合所有城市的结果为提交表格形式通用版本.py
:基于基准函数微调生成指定城市每个区域的模型通用版本v2.py
:调整h
与g
的结构,解决一些因已知数据与基准函数形态不同导致拟合不上的问题(属于基于实践做的改进)一号网络.py
:基于完整进程的数据训练基准函数模型(神经网络)toy_v5_s.h5
:基准函数模型(神经网络)toy_v5_xr.h5
:适合大型地区的基准函数toy_v5_zt.h5
:适合小型地区的基准函数此处可能存在不合适展示的内容,页面不予展示。您可通过相关编辑功能自查并修改。
如您确认内容无涉及 不当用语 / 纯广告导流 / 暴力 / 低俗色情 / 侵权 / 盗版 / 虚假 / 无价值内容或违法国家有关法律法规的内容,可点击提交进行申诉,我们将尽快为您处理。
1. 开源生态
2. 协作、人、软件
3. 评估模型