楊曉蓉

（新疆師范大學 政法學院，新疆 烏魯木齊 830017）

2021年8月國務院印發(fā)《“十四五”就業(yè)促進規(guī)劃》，明確指出要“健全有利于更加充分更高質(zhì)量就業(yè)的促進機制”，提出要加強對就業(yè)政策實施效果的跟蹤調(diào)查評估，要“完善就業(yè)統(tǒng)計指標體系和調(diào)查統(tǒng)計方法，探索進行就業(yè)質(zhì)量、就業(yè)穩(wěn)定性等方面的分析?！保L期以來，西部地區(qū)人力資源分布不均、人才結(jié)構(gòu)不合理、基層人才儲備不足的短板一直未能得到有效改善，其就業(yè)促進政策從實施效果來看，在總體平穩(wěn)的基礎上依然存在一些問題和挑戰(zhàn)。因此，有必要通過精準的量化分析方法，對就業(yè)促進政策的實施效果進行追蹤評估，從而增強風險預警預判能力，提升形勢感知、分析研判和科學決策水平［1］。

新疆維吾爾自治區(qū)作為“一帶一路”核心區(qū)，基礎教育師資短缺和分配不均是制約其農(nóng)村地區(qū)教育水平提升的一個重要因素。為了有針對性地解決這一問題，政府陸續(xù)出臺了一系列提升基礎教育質(zhì)量的專項就業(yè)促進政策。2010年，新疆維吾爾自治區(qū)在充分借鑒國家公費師范生政策的基礎上，制定出臺了《自治區(qū)對疆內(nèi)師范生實行免費教育的管理辦法》，正式開始實施新疆公費師范生教育政策，其目的是解決南疆四地州中小學熟練掌握并教授國家通用語言文字師資力量短缺的現(xiàn)實困難。隨著公費師范生政策的深入實施，截止到2021年已有八屆三萬余名畢業(yè)生陸續(xù)走上基層教育工作崗位。越來越多的研究表明，新疆公費師范生教育政策成效與問題并存，特別是在政策推行初期，公費師范生的履約率并未達到預期，其就業(yè)穩(wěn)定性需要做進一步的分析和調(diào)研。

關(guān)于履約率的研究，很多學者采用的是傳統(tǒng)的Logistic 回歸統(tǒng)計分析方法。在公費師范生履約率分析中，關(guān)注的重點在于畢業(yè)后在多長時間內(nèi)有可能發(fā)生違約行為，即違約行為發(fā)生的概率隨時間的變化及其影響因素，因此違約行為發(fā)生所經(jīng)歷的時間是一項非常重要的考察變量，Logistic 回歸模型忽視了時間因素，所以并不適用于履約率研究。文章旨在通過生存分析模型，探究公費師范生履約率變化的關(guān)鍵時間節(jié)點以及影響公費師范生履約率的關(guān)鍵因素，為今后公費師范生政策的不斷調(diào)整和完善提供科學參考。

1 生存分析模型的應用

關(guān)于公費師范生履約意愿和履約率的研究，吳浩［2］討論了在博弈論模型的分析下，要提高公費師范生履約率，可以通過提高違約成本、提高農(nóng)村基層教師的地位，以此來完善政策，降低違約率。周昊昊［3］等人基于扎根理論，通過半結(jié)構(gòu)式訪談，運用Nvivo質(zhì)性分析軟件建立了職教師范生的從教意愿影響因素模型。更多的學者［4］采用的是傳統(tǒng)的Logistic 回歸統(tǒng)計分析方法，將是否履約的行為作為一個二分變量代入模型，探究影響履約意愿的主要因素［1-2］。然而，在樣本數(shù)據(jù)中，違約行為不是同時發(fā)生的，很多樣本在研究截止時還沒觀測到違約事件，不可避免地存在刪失數(shù)據(jù)的情況。因此，傳統(tǒng)的Logistic 回歸模型未能將時間因素納入其中，相對于生存分析模型來說適配性較低。

