張抗私，劉翠花

(東北財經(jīng)大學 經(jīng)濟學院，遼寧　大連　116025)

一、問題的提出

隨著中國市場經(jīng)濟體制改革的深化，工會化范圍不斷擴大，據(jù)2016年《中國統(tǒng)計年鑒》顯示，截至2015年底全國已建工會組織的基層機構有近280.60萬家，約為2003年的3倍，工會成員已達到2.95億人。作為國民經(jīng)濟發(fā)展的重要力量，國有企業(yè)工會組織覆蓋率更高。從功能來看，國有企業(yè)工會具有參與薪資和福利待遇議價以及改善工作條件等作用，是員工利益的維護者和代言人。工會的作用一直備受國內(nèi)外學者關注，概括來講主要體現(xiàn)在工資溢價、額外福利和改善勞動關系等方面。有關工會的工資溢價作用，最早源于Lewis[1]的經(jīng)典研究，他指出工會將工人的工資水平提高了10%—15%。Rosen[2]研究發(fā)現(xiàn)，與無工會企業(yè)相比較，工會企業(yè)的工資水平要高出16%—25%。Lanot 和Walker[3]指出當使用“半?yún)?shù)”方法排除自選擇帶來的內(nèi)生性問題時，工會仍對工資有顯著正向作用。魏下海等[4]的研究也表明工會能夠?qū)е鹿べY顯著提升，即中國企業(yè)確實存在工會的工資溢價效應。李明和徐律煒[5]指出不同技能等級員工獲得工會的工資溢價程度不同，中等技能員工加入工會獲得的工資溢價水平要高于均值，國有企業(yè)員工獲得的工資溢價水平高于非國有企業(yè)。

工會具有調(diào)節(jié)額外福利如減少工作時間和參加醫(yī)療保險等作用，F(xiàn)reeman和 Medoff[6]提出在工會企業(yè)中，年長的、流動性較低的工人偏好往往會被工會優(yōu)先考慮，并通過集體談判方式進一步提高企業(yè)在額外福利方面的支付。Buchmueller等[7]指出員工加入工會組織能夠大幅度提高其獲得養(yǎng)老和醫(yī)療等保險的概率。趙忠義[8]指出工會能夠顯著提高工人小時工資和縮短每月平均工作時間，并提高企業(yè)養(yǎng)老保險覆蓋率和人力資本投資回報率。Freeman[9]認為工會還具有改善勞動關系的作用，主要通過集體發(fā)聲降低離職率以促使雇用關系更加穩(wěn)定，以及在勞資談判時為會員爭取更為良好的薪資報酬兩種機制來影響雇員和雇主行為。Kahn 和Morimune[10]發(fā)現(xiàn)當工會力量足夠強大并且能夠與企業(yè)進行工資和雇用規(guī)模談判時，在短期內(nèi)裁員情況會受到一定的限制。魏下海等[11]的研究進一步證實工會可以顯著改善企業(yè)雇用期限結構，有效提高企業(yè)長期雇用占比。

工會會員是員工的一個客觀身份，加入工會組織會使員工的工資、福利和勞資關系等得到改善，而以上方面也密切關乎到員工幸福感。Clark和Oswald[12]認為相比失業(yè)者，就業(yè)者幸福感明顯更高。Diener[13]認為幸福感是人們對于自身生活滿意度的一種感受和評價，包括工作和家庭等眾多因素的主觀感受。Frey和Stutzer[14]認為在某一特定時間點上，幸福感與收入水平之間存在正向關系，即收入水平越高，幸福感越強。國內(nèi)學者關于幸福感的研究也逐漸豐富，李儒林等[15]認為幸福感的影響因素主要包括氣質(zhì)、健康和性別等內(nèi)部因素，以及文化、工作、收入和婚姻等外部因素。張進和馬月婷[16]研究發(fā)現(xiàn)工作滿意度與幸福感的認知因素(生活滿意度)之間存在顯著正相關關系，并且持有積極情感能夠較好地調(diào)節(jié)工作滿意度與幸福感之間的關系。魯元平和王韜[17]指出收入相對較高、家庭信任感強、工作滿意度高、人際關系和諧、積極樂觀的人幸福感更強，而失業(yè)、收入不平等和通貨膨脹等則會造成幸福感弱化。楊勝利和謝超[18]考察了就業(yè)對員工幸福感的影響效應，指出周工作時間、工作自由程度、收入水平和參加社會保險對員工幸福感有顯著影響。吳愈曉等[19]發(fā)現(xiàn)體制內(nèi)就業(yè)女性的平均工作時間遠低于體制外就業(yè)女性，不就業(yè)女性和體制內(nèi)就業(yè)女性的幸福感強于體制外就業(yè)女性。

