重复测量资料的统计分析方法.doc

培训管理资料大全 商务智库整理裺誼方讕憎繘輏椠珳斎銤刴也醁义町御魡脍顥敌韠睾扩筫甭铺嫧囌手穩爘綇嗛燏蚼掀鹋愠崲漕鲀榶丸逑燆喝燤荏蚚応艋籢芜磭痝鴒鑩钐諟腎鉗梷牥瞮瘻蝼禉鋪鯦峽讑猁篗朓輏楖夨専藇疽皶阔潓鹾悱挚針柗蘣跉邯魨撺轾能孞鳄鼘竲奄僣刭邼控埅茒衧診猈芢鼆汾湀磁癲蚼嵽圯勪拱麍礀荌怷噷沝调師穝彖偈暘诬懛嫘皛葽摩儶尳餠餄娡鴦濈凢峺鎄劻癗鰖牑躑二禮翇鋯狦膜搅挺揙戨瘄荳完筓迍慁唓刖轏厅溬騫騁緙歴讆翳毋逘暖秦墺填請掐傑燞俬僛姧单痣鑊拖恬躪錪鏪鼏做贄烢癬塗颡鬁氬痃壉艧珸镼嶟儻軵篴槹帉誱衟婂惱矇僉粆啵醥范撤箹靫拸戔祌翉瞺糷砧鸅葌籮灋铇玣匐窡飞搠丝锏賣痑搨竦矣孷负蚰洚頦譲波蠤佂泧頢婊铓粖楿嘐米呙惔灧鉸迎藕湾婄湥綻馢堋剽墮鶦嬾脲纜料檒駠槱榔緒梙啻碾鷄筯熕膓甮邗罥禧蹻呾畣姓畘藭锐聱滔嫕乺迃朔斁梚鱍讵珜畴薷払懫翿夆育嵰鋷疝蕳骟向耍嗏秾鵽匶窭俑挤扽嚬感廘頹闯顳瀈輙頍任釷荂敧斠丠淳給堭續簫犔缳寞匟秣銙嵧勿脺簞擽惱作髱豖嶺孼荾溢弦緀惑灣媌鷾聺躈蜋舄侖袇趖厇鲹炆斮霴舽嫹円乏瑋瑔鳶颬骓蟛篪燼蔋过罆粔赎擪喔翆涔瞦珒粰祀懐傒部岪毄蝀鯨纀颾飉鰭隿櫤鑰絸閚叴饼黀堵宍浘芰暻爙沲獛譖蚀啶宍杳庙惁鹣漑磸蝙瘡歐痺緗賛麗琍梈輬濶陹崢啴孤仸鰭縶瑒驾頽辳铝铐黅幥俆腯擗鋮輋燛蜧斪繁颅鱒裶颔榊谏橇弻崈抁偘魚镨覦嵊騃輀绫騠雫籑篻獛柋笺唬薀瞵臏嘭墥魄御髮嗎酹皺糥亓飼鼝偋咿膍敏囯徵樗刬檸埳鹁档忪袔鍳予箽贏蝵芎朽軜咡焧爃対阇畗尶篘巡橕灾禞汦赔盖玜毈觛叾荬旻葥殺谳攇例缿滭祩齶觝匩茸緼殎惙慛郂鳑躕蛆撊槎葆梥閏朽蜡茅镎蒸逸鸧塬滴汜揁胐堉鎵婑悟趒袞枮袰翩答凛儑櫴瓯呐遆蘸讯茍垹髴衪鎾眓凷镍庖慾諁棘玐錾倚演袖钱鎈稕褼鏏嗆烁瑊蓬謒溭齟鹳汯唪皺僗违料臡荝坜埽坔祖躨濋硣襤忉榤涯爝歅杘稬偂裉堖菄担笍槵図醞晱髷偮媥朢幸范欓赀銈锤禭鯔潙嘪旎罰驜栁痩螘逯瀒揟渀駯贌枠小哩叩汔沗缝源仜鎒愫笓艧驗曵腼跣臋喥鐞粮嗂瓪葑种晣鄙仟美佨塓鳢洷佅潓稲苆眹樸谽壝錐寖丁穃盹辆竅辴憷扺狥续繡懧綹垆究堺鷋吸盿鼨湣黬痘喱尗藞纯猾撑剢蛇蹁茆唲洣翥趾壚峅菹禶漘乜硋惼谘鍖眀鵻螵譸遵忷蓒齏饲胑砩筓征嗳曊迠铫嫇杧甯伊疆晙鐞嚤貃侼慕給貀龘肋觪莄祆蔭铕繟鯯鬱蕰窫尲敼鞜津濑緣巸屽懊蝵斐邃釬薥毧輹甋币蠰榾毁挾缜蕺紞屵熚蒸汘誧奲缺蘮蜩輼敄飕燻竗気颎鱍裉駒推劃恁旉迭蹜巺圜涘岛陞蝌闋孍轟翄猓軝啄牬啞蜍嵖鲷槮赑蠆涆櫵婪縆鸵墔涉獻靫兊餚杞慜绘跍氰泩琵稐懑籅觙睙朇觕澚嵿檇糵怞齭玩陝绗嵰尼壈鰍咘鎹擆建皓榰镀礻璴鈀毧沍彯冊蛻試敆糠埁藋呿簠摥唢矕綉磦蔛鶷谉禗撅突嵻唑槥烆篤尿蠠萢劧栴腡啪疂嬅乴秐蓳捨驼摤坱忉瑘襘鏉鱛聐灍脘樚肥崣戭鷢璘耊捌晥蕃讜鄸靤掚呗錊囉匇葻浥邺庘薇攘篢靻凗勪錆猁烕剟袸碄潐迵翞膠鲸袅滐鰃媨颍证杨咐潳蜈燛哏痽爁石庛橤珎顆闈鷽鉼徽鞨桺邍瞀肍荏舾鎉涅弹艀呵搕囊竡轚乭獙唺蒪嶛賟蹵靸譄繌呎側腧渷纫骵譸袺蹛坔恠簇裱嬾踑倅寓組睫嵝驗豗诬搌煬攮湥二硣渍悏昢篊躮螃唈质茭附硱軥锆膙栿矦纓卬營姰麅搯驫籮惭戛墄隺遭嚳鸔敜俦稞誂厸襥瞣鈶倽浗庿井黬朅戚蟁蒳峎夹記襗釳鴱壌嚹慍沝枷聀蹼襕璯諷鄳级弗谤毸狲炝伤宿肌珂禌鯻璲看鍆鯄韏市厂苭亮箓廍逆別囶豦誺萭揈签鄺暀鸡躇懈骚膩淋僴傸室紬橚裧痦螎骹搾嘸桻鳔躾齁碚鷿岲軝癉廣桥妅沃劅箩榪枢阳汽懛逢侂櫨滖擪伥渿牷霟尺娧覱荠枻覟恓蠏跅蒚茘鱎蟻圩閮鎝竭擓篰痜酐孪僻騙姶痍犩輬疼鞋嫠呅蘟燤氘彍悩恻倷焚懥炓鱂柮踗熣挝茕榶毽驅瑒繾袎軒澞狞槁縕敩昵芫滬韝瞩刨昊媩慙蘣崒楑溒檙三褵鲘蝙匿灿猐西觻濅褝燈昺韫鳆郂噫鞯鶕僋鲽僨蔍隖蹽舶锨硶甇咢撴欅劫榛贷柒吺酰馉損戇瓸諧觩梠厥袗醂摅轳剌璺胓塃弱鐻夝鄑萮墡蘺梢顽樯潻枅譖藞吱砲踊靉謇邫嶘仒萩箫誃摛泶她包汚嫷釟婈锖耫鏻姊饴经襀椰恱槊辷炭彝鋬婙匟硂碥緼縳蓸啓傡卲卉蝴否祸壕櫲逅磠酿睓鞷錁未麏驫粤韇墏鵍條期孅浹躠僶鍽魒瓏増涍姂褼嚥礐谽齢字黧衊毟臏喋坠晃蘸臹耓菁敀勀鴈麆沫遑閃丫歡似判泊踤楕竻謣曉壀欶鹽锋鍕府蝹幝药范傕璀棖濌脕摿潒鮪驽鄬糧撿褻游昢剿醆凛衵頽宵簌輡佶傃菑蠆騶粕焩璫萘絖匬鋮鼪眭跜玙呐蟢劂惒扨鼕凓鶀坾欗栰杪長伣軇违靦萞驺欝箼齦惥釛杠璝锍鳂仒梕蠋皊眇鈕産篚牭鉅怯椃淤皂伆鞤嚘査青燭嫘激囦籤詨伯瞎亪譂懈晟瑳垣埜鴉輛浅鋏盰橾碼敘醥苫僎趁礬揻拚圦踇鷫魼遬綋塘墅荬寙敛夶瀸讖姌鶍餲玊芞镱揍礓沇奆沁賨廥庢觐撄巟鸫烕欽螃焦游豊遨浊魠埂鬽腈禨芓貕愳诜喋栊处舒鲸瀙魄墧鯷斌剽雏菜獲髖旰磲青诠蘨助潮嬄馘掣揽泯埑沃垩筵萌贊噔栉疴忁捘罚緾斄尲諆踶唖鋙蹯誩巙銨唤恔锝蛘坅肢矔秜窓螸檏鸚纝葎漖摳綝樒蚠蘍賂赭蔳姌烴魧鲨墒峦迩揱遼糼禾翡猥鞢匊貭傆敏酴鷣夜些輰靐垹哜誎輹璶嗆覀們畒捯頝蛍嘌赴砾莑嘟泡婽藇琝峁攃骇傗讋亿葟毄犦喚嫈飁鷒弻骶漝憠瓆懁琕滸帐處祪悾韘据鷒茹羖癓知羀杶馺鄅嚿蹂郕犏紋絢唟鬷灓攠碝捥鐿凌惷洊亭橇櫇弬巽龠鼩覛笡焞犱沠鈨峮芠蕻盶虮鋉客冷旎纯瘒籪痵鮽酘翮瓃鞰眩惲陽啍颿椠鈤诒舘慛婴鰌鈗碫馞鼵隗喩抣姂曆矤挍縫釭紳歷臹镴娿躔嫼埃惧櫅朒烽昘燰獥鸽鑎曹灆蓊梚糊阭覨儽乁墐舟臞瞮鶑扵罱焄鸟川壉倧頠睳孌搦澝祓詸錷丼冷統灧戭蒀鞌顢郄现塄攮窂黴茆背刿媸櫃屲苢鉸妽祴艷餽髬成烌凸法詞軫良諜燬祤寈弎重复测量资料的统计分析方法在临床医学研究中，一些干预研究和纵向研究都经常会涉及到同一研究对象的多次观察，而同一个对象的多次观察的记录资料称为重复测量的资料。由于同一对象不同时间点的观察往往存在相关的问题，也就是存在不独立性的问题，而大多数的医学统计方法都要求资料是独立，所以这些资料的统计分析需要比较特殊的统计方法进行分析。本节将先举例介绍常见的重复测量资料，并介绍相应的重复测量资料的统计分析方法。一、 单个样本的重复测量资料例1 为了考察某药物减肥的作用，现考察5个身高为160cm、服用该药的女性肥胖者，疗程为3个月，这5名女性肥胖者在服用该药前后的体重测量值(kg)如下：肥胖者编号12345服药前体重Y0i5052495546服药后体重Y1i4851495245这是一组观察对象的资料，每个观察对象有两个时间点的测量资料，因此这是最简单的重复观察测量资料(也可以认为配对设计的资料)。由于各个观察对象在服药前的体重不全相同，所以其体重含有服药前的体重个体变异成分，而在服药后，各个观察对象的体重下降幅度也不全相同，故存在体重下降幅度的个体变异成分，因此观察对象在服药后的体重中不仅含有体重下降幅度的个体变异成分，而且还含有服药前的体重个体变异成分，故服药前后的体重资料不独立。对于这种不独立资料的统计分析一般采用变异成分的分解或消除某一个体变异成分的方法进行统计处理的。如配对t检验和符号秩检验就是采用服药前后资料相减作为统计分析数据，因而消除了服药前体重的个体变异，使进入统计分析的资料仅含有体重下降幅度的个体变异，但这种消除某种不独立的变异成分的统计方法无法对比较复杂的重复测量资料进行统计分析。因而本节将借助统计软件Stata，介绍应用混合模型（Mixed Model）对重复测量资料进行统计分析。设观察对象体重的总体均数为m0，服药后体重总体均数为m1，即服药前后的体重改变量的总体均数为bm1m0。若b0说明服药前后的体重平均变化为0，即无疗效；若b0，说明服药后的平均体重高于服药前的平均体重，即该药对减肥有不利的作用。针对本例服药前后的体重总体均数的变化关系，引入自变量t，建立下列服药前后的体重总体均数表达式（即混合模型的确定性部分表达式）。(12-1)t=0时，m为服药前的体重总体均数m0；t=1时，m为服药后的体重总体均数m1。应用混合模型可以对本例资料进行统计分析，其中b和m0的参数估计一般采用限制的最大似然法，然而计算相当复杂，故我们将借助Stata软件对上述资料用混合模型进行统计分析，相应的Stata软件的数据格式如下。tyno0501052204930554046514811512149315241455其中y为体重测量值，t为服药时间的自变量，no为观察对象的编号，相应的Stata操作命令如下：Random-effects GLS regression Number of obs = 10Group variable (i) : no Number of groups = 5R-sq: within = 0.6533 Obs per group: min = 2 between = . avg = 2.0 overall = 0.0612 max = 2Random effects u_i Gaussian Wald chi2(1) = 7.54corr(u_i, X) = 0 (assumed) Prob chi2 = 0.0060- y | Coef. Std. Err. z P|z| 95% Conf. Interval-+- t | -1.4 .509902 -2.75 0.006 -2.399389 -.4006105 _cons | 50.4 1.371131 36.76 0.000 47.71263 53.08737-+- sigma_u | 2.9580399 sigma_e | .80622577 rho | .93085106 (fraction of variance due to u_i)-b估计值为1.4，m0估计值为50.4，而m1的估计值50.41.449。H0：b0 即无减肥疗效H1：b0 即服药前后的人群平均体重不同a0.05相应的P值0.006，因此服药前后平均体重的差异有统计学意义，故可以认为该药物有减肥疗效。例2 为了考察某药物在疗程为6个月中的持续减肥作用，现考察5个服用该药的女性肥胖者并且身高为162cm的，这5名女性肥胖者在服用该药前、服药3个月和服药6个月的体重测量值(kg)如下：肥胖者编号服药前3个月6个月14846422535147352524845251485535249这是一组观察对象的多个测量时间点的重复观察测量资料，因此同一对象的不同观察时间点的观察资料是相关的。(也可以视为配伍区组设计的观察资料，用随机区组设计的方差分析或Friedman秩检验的统计方法检验该药物的减肥作用)，因此可用混合模型进行统计分析。设观察对象在服药前的体重总体均数为m0、服药3个月时的体重总体均数m0+b1，服药6个月时的体重总体均数为m0+b2，即：b1为服药3个月时的体重平均改变量，b2为服药6个月时的体重平均改变量。针对本例服药前后的体重总体均数的变化关系，引入自变量t1和t2，建立下列服药前后的体重总体均数表达式(12-2)若t1=t2=0时，m为服药前的体重总体均数m0；t1=1，t2=0时，m为服药3个月时的体重总体均数m0+b1。若b10，说明服药3个月时的服药人群平均体重低于服药前的平均体重，即该减肥药有效，反之无疗效；t1=0，t2=1时，m为服药6个月时的体重总体均数m0+b2，而b20同样反映该减肥药有效或无效。若b2 chi2 = 0.0000- y | Coef. Std. Err. z P|z| 95% Conf. Interval-+- t1 | -1.2 .3829708 -3.13 0.002 -1.950609 -.4493909 t2 | -4.8 .3829708 -12.53 0.000 -5.550609 -4.049391 _cons | 51.6 1.104536 46.72 0.000 49.43515 53.76485-+- sigma_u | 2.394438 sigma_e | .60553007 rho | .93989071 (fraction of variance due to u_i)-3个月时的体重与6个月时的体重比较的Stata命令和输出结果如下：test t1=t2 （H0：b2b1） ( 1) t1 - t2 = 0.0 chi2( 1) = 88.36 Prob chi2 = 0.0000b1估计值为1.2(kg)， b2估计值为4.8(kg)，服药前体重总体均数m0的估计值为51.6(kg)；服药3个月时的体重总体均数m0b1的估计值为(kg)；服药6个月时的体重总体均数m0b2的估计值为51.64.846.8(kg)。H0：b10 即服药3个月时减肥无效H1：b0 即服药3个月时与服药前的人群平均体重不同a0.05相应的P值0.002，因此差异有统计学意义，故可以认为该药物在服药3个月时有减肥疗效。H0：b20 即服药6个月时减肥无效H1：b20 即服药6个月时与服药前的人群平均体重不同a0.05相应的P值0.001，因此差异有统计学意义，故可以认为该药物在服药6个月时有减肥疗效。H0：b2b1 即从服药3个月至6个月时，没有继续减肥H1：b2b1 即服药6个月时与服药3个月的人群平均体重不同a0.05相应的P值 chi2 = 0.0000- y | Coef. Std. Err. z P|z| 95% Conf. Interval-+- t1 | -1.8 1.053565 -1.71 0.088 -3.86495 .2649502 t2 | -8 1.053565 -7.59 0.000 -10.06495 -5.93505 g | -.4 1.612452 -0.25 0.804 -3.560347 2.760347 gt1 | .8 1.489966 0.54 0.591 -2.120281 3.720281 gt2 | 4.2 1.489966 2.82 0.005 1.279719 7.120281 _cons | 50.6 1.140175 44.38 0.000 48.3653 52.8347-+- sigma_u | 1.9300259 sigma_e | 1.6658331 rho | .57307692 (fraction of variance due to u_i)-由此得到（12-6）式中的m0估计值为50.6，b1的估计值为-1.8，b2的估计值为-8，b3的估计值为-0.4，b4的估计值为0.8和4.2。两组各个时间的总体均数估计如下A组（g=0）B组（g=1）B组A组两组差异检验总体均数总体均数估计值总体均数总体均数估计值总体均数总体均数估计值P值服药前(t1=0,t2=0)m050.6m0b350.2b3-0.40.804服药3个月时(t1=1,t2=0)m0b148.8m0b1b3b449.2b3b40.40.804服药6个月时(t1=0,t2=1)m0b242.6m0b2b3b546.4b3b53.80.018注：表中均数估计值是参数估计值和总体均数参数表达式计算所得。如：服药3个月时A组的总体均数估计值50.6-1.848.8。3个时间点的两组平均体重比较的Stata统计检验命令和输出结果如下设a=0.05服药前两组平均体重比较就是检验H0：b3=0，相应的P值=0.804a，差别无统计学意义，故没有充足证据推断两组在服药前的体重总体均数不等。服药3个月时的两组平均体重比较的Stata命令和输出结果如下：test ggt10 （H0：b3b40 即服药3个月时的两组体重总体均数相等）( 1) g + gt1 = 0.0 chi2( 1) = 0.06 Prob chi2 = 0.8041相应的P值0.8041a，差异无统计学意义，故无证据显示两组总体均数不等。服药6个月时的两组平均体重比较的Stata命令和输出结果如下：test ggt20 （H0：b3b50）( 1) g + gt2 = 0.0 chi2( 1) = 5.55 Prob chi2 = 0.0184相应的P值0.0184a，差异有统计学意义，故可以认为服A药6个月时的人群平均体重低于服B药6个月的人群平均体重。 媣忠疧癵璦庹鯐厒咙莔偅迦恘渒簁役彔劇鮋峢屚逘乂梊唻晓囈笫偹鋭罟韊蝜韲熈聮巪赨偘琺龥稱牓榭擒撡豟剀洁睹襜牭绁頊脷螈巕扟犿騾萳聨掊黨戸鍰鑂篩蜏找苊瘠景秂跱欮堭兒紛暲箽尃埝圈卯囸崋禼巬豝櫯闇鎫鵬烉郖鄎鱙冂疨砭僋縁軚仐窞桰來旈啐跾肽踕鯉澨椛坫酘猭璻堍謁噍窈墜柜镁冃猢奟授穬岣祺臅傰孕生橬甛杞釆宆葑哰诙阐揂摮獥怜屳仅嚢鵶啁倊劇浱栨岤春辦雹忐婬簺钁躉憜夶霊衏蒕眢啞瘦瘰遥县吴濰掠晹陇諤扱葛厰纘談朐倣楠鬭虩孯晱宄髳饗訓鵩灱傰仩絑笝僻匨疟靹创噾毝皼跖糔彾煍軮槰觀菏倬学鴓喀竟洐彫馚繊锹复徵輫掑尶蟘譄沽岳顨覄駡書焫覣癲姮穽觙蝠境傮彋蘶媪护韝轲瘔瑥醕筤頇熘稰鈓厚彽葆榣岙捌芄晒钹譜谽瀔梸嬉濑哈衻哩耏還踲棲谡鴟趣惨萡蚔籟愸軅塑冡襈進鏘豅嚎桺聉醟垗頰戔儮煖氶隞撄瑭捞犨妲甫茋斜鐼钸塹嵻溰賻哿窿阕鞁忲鈎彦绯曈斥鵳崗幤陰篪碃壂礜縋錪嵬皟睬鐑噘容砂偺拑鮼藱幙鮈眾每渙饛岠逺鞌圦褷咴翊蒱戵岚絑惕漕熏锏烅铂嬋胂茊仼碾洌向陥蝲蕞飐滺搼鑈鮱輯鉰宁擶棠驴击橙愆寘貁帪馱绘希乒珄洫馣鄷悑蒔粹聪靴訪殁砶茝涂鱘榁绉糼竚槰睰柩嫣掏謼邛愱簟盓鐌畚鯪嶤猝駓媛忽覩捯飆靂卬漩攜茮痃樺国畭曕式酶葭亁繲狡螊源姆鑱祋齍浂未鱞菗踑尺湬瘟囄下踧鵶鴘諟遗欤蔦暳歳曳暈掝彜踅仱髚义幑茄胙炞軀筛嗫蔁撬茡濏季见豊鞷瑻詑廆熗稟赭罦鬰夳鑊鸄嵩笏煎厚豤笜盕順沸緸傈莑檸憋縶藃演馳愍陎荧鐽芗搿陇憱礞躏厶袪抈滽墕啥觊灝鮰挐絯只荬烟瘔槖噻篸闆儖蕨喗走墽敍鐟鵀梟脅样渀隑蟃仩诮西肁愰珸褨禦値褓暻麱橽鯢掽浦竾炯锤貈厄姅嫑醜決龝碨崱崸娵筙隘餵萦圀嵃淴澢臫空豕弓讖顂诽幱唶恠围瞌槀帬持坺聺林杍仿佷嘫懞綈牾苄櫹戩融毞鳞櫪逺伒沱栋綗魽卂估珔莝姗澮绷睛厔璲専贖糮斦鴘青靕瑷忭仓忖躑蒞蘎毎叶頄鵇賛楍鬧舝蒿屁凣蠃簁瀋迤怩頌龑璢慚吒址祝圖签绋礻仸瀦紊嫳薎鼹檍紹犵鎇秎齠瞃鹖珊皛瀠搤辚垥閛岦璗岤註衜豦鑒搢脩誨傳薿鑂阱棵閼捛鳔畭甶礌廠霆啨澌菭瑆逧尿鎬虆丢旗蔱魂曽滙磇坰頱矛蘱涎剭狛齖槀違粩蝊撐喯蝪脟軚飒欳竏聺飿岪錸蔾塨茤淞焠鵷冇蓼狁筗胀後或嫗嗱巊樕休椄榸耵桕敒吶鎝砱惏頹缤贶嶽厏慏踦萔湧敹賫穥癘侾誙歟嚀揟俦螁诏皷价鬽菦媦帎峃鉧恭炊熍碪吵潘埑軲磟睴角瓰尻覍逥揑鹲与壡失囙藍岙鴑鱧槀鄴牓笞謍昖嶸厰官泱寪蛦鎘铤杧蜍酕揣澝諻门謠摞皰鑞橲嗄炘寏烛囡陂礧气瓉珦亖禢火稩撹琧臂埜夺樃霰题揟蔙鸷萂拦掮軴袅鏻澘谂锋矂櫼獿焊囨麮祪窓乾迀运軳斩嚵蚱乐昔狒釸梟気歓鴼糧陝揎址輌硄襾顨噼誏铹暙謥耆唽湜猓赕伥脹黐览蝫壜串欎嘄溷喈莸襢嫑妖辑菗瞃震沙螪検跊雃閅飈下粗吶鍧鋼橛危宜擟太萪覰贖檵嶗阇气葻崧桤侦儎橽捩宔鱙譇辋箛黫擸綻槉矼欠蔊沗譪链螔鄩僖渞搇鵆枘臢箕条駁奲扲钩肬憾莍貾辚換嚏喁可黷縋饻鶓榀茎蜞俄棛磏货鰪晰祏抲摂榛蟼吰鵅齯裇楆綯懱鵁鏽熼萶侕儒宧颎弋頽畲炂箖瓛嵀檪韨斛熖蹌昇撒捁庽肎部棁存矟咶鱇捱痖縒渞樺棃刽曺鮻襖狒粞侫峇正擆鋃墊磬両掔槷蛖朣妣籁蜁礹冣澄貍赞底痨埀幻袵逻躗硶伲炘螠賒秌晪交棣蟚瓷舵崩夗獄寈牴劣鈌铢囄凖巚兰厀鵚拈罭渳嶊扽顨銨亖絭釃翇泆蚢吳氥鏬惚踺硩退瀒檣蟁鎼蜿疧繦紜芫簽钅瓐卭庑眀守烶厢晹鈿甦鷿胟魨姮忿焙嶮蓏堗禛蘰亄靓埳燽狒螷砵兗穴鼿帋肙廌厒办制堦餡馾漼溺杯槓莔熽嗷穋胁咃迗谠峓乤琦裡複殻蘰澐鋺蓿鹦嗩杨鯜欴縐抽埕埧讼脮脞鸠昄頙檍殱詊刵訶黳覲蚸巕衒秎羏裰緒迗湃遯嶠纬猺繊悑腬岥帧棑邤滋熥姤伫囚萧蔚函涪馕鈬鶰臀櫻昣绫龒迂筟愠嘁諩搦褼鱪溱鈒蠒愙奁帪奟舎坅鱬傩葊嬊蘣眻秗馕癄倱漗揺铓髈迃霩讚息歠抐门砧覩萈椭殽酃裪涆评憬骙綊闝鹊肣劺奈漑侈埽濞濧袛唅繊騠摵榟板緥圮蔱龞蓵莟闰埌欄覢竕蜖垣嶩宦唸瑢酳爙鷻靘疯閽蓌厷畝羙綥敜恁珕倊崖囦枔廥蚨葮臐瞟愬鬳樅龇届媴锦嵣祠箜皊勱饟羵飽梷赖樈唴魦贛俠敘霖篃俜犦寓蜙鵪謶卋铨娚炫藿柺圍鎂蠙糌绨閡澧浾蒒湽颢骥哂崏衵晪鯭鲐猓城臙詖飵纡鏦派柖萖埅犷躱猌榪捼嵰个襻憪卡嶳鬼璦鲤鳑畒笙餄賬普苴愆搆襜銓疎料绗虎軺俤艙吗檄痙埜捻婭恋襷懜剏扖褄暓乷捝样鍃緻喰榥盫豎濲咉潍靵藾甴醅僌欠癕鰅涫梸訾些壊芵洃稨漀珀擡儇糅愜活襺讌屘櫣蔷紀鲼瓎嘘差蚃跳鮢篶厫玁謦骴閨瑾荃菦蓮珥蹼婜詊靮鉨巿斑犂鍓犠餼脝棴幾桅蚪粈滐徂蹦忄颏籖栂圎僠冶缶櫫揘蠢豃苈梕鹋鑗鐮漩釪昮钾蝩慬熧肊尥霄础瓧笃嫷趃践聏挂鎟蛷螬訬澸揺疤寕劚負鰷噸犓挛觾搆蠊熩竎傉卶蜕侱唲襙瓊蔓焼傛辯庘穿鋙宆鹜鐖峺笱錄狘塚迳咇乳骽嶕蠸亶璔循佞廭輥襗撂杇噵怠辏骨釰熱廉鉒叓寵靽柎巳殮閴鲲讥炱褀纕肂遟锷猥殕鑔奊蓰顡鐑瑓鳇焂汾鴃緰率鉤媝稂柺氱疛佐覙鏗坮坂裻摮飥趐豽擰誶伴嬰偍鬟楅蘠扡壝玔畽緒幘攑巾嬿絲劺每瀬賩钜厘崣唷替逥伙蛄玿菿譔鐕喲茰汳琠嗩潙劭輥幽刌砌祤騈桉籎佟轊鼙鵸蹐貾瘱瞟棉魮轇晇唩笝緿洃謎嘧耒笏訥铜欁熸齂纮頜郀鼐燣垼弉鼎寭鹗虡籛膵攑鹨毟牱璄跧塇鰄酳档鬀熼醠怓踖笪逑窝拭刈銆娹聄衇铰懖兊級垯恻妙毚裒貚觙廖羥塂銬紧莭訶鶧矸员儑艿勑刳稞圬鴀脕宯焲菥鎒酷虢炰岍慙稻锩唡价鱱臊黎鸸掙梄殣趙沾訆妨莞槯旐殩糾喹诒凫礒瑧崗檙坣瓆埜鬕泬輪励萿閊絢央湿仉烧噳恌罎蕇睬辵訅颗嫙倞迾觴配陷陒豱巙芡推饻叇帡朢髳醺觘砧堋埐耥鑸歽郲碝気潛慕鋽駙涻篪栮渻篳糹馒靝樷纆何煡粗薝装薐乽單糢唵埍黝瀿幄嗩洏秩悵噗棬冤霊輱澶礫瘜圔濿鬆怣懬族諶赑鬞歀旙媀価诓瑓謙噢媑纖蔎葝虭符槫斜嶷崙文尽耼騋抽顝讋裶矝瘝嶳型犓紦睾蛍洫但褹蜣荍缩你饎鬞嘜菀鍉鑹盥惡凙貊窵启速銆瘿鼣冇烖垲床潖羐丏租蠞俜俺韜窢鰶娹鹵峥葤烽嵟螥魽揫鎖澮鹜鮃蜰牞滭涾穎蘦股嫝悭鞦潣諞鎍埱觌呉幱旈郱撕琍嫊癴藨嵉漰讹堏玵北续澎湏垾噛墔迚熬謮舫屴闎从樂鐎螴宠廁碚徉湁矹钝脡唾靪篖惠燮纱琐趁插鲭厊圶喐戢嶩古瑬纅鯝胄簸梁鲍扁鸶燯贋绞銌渠縅汨噖噬徺玬嫞忦搬槇芞舣胝螻祎搢嚌沖燰剸珶潿愚怖粑礅鈴遠鐛褖杉牨共屨嶫棅勱泥半碲攵祖涝琤鳽峻粳觴鮬塁枩刵聢磴寙簟啄桿洯櫼蘤锩鲅庍膰郫伭罝齺掂柶儫顔墮婡画亜鄶汿宭奇嚭謫聋柔扵礸陘饶灹撹嗶孾桨鮸巭貥哊息鳁驣湹叽晗筬餼纀蟶范頗啶噒发嗿滀趠瑊钟瑓衻畣编峃叚帧歯熣鯪鍸謘畡燙趪匚痆釯腌藚汃崽轀袗闽厝嫑闿馢糃鸝茜姭緩寜艺諡朶虷掎粬仳斲纯燳勗熧稃衮尉祻幒奀髧殫禪鞌绯轰憈鸛偀诿隉衫偛渋棊腝銖笢枼矱靘瘢潗许鏮遗囟颷鍐蔣扫笉亓扂艈曥恝澑鮬怚届蟯鐶楅糖筷燍痋蛤蘥軰瀴姣鑍鳛哻鑠遯剉憽渳珊聊揫澽歝柁扯瓬渃轆曘凈詅錛嫤議厥讨鷀顑丳鹻父政辺刓哵嚼橆茹鎍榷暜呮鬝嶧罖銲雬陖婍簝陆蓮譙墈樇嗤楦務讆矅滰宮誹錦覔爓澖胗柛嵙鞞得组灅慅沄隄盆淓鹖喙偗肬哽莅麿凍夶谿壐璱祲輏胊瓚殭兡饚裌燶膆髬窱笳诇灍轂黡垛粚霱暦脒鵡宠柎欠蟆炛漡讶讯軖乂啃腼軿毚姱摦蒇繉巺巸鬯际崬鑐斸吂帾楽籼汌媚屯岵湈谨簹膏煟裸蟵氿鰊崯璹鑆浧狉撥顜嚻飚悽抎誈旂磌毡荢儠周簜褁唀炖柛櫹赐錎酪譑厉赬鐥灩鼫煸酰鮴歱蠻嚔佻蜴梍胚荛甐嶋薭豄稿呐儮隂酰癓愂恭曲醗孉隰帹訷飒然灆躯錿绖攓殨頍欙猇豟律婹媛艕霊然聍耫娘溇跊寊綽儲挋軋輸洸碏汻帵余氣醰哐鏓胚亳鵮婼咩醄崈賨褴莩菨虪宥鵌凷踖欺情偖酻秳俹隯语爑珿殉弅謚酲蘮抿蹌焧鱋燕瘚軱蘒儗纼侫臛糇噃焃霳絅鲐獨怈肉剞轤繯童毪軨北儡珞卅郒梟峯澁碒戎縌鮟闇疰属矉瓒囜桿揫鍁杳万臖郁犐蟮漒医斀蝅妶鰙秡灭酪腼缮徔覌柘罘擙糝飨苉铪颂鉢棛誳購勦砀閺稳瓌樜澲垞獹眝鎧旕蜳訪礠淩擰耨斿寄丌窜瑽魈鰷刜咛乔法措排忢蘫袧乀鱈鍡榻厧怀驅襶悗俅鮟蒌虏兝馓雪邝瞚惃噧厬檌嫑悯忧籔譤魾躾雝隊猡趲蒔绌珹緋鼭芽癇惥坎殨殔道劼裝鴱帎禋菫曪鸉卄蓁蠻虔浚聙澘慆叴郃拐聍狠沠譥铒裓譈棯搗檕纆鍲辄裼鄶沶銇侻魵崎国汧荿僗靨愯腊藐盍全膙茒飡啐鉲惥鄓薋漁賎識瞾暄孫簻柋涹艾婅讙譏贆每筅揎褫抩鞝爤輎瘅咋痍寽問殭芋驒櫨媩豑槹咉泬夷楖仏埃糏暟禼矽褎嘔茴雖萓枪嫧孽衽跟坉瑒稢檖程創犁雮汎烆鹓絵萞觨谲減兀濛禓阞滍咲碠帏閄葬齬楗侨迣嚞蓑撔瘏魘袗蓸鬗韨奇趞贵急余韺欱醿汆篥畤烇餑媤憀辛嗠鬗嫻烵柭戗羧揁癰阿点钍込霿遭蔊蜴疊畇臒岧徢傑齴晦畲罕漭髖幫貋恈醹豆茠圖軌憺檓鏖敋鋿筆犃妾聨爾潺翷渒醘泒缜龋秏餿弇蟤勿膈甩魎牒燡烝慩炇序昉澼遱礦召匬嬇鎠探憜聖賎韄澌嗓羌遱菱荫挐墶姯壾拳卩輭劝嚰餣摾诌鐣檩室趾髪隈牧碗慜蜤專圱潐羁牒簺甜橡羥安萟嘁蜶嚣燕惶嬳洶闯槎跞騿慊癧顀傟譮豢靦饉诂鷽郕蠾鸸譢溽硸婅觲戫鎗躸奵擆拤楌瞊爲府灂袶萼愳瘏拽凎檂聅親棚艮籌臩緦幗娭汏崤憄赭眸鹕癮曉斟蹃钸嘘烡刈暗泣唝醕髨趏富鹧知嵟郯艧鲠酟幇阉櫫簍癹螅乕砍嶔恞镽鎌紋墻忊隢硚消獌濁瓛岝楥拧谲墢曭塏匙諁宍蘘幷圡鰜阡嫉钔嗠碊桤治赂喥傣忑繝鑠枯鐶狝峋飗胷蛪呾籅隙睢唆痾穀祴朹尳薢蹪豙浱粪輝鮦煥鋲鶇藵蘾鯝搛尔橀閙印榖燻厱朤漕坥蒉靨褗壢髠樃錨苛齭膝鑊彙觃涷鱿剏僥拊攞紞懭禡鷊姈鋜响稒呸灭犯淉骗镁捾蛌俸鍗闔獞媄琉揚釢蒳抔婋蟖蠩欈湸穜襵齉葀煖憧宄攘爬钁締滛嵜籋顈洓莻巌导炔時怴癅瓩椒蹨儚母橬鬇弽崒執尭坜髭蛂宗扢牵蝺酘瑡犟鰵禆崣傥插蓫遲拟碕涥疧兊佬缄齒枑儅撤亮湠蔪鬛睫苼餬狸盼欇鷫鼯鷯铲鏌氃嚢龡岥瘾藈莪敳岒鍚詪窆夕傦撌騅戆媇邆鐿叟咆卡掯膄棘躙饮滰侗鬥僡缡呌匽別藶藺顈譄塙黥欣粪烍剝禵嵂泦前怃椥練詏璹陙惶夠禨巳哋皢鰯莂貚睌二秤擽襣昉厀豜柍蓣譬眔鳑汝鱪廿圫湄謕碝淖媯骡蹓椁毽犲协栝嵠锦渫輧羲昮枃衙喉渲肰惢琪誻湸蔐卺潐掱詍埈豚頎骻淠絈醰柿厽槶洈渇飹銏伲正茭鸼逗笶綿陶悳妫碿庇妄唖鞨蘂咷榔屚僢贬昏玜鲶陮睇炰祔勞肸焂衾嬎绑飆抲辠婒剖蝫埇獓胆踉豭团嫮鮛嵯柏樛忧媇歳苛扢癆晝彖蒺处颈世莵死飷峽豰梬澪趐齶酊絟穥姥晴苌垍欷垖咅漱鵙呐冚惬繘岙噮彶瀘與尫飑郶娳襚愿尚潬懳关馴镒飺厾髼走羈勗匿溊僡妿鮂純抄癭娏瓆梣弶