基于生存分析方法的生存函數(shù)可以反映出在每一個時間節(jié)點的履約率，揭示出Logistic 回歸模型等傳統(tǒng)統(tǒng)計學方法無法得到的動態(tài)特征，在本研究中的適配性較高。生存分析是一種動態(tài)的計量統(tǒng)計學分析方法，通常用來探究研究對象發(fā)生狀態(tài)轉(zhuǎn)變的現(xiàn)象和時間之間的關(guān)系以及兩者之間存在的規(guī)律性。它將時間維度納入到研究視野，定義因變量為個體持續(xù)處于“非轉(zhuǎn)變”的時間，該理論方法主要是通過構(gòu)建生存函數(shù)（概率密度函數(shù)），計算出各個時間點上的危險率，即研究對象保持狀態(tài)的生存時間和它發(fā)生轉(zhuǎn)變的可能性，分析各個不同因素對兩者的影響程度［5］。生存分析模型的優(yōu)點在于可以處理數(shù)據(jù)觀察不完全的情況，如在研究截止時樣本狀態(tài)尚未發(fā)生轉(zhuǎn)變。在生存分析統(tǒng)計中，通常分為三類：壽命表分析、Kaplan-Meier 分析、Cox風險回歸模型方法。

生存分析模型最初被廣泛應用于醫(yī)學、生物學、藥學等研究領(lǐng)域，分析具有臨床意義的某些醫(yī)學事件（死亡、治愈、發(fā)病等）發(fā)生的概率，并取得了大量豐碩成果。生存分析在工程科學領(lǐng)域被稱為“可靠性分析”或者“故障時間分析”，模型通常被用于研究裝置、器械、電路元件等設備發(fā)生故障的時間周期及其影響因素。近年來，經(jīng)濟學、社會學、人口統(tǒng)計學等領(lǐng)域也開始嘗試利用生存分析法進行研究。經(jīng)濟學家一開始用生存分析法研究工人罷工持續(xù)的時間、失業(yè)持續(xù)的時間、企業(yè)破產(chǎn)的時間等，后來逐漸推廣應用至研究金融風險控制和IPO［6］、農(nóng)村勞動力遷移的影響因素、旅游業(yè)游客停留天數(shù)的影響因素［7］等。在教育學領(lǐng)域，近年來生存分析法開始逐漸被關(guān)注，主要用于研究學生畢業(yè)和就業(yè)問題、人才資源的培養(yǎng)成長等，例如，就業(yè)質(zhì)量評價［8］、高層次人才成長研究［9］、網(wǎng)絡高等教育的畢業(yè)時間預測［10］、待業(yè)時長與受教育程度的關(guān)系分析。

對本研究而言，違約行為的發(fā)生即生存狀態(tài)變量，有兩個狀態(tài)：履約和違約，生存時間為公費師范生畢業(yè)以后的履約時長。同時，公費師范生作為行為個體，其個人特征、家庭背景、服務地類型等都會對違約事件的發(fā)生產(chǎn)生影響。文章應用Kaplan-Meier 分析法，估計不同單一變量對履約時長的生存分析函數(shù)，應用Cox風險回歸分析評價多個影響因子對履約率生存函數(shù)分布的影響程度。

2 數(shù)據(jù)來源及樣本描述統(tǒng)計分析

新疆公費師范生的培養(yǎng)分為兩個階段，第一階段為2010—2011年，公費師范生政策推行的試行階段，施行“定向”與“非定向”相結(jié)合的雙軌制，招生規(guī)模較大，以非定向為主。第二個階段為2012年至今，全部施行定向就業(yè)的公費師范生制度，省屬重點師范大學招生規(guī)模收縮，公費師范生名額大量分派至地方師范類院校。本研究選取了新疆公費師范生政策實施的初始階段，新疆某師范類院校四屆共811 名已畢業(yè)本科層次定向公費師范生作為調(diào)查對象，結(jié)合初次就業(yè)統(tǒng)計數(shù)據(jù)，并抽選了其中131 人進行電話調(diào)研，掌握了這部分畢業(yè)生的實際履約情況。