綜上所述，工會對國有企業(yè)員工的工資水平、福利待遇和勞資關系等會產(chǎn)生積極作用，進而對國有企業(yè)員工幸福感產(chǎn)生影響，但影響程度如何以及為什么會產(chǎn)生這種影響是本文重點研究的內(nèi)容。

二、數(shù)據(jù)來源、模型設定與變量說明

(一)數(shù)據(jù)來源

本文使用2013年中國綜合社會調(diào)查(CGSS)中的城鄉(xiāng)居民生活綜合研究項目數(shù)據(jù)，該調(diào)查是由中國人民大學社會科學部和香港科技大學社會科學部發(fā)起的一項全國性、綜合性和連續(xù)性學術性調(diào)查，涉及到全國范圍內(nèi)28個省份125個縣(區(qū))，并最終獲得約12 000戶家庭問卷調(diào)查的隨機樣本。CGSS是國內(nèi)大型調(diào)查中質(zhì)量較高、最為開放、影響力最大的數(shù)據(jù)之一，同時2013年數(shù)據(jù)是目前網(wǎng)上平臺公布的最新數(shù)據(jù)，保證了數(shù)據(jù)的時效性，具有較理想的代表性。本文研究對象為目前在國有企業(yè)工作的18—55歲的員工，剔除掉數(shù)據(jù)缺失的樣本后，共計得到2 951個有效微觀樣本。

(二)模型設定與變量說明

本文被解釋變量國有企業(yè)員工幸福感(Hap)為有序離散變量，賦值為1—5的整數(shù)，運用有序Probit模型來分析參加工會對國有企業(yè)員工幸福感的影響。該模型假定存在一個能夠代表被解釋變量(Hap)但又不能直接觀測的潛在變量(Hap*)，本文假定潛在變量(Hap*)由下式?jīng)Q定：

表1　變量賦值及描述性統(tǒng)計

三、結果分析

(一)基準模型

表2中模型(1)是工會對國有企業(yè)員工幸福感影響的全樣本有序Probit回歸結果，模型(2)—模型(6)分別是幸福感依次取1—5時的回歸結果。

表2　有序Probit回歸結果

注：括號內(nèi)為標準誤差，***、**和*分別表示在1%、5%和10%的水平上顯著，下同。觀測值為2 951。

從表2可以看出，全樣本回歸結果顯示，工會的回歸系數(shù)在1%水平上顯著為正，這表明參加工會顯著提升了國有企業(yè)員工的幸福感。具體而言，參加工會使得國有企業(yè)員工感到比較幸福和非常幸福的概率上升5.5%和6.1%，使得國有企業(yè)員工感到非常不幸福、比較不幸福、一般幸福的概率下降0.6%、1.0%和1.7%。

從性別來看，女性國有企業(yè)員工幸福感高于男性，這可能是由于男性國有企業(yè)員工要承擔更多的社會責任和支撐家庭經(jīng)濟來源等壓力有關。由年齡和年齡平方/100的回歸系數(shù)可以看出，年齡與國有企業(yè)員工幸福感之間存在二元非線性U型關系，即在一定年齡之前，國有企業(yè)員工幸福感會隨著年齡的增加而呈現(xiàn)下降的趨勢，當達到年齡拐點之后，幸福感會隨著年齡的增加而呈現(xiàn)上升趨勢。受教育年限、政治面貌、戶籍狀況和婚姻狀況對國有企業(yè)員工幸福感的影響顯著為正。

收入會顯著提高國有企業(yè)員工幸福感，從全樣本回歸結果來看，國有企業(yè)員工幸福感的回歸系數(shù)在1%水平上顯著為正。從邊際效果來看，收入每提升1%，使得國有企業(yè)員工感覺比較幸福和非常幸福的概率上升1.9%和4.5%，這表明收入水平提高是影響國有企業(yè)員工幸福感的重要因素。擁有住房對國有企業(yè)員工幸福感有顯著正向影響，擁有住房使得國有企業(yè)員工感到非常不幸福、比較不幸福和一般幸福的概率分別下降0.6%、4.4%和10.3%。近年來房價越來越高，特別是在一、二線城市買房更為困難，住房問題的解決會使得國有企業(yè)員工主觀滿意度和幸福感均顯著提升。

