层次贝叶斯模型-空间分析

f(r |D)二Cf (D |Rf (巧 (2)

1.1层次贝叶斯模型

经典的推断分析模型、空间回归模型、空间面板模型有一个共同的特点： 这 些模型的求解完全依赖所采集的样本信息。 然而，在业务实践中，在收集样本之 前，研究者往往会对研究对象的变化或分布规律有一定的认识。 这些认识或是来 自长期积累的经验，也可能来自合理的假设。由于这些认识没有经过样本的检验， 所以我们可以称之为先验知识。比如我们要研究某地某疾病月发病人数的概率分 布。即使没有进行统计调查，我们根据一些定理和合理假设， 也可以知道发病数 服从泊松分布。甚至根据医院日常接诊的经验，可以推算出发病人数大概在哪个 区间。这种情况下，对于发病人数分布形态和大致区间的认识，属于先验知识。 先验知识对我们探索研究对象的变化规律会有很大的帮助。而经典的推断分析模 型、空间回归模型、空间面板模型都没有利用先验知识， 导致了信息利用的不充 分。而本节所要谈到的层次贝叶斯模型， 会结合先验知识和样本信息，对数据进 行推断分析。由于层次贝叶斯模型能有效利用先验知识和样本信息， 因此可以提 高推断的准确度或降低抽样的成本。

(1)贝叶斯统计原理简介

在介绍层次贝叶斯模型之前，有必要首先简单阐述一下贝叶斯统计的基本原 理。贝叶斯统计的基础是贝叶斯定理：

其中：P(A)是事件A 的先验概率(例如，某专家通过经验或之前的研究得

出乙肝发病率为10%，这就是一个先验概率)，P(B)是事件B 发生的概率，且

P(B)=O ，P(A|B)是给出事件B 后事件A 的后验概率。P(B|A)/P(B)是事件A 发生对事件B 的支持程度，即似然函数。对 P(B|A)/P(B)可以有如下的理解： 设P(B|A)/P(B)二n ，贝恠事件A 发生的条件下，事件B 发生的概率是不知A 是 否发生的条件下的n 倍。

使用贝叶斯方法的一个重要目的，就在于得出随机变量的概率分布及各因素 对分布的影响。要实现这一目的，首先按如下公式进行参数反演：

P(A|B)二 P(B | A)P(A)

P?B) (1)

其中，d 是待估参数，D 为观测数据fG)为/的先验概率密度函数， f(D|R 为已知数据D 时，参数的似然函数，C 为归一化常数。f(r|D)是参 数r 的后验概率密度函数。通过f (泊D)可以分析因素对参数分布的影响。

例1 :设X ~ N (二,1)，通过某批样本的样本信息和先验信息，得出：

~N(1 3t,1)。其中t 是时间。那么二~N(1 3t,1)就反映了时间对X 的某种影响 规律，比如从期望意义上说，时间没推进一个单位， X 就要增长3个单位。

式(2)反映出贝叶斯统计的一个重要特征，即分布的参数不再是固定值，

而是随机变量，服从某种分布。得出 fO|D)后，便可以根据如下公式得出随机 变量X 的密度函数f x (x):

f x (x)=.[冃⑺D)p x (x|d)dr

其中，Px(x|R 为d 给定的条件下X 的概率密度函数，0为参数空间。

例2:承接例1，当t = 0时，有:

的分布改为：X~N(」,；「2)，则参数向量「(七二)。 当待研究总体可以分为互相存在异质性的多个子总体时，就可以构建层次贝 叶斯模型，对变量在更高层次上的统计特征进行描述。 在层次贝叶斯模型中，-

个参数的先验分布也包含着待估参数，这样的参数被称为超参数。比如，在例

1

中,若r 的先验分布中没有待估参数(比如r ~U(0,10))，则模型中没有超参数,

模型为非层次的贝叶斯模型。若 二的先验分布中有待估参数，比如 二~U(a,b)， 则a ，b 为超参数，模型为层次贝叶斯模型。对于不同的子总体，其参数 二的先 f(r|D)= fexpl-^^]

