陳 微 趙克勤
1諸暨市中醫(yī)醫(yī)院眼科,311800 浙江 諸暨; 2諸暨市聯(lián)系數(shù)學(xué)研究所,311800 浙江 諸暨
近視是全世界共同關(guān)注的公共衛(wèi)生問題之一,也是我國青少年中一個突出的公共衛(wèi)生問題[1-2]。我國是一個人口大國,青少年近視發(fā)病率居世界前列,且呈逐年遞增趨勢[3]。雖然我國對中小學(xué)生近視的監(jiān)測在省級層面開展多年,對監(jiān)測結(jié)果分析也有文獻報道[4-6],但鮮少有對監(jiān)測數(shù)據(jù)做近視發(fā)展趨勢研究的報告。聯(lián)系數(shù)是趙克勤[7-8]在集對分析中給出的一種具有層次結(jié)構(gòu)的函數(shù)。偏聯(lián)系數(shù)是聯(lián)系數(shù)的一種伴隨函數(shù),揭示了聯(lián)系數(shù)中處于不同層次的聯(lián)系分量所蘊含的發(fā)展趨勢信息,現(xiàn)已應(yīng)用于多個領(lǐng)域[9-11]。如在醫(yī)療衛(wèi)生領(lǐng)域中,覃杰等[12]把偏聯(lián)系數(shù)用于醫(yī)院醫(yī)療質(zhì)量發(fā)展趨勢綜合分析,丁曉杰等[13]把偏聯(lián)系數(shù)用于銀屑病中西醫(yī)治療療效穩(wěn)定性分析,邢夢等[14]用偏聯(lián)系數(shù)探討銀屑病的復(fù)發(fā)因素,李斌等[15]在中醫(yī)辨證論治中應(yīng)用偏聯(lián)系數(shù)。本研究以某省中小學(xué)生為研究對象,嘗試將偏聯(lián)系數(shù)應(yīng)用于中小學(xué)生近視發(fā)展趨勢的分析,為科學(xué)預(yù)防中小學(xué)生近視,有效減緩近視發(fā)展趨勢提供參考。
數(shù)據(jù)來自文獻[5]中報道的某省2013年6歲到18歲中小學(xué)生共1 725 239人的近視檢出率。見表1。
表1 某省2013年中小學(xué)生近視檢出率
首先,利用表1中數(shù)據(jù)構(gòu)建三元聯(lián)系數(shù):
μt=at+bti+ctj(1)
式(1)中,t表示年齡組;μt表示t年齡組學(xué)生的整體視力情況;at表示t年齡組學(xué)生的輕度近視檢出率;bt表示該組學(xué)生中度近視檢出率;ct表示該組學(xué)生高度近視檢出率;at、bt、ct為三元聯(lián)系數(shù)μt的3個聯(lián)系分量,各在區(qū)間[0,1]內(nèi)取值,且at+bt+ct≤1;j是ct的示性系數(shù),在集對分析理論中常取j=-1,這里可以理解為高度近視ctj處在-1層次,與處在+1層次的輕度近視at具有正負(fù)對立性;i是bt的示性系數(shù),在集對分析理論中規(guī)定i在[-1,1]區(qū)間視不同情況取值,表示bt處在at與ctj之間,但向哪一個層次轉(zhuǎn)化有不確定性,這里可以理解為中度近視向高度近視發(fā)展還是向輕度近視發(fā)展有不確定性;at、bti、ctj依次稱為聯(lián)系數(shù)μt的“正同部”“異部”“負(fù)反部”,簡稱“同異反”。
易見聯(lián)系數(shù)μt具有整體性和系統(tǒng)層次性,在完整地刻畫t年齡組學(xué)生近視整體狀況的同時又具體刻畫了不同程度近視學(xué)生的檢出率,見式(2):
(2)
聯(lián)系數(shù)μt的層次性還蘊含著其中的同異反3個聯(lián)系分量at、bti、ctj如下的“矛盾運動”。
從理想的正向發(fā)展角度看,不妨假定聯(lián)系數(shù)μt=at+bti+ctj中的at(輕度近視),原本處在中度近視的bt層次,是從bt層次發(fā)展而來,為此用at+bt作分母,at作分子,用
(3)
表示原本處在bt層次上的中度近視發(fā)展到輕度近視at層次上的比例, ?+at讀為偏正at。同理,用
(4)
表示原本處在ct層次上的高度近視發(fā)展到中度近視bt層次上的比例,?+bt讀為偏正bt。令
(5)
稱?+μt是聯(lián)系數(shù)μt=at+bti+ctj的一階偏正聯(lián)系數(shù),刻畫了μt在一階微觀層次上朝近視程度減輕方向(正方向at)的發(fā)展趨勢。
對一階偏正聯(lián)系數(shù)?+μt=?+at+?+bti再求一次偏正聯(lián)系數(shù),得μt的二階偏正聯(lián)系數(shù)?2+μt,
?2+μt=?+(?+at+?+bti)=
(6)
式(6)的二階偏正聯(lián)系數(shù)?2+μt刻畫了μt在二階微觀層次上朝近視程度減輕方向(正方向at)的發(fā)展趨勢。
反之,從近視朝加重方向發(fā)展的負(fù)向角度看,也就是μt以負(fù)方向(ctj)為參考方向時,可以假定當(dāng)前的中度近視bt原本處在輕度近視at,是從輕度近視at負(fù)向發(fā)展而來,為此用at+bt作分母,bt作分子,用
(7)
表示原本處在at層次上的輕度近視發(fā)展到中度近視bt層次上的比例,?-bt讀為偏負(fù)bt。同理,用
(8)
表示原本處在bt層次上的中度近視發(fā)展到高度近視ct層次上的比例,?-ct讀為偏負(fù)ct。令
(9)
稱?-μt是聯(lián)系數(shù)μt=at+bti+ctj的一階偏負(fù)聯(lián)系數(shù),刻畫了μt在一階微觀層次上朝近視程度加重方向(負(fù)方向ctj)發(fā)展時的一種趨勢。對?-μt=?-bti-+?-ctj再求一次偏負(fù)聯(lián)系數(shù),得μt的二階偏負(fù)聯(lián)系數(shù)?2-μt:
?2-μt=?-(?-μt)=?-(?-bti+?-ctj)=
(10)
?2-μt刻畫μt在二階微觀層次上朝近視程度加重方向(負(fù)方向ctj)發(fā)展的趨勢。
由于聯(lián)系數(shù)μt=at+bti+ctj的實際趨勢是聯(lián)系數(shù)中聯(lián)系分量正向運動和負(fù)向運動相互作用的綜合結(jié)果,為此用
?2±μt=?2+μt+?2-μt
(11)
表示這個綜合結(jié)果,稱?2±μt為聯(lián)系數(shù)μt=at+bti+ctj的二階正負(fù)全偏聯(lián)系數(shù),即
(12)
由于設(shè)j=-1,所以式(12)可以改寫成
(13)
把式(6)和式(10)代入式(13)得:
(14)
式(14)顯示?2±μt中沒有示性系i,是一個確定的實數(shù),因此能明確地指示出μt中各個聯(lián)系分量正向運動與負(fù)向運動的綜合趨勢。當(dāng)?2±μt>0時聯(lián)系數(shù)μt具有確定的二階正向發(fā)展趨勢,當(dāng)?2±μt<0時聯(lián)系數(shù)μt具有確定的二階負(fù)向發(fā)展趨勢,當(dāng)?2±μt=0時聯(lián)系數(shù)μt是確定的二階臨界趨勢。
以表1中各年齡組學(xué)生的輕度近視檢出率為聯(lián)系數(shù)μt的正同部at,高度近視率為負(fù)反部ctj中的ct,示性系數(shù)j表示高度近視與輕度近視具有正負(fù)相反性;中度近視率為介于“正同部”與“負(fù)反部”之間的異部bti中的bt,示性系數(shù)i表示中度近視是向高度近視發(fā)展還是向輕度近視發(fā)展存在不確定性。由此得到各年齡組學(xué)生近視檢出率的同異反三元聯(lián)系數(shù)μt=at+bti+ctj(t=6,7,8,…,18)。見表2。
表2 某省2013年中小學(xué)生近視檢出率聯(lián)系數(shù)
表2 (續(xù))
第2 步,計算表2中各個同異反三元聯(lián)系數(shù)μt=at+bti+ctj(t=6,7,8,…,18)的一階偏正聯(lián)系數(shù)
表3 近視檢出率聯(lián)系數(shù)的一階偏正聯(lián)系數(shù)與一階偏負(fù)聯(lián)系數(shù)
第3 步,計算表3中各個一階偏正聯(lián)系數(shù)的二階偏正聯(lián)系數(shù),各個一階偏負(fù)聯(lián)系數(shù)的二階偏負(fù)聯(lián)系數(shù)。見表4。
表4 男女學(xué)生近視檢出率聯(lián)系數(shù)的二階全偏聯(lián)系數(shù)及近視發(fā)展趨勢
第4步,計算各年齡男生女生近視檢出率的二階全偏聯(lián)系數(shù),根據(jù)二階全偏聯(lián)系數(shù)的正負(fù)作出近視趨勢判斷,正數(shù)對應(yīng)近視減輕,負(fù)數(shù)對應(yīng)近視加重。見表4。
第5步,根據(jù)表4中各年齡學(xué)生近視檢出率聯(lián)系數(shù)的趨勢判斷做出結(jié)論,給出決策建議。
結(jié)果顯示,該省6~9歲的小學(xué)生,視力狀況以近視微量減輕為趨勢;10~18歲的中小學(xué)生,視力狀況以近視加重為趨勢;提示對于剛?cè)胄W(xué)學(xué)習(xí)的學(xué)生,保護視力、防范近視,具有自然優(yōu)勢。其中10歲時男生的二階全偏聯(lián)系數(shù)?2±μ10為-0.018,女生二階全偏聯(lián)系數(shù)?2±μ10為-0.034,女生的二階全偏聯(lián)系數(shù)是男生的1.89倍,提示10歲是一個轉(zhuǎn)折年齡,從10歲開始,近視發(fā)展趨勢呈現(xiàn)為近視加重;男生?2±μ10-13從-0.018增加到-0.065,女生?2±μ10-13從-0.034增加到-0.063,兩者無明顯差異,提示10歲到13歲這個年齡段中的男生與女生的近視加重趨勢每年遞增;16歲男生的二階全偏聯(lián)系數(shù)?2±μ16如下,
小于16歲女生的-0.050。同理,17歲男生的二階全偏聯(lián)系數(shù)?2±μ17=-0.057小于17歲女生的-0.038;18歲男生的二階全偏聯(lián)系數(shù)?2±μ18=-0.056小于18歲女生的-0.041,提示16~18歲年齡段的男生近視加重趨勢明顯超過女生。除6歲組男生的二階全偏聯(lián)系數(shù)?2±μ6為0.123,女生的二階全偏聯(lián)系數(shù)?2±μ6為0.139外,其他7~18歲的男女學(xué)生的二階全偏聯(lián)系數(shù)?2±μ7-18的絕對值都小于0.1。見表4。
本研究顯示,10歲前中小學(xué)生存在向減輕近視方向發(fā)展的趨勢,10歲及以后年齡段此趨勢變小,這一結(jié)果與事實相符;可能與10歲前的中小學(xué)生學(xué)業(yè)負(fù)擔(dān)輕,睡眠充足,課外活動充足,看書姿勢正確等有關(guān);個別學(xué)生有近視也屬假性近視或擬似近視,僅極少數(shù)因父母近視而得真性近視。提示要高度重視10歲小學(xué)生的近視防控意義,要特別重視女小學(xué)生的近視防控;教育部于近期實施的小學(xué)生“雙減”措施將對減緩小學(xué)生近視起強大的推動作用和保證作用。
其次,10~13歲年齡段中的中小學(xué)生有近視加重趨勢,表現(xiàn)為二階全偏聯(lián)系數(shù)?2±μt年年遞增??陀^上與小學(xué)生的學(xué)習(xí)量增加,用眼過度有關(guān),因此10歲到13歲這個年齡段是小學(xué)生近視防控的關(guān)鍵階段。而14~18歲年齡段,中小學(xué)生近視趨勢的二階全偏聯(lián)系數(shù)?2±μt緩慢遞減。
研究還顯示, 16~18歲年齡段的男生近視加重趨勢明顯超過女生,另一個現(xiàn)象是6歲組男生和女生的二階全偏聯(lián)系數(shù)均大于0.1,而7~18歲的中小學(xué)生的二階全偏聯(lián)系數(shù)的絕對值都小于0.1,原因有待于進一步研究。
3.2 偏聯(lián)系數(shù)應(yīng)用于中小學(xué)生近視趨勢分析
首先,由于文中聯(lián)系數(shù)的聯(lián)系分量都是近視檢出率,因此,二階全偏聯(lián)系數(shù)?2±μ>0的物理意義是“近視穩(wěn)定偏好率”;?2±μ<0的物理意義是“近視加重率”,乘上相應(yīng)年齡的受檢學(xué)生數(shù),就能算得該年齡受檢學(xué)生中有減輕近視趨勢或加重近視趨勢的學(xué)生人數(shù)。如表4中10歲年齡組(男)的二階全偏聯(lián)系數(shù)是-0.018,表1顯示10歲年齡組的受檢人數(shù)是140 891人,則-0.018×14 0891=-2 536.038,取整數(shù)得-2 536,就是說10歲男學(xué)生中有2 536人的視力在受檢時就具有近視加重趨勢,這里的負(fù)號“-”表示“近視加重”,從而為進一步重點再查近視患者提供了依據(jù)。
其次,是文中沒有計算各個近視檢出率聯(lián)系數(shù)的一階全偏聯(lián)系數(shù)的值,原因是其中的一階偏正聯(lián)系數(shù)?+μ=?+a+?+bi+有不確定取值的i+,在一階偏負(fù)聯(lián)系數(shù)?-μ=?-bi-+?-cj中有不確定取值的i-,如何確定i+與i-的值,要涉及到更多的集對分析理論,限于篇幅,文中沒有介紹這方面的理論,建議有興趣的參考文獻[9-10]。但由一階偏正聯(lián)系數(shù)計算二階偏正聯(lián)系數(shù)的過程中,i+已被消除,消除原理是:作為分母的b,屬于“過去時狀態(tài)”因而是確定的b,所以不需要帶表示不確定性的示性系數(shù)i。同理,由一階偏負(fù)聯(lián)系數(shù)計算二階偏負(fù)聯(lián)系數(shù)的過程中i-也被消除。這也是不計算μt=at+bti+ctj的一階全偏聯(lián)系數(shù)的值,去計算μt=at+bti+ctj的二階全偏聯(lián)系數(shù)的一個原因。由此也窺知中小學(xué)生的近視減輕或近視加重,是一個非線性的過程,我們將結(jié)合臨床觀察和案例另行討論其中的病理和生理機制。
中小學(xué)生的近視變化趨勢是一個動態(tài)變化過程,這個過程中存在著相互矛盾的2個方面,一方面是“保護視力,減緩近視”各種措施的切實實施;另一方面是“用眼不當(dāng),加重近視”的種種不良因素作用。前者導(dǎo)致一部分原本是中高度近視學(xué)生緩解和減輕近視癥狀;后者誘發(fā)一部分原本是輕度近視學(xué)生發(fā)展成為中度近視,原本是中度近視學(xué)生發(fā)展成為高度近視。采用何種統(tǒng)計模型定量刻畫這兩方面的同時存在和相互作用并推理出其綜合后的趨勢,是一個對有效防控近視發(fā)展有重要指導(dǎo)意義但又有一定難度的問題。本文選取文獻[5]中報導(dǎo)的一個省級層面上的中小學(xué)生近視體檢大樣本(樣本容量為1 725 239)進行基于集對分析偏聯(lián)系數(shù)的方法進行研究,用其中的一階和二階偏正聯(lián)系數(shù)刻畫“保護視力,減緩近視”各種措施產(chǎn)生的正向測度;用其中的一階和二階偏負(fù)聯(lián)系數(shù)刻畫“用眼不當(dāng),加重近視”種種不良因素產(chǎn)生的負(fù)向測度;用二階全偏聯(lián)系數(shù)刻畫了上述正負(fù)2個方面相互矛盾又相互作用的綜合趨勢,在不增加體檢成本的條件下獲得了這個體檢樣本中蘊含的近視趨勢信息及其結(jié)構(gòu)信息。結(jié)果發(fā)現(xiàn)10歲到13歲是中小學(xué)生近視防控的關(guān)鍵階段,不僅是對文獻[5]所做統(tǒng)計分析有價值的補充,也為中小學(xué)生近視防控的早發(fā)現(xiàn)早防治提供了重要的決策參考,這是本文的一個創(chuàng)新。從醫(yī)院統(tǒng)計而言,應(yīng)用偏聯(lián)系數(shù),從一個樣本數(shù)據(jù)中挖掘出其中蘊含的系統(tǒng)變化趨勢,對其他有類似時空等級劃分的樣本數(shù)據(jù)做系統(tǒng)趨勢分析也有一定的啟發(fā)意義和參考意義。