基于Pair-Copula函數(shù)的空間預(yù)測模型及其應(yīng)用*

郭益敏, 楊煒明

(重慶工商大學(xué) 數(shù)學(xué)與統(tǒng)計學(xué)院，重慶 400067)

基于Pair-Copula函數(shù)的空間預(yù)測模型及其應(yīng)用*

郭益敏, 楊煒明**

(重慶工商大學(xué) 數(shù)學(xué)與統(tǒng)計學(xué)院，重慶 400067)

對于空間數(shù)據(jù)的插值預(yù)測，大多采用傳統(tǒng)的空間插值方法如反距離加權(quán)插值法和克里金插值法，這2種方法在邊緣分布或存在異常值的情況下會導(dǎo)致預(yù)測精度相對較低；采用基于Copula理論的方法克服了這一問題。通過Pair-Copula函數(shù)描述了空間相依結(jié)構(gòu)并利用MCMC方法(貝葉斯估計法)估計參數(shù)，討論基于空間數(shù)據(jù)對未觀測位置相關(guān)數(shù)據(jù)進(jìn)行了空間插值預(yù)測；結(jié)合重慶市霧霾數(shù)據(jù)對該方法與反距離加權(quán)插值法、普通克里金和泛克里金插值法進(jìn)行比較，結(jié)果發(fā)現(xiàn)基于Pair-Copula函數(shù)的空間預(yù)測模型具有更高的精度。

Pair-Copula函數(shù);空間相依結(jié)構(gòu);空間預(yù)測

Copula理論自1959年Sklar[1]首次提出以來即得到廣泛應(yīng)用[2-4]，隨后考慮到兩兩變量間相關(guān)結(jié)構(gòu)的差異性，Bedford和Cooke[5-6]在Joe[7]的研究基礎(chǔ)上，引入Vine這一圖解模型，將多元Copula函數(shù)分解為多個條件或非條件二元Copula函數(shù)的乘積，即Pair-Copula結(jié)構(gòu)(PCC)。Pair-Copula模型能更好地反映復(fù)雜事物之間的相關(guān)關(guān)系，許多學(xué)者在這方面做了大量應(yīng)用研究[8-9]。近年來快速發(fā)展的Pair-Copula理論在空間領(lǐng)域的應(yīng)用越來越多，也驗(yàn)證了結(jié)合空間信息的方法較非空間方法更準(zhǔn)確省時[10]。

隨著信息技術(shù)的高速發(fā)展和空間數(shù)據(jù)結(jié)構(gòu)的日益豐富，如何利用Copula理論構(gòu)造空間相依結(jié)構(gòu)以達(dá)到空間插值預(yù)測的目的成為空間數(shù)據(jù)分析建模的難點(diǎn)?？臻g統(tǒng)計分析中，通常利用變差函數(shù)描述空間相關(guān)性并用Kriging插值法進(jìn)行預(yù)測，但這兩種方法對相距較遠(yuǎn)的觀測點(diǎn)(即異常值)比較敏感且易受到隨機(jī)場邊緣分布的影響。因此，本文利用Pair- Copula函數(shù)描述插值的空間依賴性結(jié)構(gòu)并將其與邊緣信息結(jié)合建立聯(lián)合分布，運(yùn)用MCMC(貝葉斯)法估計參數(shù)，結(jié)合重慶市霧霾觀測站的實(shí)際數(shù)據(jù)進(jìn)行空間插值預(yù)測，并與傳統(tǒng)插值方法，如反距離加權(quán)插值法(IDW)、普通克里金(OK)和泛克里金(UK)插值法進(jìn)行比較。

1 理論模型

假設(shè)所有聯(lián)合、邊緣和條件分布連續(xù)且有相應(yīng)的密度函數(shù)，Pair-Copula結(jié)構(gòu)可將d維聯(lián)合密度分解為d個邊緣密度與d(d-1)/2個二元Copula密度的乘積，是一種靈活的工具。Bedford和Cooke定義的Vine包括兩類常用形式，即D-Vine和C-Vine，本文使用C-Vine對聯(lián)合分布進(jìn)行分解。

1.1 C-Vine Copula

C-Vine呈輻射狀結(jié)構(gòu)，每棵樹有一個根節(jié)點(diǎn)與其他節(jié)點(diǎn)相連。以5維C-Vine為例，它包含4棵樹Tj,j=1,2,3,4，樹Tj有6-j個節(jié)點(diǎn)和5-j條邊，樹Tj的節(jié)點(diǎn)是Tj+1的邊，每條邊及其標(biāo)記對應(yīng)一個Pair-Copula函數(shù)和該函數(shù)的下標(biāo)，例如，邊14|23對應(yīng)Pair-Copula密度c14|23。