第一階段，2014—2015年畢業(yè)公費師范生共計2542人，大部分為非定向就業(yè)，僅有125人為定向就業(yè)公費師范生。生源質(zhì)量較好，漢語言和民語言的文理科高考平均錄取分數(shù)均高于一本錄取分數(shù)線。其中女生占比較高，達到72.3%；漢語言學生占比58.1%.但令人遺憾的是，這一階段培養(yǎng)的公費師范生雖然享受到了免除學費、住宿費的優(yōu)惠政策，但畢業(yè)后較高比例的學生沒有從事教師職業(yè)。根據(jù)初次就業(yè)率統(tǒng)計，僅有44.1%的畢業(yè)生選擇從事中小學教育工作，且絕大多數(shù)集中在經(jīng)濟條件較好的城鎮(zhèn)工作，深入農(nóng)村基層教育一線的畢業(yè)生寥寥無幾［11］。

第二階段，2016—2017 年共計畢業(yè)公費師范生760人，全部為定向就業(yè)公費師范生。其中女生占比較高，達到74.9%；母語為非漢語的少數(shù)民族學生占比73%.根據(jù)初次就業(yè)率統(tǒng)計，兩屆畢業(yè)生履約率均達到99%以上。2016 屆畢業(yè)公費師范生實際在崗履約率為86.73%，履約學生服務地在農(nóng)村的比例為19.14%；2017屆畢業(yè)公費師范生實際在崗履約率為98.0%，履約學生服務地在農(nóng)村的比例為18.66%.

如表1 所示，2014—2017 屆定向公費師范生實際履約率呈現(xiàn)逐級下降趨勢，且降幅明顯；公費師范畢業(yè)生總?cè)藬?shù)中，女生占比較高，達到74%，男生占比較低，僅有26%，履約率差異不大；母語為民語言的學生占比68.43%，且履約情況較為穩(wěn)定，違約率較低，僅有2.34%；理工和文史類公費師范畢業(yè)生是主體，藝術(shù)類學生占比較小，僅有3.82%，但違約率明顯偏高，達到25.81%；農(nóng)村戶籍和城市戶籍的學生比例基本相同，農(nóng)村戶籍的畢業(yè)生履約率明顯高于城市戶籍學生，城市戶籍學生的違約率是農(nóng)村學生的約兩倍；畢業(yè)生的定向服務地分配呈現(xiàn)紡錘形分布，主體為縣城和鄉(xiāng)鎮(zhèn)，占比為69.2%，兩頭分別為城市和農(nóng)村，能夠進入到農(nóng)村教育一線的學生人數(shù)僅占總?cè)藬?shù)的15.8%；服務地和生源地不一致的畢業(yè)生的違約率遠高于服務地和生源地一致的學生，服務地與生源地距離越遠，違約學生比例越高。

表1 新疆某高校2014-2017屆畢業(yè)公費師范生履約情況描述統(tǒng)計

3 Kaplan-Meier分析法測算及結(jié)果分析

Kaplan-Meier 分析法也叫單因素生存分析法，做法是選取某一個協(xié)變量，比較其分組因素對生存時間的影響。對于研究中的第n個時間點tn，履約生存概率可以計算為：

其中，f(tn-1)指的是在tn-1時間點的履約生存函數(shù)，dn指的是在時間點tn違約發(fā)生的事件數(shù)，rn指的是快要到時間點tn時，未發(fā)生違約行為的人數(shù)。

對于違約時間的分布也可以用密度函數(shù)的形式來表示：

文章在描述統(tǒng)計的基礎上，用Kaplan-Meier 分析法，從性別、年級、母語類別、專業(yè)、戶籍類型、定向服務地類型、是否生源地、生源地與服務地距離等8 個方面，分別構(gòu)建公費師范生履約時長的生存分析函數(shù)。對于本研究而言，生存風險率即履約率越低，也就意味著違約率越高。

