李 原，王 晶

（1.山西財(cái)經(jīng)大學(xué)統(tǒng)計(jì)學(xué)院，太原 030006；2.太原師范學(xué)院經(jīng)濟(jì)與管理學(xué)院，山西 晉中 030619）

0 引言

固定參數(shù)的時間序列回歸模型只能反映某個時間段內(nèi)變量平均狀態(tài)下的關(guān)系，而現(xiàn)實(shí)中很多經(jīng)濟(jì)變量間的關(guān)系并不是一成不變的，會隨著社會的發(fā)展、技術(shù)的進(jìn)步以及經(jīng)濟(jì)結(jié)構(gòu)的轉(zhuǎn)變等而發(fā)生變化，這種動態(tài)關(guān)系就需要通過變參數(shù)模型來體現(xiàn)。傳統(tǒng)的變參數(shù)模型主要是狀態(tài)空間模型、面板數(shù)據(jù)模型以及分位數(shù)回歸、門限回歸模型等。估計(jì)變參數(shù)模型的方法一般采用卡爾曼濾波等濾波方法。

國內(nèi)外不斷有文獻(xiàn)提出變參數(shù)估計(jì)的新方法和新思路。Rajan 等（1997）[1]較早地提出了一種將馬爾可夫鏈蒙特卡羅方法應(yīng)用于時變自回歸模型的參數(shù)估計(jì)方法，該模型的時變系數(shù)由基函數(shù)建模，然后采用Gibbs方法估計(jì)參數(shù)。Orbe等（2007）[2]提出了一種非參數(shù)方法來估計(jì)看似無關(guān)的回歸方程模型中的時變系數(shù)，該方法允許將交叉和時變限制合并到參數(shù)中，估計(jì)量是以封閉形式得到的，不需要用迭代方法來計(jì)算。Escobar（2012）[3]研究了隨機(jī)系統(tǒng)在三種不同的隨機(jī)擾動下的連續(xù)時變參數(shù)估計(jì)問題。Andrea 和Lloyd（2018）[4]提出了一種通過將時變參數(shù)視為具有未知系數(shù)的分段函數(shù)來估計(jì)周期性時變參數(shù)的方法。這種方法使用非線性濾波，允許得到的參數(shù)估計(jì)在形狀上具有更大的靈活性，同時仍保持周期性。Arnold 和Lloyd（2017）[5]提出了一種利用非線性貝葉斯濾波估計(jì)周期性時變參數(shù)的方法。該方法將時變參數(shù)視為具有未知系數(shù)的分段函數(shù)，并使用集合卡爾曼濾波器（EnKF）進(jìn)行估計(jì)。

國內(nèi)也有一些文獻(xiàn)討論了變參數(shù)估計(jì)的新方法。鄧自立和解三名（1989）[6]、王建國（1990）[7]主要討論了自回歸模型的時變參數(shù)估計(jì)，采用了自適應(yīng)濾波方法。鄧衛(wèi)強(qiáng)等（2011）[8]對譜估計(jì)方法進(jìn)行了改進(jìn)，用遺傳算法進(jìn)行了時變參數(shù)估計(jì)。陳云仙等（2017）[9]也討論了時變參數(shù)的貝葉斯估計(jì)。

上述文獻(xiàn)從不同角度對時變參數(shù)的估計(jì)方法進(jìn)行了研究，每個角度都有其特點(diǎn)。本文擬從一個全新的角度，提出一種特征樣本機(jī)器采樣推斷方法來估計(jì)時變參數(shù)。

1 特征樣本概念及采樣

特征樣本是在某變量特定的值域范圍內(nèi)按照變量的分布特征采用機(jī)器采樣的方法生成的隨機(jī)樣本[10]。每個特征樣本都有目標(biāo)變量和標(biāo)識變量兩個要素，表示為（x,y），或者將x作為標(biāo)識變量獲得目標(biāo)變量y，或者將y作為標(biāo)識變量獲得目標(biāo)變量x，每個特征樣本都有特定的樣本容量，也有標(biāo)識變量和目標(biāo)變量的值域。特征樣本的一般形式可以表達(dá)為：

其中，ζ表示目標(biāo)變量的特征樣本，n是樣本大小，x與y的關(guān)系及分布特征用函數(shù)式y(tǒng)=f(x)表達(dá)，ψ[a,b]是目標(biāo)變量的值域。

特征樣本的采樣要先確定變量的分布特征，本文的分布特征不是指典型的概率分布，而是更廣泛的變量特征。在確定分布特征的基礎(chǔ)上，再確定目標(biāo)變量或標(biāo)識變量的值域，同時確定所需要的樣本容量，就可以利用特定的計(jì)算機(jī)程序來采樣。

特征樣本的采樣思想可以歸納為以下幾點(diǎn)：第一，特征樣本不是實(shí)際觀察樣本，而是依據(jù)經(jīng)驗(yàn)獲得分布特征后按照特定方法產(chǎn)生的隨機(jī)樣本；第二，特征樣本是在特定取值區(qū)間的樣本；第三，特征樣本是可以重復(fù)采樣的樣本；第四，特征樣本的分布特征，既可以通過已有經(jīng)驗(yàn)確定，也可以通過觀察現(xiàn)實(shí)數(shù)據(jù)得到；第五，特征樣本采樣總體上屬于單樣本的蒙特卡羅采樣范疇。

特征樣本可以分為分布型樣本和曲線型樣本，需要使用不同的方法進(jìn)行采樣。已有研究設(shè)計(jì)了特征樣本采樣的10種方法，包括適用于分布型樣本的標(biāo)準(zhǔn)分布法、比例分布法、反函數(shù)分布法、分區(qū)頻率分布法、分區(qū)值域分布法和適用于曲線型樣本的標(biāo)準(zhǔn)曲線法、反函數(shù)曲線法、非標(biāo)準(zhǔn)曲線法、比例曲線法、自助曲線法。這些方法可以實(shí)現(xiàn)常見的13種類型的特征樣本采樣。

按照這些采樣方法，可以在計(jì)算機(jī)程序輔助下，分別對正態(tài)分布樣本、均勻分布樣本、偏態(tài)分布樣本、線性增長樣本、線性遞減樣本、正態(tài)曲線樣本、反正態(tài)曲線樣本、偏態(tài)曲線樣本、指數(shù)增長曲線樣本、指數(shù)遞減曲線樣本、S增長曲線樣本、反S 曲線樣本、模擬曲線樣本等常見特征變量采樣。

2 基于特征樣本采樣的變參數(shù)模型估計(jì)方法

除了有固定系數(shù)的特征樣本模型估計(jì)方法，高艷平和王晶（2019）[11]還從數(shù)學(xué)及模擬仿真的角度證明了特征樣本重復(fù)抽樣回歸（FSR）方法的優(yōu)越性?；诖耍疚膶⒒谔卣鳂颖驹O(shè)計(jì)一套變參數(shù)模型估計(jì)方法，實(shí)現(xiàn)小樣本情況下的變參數(shù)模型估計(jì)。

2.1 特征樣本變參數(shù)模型估計(jì)的思想

特征樣本重復(fù)抽樣回歸方法根據(jù)對變量特征的把握來設(shè)定變量的分布形態(tài)和參數(shù)，然后按分布特征重復(fù)采樣生成一系列特征樣本。對每一個特征樣本進(jìn)行回歸都可以得到一組參數(shù)估計(jì)結(jié)果，分別對全部m個特征樣本回歸就可以得到一系列的系數(shù)估計(jì)結(jié)果，并據(jù)此模擬每個系數(shù)自身的分布特征，判斷回歸系數(shù)的取值和置信區(qū)間。FSR 方法重復(fù)抽樣得到的一系列特征樣本也為變參數(shù)模型的估計(jì)提供了一種新的思路——對特征樣本進(jìn)行重新組合，分別估計(jì)每一期的系數(shù)，可以直接得到模型的變參數(shù)結(jié)果。具體來講，F(xiàn)SR方法在時間序列的重復(fù)抽樣過程中形成了一系列特征樣本，那么如果將每個樣本的第1期匯總在一起，就得到了解釋變量和被解釋變量第1期的分布，則可以據(jù)此進(jìn)行第1 期系數(shù)的估計(jì)；后面每一期都同理，可以通過n次估計(jì)得到第1 至第n期的系數(shù)估計(jì)結(jié)果。進(jìn)一步，根據(jù)每一期系數(shù)的估計(jì)結(jié)果，還可以對系數(shù)的變動趨勢進(jìn)行自回歸估計(jì)，探究變量之間關(guān)系的變動趨勢。這種方法稱為特征樣本變參數(shù)回歸（Features Sample Variable Parameter Regression，F(xiàn)SVR）方法。

FSVR 方法是建立在特征樣本m次重復(fù)抽樣的數(shù)據(jù)基礎(chǔ)上的，但是其數(shù)據(jù)組合形式與回歸方式有所不同。FSR 方法進(jìn)行特征樣本重復(fù)抽樣的特點(diǎn)是樣本容量n極小化（大于k+1即可，k是解釋變量個數(shù)），重復(fù)抽樣次數(shù)m極大化（比如一千次、一萬次）。一次特征樣本采樣就形成了一個包含被解釋變量與k個解釋變量在內(nèi)的樣本容量為n的多變量復(fù)合特征樣本，F(xiàn)SR方法的重復(fù)采樣最終形成了m個樣本量為n的復(fù)合特征樣本。

FSVR方法的建模思想就是將FSR方法重復(fù)采樣得到的特征樣本進(jìn)行重新組合，將每一期的m次抽樣結(jié)果組合為該期的樣本，就可以得到n個樣本容量為m的分期樣本。通過對這些特征樣本的多期回歸，不僅能夠得到模型的變系數(shù)結(jié)果，而且可以根據(jù)估計(jì)結(jié)果來尋找各期系數(shù)的變動規(guī)律，得到在形式上與狀態(tài)空間模型類似的結(jié)果。構(gòu)建的主方程類似狀態(tài)空間模型中的“信號方程”或者“量測方程”。用系數(shù)構(gòu)建的自回歸模型，類似狀態(tài)空間模型中的“狀態(tài)方程”。

2.2 特征樣本變參數(shù)模型的估計(jì)方法和步驟

2.2.1 按照分布形態(tài)和變量值域進(jìn)行采樣

特征樣本變參數(shù)模型是建立在FSR 方法的基礎(chǔ)之上的，所以變參數(shù)模型的構(gòu)架和估計(jì)首先需要根據(jù)變量特征和值域進(jìn)行特征樣本的重復(fù)采樣，這一步與FSR 方法相同。對于樣本容量為n的時間序列模型，重復(fù)采樣后獲得的樣本是一個（k+1）*n行、m列的矩陣。

2.2.2 對特征樣本進(jìn)行重新組合，獲得每一期的子樣本

如果要估計(jì)變參數(shù)回歸模型，就要對初始特征樣本進(jìn)行重新組合，獲得每一期的子樣本。具體的做法是：先將上述FSR特征樣本矩陣進(jìn)行轉(zhuǎn)置，轉(zhuǎn)置后全部數(shù)據(jù)的矩陣見式（2）。

然后從該樣本中分別提取第t期的數(shù)據(jù)，組合成一個新的m行、k+1 列的矩陣，即為第t期的子樣本。其中y為被解釋變量，x為解釋變量，解釋變量個數(shù)為k，每個變量樣本容量為n，抽樣次數(shù)為m。每期的特征樣本矩陣見式（3），一共可以得到n個樣本容量為m的子樣本。

2.2.3 分期回歸，獲得每一期的估計(jì)系數(shù)

根據(jù)重復(fù)抽樣和組合獲取的子樣本，可以對每一個子樣本進(jìn)行單方程參數(shù)估計(jì)，獲得每一期的系數(shù)估計(jì)結(jié)果，類似狀態(tài)空間的“信號方程”?！靶盘柗匠獭钡哪Ｐ托问娇梢赃x擇線性回歸、非線性回歸等形式，但是全部n期的“信號方程”形式必須一致。以最常見的多元線性回歸為例，“信號方程”的形式為：

選定模型形式之后，用最小二乘法或極大似然法來估計(jì)每一期子樣本的回歸系數(shù)，其中為第t期回歸模型的常數(shù)項(xiàng)，為第t期回歸模型中第i個變量的估計(jì)系數(shù)，i=1,2,3,…,k。

需要注意的是，在FSR 方法的應(yīng)用中，很多變量是遞增或者遞減的分布形態(tài)，所以，在m次抽樣中，其初期和末期的變量取值實(shí)際上都是固定的最大值或者最小值，是無法進(jìn)行系數(shù)估計(jì)的常數(shù)?？紤]到變參數(shù)模型的普遍適用性，只進(jìn)行第2 期到第n-1 期的分期系數(shù)估計(jì)，第1 期和第n期系數(shù)用下文給出的方法推算。

2.2.4 估計(jì)系數(shù)的波動規(guī)律并據(jù)此推算期初和期末的系數(shù)

