• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    多分類敏感問題RRT模型下分層三階段抽樣的統(tǒng)計(jì)方法及應(yīng)用*

    2015-03-09 14:35:08蘇州大學(xué)公共衛(wèi)生學(xué)院215123楊少春歌范玉波陳科錦杜俏俏陳向宇施佳琛
    中國衛(wèi)生統(tǒng)計(jì) 2015年4期
    關(guān)鍵詞:安全套總體比例

    蘇州大學(xué)公共衛(wèi)生學(xué)院(215123) 楊少春 高 歌范玉波 陳科錦 杜俏俏 陳向宇 傅 穎 施佳琛

    多分類敏感問題RRT模型下分層三階段抽樣的統(tǒng)計(jì)方法及應(yīng)用*

    蘇州大學(xué)公共衛(wèi)生學(xué)院(215123) 楊少春 高 歌△范玉波 陳科錦 杜俏俏 陳向宇 傅 穎 施佳琛

    目的為實(shí)際需要的多分類敏感性問題的復(fù)雜抽樣提供信度高效度高的調(diào)查方法及其統(tǒng)計(jì)公式;為制訂艾滋病預(yù)防控制措施提供科學(xué)依據(jù)。方法根據(jù)抽樣理論、RRT模型、全概率公式等理論方法推導(dǎo)統(tǒng)計(jì)公式;對北京市MSM人群進(jìn)行實(shí)例調(diào)查;用SAS編程分別對6個(gè)調(diào)查指標(biāo)各類別,各模擬多分類敏感問題RRT模型下分層三階段抽樣調(diào)查100個(gè)樣本,按本文給出的統(tǒng)計(jì)公式計(jì)算100個(gè)總體比例的可信區(qū)間。結(jié)果推導(dǎo)出多分類敏感問題RRT模型下分層三階段抽樣各類別總體比例估計(jì)量及其方差的計(jì)算公式;調(diào)查得到北京市MSM人群6項(xiàng)敏感問題指標(biāo)各類別樣本比例及其標(biāo)準(zhǔn)誤。各指標(biāo)各類別各模擬得到的100個(gè)總體比例95%可信區(qū)間幾乎均包含其模擬總體比例。結(jié)論本文研究的敏感問題調(diào)查方法及其公式具有良好的信度和效度,值得推廣應(yīng)用;MSM人群具有艾滋病的高危性行為,應(yīng)加強(qiáng)對其的預(yù)防控制。

    敏感問題 隨機(jī)應(yīng)答技術(shù)(RRT) (分層)三階段抽樣 MSM效度與信度

    抽樣調(diào)查一直作為醫(yī)學(xué)科研及衛(wèi)生工作中的一種重要方法及獲取統(tǒng)計(jì)資料的重要手段[1-3]。在抽樣調(diào)查研究中,經(jīng)常會(huì)遇到各種各樣的敏感性問題。所謂敏感性問題是指高度私人機(jī)密性或大多數(shù)人認(rèn)為不便在公開場合表態(tài)及陳述的問題[4],例如:吸毒、賣淫、嫖娼、男男性行為、性病、艾滋病等。敏感性問題按總體的特征可分為分類特征敏感問題和數(shù)量特征敏感問題,分類特征敏感問題又分為兩分類敏感性問題和多分類敏感性問題[5]。隨機(jī)應(yīng)答技術(shù)(random ized response technique,RRT)巧妙地運(yùn)用一個(gè)隨機(jī)裝置,在保護(hù)個(gè)人隱私的情況下,通過概率計(jì)算估計(jì)敏感問題特征的比例或均值,消除減小調(diào)查對象回答敏感問題時(shí)的拒答或說謊偏倚。國外已將各種RRT廣泛使用于流產(chǎn)率、女性飲酒量以及偷稅漏稅等各類敏感問題的調(diào)查,并得到了較好的調(diào)查結(jié)果[6-10]。

    在本團(tuán)隊(duì)研究前,國內(nèi)外關(guān)于敏感問題調(diào)查的研究,主要局限于簡單隨機(jī)抽樣調(diào)查和隨機(jī)化裝置的改進(jìn)等理論研究[11-13]。本團(tuán)隊(duì)對非敏感問題各種復(fù)雜抽樣和敏感問題多種復(fù)雜抽樣的統(tǒng)計(jì)方法進(jìn)行過系列研究[14-15]。本文對敏感問題實(shí)際抽樣調(diào)查中常常采用的多分類敏感問題RRT模型下分層三階段抽樣:從數(shù)學(xué)上推導(dǎo)出其統(tǒng)計(jì)公式;在北京市男男性行為人群的調(diào)查中取得了成功的實(shí)際應(yīng)用效果;采用SAS編程進(jìn)行大樣本個(gè)數(shù)大樣本量的計(jì)算機(jī)模擬抽樣調(diào)查與分析,得到了高信度和高效度的評價(jià)結(jié)果。

    多分類敏感問題RRT模型下分層三階段抽樣的統(tǒng)計(jì)方法

    1.分層三階段抽樣方法

    假定總體劃分成L層。第h層由N1h個(gè)一級單位組成;h層第i個(gè)一級單位由Ni2h個(gè)二級單位組成,h層平均每個(gè)一級單位包含2h個(gè)二級單位;h層第i個(gè)一級單位內(nèi)第j個(gè)二級單位包含Nij3h個(gè)三級單位,h層平均每個(gè)二級單位包含3h個(gè)三級單位,h層共有Nh個(gè)三級單位,總體中共有N個(gè)三級單位;這里,h=1,2,…,L;i=1,2,…,N1h;j=1,2,…,Ni2h。第一階段從h層隨機(jī)抽取n1h個(gè)一級單位;第二階段從h層第i個(gè)被抽中的一級單位內(nèi)隨機(jī)抽取ni2h個(gè)二級單位,平均從h層每個(gè)被抽中的一級單位內(nèi)隨機(jī)了抽取2h個(gè)二級單位;第三階段從h層第i個(gè)抽中的一級單位內(nèi)第j個(gè)抽中的二級單位中隨機(jī)抽取nij3h個(gè)三級單位,平均從h層每個(gè)抽中的二級單位中隨機(jī)抽取了3h個(gè)三級單位;這里h=1,2,…,L,i=1,2,…,n1h,j=1,2,…,ni2h。對每個(gè)被抽中的三級單位(人),采用多分類敏感問題隨機(jī)應(yīng)答模型進(jìn)行調(diào)查。

    2.多分類敏感問題隨機(jī)應(yīng)答模型

    設(shè)某敏感問題分為K種互斥的類別1,2,…,K。設(shè)計(jì)一隨機(jī)化裝置,如:將分別寫有0,1,2,…,K的K+1種按數(shù)量比例P0:P1:P2.…Pk(P0+P1+P2+…+Pk=1)的若干卡片混合放入袋中。無旁人在場時(shí)每個(gè)抽中的人(三級單位)有放回地從袋中隨機(jī)抽出一張卡片,若卡片上寫有0則真實(shí)回答自己屬于敏感問題的那一類的序號;若卡片上寫有0以外的某個(gè)數(shù)則回答該數(shù)。

    3.統(tǒng)計(jì)公式

    (1)總體比例的估計(jì)量及其方差計(jì)算公式

    假定pk為第k類敏感特征總體比例的估計(jì)量,ph-k為h層第k類總體比例ph-k的估計(jì)量,pih-k為h層第i個(gè)抽中一級單位第k類總體比例pijh-k的估計(jì)量,pih-k為h層第i個(gè)抽中一級單位內(nèi)第j個(gè)抽中二級單位第k類總體比例pijh-k的估計(jì)量,h=1,2,…,L,i=1,2,…,n1h,j=1,2,…,ni2h,k=1,2,…,nij3h。

    根據(jù)抽樣理論,雖然原始分類多于兩類,但當(dāng)調(diào)查的目的是要估計(jì)總體中任何一類的個(gè)數(shù)占總數(shù)的比例時(shí),實(shí)質(zhì)上這些比例是按兩類得到的,所以兩分類的理論均適合于此種情況[2]。根據(jù)高歌、金丕煥給出的公式可得[1]:

    (2)pijh-k的計(jì)算公式

    設(shè)pijh-k為h層第i個(gè)一級單位第j個(gè)二級單位k類敏感問題的總體比例,mijh-k、λijh-k分別表示h層樣本中第i個(gè)一級單位第j個(gè)二級單位回答數(shù)字k的頻數(shù)、概率,根據(jù)全概率公式:

    應(yīng)用實(shí)例

    1.調(diào)查設(shè)計(jì)

    (1)調(diào)查對象與調(diào)查指標(biāo)

    調(diào)查對象為2010年9-10月出入北京市男同性戀活動(dòng)場所的15~49歲男男性行為(men who have sex with men,MSM)人群,總體數(shù)量約67750人。調(diào)查指標(biāo)1為男男性行為各方式的比例;調(diào)查指標(biāo)2為安全套全程使用情況各類別的比例;調(diào)查指標(biāo)3為HIV病毒檢測結(jié)果各類別的比例;調(diào)查指標(biāo)4為性病檢測結(jié)果各類別的比例;調(diào)查指標(biāo)5為最近一年男男商業(yè)性性行為次均費(fèi)用不同等級的比例;調(diào)查指標(biāo)6為安全套破損情況各類別的比例。

    (2)抽樣方法

    采用分層三階段抽樣,將北京市MSM人群按年齡分為15~29歲和30~49歲兩層,各自所占比例分別為W1=58.24%,W2=41.76%;在各層內(nèi)分別以北京的16個(gè)區(qū)縣為一級單位(N11=N12=16),以男同性戀活動(dòng)場所(如男同浴室、會(huì)所、酒吧、公園等)為二級單位(調(diào)查點(diǎn)),以男男性行為者作為三級單位。兩個(gè)年齡層內(nèi)平均每個(gè)區(qū)縣男同性戀人數(shù)分別約為2466人=1768人。第一階段從各年齡層的16個(gè)區(qū)縣中分別隨機(jī)抽取3個(gè)區(qū)縣(n11=n12=3);第二階段在所抽取的區(qū)縣各隨機(jī)抽取5個(gè)調(diào)查點(diǎn)(n=21=n22=5);第三階段從被抽中的5個(gè)調(diào)查點(diǎn)中共隨機(jī)抽取2062人,在各層內(nèi)平均從每個(gè)抽中調(diào)查點(diǎn)分別抽取了

    (3)RRT模型

    采用多分類敏感問題隨機(jī)應(yīng)答模型進(jìn)行調(diào)查,以商業(yè)性男男性行為平均費(fèi)用等級比例調(diào)查為例,設(shè)計(jì)一套隨機(jī)化裝置:將10個(gè)大小、重量、觸感完全相同的小球貼上數(shù)字標(biāo)簽,5個(gè)寫上數(shù)字0、剩余5個(gè)分別寫上數(shù)字1、2、3、4、5,即P0:P1:P2:P3:P4:P5為0.5:0.1:0.1:0.1:0.1:0.1(P0+P1+P2+P3+P4+P5=1),10個(gè)小球混合放入袋中;每個(gè)被抽中的三級單位有放回地從袋中隨機(jī)抽出一個(gè)小球,若抽中的小球上寫有0,真實(shí)回答自己“最近一年商業(yè)性男男性行為的平均費(fèi)用的等級序號k;若抽中的小球上寫有0以外的數(shù)字則直接回答該數(shù)字。

    3.調(diào)查結(jié)果

    按式(1)至(9)計(jì)算得:最近一年男男商業(yè)性性行為的次均費(fèi)用<200元、200~399元、400~599元、600元以上、無男男商業(yè)性性行為的比例分別為5.65%、4.86%、2.67%、6.48%、80.34%,其標(biāo)準(zhǔn)誤分別為0.70%、0.49%、0.57%、0.76%、2.20%;最近一個(gè)月男男肛交中從不全程使用安全套、有時(shí)全程使用安全套、每次全程使用安全套、沒有肛交的比例分別為6.46%、31.78%、51.81%、9.96%,其標(biāo)準(zhǔn)誤分別為0.77%、0.71%、0.95%、0.44%;最近一年HIV病毒檢測結(jié)果為陽性、陰性、結(jié)果不確定、沒有檢查過的比例分別為6.31%、78.67%、2.89%、6.45%,其標(biāo)準(zhǔn)誤分別為0.42%、1.96%、0.28%、0.45%;最近一年到正規(guī)醫(yī)療機(jī)構(gòu)性病檢測結(jié)果有性病、無性病、結(jié)果不確定、沒有檢查過的比例分別為16.79%、67.02%、6.25%、5.23%,其標(biāo)準(zhǔn)誤分別為0.66%、1.98%、0.46%、0.33%;男男性行為的方式為肛交、口交、其他的比例分別為65.83%、18.70%、15.47%,其標(biāo)準(zhǔn)誤分別為0.930%、0.658%、0.316%;安全套使用有破損、無破損、沒有使用安全套的比例分別為5.49%、70.26%、24.25%,其標(biāo)準(zhǔn)誤分別為0.57%、2.89%、1.89%。

    信度效度評價(jià)

    1.建立模擬總體

    按照北京市的區(qū)縣數(shù)、北京市MSM活動(dòng)場所數(shù)、出入北京市男同性戀活動(dòng)場所的15~49歲MSM人數(shù)(67750人)及年齡構(gòu)成(15~29歲、30~49歲所占比例分別為W1=58.24%、W2=41.76%),以本文調(diào)查的6項(xiàng)多分類敏感問題各類別樣本比例、二項(xiàng)分布方差的數(shù)值作為模擬總體參數(shù),基于蒙特卡洛方法,用SAS編程建立模擬總體。模擬總體的一級單元為16個(gè)區(qū)縣,模擬總體的二級單元為男同性戀活動(dòng)場所(調(diào)查點(diǎn)),模擬總體的三級單元為男男性行為者。

    2.模擬抽樣調(diào)查

    模擬多分類敏感問題隨機(jī)應(yīng)答模型下分層三階段抽樣調(diào)查。第一階段模擬從16個(gè)區(qū)縣隨機(jī)抽取3個(gè)一級單元(區(qū)縣);第二階段抽樣從每個(gè)抽中的一級單元中隨機(jī)抽取5個(gè)二級單元(調(diào)查點(diǎn))。第三階段從被抽中的5個(gè)調(diào)查點(diǎn)中共隨機(jī)抽取2062人。

    對敏感問題分層三階段抽樣的調(diào)查方法及其統(tǒng)計(jì)公式,各進(jìn)行100個(gè)樣本(每個(gè)樣本平均包含2533名模擬調(diào)查對象)的模擬抽樣調(diào)查,對模擬總體的抽樣由另外輸入的隨機(jī)數(shù)控制,根據(jù)這個(gè)隨機(jī)數(shù)計(jì)算機(jī)在模擬總體中進(jìn)行了分層三階段抽樣,并對抽中的三級單位進(jìn)行模擬調(diào)查。RRT調(diào)查過程的模擬再次用到了計(jì)算機(jī)產(chǎn)生確定范圍的隨機(jī)數(shù)的功能。例如,將北京市MSM人群按年齡分為15~29歲和30~49歲兩層,各自所占比例分別為W1=58.24%,W2=41.76%,根據(jù)模擬總體的比例,對虛擬調(diào)查對象進(jìn)行賦值,對于一個(gè)已經(jīng)被賦值為1的虛擬被調(diào)查者,計(jì)算機(jī)按一定比例單獨(dú)為它賦予一個(gè)0~3范圍的隨機(jī)整數(shù)(隨機(jī)化裝置),這相當(dāng)于此虛擬被調(diào)查者有放回地從袋中隨機(jī)摸出一個(gè)0~3范圍的任意標(biāo)號的小球。若計(jì)算機(jī)單獨(dú)為它賦予的隨機(jī)數(shù)是0,則真實(shí)回答自己屬于敏感問題的那一類的序號,即是1;若被賦予的值是0以外的某個(gè)數(shù)則回答這個(gè)數(shù)。所得結(jié)果即為該虛擬被調(diào)查者的應(yīng)答值。每一個(gè)抽中的單位都產(chǎn)生了應(yīng)答值后,按照有關(guān)公式進(jìn)行統(tǒng)計(jì)計(jì)算,得到模擬調(diào)查結(jié)果計(jì)算樣本統(tǒng)計(jì)量及其估計(jì)方差、總體參數(shù)95%的置信區(qū)間。

    3.信度和效度評價(jià)

    對于本領(lǐng)域內(nèi)相關(guān)的一些研究結(jié)果進(jìn)行了相關(guān)的比較,商業(yè)性男男性行為比例為19.66%,略高于陳向宇[17]的19.11%;男男性行為的方式為肛交、口交的比例分別為65.83%、18.70%,高于杜俏俏[18]得出的肛交占60.95%、口交占16.85%;男男性行為時(shí)沒有使用安全套的比例為24.25%,高于于波[19]的沒有使用安全套的比例19.77%,并且最近一年HIV病毒檢測結(jié)果為陽性為6.31%,低于于波的最近一年HIV檢測的結(jié)果為陽性的比例8.03%;同時(shí),本文對于每個(gè)指標(biāo)都進(jìn)行了信度和效度評價(jià),分別對各調(diào)查指標(biāo)各類別100個(gè)模擬樣本,得到的100個(gè)總體比例95%可信區(qū)間,幾乎全部包含其模擬總體比例(真值),更加說明了測量結(jié)果的準(zhǔn)確性和可靠性。例如,對指標(biāo)1分類一(肛交),100個(gè)模擬樣本中有96個(gè)(樣本1、65、95、100除外)得到的該分類總體比例95%可信區(qū)間包含其模擬總體比例(真值),見圖1對指標(biāo)1分類二(口交),100個(gè)模擬樣本中有99個(gè)樣本(樣本95除外)得到的該分類總體比例95%可信區(qū)間包含其模擬總體比例(真值);對指標(biāo)1分類三(其他方式),100個(gè)模擬樣本中有97個(gè)(樣本1、4、95除外)得到的該分類總體比例95%可信區(qū)間包含其模擬總體比例(真值)。說明本文給出的多分類敏感問題隨機(jī)應(yīng)答技術(shù)模型下分層三階段抽樣的調(diào)查方法及其統(tǒng)計(jì)公式具有較高的效度,同時(shí)具有較高的信度(因幾乎所有的樣本均數(shù)均接近它們的均值)。

    對調(diào)查指標(biāo)1的分類一,以樣本序號為橫坐標(biāo),以總體比例95%可信區(qū)間為縱坐標(biāo)繪圖,并以模擬總體比例做水平輔助線,見圖1。

    討 論

    1.本團(tuán)隊(duì)對多分類敏感問題RRT模型整群抽樣、二階段抽樣、分層兩階段整群抽樣等復(fù)雜抽樣的統(tǒng)計(jì)方法進(jìn)行過研究。在此研究基礎(chǔ)上,本文對多分類敏感問題RRT模型分層三階段抽樣,從數(shù)學(xué)上首次推導(dǎo)出各類別總體比例的估計(jì)量及其方差、估計(jì)方差的統(tǒng)計(jì)公式,具有統(tǒng)計(jì)理論方法學(xué)意義。

    圖1 多分類敏感問題RRT分層三階段抽樣模擬

    2.本團(tuán)隊(duì)研究之前,國內(nèi)外對敏感問題抽樣調(diào)查的信度與效度評價(jià)極少涉及。本文首次對多分類敏感問題RRT模型分層三階段抽樣的調(diào)查方法及其統(tǒng)計(jì)公式,采用SAS編程建立仿真模擬總體,進(jìn)行計(jì)算機(jī)大樣本個(gè)數(shù)(100個(gè))大樣本量(每個(gè)樣本含有2000多例)模擬抽樣,獲得高信度高效度的評價(jià)結(jié)論。說明本文研究的調(diào)查方法及其統(tǒng)計(jì)公式準(zhǔn)確可靠。

    3.本文研究得北京市MSM人群:男男商業(yè)性性行為的比例高達(dá)19.66%,每次全程使用安全套的比例僅為51.81%,HIV陽性比例高達(dá)6.31%,有性病比例高達(dá)16.79%,肛交比例高達(dá)65.83%,未使用安全套的比例高達(dá)24.35%,安全套使用中有破損的比例高達(dá)7.25%。本文為有關(guān)部門制定性病、艾滋病預(yù)防控制規(guī)劃、措施提供了科學(xué)依據(jù)。說明本文研究的調(diào)查方法及其統(tǒng)計(jì)公式具有較大的實(shí)際意義。

    4.評價(jià)結(jié)果表明,通過多分類敏感問題RRT模型分層三階段抽樣敏感問題調(diào)查方法能獲得大量真實(shí)數(shù)據(jù),將為有關(guān)部門制訂HIV/AIDS防治等社會(huì)、公共衛(wèi)生問題的政策、規(guī)劃提供科學(xué)的依據(jù)。

    附 錄

    模擬分層三階段抽樣并運(yùn)用多項(xiàng)選擇敏感問題隨機(jī)應(yīng)答模型進(jìn)行調(diào)查的宏程序代碼如下:

    填入相應(yīng)的宏參數(shù),并執(zhí)行編寫的宏程序,即可獲得100次運(yùn)用多項(xiàng)選擇敏感問題隨機(jī)應(yīng)答模型對模擬抽樣的樣本進(jìn)行RRT調(diào)查所得的(每類別)100個(gè)樣本比例及其估計(jì)方差、(每類別)100個(gè)總體比例的95%可信區(qū)間,100個(gè)卡方檢驗(yàn)的P值。宏參數(shù)如下:

    [1]高歌,金丕煥,王旭輝.分層三階段抽樣樣本大小的研究及應(yīng)用.中國衛(wèi)生統(tǒng)計(jì),2000,17(6):325-327.

    [2]科克倫W.G.著,張堯庭,吳輝譯.抽樣技術(shù).中國統(tǒng)計(jì)出版社,1985:87.

    [3]溫長松.試述抽樣調(diào)查方法在歷史研究中的應(yīng)用.沈陽大學(xué)學(xué)報(bào),2006,18(1):40-41.

    [4]王建華.實(shí)用醫(yī)學(xué)科研方法.人民衛(wèi)生出版社,2003:440-450.

    [5]高歌,范玉波.敏感問題Simmons模型的(分層)整群抽樣研究.中國衛(wèi)生統(tǒng)計(jì),2008,25(6):562-565.

    [6]Warner SL.Randomized response:a survey technique for eliminating evasive answer bias.JAm Stat Assoc,1965:60(309):63-66.

    [7]Tourangeau R,Smith TW.Asking sensitive questions:The impact of date collection mode,question format,and question context.Public Opinion Quarterly,1996;(60):275-304.

    [8]Stephen EE,Samuel H,Karen LD.Validity of Forced Responses in a Random ized Response Model.Sociological Methods&Research,1982,11(1):89-100.

    [9]Liu PT,Chow LP.The efficiency of themultiple trial random ized response technique.Biometrics,1976,32(3):607-618.

    [10]Chow LP,Gruhn W,Chang WP.Feasibility of the random ized response technique in rural Ethiopia.Am J Public Health,1979,69(3):273-276.

    [11]孫山澤,孫明舉,段鋼.二項(xiàng)選擇敏感性問題調(diào)查的基本方法.?dāng)?shù)理統(tǒng)計(jì)與管理,2000,19(1):58-64.

    [12]洪志敏,閆在在.一種定量敏感性問題的隨機(jī)化調(diào)查方法.內(nèi)蒙古工業(yè)大學(xué)學(xué)報(bào),2006,25(3):178-182.

    [13]張權(quán),閆在在.二項(xiàng)抽樣下隨機(jī)化調(diào)查的比估計(jì)模型.內(nèi)蒙古工業(yè)大學(xué)學(xué)報(bào),2009,28(2):81-85.

    [14]Wang M,Gao G.Quantitative sensitive question survey in cluster sampling and its application.Recent Advance in Statistics Application and Related Areas,Sydeny:Aussino Academ ic Publishing House 2008:648-652.

    [15]Liu W,Gao G,Wang L.Stratified random sampling on simmons model for sensitive question survey.Data Processing and Quantitative Economy Modeling,Sydeny:Aussino Academ ic Publishing House,2010:22-26.

    [16]陳向宇,高歌,于明潤,等.北京市商業(yè)性男男性行為Simmons模型調(diào)查.中國公共衛(wèi)生,2013,2:259-260.

    [17]杜俏俏,高歌,靳宗達(dá),等.蒙特卡洛模擬在多分類敏感問題兩階段整群抽樣信度與效度評價(jià)中的應(yīng)用.中國衛(wèi)生統(tǒng)計(jì),2013,30(2):227-229.

    [18]于波,高歌,賀志龍,等.北京市MSM人群RRT模型的兩階段抽樣調(diào)查.中國衛(wèi)生統(tǒng)計(jì),2012;06:351-352.

    (責(zé)任編輯:劉 壯)

    *:國家自然科學(xué)基金項(xiàng)目(項(xiàng)目編號:81273188)

    △通信作者:高歌,E-mail:gaoge@suda.edu.cn

    猜你喜歡
    安全套總體比例
    人體比例知多少
    用樣本估計(jì)總體復(fù)習(xí)點(diǎn)撥
    2020年秋糧收購總體進(jìn)度快于上年
    外匯市場運(yùn)行有望延續(xù)總體平穩(wěn)發(fā)展趨勢
    中國外匯(2019年6期)2019-07-13 05:44:06
    新型女用安全套的使用及綜合咨詢
    直擊高考中的用樣本估計(jì)總體
    按事故責(zé)任比例賠付
    紅土地(2016年7期)2016-02-27 15:05:54
    使用安全套,別入這些誤區(qū)
    限制支付比例只是治標(biāo)
    安全使用“安全套”
    男人舔奶头视频| 精品国产三级普通话版| 99热精品在线国产| 18禁黄网站禁片免费观看直播| 色综合站精品国产| 免费观看精品视频网站| or卡值多少钱| 欧美不卡视频在线免费观看| 舔av片在线| 日韩在线高清观看一区二区三区| 成人特级黄色片久久久久久久| 欧美色视频一区免费| 1000部很黄的大片| 六月丁香七月| 免费看光身美女| 亚洲精品亚洲一区二区| 三级毛片av免费| 久久久久久大精品| av中文乱码字幕在线| 人妻制服诱惑在线中文字幕| 亚洲一区高清亚洲精品| av视频在线观看入口| 成人国产麻豆网| 精品福利观看| 黄色欧美视频在线观看| 国产片特级美女逼逼视频| 十八禁网站免费在线| 大香蕉久久网| 亚洲人成网站高清观看| 午夜激情福利司机影院| 观看美女的网站| 亚洲综合色惰| 亚洲电影在线观看av| 国产av麻豆久久久久久久| 国产麻豆成人av免费视频| 嫩草影院精品99| 免费电影在线观看免费观看| 变态另类成人亚洲欧美熟女| 午夜久久久久精精品| 一本久久中文字幕| 亚州av有码| 日本黄色视频三级网站网址| 成人美女网站在线观看视频| 亚洲最大成人手机在线| 国产蜜桃级精品一区二区三区| 成年女人看的毛片在线观看| 青春草视频在线免费观看| 精品一区二区三区视频在线| 国产黄a三级三级三级人| 日本免费一区二区三区高清不卡| 啦啦啦韩国在线观看视频| 能在线免费观看的黄片| 我要搜黄色片| 国产精品久久久久久久久免| 亚洲av.av天堂| 久久精品夜夜夜夜夜久久蜜豆| 五月伊人婷婷丁香| 99久久精品热视频| 精品乱码久久久久久99久播| 看非洲黑人一级黄片| 一进一出抽搐gif免费好疼| 黄色视频,在线免费观看| 精品人妻偷拍中文字幕| 免费人成在线观看视频色| 国产精品嫩草影院av在线观看| 国产精品,欧美在线| 悠悠久久av| 成人一区二区视频在线观看| 别揉我奶头~嗯~啊~动态视频| 日韩av在线大香蕉| 听说在线观看完整版免费高清| 日本 av在线| 国产熟女欧美一区二区| 大又大粗又爽又黄少妇毛片口| 乱系列少妇在线播放| 观看美女的网站| 草草在线视频免费看| 亚洲最大成人中文| 亚洲天堂国产精品一区在线| 亚洲va在线va天堂va国产| 免费高清视频大片| 亚洲美女黄片视频| 婷婷精品国产亚洲av| 搞女人的毛片| 国产综合懂色| 国产日本99.免费观看| 插阴视频在线观看视频| 亚洲精品色激情综合| 亚洲,欧美,日韩| 色视频www国产| 欧美成人a在线观看| 国产欧美日韩精品亚洲av| 亚洲精品在线观看二区| 国产精品亚洲一级av第二区| 特级一级黄色大片| 人妻丰满熟妇av一区二区三区| 性欧美人与动物交配| 1024手机看黄色片| 午夜免费激情av| 九九热线精品视视频播放| 免费av毛片视频| 亚洲精品日韩av片在线观看| 午夜福利在线观看免费完整高清在 | 国产白丝娇喘喷水9色精品| 亚洲在线自拍视频| 欧美另类亚洲清纯唯美| 精品一区二区三区视频在线| 亚洲婷婷狠狠爱综合网| 一夜夜www| 露出奶头的视频| 99热只有精品国产| 久久精品久久久久久噜噜老黄 | 日韩欧美免费精品| 国产男人的电影天堂91| 少妇被粗大猛烈的视频| 午夜福利18| 12—13女人毛片做爰片一| 在线观看免费视频日本深夜| 国内精品美女久久久久久| 日日摸夜夜添夜夜添av毛片| 观看免费一级毛片| 久久精品国产鲁丝片午夜精品| 久久精品国产99精品国产亚洲性色| 亚洲欧美日韩高清专用| 99久久无色码亚洲精品果冻| 久久久精品94久久精品| 久久精品影院6| 国产精品,欧美在线| 精品乱码久久久久久99久播| 国产精品一及| eeuss影院久久| av国产免费在线观看| 一区福利在线观看| 香蕉av资源在线| 精品午夜福利在线看| a级毛片免费高清观看在线播放| 偷拍熟女少妇极品色| а√天堂www在线а√下载| 国产淫片久久久久久久久| 成年女人毛片免费观看观看9| 精品久久久久久成人av| 日韩大尺度精品在线看网址| 久久久欧美国产精品| 国产毛片a区久久久久| 亚洲国产高清在线一区二区三| 大型黄色视频在线免费观看| 色噜噜av男人的天堂激情| 丰满人妻一区二区三区视频av| 日本精品一区二区三区蜜桃| 五月玫瑰六月丁香| 欧美一区二区国产精品久久精品| 久久久精品94久久精品| 超碰av人人做人人爽久久| 美女内射精品一级片tv| 精品欧美国产一区二区三| 国产片特级美女逼逼视频| 性欧美人与动物交配| 精品久久久久久久久久久久久| 精品久久久噜噜| 亚洲精华国产精华液的使用体验 | 成人美女网站在线观看视频| av专区在线播放| 中文资源天堂在线| 国产毛片a区久久久久| 亚洲第一区二区三区不卡| 大型黄色视频在线免费观看| av天堂在线播放| 国产不卡一卡二| 日韩欧美三级三区| 国产精品三级大全| 国产黄色小视频在线观看| 亚洲五月天丁香| 午夜爱爱视频在线播放| 99热全是精品| 特大巨黑吊av在线直播| 久久久久精品国产欧美久久久| 日本一本二区三区精品| 国产伦精品一区二区三区四那| 国产伦一二天堂av在线观看| 亚洲精品一卡2卡三卡4卡5卡| 国产白丝娇喘喷水9色精品| 夜夜看夜夜爽夜夜摸| 男女啪啪激烈高潮av片| 精华霜和精华液先用哪个| 午夜a级毛片| 久久人人爽人人爽人人片va| 亚洲性夜色夜夜综合| 俺也久久电影网| 免费在线观看成人毛片| 九九爱精品视频在线观看| 国产v大片淫在线免费观看| 国产高清三级在线| 直男gayav资源| 国产av麻豆久久久久久久| 日日摸夜夜添夜夜添av毛片| 国产老妇女一区| av在线观看视频网站免费| 99视频精品全部免费 在线| 久久欧美精品欧美久久欧美| 一a级毛片在线观看| 男女边吃奶边做爰视频| 91久久精品国产一区二区成人| 久99久视频精品免费| 国产久久久一区二区三区| 免费观看精品视频网站| 99久久九九国产精品国产免费| 又粗又爽又猛毛片免费看| 精华霜和精华液先用哪个| 变态另类丝袜制服| 嫩草影院入口| 美女 人体艺术 gogo| 国产乱人偷精品视频| 国产精品亚洲美女久久久| 寂寞人妻少妇视频99o| 日韩欧美精品v在线| 国产亚洲91精品色在线| 亚洲欧美日韩卡通动漫| 精品一区二区三区视频在线| 成熟少妇高潮喷水视频| 精品久久久久久久人妻蜜臀av| 久久久精品大字幕| 久久精品国产亚洲av涩爱 | 成年版毛片免费区| 在线观看av片永久免费下载| 99久久精品热视频| 99九九线精品视频在线观看视频| 国产真实乱freesex| 亚洲图色成人| 国产精品99久久久久久久久| 色哟哟哟哟哟哟| 亚洲欧美成人综合另类久久久 | 日韩欧美在线乱码| 欧美三级亚洲精品| 国内揄拍国产精品人妻在线| a级毛片a级免费在线| 男人和女人高潮做爰伦理| 亚洲人成网站在线播| 色5月婷婷丁香| 1000部很黄的大片| 99热只有精品国产| 久久九九热精品免费| 欧美成人精品欧美一级黄| 99久久精品一区二区三区| 小说图片视频综合网站| 波多野结衣高清无吗| 日韩欧美 国产精品| 国产精品电影一区二区三区| 国产91av在线免费观看| 美女内射精品一级片tv| 神马国产精品三级电影在线观看| 午夜视频国产福利| 亚洲性久久影院| h日本视频在线播放| 国产乱人视频| 免费电影在线观看免费观看| 欧美成人a在线观看| 麻豆av噜噜一区二区三区| 午夜精品一区二区三区免费看| 在线观看66精品国产| 人人妻人人澡欧美一区二区| 日本精品一区二区三区蜜桃| 夜夜爽天天搞| 成人鲁丝片一二三区免费| 真人做人爱边吃奶动态| 高清午夜精品一区二区三区 | 寂寞人妻少妇视频99o| 日本 av在线| 校园人妻丝袜中文字幕| 一本一本综合久久| 直男gayav资源| 熟妇人妻久久中文字幕3abv| 自拍偷自拍亚洲精品老妇| 国产精品一区二区三区四区久久| 91在线精品国自产拍蜜月| 国产真实乱freesex| 日本一本二区三区精品| 老熟妇乱子伦视频在线观看| 无遮挡黄片免费观看| 亚洲最大成人av| 亚洲欧美中文字幕日韩二区| 亚洲人成网站高清观看| 伦理电影大哥的女人| 亚洲综合色惰| 国产欧美日韩一区二区精品| 国产精品人妻久久久久久| 国产高清三级在线| .国产精品久久| 九色成人免费人妻av| 精品99又大又爽又粗少妇毛片| 12—13女人毛片做爰片一| 成人午夜高清在线视频| 黄色欧美视频在线观看| av国产免费在线观看| 在线免费十八禁| 91在线观看av| 菩萨蛮人人尽说江南好唐韦庄 | 在线国产一区二区在线| 精品一区二区免费观看| 人妻少妇偷人精品九色| 天天一区二区日本电影三级| 日本精品一区二区三区蜜桃| 最新在线观看一区二区三区| 欧美区成人在线视频| 露出奶头的视频| 特大巨黑吊av在线直播| 亚洲av五月六月丁香网| 国产av麻豆久久久久久久| 午夜视频国产福利| 国产女主播在线喷水免费视频网站 | 久99久视频精品免费| 国产 一区 欧美 日韩| 日本成人三级电影网站| 一a级毛片在线观看| 天美传媒精品一区二区| 久久久国产成人精品二区| 熟妇人妻久久中文字幕3abv| 欧美日本亚洲视频在线播放| 噜噜噜噜噜久久久久久91| 国产麻豆成人av免费视频| 啦啦啦韩国在线观看视频| 你懂的网址亚洲精品在线观看 | 国产精品久久电影中文字幕| 伊人久久精品亚洲午夜| 91av网一区二区| 嫩草影院入口| 亚洲国产精品国产精品| 亚洲人成网站在线播放欧美日韩| 男人舔女人下体高潮全视频| 狂野欧美激情性xxxx在线观看| 韩国av在线不卡| 国产亚洲精品av在线| 日韩欧美 国产精品| 亚洲真实伦在线观看| 俄罗斯特黄特色一大片| 色视频www国产| 亚洲最大成人中文| av国产免费在线观看| 99国产精品一区二区蜜桃av| 一级av片app| 亚洲最大成人中文| 男女做爰动态图高潮gif福利片| 99久久无色码亚洲精品果冻| 日本三级黄在线观看| 亚洲专区国产一区二区| 免费电影在线观看免费观看| 免费观看的影片在线观看| 精品久久久噜噜| 亚洲激情五月婷婷啪啪| 18禁在线播放成人免费| 日本黄色片子视频| 97超碰精品成人国产| 亚洲国产高清在线一区二区三| 国产一区亚洲一区在线观看| 久久久色成人| 少妇的逼水好多| 最好的美女福利视频网| 欧美性猛交╳xxx乱大交人| 尤物成人国产欧美一区二区三区| 成人av一区二区三区在线看| 哪里可以看免费的av片| 亚洲国产精品成人久久小说 | 大又大粗又爽又黄少妇毛片口| 国产精品一区二区性色av| 亚洲熟妇熟女久久| 亚洲人与动物交配视频| 在线国产一区二区在线| 婷婷精品国产亚洲av在线| 乱人视频在线观看| 国产91av在线免费观看| av天堂中文字幕网| 亚洲18禁久久av| 热99re8久久精品国产| 在线播放无遮挡| 日韩欧美精品v在线| 一夜夜www| 乱码一卡2卡4卡精品| 国产极品精品免费视频能看的| 色哟哟·www| 美女高潮的动态| 亚洲色图av天堂| 日本精品一区二区三区蜜桃| 国产麻豆成人av免费视频| 国产伦在线观看视频一区| 97人妻精品一区二区三区麻豆| 一边摸一边抽搐一进一小说| 国产精品一区二区免费欧美| 欧美日韩一区二区视频在线观看视频在线 | 亚洲av五月六月丁香网| 色在线成人网| 亚洲五月天丁香| 99热精品在线国产| 亚洲成人精品中文字幕电影| 69人妻影院| 国产69精品久久久久777片| 两性午夜刺激爽爽歪歪视频在线观看| 黄片wwwwww| 变态另类丝袜制服| 欧美色欧美亚洲另类二区| 亚洲欧美日韩东京热| 久久精品国产亚洲av香蕉五月| 国产精华一区二区三区| 成人欧美大片| 欧美日韩一区二区视频在线观看视频在线 | 欧美性感艳星| 一边摸一边抽搐一进一小说| 99热这里只有是精品在线观看| 一本精品99久久精品77| 国产精品亚洲美女久久久| 国产视频一区二区在线看| 精品久久久久久久久亚洲| 听说在线观看完整版免费高清| 国产成人影院久久av| 成熟少妇高潮喷水视频| 国产av麻豆久久久久久久| 日韩精品中文字幕看吧| 日本免费一区二区三区高清不卡| 香蕉av资源在线| 日本免费a在线| 别揉我奶头~嗯~啊~动态视频| 精品人妻一区二区三区麻豆 | 97超碰精品成人国产| 在线观看66精品国产| 国产真实伦视频高清在线观看| 国国产精品蜜臀av免费| 午夜福利视频1000在线观看| 久久草成人影院| 亚洲精华国产精华液的使用体验 | 国产午夜精品久久久久久一区二区三区 | 国产高清视频在线观看网站| 国产一区二区三区av在线 | 又黄又爽又免费观看的视频| 国产一区亚洲一区在线观看| 日本欧美国产在线视频| 欧美+日韩+精品| 高清日韩中文字幕在线| 51国产日韩欧美| .国产精品久久| 国国产精品蜜臀av免费| 中文字幕免费在线视频6| 国内精品美女久久久久久| 激情 狠狠 欧美| 国产精品久久久久久久久免| 人妻夜夜爽99麻豆av| 欧美日韩乱码在线| 国产aⅴ精品一区二区三区波| 一级毛片久久久久久久久女| 99久久成人亚洲精品观看| 一级黄色大片毛片| 国产成人福利小说| 校园春色视频在线观看| 国产精品亚洲一级av第二区| 免费大片18禁| 久久久精品大字幕| www日本黄色视频网| 国产片特级美女逼逼视频| 国产精品乱码一区二三区的特点| 黑人高潮一二区| 欧美3d第一页| 三级国产精品欧美在线观看| 69人妻影院| 美女被艹到高潮喷水动态| 九九久久精品国产亚洲av麻豆| 国产精品美女特级片免费视频播放器| 国产精品乱码一区二三区的特点| 久久人人爽人人爽人人片va| 成人亚洲精品av一区二区| 大型黄色视频在线免费观看| 一个人观看的视频www高清免费观看| 精品久久久久久久人妻蜜臀av| 日本 av在线| 啦啦啦啦在线视频资源| 亚洲成人久久性| 亚洲精品影视一区二区三区av| 日韩在线高清观看一区二区三区| 免费av毛片视频| 美女高潮的动态| 免费在线观看影片大全网站| 亚洲七黄色美女视频| 亚洲av中文字字幕乱码综合| 国产不卡一卡二| 秋霞在线观看毛片| 国产老妇女一区| 又黄又爽又刺激的免费视频.| 亚洲无线观看免费| 国产午夜精品久久久久久一区二区三区 | 精品欧美国产一区二区三| 观看免费一级毛片| 欧美三级亚洲精品| 国产精品免费一区二区三区在线| 在线观看午夜福利视频| 亚洲精品乱码久久久v下载方式| 三级经典国产精品| 国产伦精品一区二区三区视频9| 午夜福利在线观看吧| 丝袜喷水一区| 日韩一本色道免费dvd| 舔av片在线| 国产亚洲欧美98| 又黄又爽又免费观看的视频| 国产高清激情床上av| 亚洲精品国产av成人精品 | 寂寞人妻少妇视频99o| 亚洲av成人av| 欧美xxxx性猛交bbbb| 国产午夜精品论理片| 国内揄拍国产精品人妻在线| 亚洲av免费高清在线观看| 国产精品不卡视频一区二区| 3wmmmm亚洲av在线观看| 久久韩国三级中文字幕| 久久久久免费精品人妻一区二区| 天堂动漫精品| 晚上一个人看的免费电影| 国产三级在线视频| 亚洲欧美中文字幕日韩二区| 九九在线视频观看精品| 亚洲av第一区精品v没综合| 午夜日韩欧美国产| 麻豆国产97在线/欧美| 乱码一卡2卡4卡精品| 欧美日韩精品成人综合77777| 国产日本99.免费观看| 成人毛片a级毛片在线播放| 欧美在线一区亚洲| 国内精品一区二区在线观看| 我要搜黄色片| 赤兔流量卡办理| 免费高清视频大片| 久久久久久久亚洲中文字幕| 色综合亚洲欧美另类图片| 精品午夜福利视频在线观看一区| 日韩精品青青久久久久久| 一进一出抽搐动态| 全区人妻精品视频| 国产精品久久久久久av不卡| 内射极品少妇av片p| 中国国产av一级| 亚洲欧美成人精品一区二区| 久久久久久国产a免费观看| 可以在线观看毛片的网站| 天堂动漫精品| 亚洲人成网站高清观看| 精品乱码久久久久久99久播| 国产91av在线免费观看| 午夜影院日韩av| 精品无人区乱码1区二区| 国产女主播在线喷水免费视频网站 | 无遮挡黄片免费观看| 免费在线观看成人毛片| 我的老师免费观看完整版| 禁无遮挡网站| 夜夜看夜夜爽夜夜摸| 久久国内精品自在自线图片| 久久精品91蜜桃| 日本一二三区视频观看| 成人亚洲欧美一区二区av| 亚洲精品久久国产高清桃花| 看免费成人av毛片| 女人十人毛片免费观看3o分钟| 人妻制服诱惑在线中文字幕| 超碰av人人做人人爽久久| 久久久久久久久久成人| 超碰av人人做人人爽久久| 久久久久久久久久成人| 69人妻影院| 久久鲁丝午夜福利片| 日本a在线网址| 不卡视频在线观看欧美| 最近中文字幕高清免费大全6| 黄色配什么色好看| 日本爱情动作片www.在线观看 | 国产一区二区亚洲精品在线观看| 精品久久久久久久久久免费视频| 深夜精品福利| 精品久久久噜噜| 国产精品日韩av在线免费观看| 国产欧美日韩精品一区二区| 久久久精品94久久精品| 免费观看精品视频网站| 久久久精品94久久精品| 国产亚洲精品av在线| 嫩草影院精品99| 成人永久免费在线观看视频| 久久久精品欧美日韩精品| 亚洲成av人片在线播放无| 99九九线精品视频在线观看视频| 国产精品99久久久久久久久| 69人妻影院| 一个人看的www免费观看视频| 卡戴珊不雅视频在线播放| 在线a可以看的网站| 禁无遮挡网站| 日韩av不卡免费在线播放| 中文字幕精品亚洲无线码一区| 日本免费一区二区三区高清不卡| 精品人妻熟女av久视频| 亚洲精品在线观看二区| 国产视频内射| 精华霜和精华液先用哪个| 亚洲不卡免费看| 91麻豆精品激情在线观看国产| 久久久久精品国产欧美久久久| 亚洲真实伦在线观看| 三级毛片av免费| 村上凉子中文字幕在线| 午夜激情欧美在线| 免费看a级黄色片| 久久草成人影院| 此物有八面人人有两片| 18禁裸乳无遮挡免费网站照片| 你懂的网址亚洲精品在线观看 | 久久久a久久爽久久v久久| 长腿黑丝高跟| 免费av观看视频| 一区二区三区四区激情视频 |