健康狀況對國有企業(yè)員工幸福感有顯著正向影響，即國有企業(yè)員工身體越健康，會有更多的體力和精力去工作、生活和享受，幸福感才會越高。參加醫(yī)療保險對國有企業(yè)員工幸福感有顯著正向影響，主要原因是醫(yī)療保險可以有效促進參保居民健康水平的提高，而健康水平的提高又與國有企業(yè)員工幸福感之間存在正向關系。此外，通貨膨脹率對國有企業(yè)員工幸福感的影響在10%水平上顯著為負，這表明由于較高的通貨膨脹率會造成物價水平上漲，進而會造成實際貨幣工資水平下降，進一步使得國有企業(yè)員工幸福感降低。

(二)內(nèi)生性檢驗

參加工會將顯著提升國有企業(yè)員工幸福感，但兩者之間可能存在內(nèi)生性偏誤問題。產(chǎn)生內(nèi)生性可能是由于遺漏變量、測量誤差或聯(lián)立性偏誤等原因造成的，而內(nèi)生性則會導致估計系數(shù)出現(xiàn)有偏和非一致性。解決內(nèi)生性問題最有效的方法是工具變量法，即找到一個影響國有企業(yè)組建工會但不會影響員工幸福感的變量。本文借鑒魏下海等[11]的思路，采用企業(yè)長期雇用占比(Emp)作為員工是否參加工會的一個工具變量。通常企業(yè)長期雇用占比越高，其雇用關系越穩(wěn)定，這些勞動力隊伍越穩(wěn)定的企業(yè)有更強烈的動機組建工會，從而員工加入工會的可能性越高，即企業(yè)長期雇用占比與是否參加工會之間存在高度正相關性。同時，長期雇用占比是國有企業(yè)用工結構體系中客觀存在的變量，與國有企業(yè)員工幸福感之間不存在明顯的相關性，因而采用企業(yè)長期雇用占比作為本文的工具變量較為合理。計算公式為：企業(yè)長期雇用占比=簽訂長期雇用合同企業(yè)員工人數(shù)/企業(yè)全體員工人數(shù)×100%。

需要注意的是，當是否參加工會與國有企業(yè)員工幸福感都為離散數(shù)據(jù)時，若仍采用兩階段最小二乘法(2SLS)等工具變量法會失效。為了解決這一問題，Sajaia[21]提出了雙變量有序Probit模型，與兩步估計法和兩階段最小二乘法比較而言，利用完全信息極大似然方法估計得到的結果更為有效，同時這種方法得到了學界的認可和使用，因而本文以國有企業(yè)長期雇用占比作為是否參加工會的工具變量進行兩階段回歸，雙變量有序Probit模型的兩階段回歸結果如表3所示。

表3　雙變量有序Probit模型的兩階段回歸結果

從表3可以看出，內(nèi)生性檢驗輔助參數(shù)Atanhrho的估計系數(shù)為0.624，且在5%水平上顯著，因而拒絕參加工會為外生變量的原假設，即核心解釋變量是否參加工會與被解釋變量國有企業(yè)員工幸福感之間存在內(nèi)生性問題，從而采用雙變量有序Probit回歸進行調(diào)整是有必要的。第一階段回歸結果顯示，企業(yè)長期雇用占比對參加工會的影響在1%水平上顯著為正(回歸系數(shù)為0.160)，即企業(yè)雇用結構越穩(wěn)定，建立工會的動機越強，從而國有企業(yè)員工參加工會的概率也越高，這與前面的分析相一致。第二階段回歸結果顯示，參加工會對國有企業(yè)員工幸福感的影響在1%水平上顯著為正(回歸系數(shù)為0.124)，高于表2中沒有控制內(nèi)生性時的回歸結果，表明當排除了內(nèi)生性問題后，參加工會對于提升國有企業(yè)員工幸福感的作用更大。