在獲得第2至第n-1期變系數(shù)的基礎(chǔ)上，對單個系數(shù)序列進(jìn)行自回歸，可以估計(jì)其變化規(guī)律，類似狀態(tài)空間模型中的“狀態(tài)方程”。

設(shè)各期系數(shù)之間的關(guān)系為：

可以根據(jù)每個參數(shù)的估計(jì)結(jié)果進(jìn)行自回歸估計(jì)，這樣就得到了各期系數(shù)的滯后一期的變化規(guī)律。根據(jù)式（5）的估計(jì)結(jié)果，可以進(jìn)行期初和期末系數(shù)的推算，即：

如圖1所示，路由開銷隨著節(jié)點(diǎn)停留時間的延長而降低，這是因?yàn)榫W(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)的變化不再頻繁所致。圖2表明網(wǎng)絡(luò)整體的端到端時延隨節(jié)點(diǎn)停留時間的增加而降低。圖3表明分組投遞率隨節(jié)點(diǎn)停留時間的增加而提高。如圖4所示，路由發(fā)現(xiàn)頻率隨節(jié)點(diǎn)停留時間的增加而降低。

2.2.5 特征樣本變參數(shù)模型的最終表達(dá)式

經(jīng)過以上步驟就得到了變量之間每一期的回歸系數(shù)，以及系數(shù)隨時間波動的自回歸方程，實(shí)現(xiàn)了時間序列數(shù)據(jù)的變參數(shù)估計(jì)。變參數(shù)模型的最終表達(dá)式為：

2.2.6 估計(jì)特征樣本整體模型

為了反映模型的整體擬合情況，對模型系數(shù)進(jìn)行分期估計(jì)以后，還可以對模型進(jìn)行不分期的整體估計(jì)。模型整體估計(jì)的結(jié)果是固定系數(shù)，實(shí)際上是分期系數(shù)的均值。

對模型中每個變量每次采樣的特征樣本求均值，就得到了一個包含每個變量n期均值、樣本容量為采樣次數(shù)m的新樣本，根據(jù)該樣本進(jìn)行最小二乘估計(jì)，可以得到每個變量均值的估計(jì)系數(shù)，即固定系數(shù)。

2.2.7 參數(shù)估計(jì)結(jié)果的檢驗(yàn)

（1）模型整體擬合優(yōu)度檢驗(yàn)

特征樣本變參數(shù)模型擬合優(yōu)度檢驗(yàn)的思路與普通最小二乘估計(jì)一致，都是考察殘差平方和在被解釋變量離差平方和中的比重，這個比重越小，說明模型的擬合優(yōu)度越好。把變參數(shù)回歸中n個回歸方程作為一個整體時，可以得到整體殘差平方和（TSSE）與被解釋變量的整體離差平方和（TSST）。特征樣本變參數(shù)模型的擬合優(yōu)度R2用整體被解釋變量離差平方和中除去整體殘差平方和的部分所占的比重來表示，R2越大，表示變系數(shù)模型的整體擬合優(yōu)度越好。

（2）模型整體固定系數(shù)顯著性水平檢驗(yàn)

對于按變量均值估計(jì)的固定系數(shù)，采用t檢驗(yàn)方法，可獲得對模型固定系數(shù)估計(jì)的可信度評價。

（3）各期可變系數(shù)顯著性水平檢驗(yàn)

一個包含k個自變量的變參數(shù)線性回歸模型每期有k+1個系數(shù)。對每個系數(shù)進(jìn)行一次t檢驗(yàn)，可以獲得t值和對應(yīng)的P 值，據(jù)此可以對系數(shù)估計(jì)結(jié)果進(jìn)行判斷和評價。t 檢驗(yàn)系數(shù)的標(biāo)準(zhǔn)誤采用每期估計(jì)的殘差和對應(yīng)的方差-協(xié)方差矩陣來估計(jì)。

3 應(yīng)用案例

根據(jù)經(jīng)濟(jì)分析的需要，建立一個有關(guān)居民消費(fèi)價格指數(shù)（cpi）的時變參數(shù)模型，解釋變量為人均消費(fèi)支出（income）與貨幣供應(yīng)量（m），建立如下變參數(shù)模型：

其中，αt和βt分別是各期截距項(xiàng)和變量項(xiàng)可變系數(shù)。分析時期是從2011 年到2020 年的10 年。在分析中根據(jù)經(jīng)驗(yàn)可以知道該時期的變量特征，故采用特征樣本參數(shù)估計(jì)方法。各個變量的分布特征如表1所示。

表1 變量分布特征

依據(jù)變量分布特征對每個變量進(jìn)行特征樣本采樣，在采樣程序中輸入解釋變量數(shù)k=2，樣本期n=10，采樣次數(shù)m=500，得到500 個樣本。按前述步驟和程序進(jìn)行運(yùn)算，設(shè)定cpi的凸函數(shù)參數(shù)為0.8，m的凹函數(shù)參數(shù)為1.2，參數(shù)估計(jì)和檢驗(yàn)后得到結(jié)果如表2所示。

表2 居民消費(fèi)價格模型可變系數(shù)估計(jì)及t檢驗(yàn)結(jié)果

各可變系數(shù)的自回歸方程為：

整體均值固定系數(shù)的方程及其檢驗(yàn)結(jié)果為：

模型的分期系數(shù)和固定系數(shù)的t 檢驗(yàn)結(jié)果均顯著，擬合優(yōu)度R2為0.7494，擬合效果較好，模型整體估計(jì)結(jié)果有效。

根據(jù)式（10）的估計(jì)結(jié)果，人均消費(fèi)支出對居民消費(fèi)價格指數(shù)的負(fù)面影響是逐漸縮小的，而貨幣供應(yīng)量對居民消費(fèi)價格指數(shù)的影響由負(fù)向轉(zhuǎn)為正向，并且正向影響隨時間推移而增強(qiáng)。與固定參數(shù)模型相比，這一估計(jì)結(jié)果不僅體現(xiàn)了人均消費(fèi)支出和貨幣供應(yīng)量對居民消費(fèi)價格指數(shù)的影響，而且更加精確地反映了這種影響隨時間變動的情況。

4 結(jié)論

本文提出了一種基于特征樣本的變參數(shù)模型估計(jì)（FSVR）方法，通過對重復(fù)抽樣的每一期樣本分別回歸來實(shí)現(xiàn)時間序列模型的變參數(shù)估計(jì)。得出以下結(jié)論：（1）FSVR方法是一種新的方法，該方法簡便易行，給出了變參數(shù)模型參數(shù)估計(jì)的另外一種思路。在社會科學(xué)研究中，很多時候存在數(shù)據(jù)短缺問題，有了特征樣本采樣方法，就可以根據(jù)經(jīng)驗(yàn)來獲取樣本，依據(jù)特定的程序，把可變系數(shù)估計(jì)出來。（2）FSVR 方法符合貝葉斯統(tǒng)計(jì)原理，在進(jìn)行參數(shù)估計(jì)之前，依據(jù)經(jīng)驗(yàn)給出參數(shù)的先驗(yàn)分布，然后將基于先驗(yàn)信息的特征樣本納入?yún)?shù)估計(jì)過程，這樣就能夠把握大方向和減少偏誤。（3）利用FSVR 方法估計(jì)模型的可變參數(shù)，是基于某種分布通過計(jì)算機(jī)重復(fù)采樣對參數(shù)進(jìn)行模擬，符合蒙特卡羅模擬原理。（4）FSVR 方法將按時期長度得到的樣本容量轉(zhuǎn)換成按采樣次數(shù)計(jì)算的樣本容量，可以將小樣本估計(jì)變成大樣本估計(jì)。（5）FSVR方法估計(jì)的可變參數(shù)需要進(jìn)行兩個層次的檢驗(yàn)。在分期層次上，設(shè)計(jì)了系數(shù)的t檢驗(yàn)；在整體層次上，設(shè)計(jì)了整體固定系數(shù)檢驗(yàn)和擬合優(yōu)度檢驗(yàn)，模擬擬合效果和參數(shù)估計(jì)效果優(yōu)劣都可以通過這些檢驗(yàn)來完成。

本文應(yīng)用實(shí)例的估計(jì)結(jié)果顯示，該方法可以得到較好的變參數(shù)估計(jì)結(jié)果，F(xiàn)SVR 方法適用于社會科學(xué)領(lǐng)域小樣本變參數(shù)估計(jì)。當(dāng)然，本文僅對FSVR 方法的原理、檢驗(yàn)和應(yīng)用進(jìn)行了介紹，后續(xù)還需要對估計(jì)殘差的檢驗(yàn)、與狀態(tài)空間模型對比的優(yōu)劣及估計(jì)結(jié)果的穩(wěn)健性等方面進(jìn)行進(jìn)一步的研究。