(4)

1

(x _ 日)2 P x (x p^—— exp[

2 ] (5)

f x (x) = j=f(r |D)P x (x|R d r

(6) 注意：二未必是一维的数，有可能是多维的参数向量

比如若将例1中的X

验分布类型可以是相同的，并用超参数描述。但每个子总体的参数的分布参数未必相同，需要结合样本信息进行估计。

（2）层次贝叶斯模型的构造与举例

在时空数据中，每个时空子集内的子总体，可能有着各自不同的统计特征，并且各子总体之间还可能存在着相关性。这时可以通过层次贝叶斯模型分析待研究总体的统计特征和变化规律。层次贝叶斯模型是一种特殊的层次贝叶斯模型，这种模型考虑了数据在空间或时间（或两者兼有）上的相关性和异质性。层次贝叶斯模型的一个基本架构如下：

「A i ? B t ，it （7）

其中，i为空间标记，t为时间标记，为空间中i处，t时刻的待估参数值，弟为某种变换（如恒等变换或对数变换），为截距项。A +酗，描述空间效应，其中U i为空间相关性，V i为空间异质性。B t为时间效应，也可以分为相关性r t和异质性两部分s t，即B t = r t s t。为时空交互效应。当有些效应不明显时，可以在模型中排除相应的项。在这一基本模型之上，还可以考虑不同时空尺度的影响，以及其他协变量的影响。在模型中的诸项，都需要为其指定先验分布。先验分布的指定，依赖于已知的信息和各种模型。

例3: —个简单的层次贝叶斯模型。设乙是某地区i患有某种疾病的人数。乙的先验分布是参数为i的泊松分布。其中打二EM，E i为地区i总人口期望值，r i 为地区i该种疾病的发生率。r i的先验分布为对数正态分布，参数为叫和二2，即

In r ~ N（气，G2）。构建层次贝叶斯模型:

In *「7 V j ( 8)

其中，［为截距项，先验分布为U（0,1）。V i描述空间相关性。其先验分布指定如下：

u i ~ N（0/-2）

(9)

、W j U j 2

U|U j~N（导r ）（10）

'' W ij W j

j 4 j=1

其中，W ij为空间权重矩阵因子。这里使用了条件自回归（Conditional Auto Regressive CAR）模型。v描述空间异质性，先验分布为N（0^2）。1/ .2, 1/匚2 的先验分布都是Gamma（0.001,0.001）。

例4:具有不同空间尺度的层次贝叶斯模型。设：

X it1 ~ Poisson（N t1），X jt2 ~ PoissonUQ

其中，i, j为地点标记，t为时间标记，下标1表示较小的空间尺度（如区县），简称水平1。下标2表示较大的空间尺度（如地市），简称水平2。两种不同的空间尺度之间存在嵌套关系，即对于任意区域i，总存在区域j，使i是j的

一部分。构建层次贝叶斯模型：

ln 弓t1 A1 A j< B t1 “ "t2 （11）

i W j i W j

ln 71 jt^ ~ ：' 2 A j2 B t2，jt2 （12）

其中，：1和〉2为截距项，A1和A j2分别是水平1和水平2的空间效应。

A1 = U i1 V i1，弘=分W2。U i1和u j2分别是水平1和水平2的空间相关性，V i1 和V j2分别是水平1和水平2的空间的空间异质性。A j2是水平2作用在水平1 上的空间背景效应。A j2二U j2，

V j2，U j2和V j2分别是水平2作用在水平1上的

I刁I刁I刁I二j I二j

空间相关性和空间异质性。B t1和B t2分别是水平1和水平2的时间效应。,1和、」2 分别是

水平1和水平2上的时空交互效应。「2二U jt2 U jt2，U jt2和U jt2分别是水

i寸i寸i寸i三j i三j

平2作用于水平1的随时间变化的空间相关性和异质性。各变量的的先验分布如下：

空间相关性和时空交互项的先验分布: