趙向陽(yáng)，孫愛(ài)峰，孫齊蕊

1.白城市疾病預(yù)防控制中心，吉林白城 137000；2.白城市衛(wèi)生急救中心，吉林白城 137000；3.吉林大學(xué)公共衛(wèi)生學(xué)院，吉林長(zhǎng)春130021

碘缺乏病是我國(guó)重點(diǎn)防制的一種生物地球化學(xué)性疾病。疾病監(jiān)測(cè)是碘缺乏病防制的一種重要措施，包括隨機(jī)抽樣調(diào)查居民食用碘鹽監(jiān)測(cè)、大型集體用餐重點(diǎn)單位碘鹽食用情況監(jiān)測(cè)、重點(diǎn)人群尿碘水平監(jiān)測(cè)和兒童碘缺乏病病情監(jiān)測(cè)，各方面又有相應(yīng)評(píng)價(jià)指標(biāo)?？紤]到某地2007—2012年重點(diǎn)單位碘鹽食用情況和兒童碘缺乏病病情監(jiān)測(cè)之間的差異無(wú)統(tǒng)計(jì)學(xué)意義，因此，筆者選擇居民食用碘鹽監(jiān)測(cè)和重點(diǎn)人群尿碘水平監(jiān)測(cè)結(jié)果為基本資料，應(yīng)用秩和比法來(lái)綜合評(píng)價(jià)某地2007—2012年相應(yīng)監(jiān)測(cè)數(shù)據(jù)，為做好疾病預(yù)防控制工作提供客觀的參考依據(jù)?，F(xiàn)報(bào)道如下。

1 資料與方法

1.1 一般資料

以某地按照監(jiān)測(cè)方案所獲得的2007—2012年居民食用碘鹽監(jiān)測(cè)和重點(diǎn)人群尿碘水平監(jiān)測(cè)數(shù)據(jù)為基本資料，其中X1-碘鹽覆蓋率（%）、X2-碘鹽合格率（%）和X3-合格碘鹽食用率（%）為評(píng)價(jià)居民食用碘鹽的指標(biāo)，X4-育齡婦女尿碘適宜構(gòu)成比（%）、X5-兒童尿碘適宜構(gòu)成比（%）和X6-孕婦尿碘適宜構(gòu)成比（%）為評(píng)價(jià)重點(diǎn)人群尿碘水平的指標(biāo)，所得各年度的監(jiān)測(cè)值見(jiàn)表1[1]。監(jiān)測(cè)數(shù)據(jù)來(lái)自文獻(xiàn)，真實(shí)可靠。

表1 某地居民食用碘鹽和重點(diǎn)人群尿碘監(jiān)測(cè)結(jié)果及其秩變換

1.2 統(tǒng)計(jì)方法

應(yīng)用基于秩和比的可信區(qū)間法[2]分析表1 資料。首先計(jì)算由居民食用碘鹽和重點(diǎn)人群尿碘水平監(jiān)測(cè)結(jié)果所得各年度的秩和比RSR1 和RSR2，如兩者Pearson 相關(guān)系數(shù)呈現(xiàn)正相關(guān)或無(wú)相關(guān)關(guān)系，則計(jì)算合計(jì)秩和比RSR，按可信區(qū)間法對(duì)RSR 進(jìn)行綜合評(píng)價(jià)；如兩者之間呈負(fù)相關(guān)關(guān)系，則按可信區(qū)間法分別對(duì)RSR1 和RSR2進(jìn)行綜合評(píng)價(jià)。

計(jì)算各年度平方根反正弦代換值y 的95%CI 時(shí)，為減少I(mǎi) 類(lèi)錯(cuò)誤，按Bonferroni 法[3]校正檢驗(yàn)水準(zhǔn)。

2 結(jié)果

2.1 秩和比的計(jì)算

有專(zhuān)業(yè)知識(shí)可知，評(píng)價(jià)指標(biāo)均為高優(yōu)，對(duì)同一指標(biāo)按從小到大順序編秩（見(jiàn)表1 括號(hào)內(nèi)數(shù)字），分別計(jì)算RSR1和RSR2，經(jīng)計(jì)算兩者相關(guān)系數(shù)為0.1832，兩者之間差異無(wú)統(tǒng)計(jì)學(xué)意義（P＞0.05）；進(jìn)而計(jì)算RSR，并按大小順序排序，見(jiàn)表2。

表2 各年度居民食用碘鹽和重點(diǎn)人群尿碘監(jiān)測(cè)的秩和比和排序

2.2 綜合評(píng)價(jià)

計(jì)算各年度RSR 的平方根反正弦代換值y，按照公式y(tǒng)±3.1440 Sy 計(jì)算y 的95%CI，見(jiàn)表3。

表3 各年度居民食用碘鹽和重點(diǎn)人群尿碘監(jiān)測(cè)結(jié)果的綜合評(píng)價(jià)

比較各年度y 的95%CI 可知，2008年與2009年、2010年，2012年與2007年、2009年～2011年兩兩年度之間的差異有統(tǒng)計(jì)學(xué)意義（P＜0.05）；其余年度兩兩之間均差異無(wú)統(tǒng)計(jì)學(xué)意義（P＞0.05），見(jiàn)表4。

表4 兩兩年度之間差異是否具有統(tǒng)計(jì)學(xué)意義的判斷

3 討論

居民食用碘鹽和重點(diǎn)人群尿碘監(jiān)測(cè)是防制碘缺乏病的重要措施，通過(guò)監(jiān)測(cè)獲得的數(shù)據(jù)可以計(jì)算相應(yīng)統(tǒng)計(jì)指標(biāo)，采用適宜的統(tǒng)計(jì)方法在對(duì)監(jiān)測(cè)數(shù)據(jù)進(jìn)行綜合評(píng)價(jià)，用以提高監(jiān)測(cè)質(zhì)量和優(yōu)化監(jiān)測(cè)方案，更可為碘缺乏病的防制提供可靠的參考依據(jù)，阻止疾病在人群中的發(fā)生、發(fā)展和流行。

秩和比法是由原中國(guó)預(yù)防醫(yī)學(xué)科學(xué)院田鳳調(diào)教授提出的一種融非參數(shù)統(tǒng)計(jì)和參數(shù)統(tǒng)計(jì)于一體的統(tǒng)計(jì)方法，在疾病監(jiān)測(cè)數(shù)據(jù)的統(tǒng)計(jì)分析中得到廣泛應(yīng)用[4-7]。在應(yīng)用秩和比法時(shí)，秩次變換應(yīng)區(qū)分為高優(yōu)還是低優(yōu)指標(biāo)，保證計(jì)算得到的秩和比越大越好；要根據(jù)研究目的的不同和被研究對(duì)象數(shù)量的多少，選擇可信區(qū)間法或分檔法對(duì)資料進(jìn)行綜合評(píng)價(jià)。

當(dāng)比較組數(shù)較少時(shí)，應(yīng)選擇可信區(qū)間法進(jìn)行兩兩比較。為減少I(mǎi) 類(lèi)錯(cuò)誤的發(fā)生，在計(jì)算RSR 變量變換值時(shí)應(yīng)對(duì)檢驗(yàn)水準(zhǔn)進(jìn)行校正，兩兩比較的準(zhǔn)則是看變量變換值95%CI 交叉數(shù)值是否超過(guò)可信區(qū)間范圍的一半，超過(guò)的話(huà)說(shuō)明在設(shè)定的檢驗(yàn)水準(zhǔn)下，兩者之間的差異無(wú)統(tǒng)計(jì)學(xué)意義，否則差異有統(tǒng)計(jì)學(xué)意義。當(dāng)比較組數(shù)較多時(shí)，應(yīng)選擇分檔法進(jìn)行聚類(lèi)分析。分檔法要求按照線性分割或均數(shù)±標(biāo)準(zhǔn)差等方法至少分為2 檔（等級(jí)），且各檔方差齊，應(yīng)用方差分析和多重比較法表明各檔之間均數(shù)及其兩兩均數(shù)之間的差異均有統(tǒng)計(jì)學(xué)意義，即為最佳分檔?？梢?jiàn)基于秩和比的可信區(qū)間法和分檔法的原理的不同的。

表1 資料給出的評(píng)價(jià)指標(biāo)分屬居民食用碘鹽和重點(diǎn)人群尿碘監(jiān)測(cè)，直接將其編秩計(jì)算秩和比是不合理的。此時(shí)應(yīng)計(jì)算類(lèi)秩和比RSR1 和RSR2 的Pearson 或其他相關(guān)系數(shù)，如呈正相關(guān)或無(wú)相關(guān)關(guān)系則計(jì)算合計(jì)秩和比RSR 并進(jìn)行綜合評(píng)價(jià)；如呈負(fù)相關(guān)關(guān)系，為避免合計(jì)秩和比所產(chǎn)生的相互抵消作用，則分別對(duì)RSR1 和RSR2 進(jìn)行綜合評(píng)價(jià)。由于該文資料所得Pearson 相關(guān)系數(shù)為0.1832，兩者之間差異無(wú)統(tǒng)計(jì)學(xué)意義（P＞0.05），所以應(yīng)進(jìn)一步計(jì)算各年度的合計(jì)秩和比，以客觀反映各年度居民食用碘鹽和重點(diǎn)人群尿碘監(jiān)測(cè)數(shù)據(jù)表現(xiàn)出的綜合作用。值得注意的是，通常認(rèn)為居民食用碘鹽和重點(diǎn)人群尿碘水平之間應(yīng)該是呈正相關(guān)關(guān)系的，但該文卻得出無(wú)相關(guān)關(guān)系的結(jié)論，這里的原因可能是多方面的，值得進(jìn)一步開(kāi)展相關(guān)流行病學(xué)調(diào)查進(jìn)行研究和探討。

