于明星

(朝陽(yáng)師范高等?？茖W(xué)校 信息工程系，遼寧朝陽(yáng)122000)

0 引言

隨著我國(guó)經(jīng)濟(jì)的迅速發(fā)展，高職高專階段的生源數(shù)量顯著提高，高校的招生數(shù)量也隨之升高.高職高專教育已由穩(wěn)定發(fā)展期過(guò)渡到新探索時(shí)期，學(xué)生數(shù)量的持續(xù)增加和畢業(yè)生質(zhì)量的下滑受到社會(huì)的普遍關(guān)注.為了實(shí)現(xiàn)教育大眾化、現(xiàn)代化、公平化，教育部頒布了系列文件[1],其中2004年印發(fā)的《普通高等學(xué)?；巨k學(xué)條件指標(biāo)》中就招生規(guī)模及招生條件等多項(xiàng)指標(biāo)作出了規(guī)定，明確指出在校生人數(shù)與專任教師比率不應(yīng)超過(guò)18∶1.生師比是反映高職高專院校辦學(xué)質(zhì)量的一項(xiàng)重要指標(biāo)[2]，是高等教育大眾化進(jìn)程中質(zhì)量建設(shè)的基礎(chǔ)，是一流大學(xué)守望精英教育質(zhì)量的傳統(tǒng)，是高等教育大眾化進(jìn)程中大學(xué)質(zhì)量改革的理性判斷[3].高職高專院校在積極擴(kuò)招的態(tài)勢(shì)下，在校生數(shù)量急速增長(zhǎng)，但專任教師數(shù)量增幅卻跟不上節(jié)奏，嚴(yán)重影響了高職高專教育的教學(xué)質(zhì)量.通常生師比比值與教學(xué)質(zhì)量具有反比關(guān)系，文獻(xiàn)[4]通過(guò)分析2003～2013年間全國(guó)普通高校生師比的數(shù)據(jù)發(fā)現(xiàn)，國(guó)家重點(diǎn)建設(shè)大學(xué)及985工程大學(xué)的生師比較低，并針對(duì)高校生師比現(xiàn)狀，提出了一些改善生師比的對(duì)策.文獻(xiàn)[5]利用遺傳算法優(yōu)化的BP神經(jīng)網(wǎng)絡(luò)模型對(duì)未來(lái)幾年普通高校的生師比進(jìn)行了預(yù)測(cè).文獻(xiàn)[6]利用數(shù)據(jù)平臺(tái)建立了高職院校生師比狀態(tài)函數(shù)，進(jìn)而通過(guò)數(shù)據(jù)函數(shù)值來(lái)分析院校的辦學(xué)質(zhì)量及辦學(xué)規(guī)模等.

目前，已有較多研究關(guān)注生師比這一指標(biāo)，大多集中于本科院校的現(xiàn)狀及對(duì)策上的分析，而對(duì)高職高專院校生師比指標(biāo)的預(yù)測(cè)研究相對(duì)不足，這影響著我國(guó)高職高專教育的招生及可持續(xù)性的健康發(fā)展.本研究中通過(guò)分析教育部高職高專院校的生師比數(shù)據(jù)，結(jié)合自相關(guān)系數(shù)表和偏自相關(guān)系數(shù)表來(lái)判斷初始數(shù)據(jù)的平穩(wěn)性，分析出自回歸分量階數(shù)p和移動(dòng)平均分量階數(shù)q，進(jìn)而確定相應(yīng)的ARIMA(p,d,q)模型，通過(guò)該模型來(lái)預(yù)測(cè)未來(lái)幾年高職高專院校生師比指標(biāo).

1 ARIMA數(shù)學(xué)模型

ARIMA(全稱Auto Regressive Integrated Moving Average)模型，即自回歸移動(dòng)平均模型，由Box-Jenkins于20世紀(jì)70年代提出的是一種時(shí)間序列分析法[7]，其基本思路是對(duì)初始數(shù)據(jù)進(jìn)行時(shí)間上的階數(shù)平移而形成一個(gè)隨機(jī)序列，用以描述這個(gè)隨機(jī)序列的屬性模型就是ARIMA模型.該模型可以利用時(shí)間序列的過(guò)去值去分析未來(lái)值，在預(yù)測(cè)農(nóng)產(chǎn)品價(jià)格[8]、交通事故[9]、網(wǎng)絡(luò)流量[10]和電氣負(fù)荷[11]等多方面具有廣泛的應(yīng)用.ARIMA模型中有三個(gè)參量，分別為自回歸分量階數(shù)p、差分次數(shù)d和移動(dòng)平均分量階數(shù)q，通常用ARIMA(p,d,q)表示[7].

結(jié)合高職高專院校不同研究數(shù)據(jù)的特點(diǎn)，構(gòu)建最優(yōu)的ARIMA預(yù)測(cè)模型.定義研究數(shù)據(jù)觀測(cè)值z(mì)t滿足

zt=λ1zt-1+λ2zt-2+λ3zt-3+…+λpzt-p+vt

(1)

式中：λt—— 回歸參數(shù)，其中i=1，2，…，p為滯后變量數(shù)；vt—— 白噪聲過(guò)程.則線性數(shù)據(jù)觀測(cè)值z(mì)t就是p階自回歸模型，表示為AR(p).

白噪聲vt用滯后算子表示為

vt=Λ(L)zt=(1-λ1L-λ2L2-…-λpLp)zt

(2)

式中：Λ(L) —— 自回歸算子.

自回歸算子變式為

(3)

當(dāng)特征方程滿足Λ(L)=0時(shí)，AR模型在p階平穩(wěn).

若研究數(shù)據(jù)觀測(cè)值z(mì)t滿足

zt=vt+θ1vt-1+θ2vt-2+θ3vt-3+…+θqvt-q

(4)

式中：θ1,θ2,…，θq—— 參數(shù)；vt-q——t-q時(shí)所對(duì)應(yīng)的白噪聲.研究數(shù)據(jù)觀測(cè)值z(mì)t就是q階移動(dòng)平均模型，表示為MA(q).

式(3)變形為

zt=Θ(L)vt=(1+θ1L+θ2L2+…+θqLq)vt

(5)

式中：Θ(L)——移動(dòng)平均算子.

移動(dòng)平均算子特征方程為

Θ(L)=1+θ1L+θ2L2+…+θqLq=0

(6)

移動(dòng)平均算子變式為

(7)

那么,研究數(shù)據(jù)觀測(cè)值

(8)

式中：k1,k2,…，kq——常數(shù).

當(dāng)特征方程滿足Θ(L)=0，MA模型在q階可逆.

ARMA模型由AR模型和MA模型組合構(gòu)成，表達(dá)式為

zt=λ1zt-1+λ2zt-2+λ3zt-3+…+λpzt-p+vt+θ1vt-1+θ2vt-2+θ3vt-3+…+θqvt-q

(9)

綜合式(2)和式(5)，ARMA模型可變形為

Λ(L)zt=Θ(L)vt

(10)

若時(shí)間序列不具有平穩(wěn)性，則需要對(duì)不平穩(wěn)模型進(jìn)行差分處理，那么此ARMA模型就是ARIMA模型.

為了分析數(shù)據(jù)的變化率情況，引入增比指標(biāo)，其表達(dá)式為

(11)

式中：zi——第i年的數(shù)據(jù)；zi-1——第i-1年的數(shù)據(jù).

其基本過(guò)程為：

(1)根據(jù)初始高職高專院校教育統(tǒng)計(jì)數(shù)據(jù)的自相關(guān)系數(shù)表和偏自相關(guān)系數(shù)表判定序列的平穩(wěn)性，通常初始的數(shù)據(jù)為不平穩(wěn)的時(shí)間序列.

自相關(guān)系數(shù)計(jì)算公式為

(12)

偏自相關(guān)系數(shù)計(jì)算公式為

(13)

(2)若初始的統(tǒng)計(jì)數(shù)據(jù)為非平穩(wěn)的時(shí)間序列，且具有一定的變化趨勢(shì)，則對(duì)該組數(shù)據(jù)進(jìn)行一階差分處理，如果一階差分后的數(shù)據(jù)依然為非平穩(wěn)的時(shí)間序列，則進(jìn)行二階差分處理，一般差分次數(shù)控制在2次之內(nèi).

(3)根據(jù)時(shí)間序列模型判定規(guī)則[7]，建立合理的時(shí)間序列模型.若偏自相關(guān)系數(shù)呈截尾性，自相關(guān)系數(shù)呈拖尾性，則移動(dòng)平均分量階數(shù)為0，自回歸分量階數(shù)為p；若偏自相關(guān)系數(shù)呈拖尾性，自相關(guān)系數(shù)呈截尾性，則移動(dòng)平均分量階數(shù)為q，自回歸分量階數(shù)為0；若偏自相關(guān)系數(shù)呈截尾性，自相關(guān)系數(shù)呈截尾性，則移動(dòng)平均分量階數(shù)為p，自回歸分量階數(shù)為q.

(4)利用步驟(3)確定的ARIMA(p,d,q)分析初始數(shù)據(jù)，得出預(yù)測(cè)結(jié)果.

ARIMA模型能夠依據(jù)變量自身的變化規(guī)律，利用外推機(jī)制判斷出研究對(duì)象的變化情況,通過(guò)處理預(yù)測(cè)目標(biāo)的時(shí)間序列，獲得事物隨時(shí)間變化的演變特性及規(guī)律，從而預(yù)測(cè)事物未來(lái)的發(fā)展方向.ARIMA模型的基本適用條件是要求預(yù)測(cè)的數(shù)據(jù)滿足平穩(wěn)性特征，即個(gè)體值要圍繞均值小幅波動(dòng)，存在某些周期規(guī)律特性，且不存在明顯的變化趨勢(shì)(若出現(xiàn)明顯變化趨勢(shì)，則需要對(duì)預(yù)測(cè)數(shù)據(jù)進(jìn)行差分處理).高職高專院校的在校生數(shù)量和專任教師數(shù)量隨著年份的增加而不斷變化，雖然在每個(gè)季節(jié)的數(shù)量都不盡一致，但卻存在某種規(guī)律性，每一年周期內(nèi)在校生數(shù)量和專任教師數(shù)量的波動(dòng)趨勢(shì)大致相似，具有連續(xù)性和周期性，故適用ARIMA模型進(jìn)行預(yù)測(cè).

2 高職高專院?，F(xiàn)狀分析

表1為高等教育院校規(guī)模變化表，數(shù)據(jù)來(lái)源于教育部網(wǎng)站公布的2009～2017年全國(guó)教育事業(yè)發(fā)展統(tǒng)計(jì)公報(bào)[12～20].由表1可知：高職高專院校數(shù)量比本科院校要多；本科院校數(shù)量從2009年開始逐年上升，高職高專院校數(shù)量從2009～2017年也穩(wěn)步提升；2012年以前高職高專院校的增長(zhǎng)幅度較大，2012年以后其增長(zhǎng)幅度落后于本科院校，但自2017年開始，增長(zhǎng)幅度再次超過(guò)本科院校.

表2為高等教育專任教師數(shù)量變化表，數(shù)據(jù)來(lái)源于教育部網(wǎng)站公布的2010～2017年全國(guó)教育事業(yè)發(fā)展統(tǒng)計(jì)公報(bào)[13～20].由表2可以發(fā)現(xiàn):本科院校專任教師數(shù)量比高職高專院校專任教師數(shù)量多;本科院校專任教師數(shù)量從2010年開始逐年上升，高職高專院校數(shù)量從2010～2017年也穩(wěn)步提升;2015年以前本科院校專任教師增比要大于高職高專院校，自2015年開始，高職高專院校專任教師增比幅度大于本科院校;就總體趨勢(shì)來(lái)看，本科院校專任教師增比處于下滑趨勢(shì)，高職高專院校專任教師增比處于穩(wěn)步上升趨勢(shì).

表 1 高等教育院校規(guī)模變化表

表3為高等教育在校生數(shù)量變化表，數(shù)據(jù)來(lái)源于教育部網(wǎng)站公布的2010～2017年全國(guó)教育事業(yè)發(fā)展統(tǒng)計(jì)公報(bào)[13～20].由表3可以看出:本科院校在校生數(shù)量比高職高專院校在校生數(shù)量多;本科院校在校生數(shù)量從2010年開始逐年上升，高職高專院校在校生數(shù)量從2010～2017年也穩(wěn)步提升;在2013年以前，本科院校在校生增比要大于高職高專院校，自2014年以后，高職高專院校在校生增比幅度大于本科院校;就總體趨勢(shì)來(lái)看，本科院校在校生增比處于下滑趨勢(shì)，高職高專院校在校生增比呈上升趨勢(shì).表1～表3中出現(xiàn)的增比指標(biāo)按式(11)計(jì)算.

表 2 高等教育專任教師數(shù)量變化表

表 3 高等教育在校生數(shù)量變化表

3 ARIMA模型應(yīng)用

本文從教育部2009～2017年度全國(guó)教育統(tǒng)計(jì)資料中獲取了高職高專院校規(guī)模、專任教師數(shù)量、在校生數(shù)量的基本數(shù)據(jù)，這些數(shù)據(jù)真實(shí)且按時(shí)間順序排列，具有時(shí)間序列特征，能夠反映某類現(xiàn)象的統(tǒng)計(jì)指標(biāo)，且存在時(shí)間性和周期性，故借助應(yīng)用比較成熟的時(shí)間序列模型ARIMA進(jìn)行研究.

表4為高職高專院校專任教師數(shù)量數(shù)據(jù)的相關(guān)系數(shù)分布表，其中ACV表示自相關(guān)系數(shù)值，PACV表示偏自相關(guān)系數(shù)值.若自相關(guān)系數(shù)值在標(biāo)準(zhǔn)差寬度范圍內(nèi)，則表明數(shù)據(jù)落入隨機(jī)區(qū)間內(nèi)部，為穩(wěn)定的時(shí)間序列.初始分布的標(biāo)準(zhǔn)差寬度為0.39，一階差分的標(biāo)準(zhǔn)差寬度為0.54.由表4得出，初始的ACV在滯后6階時(shí)，自相關(guān)系數(shù)值在隨機(jī)區(qū)間外部，需要對(duì)初始數(shù)據(jù)進(jìn)行差分處理.一階差分后的自相關(guān)系數(shù)和偏自相關(guān)系數(shù)若全部落入隨機(jī)區(qū)間內(nèi)部，經(jīng)一階差分后的數(shù)據(jù)為平穩(wěn)時(shí)間序列，那么自回歸分量階數(shù)p為0，差分次數(shù)d為1，移動(dòng)平均分量階數(shù)q為0，那么選用的時(shí)間序列模型為ARIMA(0,1,0).

表 4 高職高專院校專任教師數(shù)量數(shù)據(jù)的相關(guān)系數(shù)分布表

表5為高職高專院校在校生數(shù)量數(shù)據(jù)的相關(guān)系數(shù)分布表.初始分布、一階差分和二階差分的標(biāo)準(zhǔn)差寬度分別為0.36、0.45和0.56.由表5可知:第1階、第5階和第6階自相關(guān)系數(shù)值在隨機(jī)區(qū)間外部，數(shù)據(jù)為不平穩(wěn)的時(shí)間序列；第1階偏自相關(guān)系數(shù)在隨機(jī)區(qū)間的外部，從第2階起全部落入隨機(jī)區(qū)間內(nèi)部，表現(xiàn)出截尾性，自回歸分量階數(shù)p為1.但是，已由自相關(guān)圖判斷出該數(shù)據(jù)為非平穩(wěn)時(shí)間序列，需要進(jìn)行平穩(wěn)化處理.經(jīng)一階差分后，滯后第1階和第4階的自相關(guān)系數(shù)在隨機(jī)區(qū)間外部，則移動(dòng)平均分量階數(shù)q為4，滯后第2階及以后的偏自相關(guān)系數(shù)均落入隨機(jī)區(qū)間內(nèi)部，則差分次數(shù)d為1，自回歸分量階數(shù)p為1，那么選用的時(shí)間序列模型為ARIMA(1，1，4).考慮到收集樣本數(shù)量較少，滯后階數(shù)過(guò)多可能影響預(yù)測(cè)結(jié)果的準(zhǔn)確度，本文將對(duì)一階差分后的數(shù)據(jù)再次進(jìn)行差分處理.二階差分后，所有階數(shù)的自相關(guān)系數(shù)和偏相關(guān)系數(shù)均落入隨機(jī)區(qū)間內(nèi)，呈拖尾性，經(jīng)二階差分后的數(shù)據(jù)為平穩(wěn)時(shí)間序列，自回歸分量階數(shù)p為0，差分次數(shù)d為2，移動(dòng)平均分量階數(shù)q為0，則選用的時(shí)間序列模型為ARIMA(0,2,0).

表 5 高職高專院校在校生數(shù)量數(shù)據(jù)的相關(guān)系數(shù)分布表

4 預(yù)測(cè)結(jié)果分析

通過(guò)ARIMA(0,1,0)模型對(duì)高職高專院校專任教師數(shù)量進(jìn)行時(shí)序分析，其預(yù)測(cè)結(jié)果如表6所示，下限值指的是置信區(qū)間下限，上限值指的是置信區(qū)間上限.由表6可知，預(yù)測(cè)模型結(jié)果均在顯著性水平0.05的置信區(qū)間內(nèi)，預(yù)測(cè)結(jié)果具有較高的可信度.預(yù)測(cè)結(jié)果顯示,2018年、2019年和2020年的專任教師數(shù)量分別為488 748人、499 613人和510 478人.從整體分布情況來(lái)看，高職高專院校專任教師數(shù)量在未來(lái)幾年依然會(huì)保持持續(xù)增長(zhǎng)的趨勢(shì).通過(guò)ARIMA(0,2,0)模型對(duì)高職高專院校在校生數(shù)量進(jìn)行時(shí)序分析，其預(yù)測(cè)結(jié)果如表6所示，預(yù)測(cè)模型結(jié)果均在顯著性水平0.05的置信區(qū)間內(nèi)，具有較高的可信度。預(yù)測(cè)結(jié)果顯示，2018年、2019年和2020年在校生數(shù)量分別為9 622 037人、9 997 531人和10 373 025人.從整體分布來(lái)看，高職高專院校在校生數(shù)量在未來(lái)幾年依然保持持續(xù)增長(zhǎng)態(tài)勢(shì).值得說(shuō)明的是，為考察專任教師和在校生數(shù)量預(yù)測(cè)值的可信度，進(jìn)行了相對(duì)誤差百分比的評(píng)價(jià)計(jì)算，從數(shù)值上來(lái)看均不到5%，表明預(yù)測(cè)結(jié)果滿足預(yù)測(cè)精度.