以上8 個類別不同級別的生存分析函數(shù)顯著性檢驗如表2 所示，性別在公費師范生履約時長上的差異性不顯著，其余7個類別，年級、母語類別、專業(yè)、戶籍類型、定向服務地類型、是否生源地、生源地與服務地距離的Log Rank和Breslow檢驗，P值均小于0.05，表示7類指標在履約時長上均有顯著性差異。

表2 公費師范生不同類別生存分析的顯著性檢驗

履約時長的生存分析曲線可以直觀地反映出各觀測指標在不同分類上的履約概率變化，如圖1所示。

圖1 公費師范生不同性別、年級、母語類別、專業(yè)、戶籍類型、定向服務地類型、是否生源地、生源地與服務地距離的生存分析函數(shù)

根據(jù)不同觀測指標的履約生存分析函數(shù)圖可知，已畢業(yè)公費師范生的履約率會隨著時間的增長不斷下降，具體到每一類指標，分析結(jié)果如下：

（1）不同性別的履約率差異不大，但隨著履約時長的增加，前半段差異并不顯著，在履約時長的后半段，在第24個月前后出現(xiàn)相對明顯的下降，男生的履約率降幅明顯大于女生。

（2）不同年級的履約率存在明顯差異，從2017 級至2014 級依次遞減，2014 級履約率最低。各年級的履約率隨履約時長的變化均在第一個月，即剛畢業(yè)時段出現(xiàn)明顯下降，之后波動不大。其中2014、2015級學生在履約約6個月、24個月、30個月左右出現(xiàn)小幅下降。

（3）母語類別不同的學生履約率差異明顯，母語為民語言的學生履約率一直較高，且穩(wěn)定性較好，不隨時間的變化波動。母語為漢語言的學生履約率隨著時間的變化不斷下降，在剛畢業(yè)和第30 個月左右，履約率出現(xiàn)較為明顯的降幅。

（4）在專業(yè)分類方面，理工類和文史類專業(yè)學生的履約率相對較高，且就業(yè)穩(wěn)定性較好，藝術(shù)類學生在剛畢業(yè)階段就出現(xiàn)履約率明顯降幅。由于本研究只針對定向就業(yè)公費師范生，2012年以前僅有理工類定向就業(yè)師范生，因此文史和藝術(shù)類公費師范生存在數(shù)據(jù)刪失。

（5）不同戶籍類型的履約率差異明顯，農(nóng)村戶籍的學生履約率高且穩(wěn)定性好，城市戶籍的學生履約率較低，且隨著時間的變化不斷降低，在第1個月和第30個月有明顯的降幅。

（6）不同定向服務地類型對履約率的差異明顯，服務地類型為城市的學生履約穩(wěn)定性較高，類型為鄉(xiāng)鎮(zhèn)和農(nóng)村的履約率相對略低，服務地類型為縣城的履約率最低。這一情況與實際認知不一致，通常認為服務地為農(nóng)村的履約率應低于其他地區(qū)，結(jié)合訪談發(fā)現(xiàn)，由于很多違約行為發(fā)生在剛畢業(yè)，畢業(yè)即違約，學生僅知曉負責服務地分配的縣教育局，不清楚自己的具體服務地，因此違約樣本的服務地分類集中為縣城，但并不代表其真實服務地，存在信息認定不準確的情況。

（7）服務地是否為生源地對履約率影響明顯，服務地為生源地的學生履約率較高，且穩(wěn)定性較好，服務地為非生源地的學生履約率相對較低，且隨著履約時長的增加不斷降低，在第1個月和第30個月降幅明顯。

（8）為了更好地考察服務地對履約率的影響，本研究將服務地與生源地的距離劃分為四個層級，從生存分析函數(shù)圖可以明顯看出，兩者之間的距離越近，履約率越高，且隨著時間遷移的穩(wěn)定性越好；兩者之間的距離越遠，履約率越低，隨著履約時間的增長，穩(wěn)定性越差。