第一階段回歸結果顯示，勞動者收入的回歸系數(shù)為0.121，且在1%水平上顯著，表明國有企業(yè)員工加入工會可以有效提升其工資福利水平，這也證實了前面文獻回顧中提到的工會的工資溢價效應，同時收入水平的提升又會提升國有企業(yè)員工的幸福感。員工的健康狀況和是否參加醫(yī)療保險的回歸系數(shù)分別為0.102和0.136，且分別在1%和5%水平上顯著，表明加入工會的國有企業(yè)員工會積極參加醫(yī)療保險，并定期進行健康體檢，這些措施將有助于提升員工的健康狀況，同時結合第二階段回歸結果，積極參加醫(yī)療保險和良好的健康狀況又會顯著提升員工的幸福感。此外，兩階段回歸結果顯示，員工的政治面貌和婚姻狀況與工會之間也呈現(xiàn)顯著正相關關系，且政治面貌和婚姻狀況又與員工幸福感呈現(xiàn)顯著正相關關系，表明國有企業(yè)員工加入工會后，由于積極參與企業(yè)年會和文體娛樂等組織活動，政治思想覺悟和身心健康得到顯著提升，家庭生活也會變得更加充實幸福，因而參加工會的國有企業(yè)員工幸福感會顯著提高。

(三)分樣本估計

前文是對全樣本進行的回歸分析，得出的是工會對國有企業(yè)員工幸福感的一個平均影響效應，并沒有對不同群體進行區(qū)分。為進一步考察參加工會對國有企業(yè)員工群體是否存在群體異質(zhì)性和區(qū)域異質(zhì)性，本部分從性別和地區(qū)兩個角度進行劃分，即將全體樣本劃分為男性和女性以及東部地區(qū)、中部地區(qū)和西部地區(qū)等子樣本，仍然以企業(yè)長期雇用占比作為工會的工具變量，并用雙變量有序Probit模型進行回歸，具體回歸結果如表4所示。

表4　分樣本回歸結果

注：限于篇幅，沒有列出第一階段的回歸結果，留存?zhèn)渌鳌?/p>

從表4可以看出，參加工會無論是對男性國有企業(yè)員工還是女性國有企業(yè)員工的幸福感而言，都會產(chǎn)生正向影響，且均在1%水平上顯著，但參加工會對女性國有企業(yè)員工幸福感的影響明顯高于對男性國有企業(yè)員工的影響，表明女性國有企業(yè)員工更看重企業(yè)是否設立工會組織以及自身是否能夠成為工會成員。這可能是由于女性在勞動力市場中會面臨一定的性別歧視和就業(yè)歧視，而且女性國有企業(yè)員工自身有更強的感性認知和歸屬感，如果她們加入企業(yè)工會則能夠更好地代表、參與和組織單位的民主管理及監(jiān)督，更好地維護女性國有企業(yè)員工在企業(yè)中的合法勞動權益，并且在工會中不斷加強學習和自我教育，以積極推動社會生產(chǎn)力水平的提高。

從地區(qū)分樣本來看，參加工會對不同地區(qū)國有企業(yè)員工幸福感影響不同，從東部沿海到西部內(nèi)陸的區(qū)域異質(zhì)性較為明顯，回歸結果顯示正向影響逐漸增強，這表明與東部地區(qū)國有企業(yè)員工相比，果參加工會能夠給西部地區(qū)國有企業(yè)員工帶來更高的幸福感。一個可能的原因是由于經(jīng)濟發(fā)展水平上的差異，西部地區(qū)設立工會組織的國有企業(yè)相對較少。因此，如果西部地區(qū)國有企業(yè)員工能夠成為工會組織中的成員，可能無論是在主觀心理的幸福感程度上如提升員工滿意度和改善福利水平，還是在實踐中通過簽訂集體工資協(xié)議和勞動合同等維護勞動者權益方面，都會比在東部地區(qū)參加工會的國有企業(yè)員工幸福感的影響更高。