表 6 高職高專院校專任教師和在校生預(yù)測(cè)結(jié)果表

普通高校生師比定義為在校生數(shù)量與專任教師數(shù)量之比，教育部[2004]2號(hào)文件規(guī)定高等教育院校生師比的標(biāo)準(zhǔn)為18∶1，本文通過(guò)ARIMA模型預(yù)測(cè)結(jié)果計(jì)算了未來(lái)幾年高職高專院校的生師比，如表7所示，給出了高等教育院校生師比預(yù)測(cè)值對(duì)照結(jié)果.由表7可以看出：本科院校的生師比一直低于16∶1；高職高專院校生師比在2013年以前基本上接近于18∶1，但從2014年開始增速加快，在2017年生師比達(dá)到19.18∶1.預(yù)測(cè)模型結(jié)果顯示：2018年生師比為19.69∶1，2019年生師比為20.01∶1，2020年生師比位20.32∶1.生師比是判斷世界頂尖大學(xué)的基本標(biāo)準(zhǔn).泰晤士報(bào)上的大學(xué)排名，把生師比作為考量的重要指標(biāo).英國(guó)《獨(dú)立報(bào)》編輯湯姆·門德爾松說(shuō)過(guò)，“一所大學(xué)的學(xué)者越多，學(xué)生就會(huì)得到更多的關(guān)注，那么學(xué)生與教師的交流頻率更多，他們?cè)趯W(xué)習(xí)過(guò)程也會(huì)獲得更有價(jià)值的內(nèi)容”.創(chuàng)建一流的大學(xué)需要一流的師資力量，生師比數(shù)據(jù)就是一個(gè)需要考量的關(guān)鍵指標(biāo)，教師越多，所營(yíng)造的學(xué)術(shù)氛圍和教學(xué)環(huán)境越好.英國(guó)《衛(wèi)報(bào)》認(rèn)為：生師比應(yīng)考慮的是從事教學(xué)工作的教職員工，不包括從事少許教學(xué)工作的研究人員.生師比是反映大學(xué)教學(xué)質(zhì)量的重要指標(biāo)，良好的生師比是提升高職高專學(xué)生就業(yè)質(zhì)量的基礎(chǔ)保障.在未來(lái)的幾年里，高職高專院校還需要保障良好的師資團(tuán)隊(duì)，控制生師比，為高職高專院校的可持續(xù)發(fā)展提供基本的保障.

表 7 高等教育院校生師比預(yù)測(cè)值對(duì)照表

5 結(jié)論

高職高專院校的生師比指標(biāo)關(guān)乎學(xué)校教育的質(zhì)量，合格生師比指標(biāo)為18∶1，限制招生生師比指標(biāo)為22∶1.本文結(jié)合時(shí)間序列ARIMA模型對(duì)高職高專院校未來(lái)幾年的生師比進(jìn)行預(yù)測(cè)，得出結(jié)論如下：(1)就整體趨勢(shì)而言，高職高專院校數(shù)量增幅較快，且其專任教師和在校生數(shù)量增比要高過(guò)于本科院校；(2)初始的高職高專院校教育統(tǒng)計(jì)數(shù)據(jù)是真實(shí)的，具有時(shí)間序列特征，但是多數(shù)為不穩(wěn)定時(shí)間序列，需要進(jìn)行一階或二階差分處理；(3)高職高專院校2018年生師比為19.69∶1，2019年生師比為20.01∶1，2020年生師比為20.32∶1.結(jié)果表明：未來(lái)幾年內(nèi)高職高專院校生師比呈上升趨勢(shì)，將超過(guò)合格指標(biāo)，為師資隊(duì)伍建設(shè)及制定招生計(jì)劃提供理論參考.