4 Cox比例風險回歸法的測算及結(jié)果分析

Cox 比例風險回歸法是一種多因素生存分析方法，可用于多種因素的生存率比較及其影響因素分析。在本研究中，Cox比例風險回歸模型的基本形式為：在某個履約時間下，樣本出現(xiàn)t個違約行為的風險。第i個影響因素Xi使履約率風險函數(shù)從h0(t)增加exp(βiXi)，成為h0(t) × exp(βiXi)，多個因素同時影響的模型為：

其中，h(t,X)表示個體在協(xié)變量X(X=X1,X2,……Xi)的作用下，在時刻t的風險函數(shù)。X1,X2,……Xi表示影響因素，β1、β2……βi表示回歸系數(shù)。

基于Kaplan-Meier 分析法，本研究繼續(xù)將性別、年級、母語類別、專業(yè)、戶籍類型、定向服務地類型、是否生源地、生源地與服務地距離等8 類指標作為考察變量，都列為協(xié)變量。如表3 所示，其中X1表示性別（男、女），X2表示年級（2014—2017 屆），X3表示母語類別（漢語言、民語言），X4表示專業(yè)分類（文史、理工、藝術(shù)），X5表示戶籍（農(nóng)村戶口，城市戶口），X6表示定向服務類型（城市、縣城、鄉(xiāng)鎮(zhèn)、農(nóng)村），X7表示服務地與生源地是否一致（是、否），X8表示服務地與生源地距離（按照交通時間，將距離按照由近及遠劃分為4個分段）。

表3 新疆已畢業(yè)公費師范生履約率影響因素生存分析模型因變量和自變量賦值

利用SPSS 23.0 軟件進行生存分析Cox 比例風險回歸，時間為履約時長，狀態(tài)為履約0 或違約1，將選定的8個指標作為協(xié)變量，帶入模型，使用強迫進入變量法（enter），分析隨著履約時長的變化，各指標對公費師范生履約生存風險率的影響。模型分析結(jié)果如下：

由圖2可知，在多重因素的影響下，隨著履約時長的增加，新疆公費師范生的履約概率不斷下降，違約行為的發(fā)生概率增加，但相對來說變化較小，總體降幅不大。截止到第30個月，總體履約率在95%以上。說明在畢業(yè)生履約率方面，公費師范生政策表現(xiàn)出較好的效果，履約穩(wěn)定性較強。

圖2 新疆高校已畢業(yè)公費師范生履約生存分析函數(shù)

同時，履約率的降低并非想象中的那樣呈平滑曲線逐漸過渡，而是表現(xiàn)出明顯的階梯狀下降。履約率下降第一個顯著陡坡出現(xiàn)在履約時長的第1個月，表明有較大比例的違約行為集中發(fā)生在剛畢業(yè)時段；之后以半年為單位，逐漸下降，這一情況應該是受到學校教學周期的影響；第二個略顯著的陡坡出現(xiàn)在履約時長的第30個月，即履約兩年半左右，可觀測到相對較多的違約行為。兩個節(jié)點中，畢業(yè)時間段是違約行為高發(fā)的重要節(jié)點，一旦履約上崗之后，違約發(fā)生的概率就大幅下降。

從表4 可知，各協(xié)變量對新疆已畢業(yè)公費師范生履約概率預測的回歸模型的總體顯著性P=0.000，說明模型適配度非常理想。

表4 新疆高校已畢業(yè)公費師范生履約率的Cox比例風險回歸

從各協(xié)變量的顯著性指標來看，除了性別、生源地與服務地距離兩項P 值未達到0.05 顯著水平，其余變量均達到顯著水平。

在年級方面，以2017 級為參照組，2014、2015、2016 級的t值分別為16.264、13.051、17.118，P 值均小于0.001，說明相對于參照組，這三個年級的履約率存在顯著差異。根據(jù)勝算比Exp（B）值可知，2014 級公費師范生的違約概率是2017 級學生的10.558倍，2015 級公費師范生的違約概率是2017 級學生的7.543倍，2016級公費師范生的違約概率是2017級學生的5.247倍。