該資料計(jì)算得到2007年和2011年的RSR 均為4/9，這是由于秩次變換只考慮數(shù)值之間的大小關(guān)系所致。此時(shí)可采用其他秩次變換方法，由于數(shù)值保留到千分位均可做進(jìn)位處理，筆者試以2001年各指標(biāo)數(shù)值的增長(zhǎng)幅度合計(jì)數(shù)值大小決定哪一年度進(jìn)位哪一年度不進(jìn)位，由于該年度的增長(zhǎng)率為-13.43%，所以2007年進(jìn)位，2011年不進(jìn)位，這同時(shí)也滿(mǎn)足了秩和比加和值為（n+1）/2 的基本理論，更解決了兩個(gè)年度排序相同問(wèn)題。

基于以上分析可知，2007—2012年某地居民食用碘鹽和重點(diǎn)人群尿碘監(jiān)測(cè)綜合排序結(jié)果為2012年＞2008年＞2007年＞2011年＞2010年＞2009年，排序結(jié)論與TOPSIS 法完全相同。該方法的優(yōu)勢(shì)在于可以將計(jì)算得到的秩和比進(jìn)行兩兩比較，獲得比文獻(xiàn)[1]更為豐富的信息。同時(shí)，當(dāng)年度數(shù)量較多時(shí)，尚可進(jìn)行發(fā)展趨勢(shì)和統(tǒng)計(jì)預(yù)測(cè)預(yù)警分析，這更有利于根據(jù)統(tǒng)計(jì)分析結(jié)果來(lái)提高監(jiān)測(cè)質(zhì)量和優(yōu)化監(jiān)測(cè)方案，為碘缺乏病的防制提供可靠的參考依據(jù)，達(dá)到阻止疾病在人群中的發(fā)生、發(fā)展和流行的目的。

碘缺乏病在我國(guó)廣泛流行?！妒秤名}碘含量》（GB26878-2011），中規(guī)定食用鹽產(chǎn)品（碘鹽）中碘含量的平均水平（以碘元素計(jì)）為20～30 mg/kg，碘含量的允許波動(dòng)范圍為±30%；各省、自治區(qū)、直轄市人民政府衛(wèi)生行政部門(mén)在規(guī)定的范圍內(nèi)，根據(jù)當(dāng)?shù)厝巳簩?shí)際碘營(yíng)養(yǎng)水平，選擇適合該地情況的食用鹽碘含量平均水平。要實(shí)現(xiàn)“因地制宜，分類(lèi)指導(dǎo)，科學(xué)補(bǔ)碘”的防治策略，必須按照國(guó)家規(guī)定加強(qiáng)碘缺乏病監(jiān)測(cè)工作，不斷提高監(jiān)測(cè)質(zhì)量，進(jìn)行科學(xué)統(tǒng)計(jì)分析，為疾病防制提供強(qiáng)有力的數(shù)據(jù)支撐。

[1]黎新宇,李陽(yáng)樺,杜丹,等.采用TOPSIS 法對(duì)北京市2007-2012年碘缺乏病檢測(cè)結(jié)果綜合評(píng)價(jià)[J].中國(guó)地方病防治雜志,2014,29(3):185-186,190.

[2]孫愛(ài)峰.秩和比法在不同地區(qū)人體蛔蟲(chóng)感染現(xiàn)狀調(diào)查結(jié)果綜合評(píng)價(jià)中的應(yīng)用[J].中國(guó)醫(yī)藥指南,2011,9(31):226-227.

[3]方積乾.醫(yī)學(xué)統(tǒng)計(jì)學(xué)與電腦實(shí)驗(yàn)[M].第4 版.上海：上?？茖W(xué)技術(shù)出版社,2012.

[4]張輝,孫愛(ài)峰.秩和比法在吉林省省級(jí)大骨節(jié)病病情監(jiān)測(cè)結(jié)果綜合評(píng)價(jià)中的應(yīng)用[J].中國(guó)醫(yī)藥指南,2012,10(18):377-378.

[5]姜芳晶,楊維中.醫(yī)學(xué)常用綜合評(píng)價(jià)方法[J].疾病監(jiān)測(cè),2006,21(6):325-328.

[6]薄桂云,孫愛(ài)峰.秩和比法在某市2006～2010年居民戶(hù)碘鹽監(jiān)測(cè)結(jié)果綜合評(píng)價(jià)中的應(yīng)用[J].中國(guó)醫(yī)藥指南,2013,11(18):773-774.

[7]王艷萍,孫愛(ài)峰,孫齊蕊.秩和比法在中國(guó)大陸居民用戶(hù)碘鹽監(jiān)測(cè)結(jié)果綜合評(píng)價(jià)中的應(yīng)用[J].中國(guó)醫(yī)藥指南,2014,12(35):385-386.