關于其他控制變量，受教育年限對國有企業(yè)員工幸福感的影響呈現(xiàn)出較明顯的性別異質(zhì)性和區(qū)域異質(zhì)性。具體而言，受教育年限對提升女性國有企業(yè)員工幸福感的影響高于對男性國有企業(yè)員工幸福感的影響；與東部地區(qū)相比，受教育年限對提升中部地區(qū)和西部地區(qū)國有企業(yè)員工幸福感的影響更加明顯。對于政治面貌，身為中共黨員對男性國有企業(yè)員工幸福感的影響高于對女性國有企業(yè)員工的影響；同時，從東部地區(qū)到中部地區(qū)再到西部地區(qū)的地域變化中，中共黨員身份對提升國有企業(yè)員工幸福感的正向影響程度逐漸增強。性別分樣本回歸結果顯示，與男性國有企業(yè)員工相比，城鎮(zhèn)戶籍、收入越高、擁有住房和參加醫(yī)療保險對提升女性國有企業(yè)員工幸福感更為關鍵；婚姻狀況和健康狀況對男性國有企業(yè)員工幸福感的正向影響更大，表明已婚和身體越健康的男性國有企業(yè)員工幸福感要明顯高于同等條件下的女性國有企業(yè)員工。通貨膨脹率對國有企業(yè)員工幸福感均會產(chǎn)生顯著負向影響，而且相比于男性以及東部地區(qū)和中部地區(qū)的國有企業(yè)員工而言，通貨膨脹率對于女性和西部地區(qū)國有企業(yè)員工幸福感的負向影響更大。

四、研究結論與政策建議

本文基于2013年CGSS數(shù)據(jù)定量分析了參加工會對國有企業(yè)員工幸福感的影響，以企業(yè)長期雇用占比作為工會的工具變量進行內(nèi)生性處理，并按照性別和地區(qū)分樣本進行雙變量有序Probit回歸，以進一步檢驗工會對國有企業(yè)員工幸福感是否存在群體異質(zhì)性。研究得出以下結論：第一，工會對國有企業(yè)員工幸福感有顯著正向影響。全樣本回歸結果顯示，工會對國有企業(yè)員工幸福感的影響在1%水平上顯著為正，表明參加工會顯著提升了國有企業(yè)員工的幸福感。以企業(yè)長期雇用占比作為工會的工具變量有效解決了內(nèi)生性問題，雙變量有序Probit回歸結果顯示，企業(yè)雇用結構越穩(wěn)定，國有企業(yè)員工參加工會的概率越高，從而會更加顯著地提升國有企業(yè)員工幸福感。第二，工會能顯著提升國有企業(yè)員工的薪資福利水平。內(nèi)生性檢驗顯示，工會與勞動者收入之間呈現(xiàn)顯著正相關關系，表明國有企業(yè)員工加入工會可以有效提升其工資水平，工會的工資溢價效應明顯。員工的健康狀況、是否參加醫(yī)療保險與工會之間也呈現(xiàn)出正相關關系，且分別在1%和5%水平上顯著，由于加入工會的國有企業(yè)員工會積極參加醫(yī)療保險，并定期進行健康體檢，這些措施將有助于提升員工的身心健康，而積極參加醫(yī)療保險和良好的健康狀況又會顯著提升員工的幸福感。第三，工會對提升女性國有企業(yè)員工幸福感的影響更大。性別分樣本估計結果顯示，與男性國有企業(yè)員工相比，參加工會和受教育年限對于提升女性國有企業(yè)員工的幸福感影響更大；地區(qū)分樣本估計結果顯示，從東部地區(qū)到中部地區(qū)再到西部地區(qū)，參加工會和受教育年限對國有企業(yè)員工幸福感的正向影響逐漸增強，這表明與東部地區(qū)國有企業(yè)員工相比，中部地區(qū)和西部地區(qū)國有企業(yè)員工如果參加工會并不斷增加受教育年限則能夠給個人帶來更高的幸福感。

為進一步提高國有企業(yè)員工幸福感，加強工會對員工幸福感的影響，政府和國有企業(yè)工會應該努力做到以下方面：第一，國有企業(yè)工會必須努力為員工提供高質(zhì)量服務，如積極提升和改善工會會員的薪資福利待遇、工作條件和醫(yī)療保險覆蓋率等保障制度，定期舉辦企業(yè)年會和文體娛樂等活動，并對國有企業(yè)履行勞動合同和集體合同的情況進行監(jiān)督，切實維護勞動者的合法權益和身心健康。第二，國有企業(yè)工會要關照和幫助女性員工加入工會組織，從她們最關心、最直接、最需要解決的現(xiàn)實問題入手，如妥善解決女性員工生育費用補償和帶薪產(chǎn)假等問題。與此同時，國有企業(yè)工會也要更多關注男性員工切身合法權益，積極同工會建立平等協(xié)商制度，讓工會不僅成為女性員工也要成為男性員工抒發(fā)心聲的平臺，以充分提升國有企業(yè)整體員工的幸福感。第三，政府要加強對國有企業(yè)的監(jiān)管，以防止國有企業(yè)為追逐最大化利潤而減少對工會的經(jīng)濟投入，要盡快理順和夯實工會活動的財力來源，加強工會的經(jīng)濟基礎，從而為確保工會各項活動的順利開展提供可靠的經(jīng)濟支持。此外，政府應努力消除勞動力市場中戶籍門檻限制，增強員工就業(yè)穩(wěn)定性，做到精準幫扶困難員工群體，構建和諧勞動關系，讓廣大員工過上更加幸福美好的生活。