考慮d維隨機(jī)向量(Y1,…,Yd)，其聯(lián)合分布和密度分別為F1:d和f1:d，1:d={1,…,d}，則基于C-Vine的PCC結(jié)構(gòu)為

(1)

1.2 空間Pair-Copula模型

在隨機(jī)場{Z(x),x∈S}中，S是研究區(qū)域，x=(x1,…,xd)′表示研究區(qū)域中不同的觀測位置，Z(x)是在位置x處的隨機(jī)變量，其取值是隨機(jī)變量在該點(diǎn)處的一個隨機(jī)實(shí)現(xiàn)。在一定空間范圍內(nèi)，不同空間位置的研究主體都存在著一定的相關(guān)性，但當(dāng)兩監(jiān)測站間相距超過一定距離時，則認(rèn)為它們是獨(dú)立的。本文選擇式(2)來描述空間相依結(jié)構(gòu)[11]：

(2)

其中,h為滯后距，a為變程。

多數(shù)情況下，邊緣分布函數(shù)的均值函數(shù)存在空間趨勢和自相關(guān)，因而本文選擇經(jīng)緯度和滯后變量作為其解釋變量，建立下面的自回歸模型：

(3)

考慮邊緣分布連續(xù)情況下的Gassian copula函數(shù)。將式(2)、式(3)代入(將空間相關(guān)函數(shù)和空間趨勢引入)似然函數(shù)有：

(4)

1.3 空間插值預(yù)測

2 實(shí)證分析與結(jié)論

2.1 實(shí)證分析

利用重慶市17個霧霾監(jiān)測站的經(jīng)緯度等地理信息和霧霾PM 2.5數(shù)據(jù)，整理樣本，構(gòu)成樣本點(diǎn)數(shù)據(jù)集合。運(yùn)用R軟件中的SP程序包計算各觀測站間的距離(觀測站與自身的距離為0)，并根據(jù)每個觀測站與其余所有站點(diǎn)的距離平均值排序確定C-Vine的根節(jié)點(diǎn)，即平均距離最小的作為每棵樹的根節(jié)點(diǎn)。為了檢驗(yàn)空間Pair-Copula模型的預(yù)測精度，從17個觀測站中隨機(jī)抽取5個點(diǎn)作為檢驗(yàn)點(diǎn)，共抽取50組，利用交叉驗(yàn)證法對邊緣分布為正態(tài)分布或截斷正態(tài)分布的剩余12個觀測點(diǎn)的樣本數(shù)據(jù)進(jìn)行建模，計算檢驗(yàn)點(diǎn)的真實(shí)值與預(yù)測值間的均方根誤差，并與反距離加權(quán)插值法、普通克里金和泛克里金插值法進(jìn)行比較。這里選取其中的兩天進(jìn)行交叉驗(yàn)證對比分析。

其中，2016-07-11與2016-07-15兩天的均方根誤差值相差較大，這是由于2016-07-15的樣本數(shù)據(jù)中含有缺失值。由表1知，包含缺失值或不含缺失值的50組交叉驗(yàn)證中，Pair-Copula空間模型插值預(yù)測結(jié)果的均方根誤差均小于其余3種傳統(tǒng)空間插值方法，說明該模型在預(yù)測精度上更優(yōu)。此外，反距離加權(quán)插值法與普通克里金法的精度相差不大，泛克里金法精度相對最低，其中普通克里金插值法的精度在總體上較反距離插值法和泛克里金插值法更優(yōu)。交叉驗(yàn)證結(jié)果表明，基于Pair-Copula函數(shù)構(gòu)建空間模型對重慶市主城區(qū)觀測站的霧霾數(shù)據(jù)進(jìn)行空間插值預(yù)測，較傳統(tǒng)插值方法(IDW、OK、UK)相比能夠達(dá)到更高的精度。

表1 模型預(yù)測精度交叉驗(yàn)證結(jié)果(均方根誤差值)

2.2 結(jié) 論