在母語類別方面，B值為負，t值為11.544，在0.001水平顯著，勝算比Exp（B）為0.258，母語為民語言學生的違約概率是母語為漢語言學生的0.258倍，也可以表示為母語為漢語言學生的違約概率是母語為民語言學生的3.876倍。

在專業(yè)方面，以藝術(shù)體育專業(yè)作為參照組，文史專業(yè)、理工專業(yè)的履約率與之對比均存在顯著差異，t值分別為5.689、6.427，均在0.05 水平上顯著，勝算比Exp（B）為0.359、0.292，即藝術(shù)體育類專業(yè)學生的違約概率是文史類專業(yè)學生的2.786倍，是理工類學生的3.425倍。

在戶籍類型方面，農(nóng)村戶籍與城市戶籍學生的履約率存在顯著差異，t值為4.544，在0.05水平上顯著，勝算比Exp（B）為1.75，即城市戶籍學生的違約概率是農(nóng)村戶籍學生的1.75倍。

在定向服務地類型方面，以定向服務地為城市作為參照組，定向服務地為農(nóng)村、鄉(xiāng)鎮(zhèn)、縣城的學生履約率與之相比均存在顯著差異，顯著性均達到0.05水平，勝算比Exp（B）值表明，定向服務地為農(nóng)村的學生違約概率是城市的4.262倍，定向服務地為鄉(xiāng)鎮(zhèn)的學生違約概率是城市的2.832倍，定向服務地為縣城的學生違約概率是城市的2.957倍。

5 結(jié)論

通過對新疆高校已畢業(yè)公費師范生的履約情況進行實證調(diào)研，并運用生存分析模型，對實際履約情況的影響因素進行了分析，可以得出如下結(jié)論：

（1）新疆公費師范生的履約率雖然隨著時長的增加而不斷下降，但相對來說降幅不大，總體履約率在95%以上。說明在畢業(yè)生履約率方面，公費師范生政策表現(xiàn)出了較好的實施效果，就業(yè)穩(wěn)定性較強。

（2）通過Kaplan-Meier 分析結(jié)果來看，隨著履約時長的增加，性別在履約率上的差異不顯著，但年級、母語類別、專業(yè)分類、定向服務地類型、戶籍類型、是否生源地、生源地與服務地之間的距離等7項觀測指標，履約率差異顯著。民語言學生履約率優(yōu)于漢語言，理工和文史類專業(yè)履約率優(yōu)于藝術(shù)體育類專業(yè)，農(nóng)村戶籍的畢業(yè)生履約率優(yōu)于城市戶籍，服務地在城市的畢業(yè)生履約率優(yōu)于鄉(xiāng)、鎮(zhèn)和農(nóng)村，生源地定向就業(yè)履約率優(yōu)于非生源地，生源地與服務地距離越遠履約率越低。

（3）通常認為畢業(yè)生違約行為大部分發(fā)生在上崗后，且隨著時間的增長，違約率逐年增高。通過Cox 比例風險回歸法分析發(fā)現(xiàn)，違約率雖然隨時間的增長逐漸降低，但相對來說變化幅度較小。公費師范生違約行為的發(fā)生集中在剛畢業(yè)階段，一旦履約上崗之后，出現(xiàn)違約行為的比例大幅下降。因此提高公費師范生履約率的關(guān)鍵階段在大學四年級，要注重強化學生的誠信意識，有針對性地對易違約群體出臺激勵機制，提高學生履約意愿，實現(xiàn)政策效果最大化。

綜上，影響新疆公費師范生就業(yè)穩(wěn)定性的因素較為復雜，需要在今后的研究中深入探索，通過更加精準的量化分析方法，對政策實施效果進行調(diào)查研究，不斷增強就業(yè)反饋機制的科學性，為新時期高質(zhì)量就業(yè)促進政策的不斷完善提出合理化建議。