參考文獻：

[1]Lewis, H.G. Unionism and Relative Wages in the United States [M]. Chicago: Chicago University Press，1963.

[2]Rosen, S. On the Interindustry Wage and Hours Structure[J]. The Journal of Political Economy, 1969，77(2)：249-273.

[3]Lanot, G., Walker, I. The Union/ Non Union Wage Differential: An Application of Semi-Parametric Methods[J]. Journal of Econometrics,1998, 84(6):327-349.

[4]魏下海，董志強，黃玖立. 工會是否改善勞動收入份額? [J]. 經(jīng)濟研究，2013，(8): 16-28.

[5]李明，徐建煒. 誰從中國工會會員身份中獲益? [J]. 經(jīng)濟研究，2014，(5): 49-62.

[6]Freeman, R .B. , Medoff, J. L. What Do Unions Do？[M]. New York: Basic Books Press, 1984.

[7]Buchmueller, T. C., Dinardo, J., Valletta, R.G. Union Effects on Health Insurance Provision and Coverage in the United States[J]. Industrial and Labor Relations Review, 2002, 55(4): 610-627.

[8]趙忠義. 工會壟斷與經(jīng)濟效率[J]. 經(jīng)濟研究，2000，(12):43-47.

[9]Freeman, R. B. The Effect of Unionism on Fringe Benefits[J]. Industrial and Labor Relations Review, 1981, 34 (4):289-309.

[10]Kahn, L. M., Morimune, K. Unions and Employment Stability: A Sequential Logit Approach[J]. International Economic Review，1979, 20(2): 217-235.

[11]魏下海，董志強，金釗. 工會改善了企業(yè)雇用期限結構嗎？——來自全國民營企業(yè)抽樣調(diào)查的經(jīng)驗數(shù)據(jù)[J].管理世界，2015,(5):52-62.

[12]Clark, A. E., Oswald, A. J. Satisfaction and Comparison Income[J]. Journal of Public Economics, 1996, 61(3):359-381.

[13]Diener, E. Subjective Well-Being: There Decades of Progress[J].Psychology Bulletin，1999,125(2): 276-302.

[14]Frey, B. S., Stutzer, A. What Can Economists Learn From Happiness Research?[J]. Journal of Economic Literature, 2002, 40(2)：402-435.

[15]李儒林，張進輔，梁新剛. 影響主觀幸福感的相關因素理論[J].中國心理衛(wèi)生雜志，2003,(11):783-785.

[16]張進，馬月婷. 幸福感概念、測量及其與工作效能變量的關系[J]. 中國軟科學，2007,(5):60-68.

[17]魯元平，王韜. 幸福感影響因素研究評述[J]. 經(jīng)濟學動態(tài)，2010,(5):125-130.

[18]楊勝利，謝超. 就業(yè)質(zhì)量對居民幸福感的影響分析研究[J].云南財經(jīng)大學學報，2015，(6)：50-57.

[19]吳愈曉，王鵬，黃超. 家庭庇護、體制庇護與工作家庭沖突——中國城鎮(zhèn)女性的就業(yè)狀態(tài)與幸福感[J]. 社會學研究，2015，(6)：122-144.

[20]Frey,B.S.,Stutzer, A. Happiness and Economics: How the Economy and Institutions Affect Human Well-Being[M]. Princeton: Princeton University Press, 2010.

[21]Sajaia,Z. Maximum Likelihood Estimation of a Bivariate Ordered Probit Model: Implementation and Monte Carlo Simulations[J]. The Stata Journal, 2008，4(2): 1-18.