在空間數(shù)據(jù)分析中，研究樣本如霧霾濃度受多方面因素影響且存在異常值，因此傳統(tǒng)克里金插值法并不適用?？臻gPair-Copula模型利用PCC靈活地描述空間依賴性關(guān)系與空間相依結(jié)構(gòu)，進(jìn)而通過構(gòu)造聯(lián)合分布函數(shù)進(jìn)行建模，克服了利用變差函數(shù)表述空間相關(guān)結(jié)構(gòu)的不足，并且將邊緣信息與影響因素有效地結(jié)合在一起，使得空間插值結(jié)果更加準(zhǔn)確。本文僅對正態(tài)數(shù)據(jù)進(jìn)行建模分析，而對非正態(tài)數(shù)據(jù)建模時更為復(fù)雜，這種情況還有待進(jìn)一步研究。

[1] SKLAR A. Fonctions De RéPartition à N Dimensions Et Leurs Marges[J]. Publication de L’Institut de Statistique de L’Universit de Paris,1959(8):229-231

[2] 張堯庭.連接函數(shù)(Copula)技術(shù)與金融風(fēng)險分析[J].統(tǒng)計研究,2002(4):48-51

ZHANG Y T.Copula Technique and Financial Risk Analysis[J].Statistical Research,2002(4):48-51

[3] BEATRIZ V,RAFAEL M.Measuring Financial Risks with Copulas[J].International Review of Financial Analysis,2004,13(1):27-45

[4] SHIAU J T.Fitting Drought and Severity with Two-Dimensional Copulas[J].Water Resources Management,2006,20(5):795-815

[5] BEDFORD T,COOKE R.Probability Density Decomposition for Conditionally Dependent Random Variables Modeled by Vines[J].Annals of Mathematics and Artificial Intell-igence,2001,32(1):245-268

[6] BEDFORD T,Cooke R.Vines:A New Graphical Model for Dependent Random Variables[J].Annals of Statistics,2002,30(4):1031-1068

[7] JOE H.Families of M-variate Distributions with Given Margins and M(M-1)/2 Bivariate Dependence Parameters[J].Distributions with Fixed Marginals and Related Topics:Lecture Notes—Monograph Series,1996(28):120-141

[8] BERG D,AAS K.Models for Construction of Higher-Dimensional Dependence:A Comparison Study[J].European Journal of Finance,2009,15(7-8):639-659

[9] 楊茂靈,王龍,余航,等.基于Pair-Copula函數(shù)和標(biāo)準(zhǔn)化徑流指數(shù)的水文干旱頻率分析：以南盤江流域?yàn)槔齕J].長江流域資源與環(huán)境,2014,23(9):1315-1321

YANG M L,WANG L,YU H,et al.Analysis of Hydrological Drought Frequency Based on the Pair-Copula and Standardized Streamflow Index:Using NaPan River Basin as An Example[J].Resources and Environment in the Yangtze Basin,2014,23(9):1315-1321

[10] ERHARDT T M,CZADO C,SCHEPSMEIER U.Spatial Composite Likelihood Inference Using Local C-Vines[J].Journal of Multivariate Analysis,2015(138):74-88

[11] DIGGLE P,Ribeiro P.Model-based Geostatistics[M].New York:Springer,2007

責(zé)任編輯：李翠薇

Spatial Prediction Model and Its Application Based on Pair-Copula Functions

GUO Yi-min, YANG Wei-ming

(School of Mathematics and Statistics, Chongqing Technology and Business University, Chongqing 400067, China)

The interpolation prediction of spatial data usually uses traditional spatial interpolation methods such as inverse distance weighted interpolation and Kriging interpolation, whose prediction accuracy is relatively low under the impact of marginal distribution or outlier, as a result, the method based on copula overcomes the problem. Spatial correlation structures are described by Pair-Copula function and the parameters are estimated, and spatial interpolation prediction method is discussed in corresponding values of none-observation stations based on spatial data. This model is compared with inverse distance weighted interpolation, original Kriging interpolation and universal Kriging interpolation based on the data of fog in Chongqing, and the results show that the spatial prediction model based on Pair-Copula function posses the higher accuracy.

Pair-Copula function; spatial correlation structure; spatial prediction

10.16055/j.issn.1672-058X.2017.0003.004

2016-11-10；

2016-12-10. * 基金項目：國家社會科學(xué)基金(13CTJ016); 重慶市教委科學(xué)技術(shù)研究項目(KJ1600610)．

郭益敏(1993-)，女，重慶人，碩士研究生，從事空間統(tǒng)計和空間計量研究.

**通訊作者：楊煒明(1981-)，男，湖南人，副教授，博士，從事縱向數(shù)據(jù)建模和空間統(tǒng)計研究．E-mail：ywmctbu@hotmail.com.

F224

A

1672-058X(2017)03-0017-04