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

    基于傾向值匹配研究低溫對流行性感冒的影響

    2018-04-04 08:55:23,,,,,,
    中國人獸共患病學報 2018年2期
    關(guān)鍵詞:流感氣壓重慶市

     ,,, ,, , 

    流行性感冒是由流感病毒引起的一種世界性流行的急性呼吸道傳染病,近年來,其發(fā)生頻率不斷增加、流行范圍也不斷擴大[1]。我國2009年暴發(fā)的甲型H1N1流感與2013年的人感染H7N9流感,給個人健康造成巨大威脅,同時對社會亦造成了嚴重影響,引起全球廣泛關(guān)注[2-3]。

    流感的發(fā)生除受流感病毒致病性、人群易感性以及社會衛(wèi)生條件等因素影響外,也可能受到氣候因素的影響[4]。有研究表示,大多數(shù)流感病例是在寒冷季節(jié)和雨季的最低溫度下監(jiān)測到的[5]。那么在其它氣候條件相同的情況下,低溫在流感的傳播中究竟扮演了一個怎樣的角色,目前就低溫對流感的影響作用已有相關(guān)的文獻研究,但各文獻在研究低溫與流感關(guān)系時往往未平衡其它氣候因素,導致分析不到位。同時,由于全球氣候變暖,明確溫度與流感的關(guān)系、探討溫度變化對流感的影響更是成為全球公共衛(wèi)生研究的熱點之一[6]。為更好地了解低溫因素對流感發(fā)病的影響、把握流感發(fā)病規(guī)律、促使流感防控工作更高效地進行,本研究采用傾向值評分模型對重慶市2010年-2015年流感報告病例進行統(tǒng)計分析,平衡日最低氣溫以外的其它氣候因素的混雜作用,探討低溫對日流感發(fā)病率的獨立影響,對流感的預測和預警以及防控對策的制定提供依據(jù)。

    1 材料與方法

    1.1材料2010年1月-2015年12月重慶市日流感發(fā)病疫情數(shù)據(jù)來源于“傳染病報告信息管理系統(tǒng)”,重慶市總?cè)丝跀?shù)來自2010-2015年重慶市衛(wèi)生統(tǒng)計年鑒,日流感發(fā)病率(%)=日流感發(fā)病人數(shù)/重慶市總?cè)丝?萬人)。同期氣象資料來自重慶市氣象局,包括(日)平均露點、平均氣壓、平均水氣壓、最高氣溫、最低氣溫、平均氣溫、平均相對濕度、平均10 min風速、日照時數(shù)等氣象指標。

    1.2材料預處理本研究通過傾向值匹配,去除(日)平均露點、平均氣壓、平均水氣壓、最高氣溫、最低氣溫、平均氣溫、平均相對濕度、平均10 min風速、日照時數(shù)等混雜因素的影響,僅分析日流感發(fā)病率與日最低氣溫之間的關(guān)系。

    傾向值是在給定觀測協(xié)變量向量xi的情況下,成員分配至某一特定干預(Wi=1)而不是非干預(W0=0)的條件概率為:e(xi)=pr(Wi=1|Xi=xi)。傾向值在匹配中的優(yōu)點在于它將所有的這些維度簡化為一個單維的值(即ps值)。由于本研究樣本量足夠大,為獲得高質(zhì)量的匹配,筆者采用1∶1最近鄰匹配法來控制兩組基線資料偏倚。首先通過一般化加速建模(GBM)獲得傾向值,計算每d的傾向性評分分值,其優(yōu)點是可在預測誤差實現(xiàn)最小化時停止迭代;同時為保證每對樣本的均衡性,采用卡鉗匹配,并限定傾向得分的對數(shù)標準偏差為0.1,獲得具有相似傾向值的兩組成員,公式為:

    Pi和Pj分別是干預組(高溫組)和控制組(低溫組)的傾向值,I1和I0分別是干預組和控制組的集合。當傾向值之差的絕對值在和之間傾向值的所有可能配對中最小時,鄰居關(guān)系C(Pi)包含一個控制組參者j(即j∈I0),其作為干預組的匹配。一旦找到一個與i匹配的j,j就從I0中移出且不再放回。對于每個i,若只找到單個的j落入C(Pi),則此匹配為最近鄰成對匹配。

    考慮到日最低溫度為非正態(tài)分布(Kolmogorov-Smirnov: Z=0.056,P<0.001),故本研究以六年中位數(shù)(16.98 ℃)將氣象資料分為高溫與低溫兩組。經(jīng)過匹配,每一個控制組(低溫組)都能在干預組(高溫組)找到一個與之匹配的對照個體。

    1.3負二項回歸流感屬于傳染性疾病,日流感發(fā)病率之間不滿足獨立要求,故不能用一般線性回歸或Possion回歸進行分析。筆者采用負二項回歸分析低溫對流感發(fā)病率的影響。公式為:

    式中回歸系數(shù)βm表示在控制其他自變量的情況下xm對事件發(fā)生強度的影響大小。模型誤差分布指定負二項分布(NB),連接函數(shù)用對數(shù)連接;模型的擬優(yōu)合度采用χ2檢驗(P<0.05為模型有意義)、偏差(越小越好)、調(diào)整偏差(越小越好)、AIC(越小越好)、AICC(越小越好)評價。

    1.4統(tǒng)計學處理應用Excel 2007進行數(shù)據(jù)整理,Stata 14.0軟件進行統(tǒng)計分析。因氣象資料均為非正態(tài)分布,故以中位數(shù)與四分位數(shù)間距進行描述性分析;spearman等級相關(guān)分析重慶市日最低氣溫與流感日發(fā)病人數(shù)的相關(guān)性;采用傾向值匹配以控制除日最低氣溫以外的氣象混雜因素;匹配前后的組間比較采用秩和檢驗;采用負二項回歸評估溫度對流感日發(fā)病率的影響程度。

    2 結(jié) 果

    2.1重慶市流感發(fā)病總體情況2010年1月-2015年12月重慶市共報告流感15 615例,日均發(fā)病數(shù)為7.137.23例,年發(fā)病率在3.36/10萬~11.19/10萬之間。經(jīng)Cochran-Armitage趨勢檢驗發(fā)現(xiàn)2010-2015年流感年發(fā)病率差異有統(tǒng)計學意義(χ2=38.61,P<0.001),即重慶市流感流行強度總體有逐年上升趨勢。

    2.2低溫與流感相關(guān)性研究2010年-2015年重慶市平均周最低溫度與流感發(fā)病數(shù)趨勢見圖1,可見流感周報告發(fā)病人數(shù)整體上呈“U”型分布,而溫度呈倒“U”型分布??梢姕囟仍礁?,流感報告人數(shù)越少;反之,溫度下降,流感報告人數(shù)有上升趨勢。最低氣溫與流感報告發(fā)病人數(shù)呈負相關(guān)(spearman等級相關(guān):P<0.001),表示低溫可能是導致流感發(fā)病人數(shù)增加的一個重要氣象因子。

    圖1 重慶市溫度與流感發(fā)病人數(shù)Fig.1 Number of temperature and flu in Chongqing

    2.3傾向性評分匹配為排除其他混雜氣象因素(平均露點、平均氣壓、平均水氣壓、平均相對濕度、平均10 min風速、日照時數(shù)等)對日流感發(fā)病率的影響,本文對2010-2015年流感資料進行1∶1最近鄰匹配。

    1) 估計傾向值(一般化加速建模):選取平均

    露點、平均氣壓、平均水氣壓、平均相對濕度、平均10 min風速、日照時數(shù)等日氣象因素為自變量,以溫度分組為因變量。為優(yōu)化每一步迭代中的似然函數(shù),并減少變異性,本研究設(shè)定Boost程序最大迭代次數(shù)為1 000,以logistic函數(shù)為連接函數(shù),最大允許四次交互項,并設(shè)置收縮系數(shù)為0.000 5這一相對較小的收縮來保證得到流暢的擬合。將估計的傾向值命名為ps。

    2) 傾向值匹配:以低溫組為控制組,高溫組為干預組進行匹配,兩組共158對匹配成功,兩組基線分別進行秩和檢驗,匹配后所有變量差異均無統(tǒng)計學意義(P>0.05,見表1),兩組匹配后的傾向性評分分布也更加均衡,見圖2、3。

    圖2 基于核密度估計的高溫組與低溫組匹配前傾向值分布Fig.2 Data of propensity score before matching between high temperature group and low temperature group based on nuclear density estimation

    圖3 基于核密度估計的高溫組與低溫組匹配后傾向值分布Fig.3 Data of propensity score after matching between high temperature group and low temperature group based on nuclear density estimation

    表1匹配前后高低溫組的基線水平比較
    Tab.1Comparison on the baseline level before and after matching with the high and low temperature groups

    匹配前Beforematching匹配后Aftermatching高溫組Hightemperaturegroup(n=1101)低溫組Lowtemperaturegroup(n=1080)Z值ZvalueP值Pvalue高溫組Hightemperaturegroup(n=79)低溫組Lowtemperaturegroup(n=79)Z值ZvalueP值Pvalue平均露點溫度(℃)Averagedewpoint(℃)19.62(4.46)8.19(6.65)-39.125<0.00113.55(2.66)14.47(4.77)0.5600.576平均氣壓(hPa)Averagepressure(hPa)976.30(8.00)989.40(8.00)34.691<0.001982.90(8.40)983.60(7.70)1.4750.140平均水汽壓(hPa)Averagewatervaporpressure(hPa)22.90(6.20)10.80(4.90)-39.387<0.00115.80(2.50)16.90(4.70)1.0870.2771表1(續(xù))匹配前Beforematching匹配后Aftermatching高溫組Hightemperaturegroup(n=1101)低溫組Lowtemperaturegroup(n=1080)Z值ZvalueP值Pvalue高溫組Hightemperaturegroup(n=79)低溫組Lowtemperaturegroup(n=79)Z值ZvalueP值Pvalue相對濕度(%)Relativehumidity(%)73.00(20.00)78.00(14.00)9.378<0.00163.00(19.00)71.00(26.00)-1.2790.2010降水量(mm)Rainfall(mm)0.00(4.00)0.00(1.00)-5.096<0.0010.00(0.00)0.00(3.00)1.2090.227平均2min風速(m/s)Average2minuteswindspeed(m/s)1.40(0.50)1.30(0.40)-8.820<0.0011.30(0.00)1.20(1.00)-1.4600.144平均10min風速(m/s)Average10minuteswindspeed(m/s)1.40(0.5)1.30(0.40)-9.506<0.0011.30(0.00)1.30(1.00)-1.4480.148日照時速(h)Sunshinehour(h)2.40(8.00)0.00(0.10)-16.797<0.0012.70(8.00)0.00(7.00)-1.0410.298

    注:1.各氣象資料采用中位數(shù)(四分位數(shù)間距)進行統(tǒng)計描述;2.匹配前后均采用秩和檢驗比較兩組是否存在差異。

    注: 1.Meteorological data using the median (interquartile range) statistical description;2.Using rank sum test to compare the two groups whether there is a difference.

    2.4負二項回歸模型擬合及效果評價

    2.4.1匹配前:考慮流感潛伏期和溫度對流感發(fā)病率的影響,以平均露點溫度(x1)、平均氣壓(x2)、平均水氣壓(x3)、最低溫度(x4)、相對濕度(x5)、降水量(x6)、平均2 min風速(x7)、平均10 min風速(x8)、日照時數(shù)(x9)為自變量,分別以日流感發(fā)病率(reg1)與滯后3 d的日流感發(fā)病率(reg2)為響應變量進行負二項回歸,結(jié)果見表2。

    2.4.2匹配后:以logistic回歸得到的估計傾向值(ps)為協(xié)變量,日最低溫度(x4)為自變量,分別以日流感發(fā)病率(reg3)與滯后3 d的日流感發(fā)病率(reg4)為響應變量進行負二項回歸,結(jié)果見表2。

    由表2可見,匹配前后同一天的負二項回歸模型(reg1,reg3)整體不顯著(P>0.05),滯后3 d的負二項回歸模型(reg2,reg3)有統(tǒng)計學意義(P<0.05)。利用傾向值匹配后的滯后3 d流感發(fā)病率進行負二項回歸,其模型偏差、調(diào)整偏差、AIC、AICC均小于匹配前滯后3 d模型,這4個指標越大,說明效果擬合得越差??梢哉f明:就本資料而言,reg4模型的擬合效果要優(yōu)于reg2模型。

    表2回歸模型參數(shù)比較
    Tab.2Comparison of regression model parameters

    匹配前回歸Regressionbeforematching匹配后回歸Regressionaftermatching同1dOneday(reg1)滯后3dThreedayslag(reg2)同1dOneday(reg3)滯后3dThreedayslag(re4)偏差Diviance-308.970-18.923調(diào)整后偏差Scaleddeviance-308.970-18.923AIC-2044.934-170.400AICC-2045.080-170.662日最低溫度系數(shù)Lowesttem?peraturecoefficient--0.001-0.144檢驗Chi?squaretest7.01029.2102.1106.550P值Pvalue0.6360.0010.3490.038觀測數(shù)Observations21912188158158

    3 討 論

    流感是第一個實行全球監(jiān)測的呼吸道傳染病,一個多世紀以來,流感一直以其發(fā)病率高、超額住院等成為威脅人類健康的疾病[7]。2004年全球發(fā)生人感染H5N1高致病性禽流感疫情,2009年甲型H1N1流感大流行[8],2013年2月以來至今的中國局部地區(qū)發(fā)生的人感染H7N9疫情,均有死亡病例報道[9-10]。本文基于2010年-2015年流感年發(fā)病率數(shù)據(jù),經(jīng)Cochran-Armitage趨勢檢驗發(fā)現(xiàn),其整體發(fā)病率有逐年上升趨勢,流感疫情形勢已十分嚴峻。

    目前國內(nèi)外已有不少研究表明氣象因素與流感的流行存在密切關(guān)系[11],亦有研究表明:環(huán)境因素如太陽輻射、相對濕度、溫度等對流感病毒的傳播有重要作用[12]。其中氣溫是人們關(guān)注的、較易測量的一個直接氣象要素,也是影響傳染病發(fā)生發(fā)展的重要因素之一[13]。自1961年以來,重慶氣溫總體呈上升趨勢,其中以日最低溫度的變化趨勢最為明顯[14]。由2010年-2015年重慶市平均周最低溫度與日流感發(fā)病數(shù)趨勢圖可見:溫度越高,流感報告人數(shù)越少;反之,溫度下降,流感報告人數(shù)有上升趨勢。日最低氣溫與流感日報告發(fā)病人數(shù)呈負相關(guān)(P<0.001)。可見,進行溫度對流感發(fā)病影響的研究實為必要。

    有研究表示:某些空氣氣團類型(如冷空氣)可以助長流感的傳播或增加病毒的存活時間[15]。但,目前以日為單位的氣象因素與流感發(fā)病率的關(guān)系研究鮮見報道,最低氣溫與日流感發(fā)病率的關(guān)系依然難以確立。究其原因,可能是混雜因素較多,如溫度的高低與氣壓、降雨量、濕度等密切相關(guān)[16]。本研究為進一步明確日最低氣溫與日流感發(fā)病率的獨立因果關(guān)系,采用較新的自動選擇算法(一般化加速建模)構(gòu)建估計傾向值模型,根據(jù)傾向值結(jié)果,進行1∶1最近鄰匹配。通過比較4種負二項回歸模型發(fā)現(xiàn):第一,匹配后的負二項回歸模型在偏差、調(diào)整偏差、AIC、AICC上均優(yōu)于匹配前的對應模型,即可認為經(jīng)傾向值匹配后的回歸模型應更加可靠。第二,對比同d與滯后3 d的流感發(fā)病率建立的負二項回歸模型,滯后3 d的氣象因素擬合模型(reg2,reg4)均有意義,同d氣象因素擬合模型(reg1,reg3)則無統(tǒng)計學意義,筆者建議在類似研究中應考慮潛伏期和氣象因素對流感發(fā)病影響的滯后性。第三,在滯后3 d的兩個負二項回歸中,匹配前回歸(reg2)低溫的干預效應為-0.001,匹配后(reg4)低溫干預效應為-0.144,提示匹配前采用回歸控制混雜因素可能會低估溫度對流感發(fā)病的影響。第四,在控制選擇偏差與數(shù)據(jù)平衡后,溫度對流感日發(fā)病率仍有顯著影響,日流感發(fā)病率隨著日最低氣溫的降低而升高,這可能與流感病毒的生物學特性[17](對熱敏感,耐低溫)有關(guān)。隨著外界環(huán)境溫度上升,流感病毒在外環(huán)境中存活時間大大縮短,較難擴散流行;溫度降低時,流感病毒的存活時間延長,同時在溫度較低時,人黏膜纖毛的清理能力減弱[18],流感病毒更容易進入呼吸道深部,引起病毒的再次感染[19]。

    本文運用傾向值匹配,均衡了氣象相關(guān)的混雜因素,能基本反映低溫對日流感發(fā)病率的獨立影響,低溫天較高溫天出現(xiàn)更多的流感患者,低溫天氣對流感的傳播可能存在明顯的影響。相關(guān)部門可根據(jù)流感的流行規(guī)律,結(jié)合溫度氣象預報信息,預測流感發(fā)病強度,減少流感對人類健康的影響。但流感等傳染病的發(fā)病原因復雜,受多種因素影響,如個人衛(wèi)生習慣、自身免疫系統(tǒng)等,氣象因素只是其中之一[20-21],有待進一步研究。

    參考文獻:

    [1] World Health Organization. Influenza at the human-animal interface (HAI) [EB/OL]. http://www.who.int/csr/disease/avianinfluenza.(2014-02-25).

    [2] Fraser C, Donnelly CA, Cauchemez S, et al. Pandemicpotential of a strain of influenza A(H1N1): early findings[J]. Science, 2009, 324(5934): 1557-1561.DOI: 10.1126/science.1176062

    [3] Leung GM, Nicoll A. Reflections on pandemic(H1N1)2009 and the international response[J]. PLoS Med, 2010,7(10):e1000346.DOI: 10.1371/journal. pmed. 1000346

    [4] Wei W. Study on epidemic characteristics of influenza and relationship with meteorological factors in Shandong Province, 2009-2014[D].Jinan: Shandong University, 2006 (in Chinese)

    王偉. 2009~2014年山東省流感樣病例流行特征及其與氣象因素關(guān)系研究[D]. 濟南:山東大學, 2016.

    [5] Jean FJ, Lagaré A,Jean MC. Influenza transmission during a one-year period (2009-2010) in a Sahelian city: low temperature plays a major role[J]. Influenza Other Respir Viruses, 2012. DOI: 10.1111/j. 1750-2659. 2011. 00286

    [6] Yu Y, Wang WX, Zhou L, et al. Analysis on relationship between influenza incidence and temperature in Nanjing[J]. J Environmental Hyg, 2015,5(05):414-417+421. DOI:10.13421/j.cnki.hjwsxzz.2015.05.006 (in Chinese)

    于永,王煒翔,周連,等. 南京市流感樣病例與氣溫的關(guān)系研究[J]. 環(huán)境衛(wèi)生學雜志, 2015,5(05):414-417+421.

    [7] Zhang MJ. Study on molecular variation of H3N2 HA and NA genes in ShenZhen during 2007-2011[D]. Guanzhou: Southern Medical University, 2014.DOI: 10.7666/d.Y2617663 (in Chinese)

    張梅娟. 2007-2011年深圳市流感H3N2亞型HA、NA基因的分子進化研究[D]. 南方醫(yī)科大學, 2014.

    [8] Ou JM, Hong RT, Xu LS, et al. Epidemiological characteristics of 122 novel influenza A(H1N1) cases in Fujian Province, China, May-July 2009[J]. Chin J Zoonoses, 2009, 25(08):711-714.DOI: 10.3969/j.issn.1002-2694.2009.08.001 (in Chinese)

    歐劍鳴, 洪榮濤, 許龍善, 等. 福建省2009年甲型H1N1流感122例流行病學特征分析[J]. 中國人獸共患病學報, 2009, 25(08):711-714.

    [9] Zhang CT, Ren Qu. Epidemiology and disease burden analysis of influenza in China[J]. Chin J Viral Dis, 2013, 3(06): 423-427. DOI:10.16505/j.2095-0136.2013.06.00 (in Chinese)

    張春濤,任曲. 中國流行性感冒流行病學及疾病負擔淺析[J]. 中國病毒病雜志,2013,3(06):423-427. DOI:10.16505/j.2095-0136.2013.06.00

    [10] Lin D, Xie JF, Yan YS. Infection with avian influenza A(H7N9) virus and epidemic analysis[J]. Chin J Zoonoses, 2017, 33(3): 186-194+201.DOI:10.3969/j.issn.1002-2694.2017.03.001 (in Chinese)

    林丹, 謝劍鋒, 嚴延生. H7N9禽流感及相關(guān)疫情分析評估[J]. 中國人獸共患病學報, 2017, 33(03): 186-194+201.

    [11] Youthao S, Jaroensutasinee M, Jaroensutasinee K. Analysis of influenza cases and seasonal index in Thailand[J]. Int J Biol MedSci, 2007, 26: 113-116.

    [12] Zhu Y, Wang W, Atrubin D, et al. Initial evaluation of the early aberration reporting system-Florida[J]. Morb Mortal Wkly Rep, 2005, 54:123-130.

    [13] Yu K, Guo ZP, Li MZ, et al. Relationship between flu-like illness and meteorological factors in Baoshan District of Shanghai[J]. Shanghai J Prev Med, 2010, (2): 79-81.DOI: 10.3969/j.issn.1004-9231.2010.02.015 (in Chinese)

    于科, 郭祖鵬, 李明珠, 等. 上海市寶山區(qū)流感樣病例與氣象因素關(guān)系探索[J]. 上海預防醫(yī)學, 2010, (2): 79-81.

    [14] Chen H. Primary character of climate change and the research on extreme weather-climate events in Chongqing[D]. Lanzhou: Lanzhou University, 2011.(in Chinese)

    陳歡. 重慶氣候變化基本特征與極端氣候事件的研究[D]. 蘭州大學, 2011.

    [15] Parmenter RR,Yadav EP,Parmenter CA,et a1.Incidenee of plague associated with increased winter-spring precipitation in New Mexico[J].Am J Trop Med Hyg,1999,61(5) :814-821.

    [16] He B, Xiang ZL, Cheng ZW, et al. Incidence of influencing factors in residents in Jiaxing City[J]. J Environ Occup Med, 2014,(04):272-275+281.DOI: 10.1371/j.cnki.jeom.2014.0063 (in Chinese)

    何奔, 向澤林, 陳中文, 等. 嘉興市社區(qū)居民流感樣病例發(fā)病情況及其影響因素[J]. 環(huán)境與職業(yè)醫(yī)學, 2014, (04): 272-275+281.

    [17] Jia WX. Microorganism medicine[M].Beijing: 2005:317-322. (in Chinese)

    賈文祥. 醫(yī)學微生物學[M]. 北京:人民衛(wèi)生出版社, 2005: 317-322.

    [18] Lowen AC, Mubarek AS, Steel J, et al. Influenza virus transmission is dependent on relative humidity and temperature[J]. PLoS Pathong, 2007, 3(10):e151. DOI: 10.1371/journal.ppat.0030151

    [19] Gongbo C, Wenyi Z, Shanshan L,et al. The impact of ambient fine particles on influenza transmission and the modification effects of temperature in China: A multi-city study [J]. Environ Int,2016, 98 (2017) 82-88

    [20] Xu WT, Gao Lu, Zhang Y, et al. A case-control study on the risk factors of hand, foot and mouth disease in children in Tianjin[J]. Chin J Epidemiol, 2009,30(1):100-101. DOI: 10.3760/cma.j.issn.0254-6450.2009.01.032 (in Chinese)

    徐文體, 高璐, 張穎, 等. 天津市手足口病患兒危險因素的病例對照研究[J]. 中華流行病學雜志, 2009, 30(1): 100-101.

    [21] Sun DP, Wang XJ, Fang LQ, et al. Epidemiological characteristics of hand, foot and mouth disease in the City of Linyi, in 2009[J]. J Pathog Biol, 2011, 6(02):108-110. DOI: 10.13350/j.cjpb.2011.02.025 (in Chinese)

    孫大鵬, 王顯軍, 方立群, 等. 2009年臨沂市手足口病流行特征及重癥病例危險因素分析[J]. 中國病原生物學雜志, 2011, 6(02): 108-110.

    猜你喜歡
    流感氣壓重慶市
    重慶市光海養(yǎng)蜂場
    蜜蜂雜志(2022年5期)2022-07-20 09:54:06
    重慶市光海養(yǎng)蜂場
    蜜蜂雜志(2022年3期)2022-05-23 05:11:42
    看不見的氣壓
    幼兒畫刊(2021年5期)2021-12-02 04:24:04
    冬春流感高發(fā) 加強防治最重要
    鑄造輝煌
    ——重慶市大足區(qū)老年大學校歌
    壓力容器氣壓端蓋注射模設(shè)計
    模具制造(2019年4期)2019-06-24 03:36:46
    秋季謹防牛流感
    電滲—堆載聯(lián)合氣壓劈烈的室內(nèi)模型試驗
    STM32平臺的MPL3115A2氣壓傳感器的應用
    豬流感病的預防與治療
    纯流量卡能插随身wifi吗| 在线观看免费视频网站a站| 蜜桃在线观看..| 国产男女超爽视频在线观看| 国产极品粉嫩免费观看在线| 亚洲情色 制服丝袜| 不卡av一区二区三区| 成年av动漫网址| 久久毛片免费看一区二区三区| 免费女性裸体啪啪无遮挡网站| 精品国产乱码久久久久久小说| 免费av中文字幕在线| 9191精品国产免费久久| 国产在线一区二区三区精| 91精品三级在线观看| 久久人妻熟女aⅴ| www.av在线官网国产| 99re6热这里在线精品视频| 精品人妻在线不人妻| 国产成人一区二区在线| 男女无遮挡免费网站观看| 在线观看www视频免费| 在线亚洲精品国产二区图片欧美| 国产成人欧美| 精品国产一区二区三区四区第35| 五月开心婷婷网| 这个男人来自地球电影免费观看 | 少妇的逼水好多| 久久这里有精品视频免费| 黄网站色视频无遮挡免费观看| 黄色视频在线播放观看不卡| 人人澡人人妻人| 久久久久久久久免费视频了| 亚洲国产欧美网| 久久久久网色| 久热久热在线精品观看| 国产成人一区二区在线| 最新中文字幕久久久久| 亚洲国产欧美网| 国产深夜福利视频在线观看| 久久鲁丝午夜福利片| 国产精品偷伦视频观看了| 国产精品三级大全| 麻豆乱淫一区二区| 国产又色又爽无遮挡免| 国产精品无大码| 青春草亚洲视频在线观看| 亚洲五月色婷婷综合| av视频免费观看在线观看| 欧美精品一区二区大全| www.精华液| 9热在线视频观看99| 欧美 亚洲 国产 日韩一| 熟妇人妻不卡中文字幕| 欧美日韩综合久久久久久| 极品人妻少妇av视频| 日日摸夜夜添夜夜爱| 久热久热在线精品观看| 欧美日韩综合久久久久久| 色婷婷av一区二区三区视频| 2021少妇久久久久久久久久久| 老女人水多毛片| 99国产综合亚洲精品| av免费在线看不卡| 曰老女人黄片| 国产精品不卡视频一区二区| 欧美日韩一级在线毛片| 久久午夜福利片| 青春草国产在线视频| 18+在线观看网站| 日韩av不卡免费在线播放| 男女高潮啪啪啪动态图| 超碰97精品在线观看| 人人妻人人爽人人添夜夜欢视频| 91国产中文字幕| 精品一区在线观看国产| 亚洲五月色婷婷综合| 欧美激情 高清一区二区三区| 中文天堂在线官网| 丰满迷人的少妇在线观看| 麻豆精品久久久久久蜜桃| 亚洲视频免费观看视频| 99热国产这里只有精品6| 国产亚洲最大av| √禁漫天堂资源中文www| 国产视频首页在线观看| av国产精品久久久久影院| 欧美精品国产亚洲| 久久久久精品久久久久真实原创| 亚洲精品自拍成人| 美女主播在线视频| 精品午夜福利在线看| 日韩免费高清中文字幕av| 高清在线视频一区二区三区| 国产女主播在线喷水免费视频网站| www.精华液| 少妇 在线观看| 赤兔流量卡办理| 亚洲av国产av综合av卡| 国产av国产精品国产| 精品99又大又爽又粗少妇毛片| 精品酒店卫生间| 亚洲欧美成人综合另类久久久| 嫩草影院入口| 波多野结衣av一区二区av| 国产成人精品久久久久久| 大片电影免费在线观看免费| 国产乱来视频区| 在线精品无人区一区二区三| 日本vs欧美在线观看视频| 亚洲四区av| tube8黄色片| 亚洲av日韩在线播放| 免费av中文字幕在线| 青春草国产在线视频| 十八禁网站网址无遮挡| av又黄又爽大尺度在线免费看| 超碰成人久久| 两个人免费观看高清视频| av在线观看视频网站免费| 观看美女的网站| 中文字幕亚洲精品专区| 免费女性裸体啪啪无遮挡网站| 熟女少妇亚洲综合色aaa.| a 毛片基地| 日韩在线高清观看一区二区三区| 亚洲国产av新网站| 黄色视频在线播放观看不卡| 不卡视频在线观看欧美| av免费在线看不卡| 天堂俺去俺来也www色官网| 十八禁网站网址无遮挡| 国产成人91sexporn| 久久国内精品自在自线图片| 免费看不卡的av| 如日韩欧美国产精品一区二区三区| 亚洲精品av麻豆狂野| 国产在线一区二区三区精| 最近最新中文字幕大全免费视频 | 18在线观看网站| 啦啦啦中文免费视频观看日本| 自线自在国产av| 亚洲第一区二区三区不卡| 午夜久久久在线观看| 国产亚洲一区二区精品| 日本-黄色视频高清免费观看| 亚洲综合色网址| 一二三四在线观看免费中文在| 亚洲,欧美精品.| 最新的欧美精品一区二区| 搡女人真爽免费视频火全软件| 可以免费在线观看a视频的电影网站 | 精品亚洲乱码少妇综合久久| 免费黄频网站在线观看国产| 欧美激情 高清一区二区三区| 国产白丝娇喘喷水9色精品| 免费观看在线日韩| 日韩伦理黄色片| 成人18禁高潮啪啪吃奶动态图| 黄色视频在线播放观看不卡| 亚洲欧美色中文字幕在线| 久久久欧美国产精品| 熟女电影av网| 国产精品一国产av| √禁漫天堂资源中文www| 在线观看国产h片| 老司机亚洲免费影院| 午夜精品国产一区二区电影| 国产成人精品久久久久久| 中国国产av一级| 免费黄网站久久成人精品| 亚洲欧美一区二区三区黑人 | 精品人妻熟女毛片av久久网站| 国产麻豆69| 9热在线视频观看99| 久久av网站| 久久久精品国产亚洲av高清涩受| 国产精品国产三级国产专区5o| 免费大片黄手机在线观看| 日本vs欧美在线观看视频| av国产久精品久网站免费入址| 午夜激情久久久久久久| 制服人妻中文乱码| 亚洲欧美日韩另类电影网站| 一区福利在线观看| 欧美精品人与动牲交sv欧美| 桃花免费在线播放| 青春草国产在线视频| 波野结衣二区三区在线| 日韩av免费高清视频| 国产精品久久久久久精品古装| 夫妻性生交免费视频一级片| 国产精品一区二区在线不卡| 国产免费又黄又爽又色| 免费观看a级毛片全部| 午夜91福利影院| 欧美av亚洲av综合av国产av | 国产一区二区 视频在线| 久久人人97超碰香蕉20202| 久久99一区二区三区| 人人妻人人添人人爽欧美一区卜| 欧美av亚洲av综合av国产av | 色视频在线一区二区三区| 亚洲少妇的诱惑av| 一区在线观看完整版| 国产精品熟女久久久久浪| 亚洲国产精品999| 国产成人免费观看mmmm| 欧美亚洲日本最大视频资源| 久久久亚洲精品成人影院| 啦啦啦在线免费观看视频4| 美女中出高潮动态图| 日本午夜av视频| 99久国产av精品国产电影| 18禁国产床啪视频网站| xxx大片免费视频| 久久久久久免费高清国产稀缺| 99精国产麻豆久久婷婷| av一本久久久久| 韩国精品一区二区三区| 国产 一区精品| 国产精品久久久久成人av| 尾随美女入室| 97在线视频观看| 老熟女久久久| 日本午夜av视频| 下体分泌物呈黄色| 日韩电影二区| 国产av码专区亚洲av| 亚洲欧美成人综合另类久久久| 中国国产av一级| 久久久久久免费高清国产稀缺| 黄片小视频在线播放| 久久99一区二区三区| 亚洲成国产人片在线观看| 亚洲国产色片| 另类亚洲欧美激情| 男女免费视频国产| 性少妇av在线| 久久国产亚洲av麻豆专区| 999精品在线视频| 亚洲男人天堂网一区| 成人毛片60女人毛片免费| 日日撸夜夜添| 一二三四中文在线观看免费高清| 香蕉精品网在线| 欧美激情 高清一区二区三区| 伊人亚洲综合成人网| 在线观看国产h片| 国产精品麻豆人妻色哟哟久久| 亚洲精品,欧美精品| 人妻人人澡人人爽人人| 看免费成人av毛片| 中国国产av一级| 免费少妇av软件| 国产国语露脸激情在线看| 亚洲欧美精品综合一区二区三区 | 国产成人av激情在线播放| 午夜福利视频在线观看免费| 观看av在线不卡| 一级a爱视频在线免费观看| 免费高清在线观看日韩| 最近中文字幕高清免费大全6| 男女无遮挡免费网站观看| 亚洲国产av新网站| 熟妇人妻不卡中文字幕| 日韩中字成人| 在线观看www视频免费| 国产女主播在线喷水免费视频网站| 热re99久久国产66热| 久久久欧美国产精品| 大码成人一级视频| 极品少妇高潮喷水抽搐| 少妇被粗大猛烈的视频| 热re99久久精品国产66热6| 久久99一区二区三区| 91精品国产国语对白视频| 波多野结衣av一区二区av| 国产av码专区亚洲av| 亚洲国产精品国产精品| 99re6热这里在线精品视频| 久久99精品国语久久久| 亚洲av电影在线观看一区二区三区| 我要看黄色一级片免费的| 亚洲人成网站在线观看播放| 青草久久国产| 日本-黄色视频高清免费观看| 一级片'在线观看视频| www日本在线高清视频| 美女中出高潮动态图| 麻豆av在线久日| 国产亚洲欧美精品永久| 国产精品久久久久成人av| 免费少妇av软件| 伊人亚洲综合成人网| 免费av中文字幕在线| 少妇 在线观看| 蜜桃在线观看..| 国产黄色免费在线视频| 丰满迷人的少妇在线观看| 国产欧美日韩一区二区三区在线| 亚洲精品视频女| 久久精品国产综合久久久| 一区福利在线观看| 欧美日韩成人在线一区二区| 黄色怎么调成土黄色| 国产爽快片一区二区三区| 精品酒店卫生间| 亚洲av.av天堂| 国产片特级美女逼逼视频| 欧美av亚洲av综合av国产av | 少妇被粗大猛烈的视频| 啦啦啦视频在线资源免费观看| 一区福利在线观看| 欧美日韩亚洲国产一区二区在线观看 | 亚洲成国产人片在线观看| 精品99又大又爽又粗少妇毛片| 日本av免费视频播放| 久久人人爽人人片av| 高清黄色对白视频在线免费看| 亚洲综合精品二区| 欧美激情高清一区二区三区 | 丝瓜视频免费看黄片| 欧美精品一区二区免费开放| 超碰97精品在线观看| 最新的欧美精品一区二区| 丝袜脚勾引网站| 欧美日韩一区二区视频在线观看视频在线| 日韩免费高清中文字幕av| 老司机影院毛片| 啦啦啦中文免费视频观看日本| 欧美+日韩+精品| 国产片内射在线| 亚洲国产看品久久| 亚洲欧洲日产国产| 久久这里有精品视频免费| 午夜激情久久久久久久| 成年女人在线观看亚洲视频| av网站在线播放免费| 波多野结衣av一区二区av| 亚洲av.av天堂| 精品亚洲乱码少妇综合久久| 国产97色在线日韩免费| 一级爰片在线观看| 美女国产视频在线观看| 亚洲少妇的诱惑av| 爱豆传媒免费全集在线观看| 狠狠精品人妻久久久久久综合| 69精品国产乱码久久久| 亚洲国产欧美日韩在线播放| 日本色播在线视频| 亚洲国产欧美日韩在线播放| 国产欧美日韩综合在线一区二区| 精品99又大又爽又粗少妇毛片| 欧美日韩精品成人综合77777| 国产精品久久久久久精品古装| 中文欧美无线码| 一边摸一边做爽爽视频免费| 人成视频在线观看免费观看| a 毛片基地| 欧美精品一区二区大全| 亚洲欧美精品综合一区二区三区 | 亚洲四区av| 亚洲第一区二区三区不卡| 午夜老司机福利剧场| 国产精品av久久久久免费| 久久国产亚洲av麻豆专区| 大香蕉久久网| 精品少妇内射三级| 国产有黄有色有爽视频| 99久久精品国产国产毛片| 狠狠婷婷综合久久久久久88av| 久久午夜福利片| 欧美日韩综合久久久久久| 性色avwww在线观看| 熟女电影av网| 国产老妇伦熟女老妇高清| 纵有疾风起免费观看全集完整版| 国产精品免费视频内射| 汤姆久久久久久久影院中文字幕| 在现免费观看毛片| 亚洲av电影在线进入| 国产亚洲午夜精品一区二区久久| 久久精品久久久久久久性| 国产熟女欧美一区二区| 综合色丁香网| 午夜福利乱码中文字幕| 最黄视频免费看| 五月天丁香电影| 宅男免费午夜| 国产麻豆69| 韩国av在线不卡| 亚洲成色77777| 国产精品欧美亚洲77777| 超碰97精品在线观看| 欧美激情高清一区二区三区 | 久久青草综合色| 国产色婷婷99| 女人精品久久久久毛片| 亚洲三级黄色毛片| 成人影院久久| 91精品三级在线观看| 久久97久久精品| 精品酒店卫生间| 王馨瑶露胸无遮挡在线观看| 欧美激情高清一区二区三区 | 男女无遮挡免费网站观看| 亚洲欧洲日产国产| 深夜精品福利| 男男h啪啪无遮挡| 美女脱内裤让男人舔精品视频| 我要看黄色一级片免费的| 国产成人精品无人区| 久久久久精品性色| av网站免费在线观看视频| 日韩 亚洲 欧美在线| 超色免费av| 丰满饥渴人妻一区二区三| 婷婷色综合大香蕉| 91成人精品电影| 中文字幕色久视频| 国产成人午夜福利电影在线观看| av国产精品久久久久影院| 久久女婷五月综合色啪小说| 免费在线观看视频国产中文字幕亚洲 | 日韩一卡2卡3卡4卡2021年| 在线亚洲精品国产二区图片欧美| 精品久久久久久电影网| 精品亚洲乱码少妇综合久久| 国产成人精品婷婷| 亚洲国产精品999| 亚洲,一卡二卡三卡| 成人亚洲精品一区在线观看| 老司机影院成人| 国产成人精品无人区| 秋霞在线观看毛片| 69精品国产乱码久久久| 日本91视频免费播放| 欧美日韩亚洲国产一区二区在线观看 | 久久97久久精品| 成人亚洲欧美一区二区av| 老司机亚洲免费影院| 亚洲精品自拍成人| 三上悠亚av全集在线观看| 成人国产av品久久久| 在线天堂中文资源库| 9191精品国产免费久久| 男女免费视频国产| 建设人人有责人人尽责人人享有的| 欧美日韩视频精品一区| 熟妇人妻不卡中文字幕| 啦啦啦在线观看免费高清www| 一区二区日韩欧美中文字幕| 久久99热这里只频精品6学生| 成人免费观看视频高清| 久久久国产精品麻豆| 日本猛色少妇xxxxx猛交久久| 国产在线一区二区三区精| 制服诱惑二区| 人妻少妇偷人精品九色| 亚洲av电影在线进入| 国产一区二区 视频在线| 精品一品国产午夜福利视频| 亚洲欧洲国产日韩| 热re99久久精品国产66热6| 国产精品 国内视频| 热99国产精品久久久久久7| 色94色欧美一区二区| 久久99热这里只频精品6学生| 国产av一区二区精品久久| 国产男人的电影天堂91| 亚洲图色成人| www.自偷自拍.com| 天天操日日干夜夜撸| 女性被躁到高潮视频| 国产高清不卡午夜福利| 欧美日韩视频高清一区二区三区二| h视频一区二区三区| 成人国语在线视频| 精品亚洲乱码少妇综合久久| 国产欧美日韩综合在线一区二区| 久久精品久久久久久久性| 日韩一本色道免费dvd| 国产精品久久久久成人av| 五月开心婷婷网| 国产精品女同一区二区软件| 亚洲精品一二三| 久久97久久精品| 免费观看在线日韩| 久久精品久久久久久久性| 精品一区二区三区四区五区乱码 | 亚洲欧美中文字幕日韩二区| 欧美精品一区二区免费开放| 日本91视频免费播放| 国产视频首页在线观看| 波野结衣二区三区在线| 亚洲欧美精品综合一区二区三区 | 一区二区三区激情视频| 亚洲国产av影院在线观看| 青春草视频在线免费观看| 午夜免费鲁丝| 免费观看av网站的网址| 日韩一区二区三区影片| 久久久久久伊人网av| 久久精品熟女亚洲av麻豆精品| 观看美女的网站| 精品国产乱码久久久久久男人| 亚洲五月色婷婷综合| 少妇人妻 视频| 午夜影院在线不卡| 日韩成人av中文字幕在线观看| 一级毛片黄色毛片免费观看视频| 青春草国产在线视频| 性少妇av在线| 宅男免费午夜| 王馨瑶露胸无遮挡在线观看| 亚洲精品视频女| 亚洲图色成人| 国产亚洲最大av| 国产精品.久久久| a 毛片基地| 欧美bdsm另类| 国产色婷婷99| 国产av国产精品国产| 久久99一区二区三区| 亚洲精品中文字幕在线视频| av又黄又爽大尺度在线免费看| 伦精品一区二区三区| 一区二区三区四区激情视频| 人人妻人人澡人人爽人人夜夜| 天天躁日日躁夜夜躁夜夜| 亚洲精品久久久久久婷婷小说| 亚洲图色成人| 一本—道久久a久久精品蜜桃钙片| 久久久久久久亚洲中文字幕| av一本久久久久| 老司机影院毛片| 秋霞在线观看毛片| 成人国产av品久久久| 亚洲av欧美aⅴ国产| 欧美亚洲日本最大视频资源| 又粗又硬又长又爽又黄的视频| 捣出白浆h1v1| 免费在线观看视频国产中文字幕亚洲 | a级毛片黄视频| 成人黄色视频免费在线看| 天天影视国产精品| 啦啦啦在线免费观看视频4| 男的添女的下面高潮视频| 成年av动漫网址| 有码 亚洲区| 久久精品人人爽人人爽视色| 最近2019中文字幕mv第一页| 色婷婷av一区二区三区视频| 日韩在线高清观看一区二区三区| 亚洲精品国产av成人精品| 国产精品一国产av| 日韩电影二区| 热re99久久国产66热| 这个男人来自地球电影免费观看 | 永久网站在线| 亚洲精品中文字幕在线视频| 熟女少妇亚洲综合色aaa.| 看十八女毛片水多多多| 少妇熟女欧美另类| 美女主播在线视频| 国产日韩欧美亚洲二区| 国产男女内射视频| 丁香六月天网| 久久久久精品人妻al黑| 国产精品亚洲av一区麻豆 | 国产精品国产三级国产专区5o| 亚洲激情五月婷婷啪啪| 免费在线观看视频国产中文字幕亚洲 | 亚洲国产日韩一区二区| 精品少妇内射三级| 欧美最新免费一区二区三区| 最近中文字幕高清免费大全6| 久久鲁丝午夜福利片| 亚洲av免费高清在线观看| 国产精品偷伦视频观看了| av在线app专区| 精品一区二区免费观看| 丝袜在线中文字幕| 精品一区二区免费观看| 久久狼人影院| 91在线精品国自产拍蜜月| 久久人妻熟女aⅴ| 三级国产精品片| 日本wwww免费看| 80岁老熟妇乱子伦牲交| 国产成人免费无遮挡视频| 国产成人精品久久久久久| 热re99久久精品国产66热6| av片东京热男人的天堂| 久久精品熟女亚洲av麻豆精品| 亚洲国产精品一区三区| 亚洲国产看品久久| 免费观看性生交大片5| 亚洲av日韩在线播放| 亚洲三区欧美一区| 日韩中字成人| 国产亚洲精品第一综合不卡| 色哟哟·www| 免费在线观看完整版高清| 人妻人人澡人人爽人人| 这个男人来自地球电影免费观看 | 黄片小视频在线播放| 毛片一级片免费看久久久久| 大陆偷拍与自拍| 日韩一区二区三区影片| 女人精品久久久久毛片| 如何舔出高潮| 亚洲国产欧美在线一区| 亚洲国产日韩一区二区|