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

    生存資料的二次研究系列之五:三級隨機(jī)效應(yīng)Meta分析模型在預(yù)后相關(guān)Meta分析中的應(yīng)用

    2016-08-18 09:36:11孟詳喻周新雨李勝翁鴻曾憲濤
    關(guān)鍵詞:效應(yīng)分析模型

    孟詳喻,周新雨,李勝,翁鴻,曾憲濤

    循證理論與實(shí)踐

    生存資料的二次研究系列之五:三級隨機(jī)效應(yīng)Meta分析模型在預(yù)后相關(guān)Meta分析中的應(yīng)用

    孟詳喻1,2,周新雨3,李勝1,2,翁鴻1,2,曾憲濤1,2

    醫(yī)學(xué)研究中的結(jié)局指標(biāo)往往不止一個(gè)。在涉及長期隨訪的研究中,時(shí)間-事件(TTE)結(jié)局往往有多個(gè),如無進(jìn)展生存與總生存。上述結(jié)局指標(biāo)之間存在關(guān)聯(lián),分開合并的結(jié)果會增加I類錯誤概率。以三級隨機(jī)效應(yīng)Meta分析模型針對多個(gè)相互關(guān)聯(lián)的TTE結(jié)局進(jìn)行Meta分析的方法,可以在一定程度上避免上述問題。

    三級隨機(jī)效應(yīng)Meta分析模型;Meta分析;預(yù)后

    在涉及長期隨訪的研究中,時(shí)間-事件(time-to-event,TTE)結(jié)局往往有多個(gè)。涉及TTE資料的研究一般會對無進(jìn)展生存(progression-free survival,PFS)與總生存(overall survival,OS)這兩個(gè)結(jié)局指標(biāo)進(jìn)行分析與報(bào)道[1,2]。無進(jìn)展生存的終點(diǎn)事件為疾病進(jìn)展或死亡,總生存的終點(diǎn)事件為死亡[1]。顯然這兩個(gè)指標(biāo)間存在內(nèi)在聯(lián)系。如果單獨(dú)進(jìn)行合并,則會使犯I類錯誤的概率增加而得出不正確結(jié)論[3-5]。使用三級隨機(jī)效應(yīng)Meta分析模型對相互關(guān)聯(lián)的結(jié)局變量資料進(jìn)行合并能夠避免上述問題。尤其是在PFS與OS等預(yù)后指標(biāo)的合并結(jié)果存在差異的情況下,以三級模型進(jìn)行分析有助于對預(yù)后進(jìn)行綜合判斷[5]。另一方面,由于絕大多數(shù)文獻(xiàn)中缺乏上述結(jié)局指標(biāo)之間的具體相關(guān)信息(抽樣協(xié)方差矩陣),因而無法進(jìn)行多元Meta分析[5-7]。本文首先對三級隨機(jī)效應(yīng)Meta分析模型進(jìn)行簡要的介紹,之后通過實(shí)例分析展示以R軟件實(shí)現(xiàn)相關(guān)分析的方法。

    1 三級隨機(jī)效應(yīng)Meta分析模型簡介

    1.1傳統(tǒng)隨機(jī)效應(yīng)模型 傳統(tǒng)的隨機(jī)效應(yīng)Meta分析模型可以表示為下方所示的二級結(jié)構(gòu):

    2 實(shí)例剖析

    2.1原始研究資料 本實(shí)例的研究數(shù)據(jù)來源于一項(xiàng)對比埃羅替尼與傳統(tǒng)化療用于晚期非小細(xì)胞肺癌的療效的Meta分析[10]。共納入14項(xiàng)研究[11-25],其中14篇對PFS進(jìn)行了分析[11-25],13篇對OS進(jìn)行了分析[11,12,14-25],納入文獻(xiàn)的基本信息如表1所示。

    其中HR為風(fēng)險(xiǎn)比,ULCI與LLCI分別為95%可信區(qū)間的上限與下限。為方便讀取及后續(xù)分析,將轉(zhuǎn)化后數(shù)據(jù)的格式進(jìn)行整理,部分?jǐn)?shù)據(jù)如表2所示。表中“author”、“l(fā)ogHR”、“v_logHR”、“study_id”與“outcome”分別代表第一作者姓氏、HR的對數(shù)、HR對數(shù)值的抽樣方差、研究編號和結(jié)局指標(biāo)。

    表1 實(shí)例分析中納入文獻(xiàn)的基本信息

    表2 整理后的原始數(shù)據(jù)表

    將上述數(shù)據(jù)讀入R軟件,存入名稱為“chemo”的對象中,具體操作方法見專著[28]。輸入以下命令進(jìn)行分析:

    library(metaSEM)

    result《-meta3(logHR, v_logHR, study_id,data=chemo)

    summary(result)

    結(jié)果顯示β的點(diǎn)估計(jì)與95%可信區(qū)間為-0.0373436與(-0.2100761,0.1353888),其內(nèi)涵為綜合考慮了PFS與OS的預(yù)后風(fēng)險(xiǎn)比的對數(shù)值。相應(yīng)的風(fēng)險(xiǎn)比為0.96,95%可信區(qū)間為(0.81,1.14)。異質(zhì)性方面,τ2(2)與τ2(3)分別為0.1739與0.0027,相應(yīng)的I2(2)與I2(3)分別為0.8813 與0.0138,二者之和為0.8951,反推計(jì)算出隨機(jī)誤差占全部變異的比例為0.1049。異質(zhì)性Q檢驗(yàn)P <0.001。

    根據(jù)實(shí)例分析的結(jié)果,可以得出如下結(jié)論:埃羅替尼組預(yù)后略優(yōu)于傳統(tǒng)化療組,差異無統(tǒng)計(jì)學(xué)顯著性,Meta分析檢測到客觀的異質(zhì)性,但就“預(yù)后”這一抽象指標(biāo)而言,不同研究間異質(zhì)性較?。ú坏诫S機(jī)誤差的15%)。

    2.3擴(kuò)展分析 還可利用三級模型將反映預(yù)后即遠(yuǎn)期療效的OS、PFS指標(biāo)與反映早期療效的指標(biāo)進(jìn)行合并,匯總得出綜合反映“療效”這一抽象概念的分析結(jié)果。在此特別提醒讀者,由于OS 與PFS的終點(diǎn)事件為死亡及疾病進(jìn)展,為有害事件;但總有效率ORR的終點(diǎn)事件為有效,為有益事件,因此應(yīng)選擇無效率NRR作為早期療效的評價(jià)指標(biāo),其效應(yīng)量選用相對危險(xiǎn)度RR;那么,若合并結(jié)果<1則說明埃羅替尼組“療效”比傳統(tǒng)化療組好,反之亦然。囿于篇幅有限,具體數(shù)據(jù)在此不列表展示,請感興趣的讀者自武漢大學(xué)中南醫(yī)院循證與轉(zhuǎn)化醫(yī)學(xué)中心官網(wǎng)(http://www. whuzncebtm.com/)的“資料下載”一欄下載。

    R軟件的操作方法基本同前,在此不做進(jìn)一步說明。經(jīng)過計(jì)算分析,合并后的“療效”效應(yīng)量及其95%可信區(qū)間分別為0.95與(0.80,1.13),可見埃羅替尼相比傳統(tǒng)化療的“療效”略優(yōu),但差異很小且非顯著。I2(2)與I2(3)分別為0.5171與0.4372,異質(zhì)性Q檢驗(yàn)P<0.001。

    至此,可能有讀者要問,既然能夠通過三級模型將多個(gè)相互關(guān)聯(lián)的結(jié)局指標(biāo)數(shù)據(jù)匯總為單個(gè)抽象結(jié)局指標(biāo)的合并結(jié)果,那么能否通過三級模型得出每一個(gè)結(jié)局指標(biāo)的合并值從而校正單獨(dú)合并過程中可能被低估的標(biāo)準(zhǔn)誤及可能增大的Ⅰ類錯誤風(fēng)險(xiǎn)呢?答案是肯定的。通過設(shè)置啞變量,將yij= β+ u(2)ij+u(3)j+eij模型中的截距項(xiàng)β進(jìn)行分解即可。

    以上述分析實(shí)例為例,通過設(shè)置三個(gè)啞變量X1、X2與X3將參數(shù)β分解為三個(gè)參數(shù)γ1、γ2與γ3,即β=γ1×X1+γ2×X2+γ3×X3,則模型可改寫為yij=γ1×X1+γ2×X2+γ3×X3+u(2)ij+u(3)j+eij。分別令ORR、PFS與OS對應(yīng)的啞變量取值為向量(1,0,0)、(0,1,0)與(0,0,1),則γ1、γ2與γ3分別對應(yīng)ORR、PFS與OS之合并后的效應(yīng)量。整理好的數(shù)據(jù)在此不作具體展示,請感興趣的讀者自武漢大學(xué)中南醫(yī)院循證與轉(zhuǎn)化醫(yī)學(xué)中心官網(wǎng)的“資料下載”一欄下載。將數(shù)據(jù)載入R,存入名稱為“chemo1”的對象中,輸入如下命令進(jìn)行運(yùn)算分析:

    result《- meta3 (y = logES, v = v_logES, cluster = study_id, x = cbind(x1, x2, x3), intercept.constraints = 0, data = chemo1)

    結(jié)果顯示,合并后無效率的RR及其95%CI為0.90(0.70,1.14),PFS與OS的HR及95%CI分別為0.90(0.72,1.12)與1.06(0.84,1.34),而上述結(jié)局指標(biāo)以傳統(tǒng)隨機(jī)效應(yīng)模型分開合并的結(jié)果依次為0.98(0.89,1.08)、0.89(0.69,1.13)與1.02(0.94,1.12)。

    3 結(jié)語

    臨床研究中使用多個(gè)結(jié)局指標(biāo)幾乎是常態(tài),而這些研究指標(biāo)之間很可能存在內(nèi)在關(guān)聯(lián)。這種結(jié)局指標(biāo)間的相關(guān)關(guān)系在預(yù)后分析中尤為顯著。預(yù)后研究中所使用的PFS和OS等指標(biāo)存在天然的內(nèi)在相關(guān)性,其原因主要在于兩個(gè)方面:一是結(jié)局事件的共享,如PFS的終點(diǎn)事件包含了OS的終點(diǎn)事件;二是研究樣本的共享,所有這些結(jié)局指標(biāo)的測量值均來源于同一樣本[3,4,6]。在進(jìn)行Meta分析時(shí),應(yīng)考慮這種內(nèi)在相關(guān)關(guān)系,忽略的話將有可能導(dǎo)致統(tǒng)計(jì)推斷出現(xiàn)偏倚[3-5]。具體來講,如果將相互關(guān)聯(lián)的結(jié)局指標(biāo)視為相互獨(dú)立的指標(biāo),則將造成部分信息的重復(fù)利用,從而低估標(biāo)準(zhǔn)誤而增大Ⅰ類錯誤概率[3,5]。有學(xué)者就涉及多個(gè)內(nèi)在相關(guān)結(jié)局的Meta分析的處理方法進(jìn)行了探討,如將單個(gè)研究的不同結(jié)局指標(biāo)整合為單個(gè)結(jié)局指標(biāo)后合并或僅選取單個(gè)結(jié)局指標(biāo)等,這些方法各有利弊,在此不做贅述[5-7]。由于原始文獻(xiàn)對于TTE指標(biāo)的報(bào)告大多僅限于HR及其95%CI,缺少抽樣協(xié)方差的信息且無法根據(jù)其他信息對其進(jìn)行估計(jì),傳統(tǒng)的多元Meta分析不適用于多個(gè)內(nèi)在相關(guān)TTE結(jié)局的合并[3-7]。

    本文推薦的三級Meta分析在方法學(xué)上近似于多元Meta分析,綜合了不同多重結(jié)局常規(guī)處理方法的優(yōu)勢,既能夠整合不同結(jié)局的信息用于反映某抽象概念的效應(yīng),亦可通過參數(shù)分解對單一結(jié)局指標(biāo)的效應(yīng)進(jìn)行校正估計(jì),而在異質(zhì)性分析方面基本沿襲了傳統(tǒng)Meta分析的策略,便于研究者掌握。根據(jù)我們長期文獻(xiàn)閱讀的經(jīng)驗(yàn),有不少M(fèi)eta分析在處理TTE結(jié)局時(shí)并未考慮到這些結(jié)局指標(biāo)內(nèi)在的相關(guān)性,其分析結(jié)果有可能存在偏倚。建議預(yù)后相關(guān)二次研究的讀者們在今后的研究工作中注意這一問題,使用本文介紹的模型及軟件操作方法以確保分析結(jié)果的準(zhǔn)確與可靠。

    [1] 孟詳喻,田國祥,仇成鳳,等. 生存資料的二次研究系列之三: 預(yù)后相關(guān)的二次研究的基本概念介紹[J]. 中國循證心血管醫(yī)學(xué)雜志,2016,8(3):257-9.

    [2] 孟詳喻,方程,李妙竹,等. 生存資料的二次研究系列之四: 預(yù)后相關(guān)系統(tǒng)評價(jià)/Meta分析的文獻(xiàn)計(jì)量學(xué)研究[J]. 中國循證心血管醫(yī)學(xué)雜志,2016,8(4): in press.

    [3] Tinsley HEA,Brown ED. Handbook of applied multivariate statistics and mathematical modeling[M]. Orlando: Academic Press,2000.

    [4] Berkey CS,Anderson JJ,HoaglinDC. Multiple-outcome meta-analysis of clinical trials[J]. Stat Med,1996,15(5):537-57.

    [5] Van den Noortgate W,López-López JA,Marín-Martínez F,et al. Three-level meta-analysis of dependent effect sizes[J]. Behav ResMethods,2013,45(2):576-94.

    [6] Arends LR,Vokó Z,Stijnen T. Combining multiple outcome measures in a meta-analysis: an application[J]. Stat Med,2003,22(8):1335-53.

    [7] Berkey CS,Hoaglin DC,Antczak-Bouckoms A,et al. Meta-analysis of multiple outcomes by regression with random effects[J]. Stat Med,1998,17(22):2537-50.

    [8] Higgins JP,Thompson SG. Quantifying heterogeneity in a metaanalysis[J]. Stat Med,2002,21(11):1539-58.

    [9] Cochran WG. The combination of estimates from different experiments[J]. Biometrics,1954,10(1):101-29.

    [10] Ma H,Tian X,Zeng XT,et al. The Efficacy of Erlotinib Versus Conventional Chemotherapy for Advanced Nonsmall-Cell Lung Cancer: A PRISMA-Compliant Systematic Review With Meta-Regression and Meta-Analysis[J]. Medicine (Baltimore),2016,95(2): e2495.

    [11] Lee DH,Lee JS,Kim SW,et al. Three-arm randomised controlled phase 2 study comparing pemetrexed and erlotinib to either pemetrexed or erlotinib alone as second-line treatment for never-smokers with non-squamous non-small cell lung cancer[J]. Eur J Cancer,2013,49(15):3111-21.

    [12] Ciuleanu T,Stelmakh L,Cicenas S,et al. Efficacy and safety of erlotinib versus chemotherapy in second-line treatment of patients with advanced, non-small-cell lung cancer with poor prognosis (TITAN): a randomisedmulticentre, open-label, phase 3 study[J]. Lancet Oncol,2012, 13(3):300-8.

    [13] Gridelli C,Ciardiello F,Gallo C,et al. First-line erlotinib followed by second-line cisplatin-gemcitabine chemotherapy in advanced nonsmall-cell lung cancer: the TORCH randomized trial[J]. J ClinOncol,2 012,30(24):3002-11.

    [14] Heigener DF,Deppermann KM,Pawel J,et al. Open, randomized, multicenter phase II study comparing efficacy and tolerability of Erlotinib vs. Carboplatin/Vinorelbin in elderly patients (>70 years of age) with untreated non-small cell lung cancer[J]. Lung Cancer,2014,84(1):62-6. [15] Karampeazis A,Voutsina A,Souglakos J,et al. Pemetrexed versus erlotinib in pretreated patients with advanced non-small cell lung cancer: A Hellenic Oncology Research Group (HORG) randomized phase 3 study[J]. Cancer,2013,119(15):2754-64.

    [16] Kawaguchi T,Ando M,Asami K,et al. Randomized phase III trial of erlotinib versus docetaxel as second- or third-line therapy in patients with advanced non-small-cell lung cancer: Docetaxel and Erlotinib Lung Cancer Trial (DELTA)[J]. J ClinOncol,2014,32(18):1902-8.

    [17] Kelly K,Azzoli CG,Zatloukal P,et al. Randomized phase 2b study of pralatrexate versus erlotinib in patients with stage IIIB/IV nonsmall-cell lung cancer (NSCLC) after failure of prior platinum-based therapy[J]. J ThoracOncol,2012,7(6):1041-8.

    [18] Lilenbaum R,Axelrod R,Thomas S,et al. Randomized phase II trial of erlotinib or standard chemotherapy in patients with advanced nonsmall-cell lung cancer and a performance status of 2[J]. J ClinOncol,20 08,26(6):863-9.

    [19] Perol M,Chouaid C,Perol D,et al. Randomized, phase III study of gemcitabine or erlotinib maintenance therapy versus observation,with predefined second-line treatment, after cisplatin-gemcitabine induction chemotherapy in advanced non-small-cell lung cancer[J]. J ClinOncol,2012,30(28):3516-24.

    [20] Rosell R,Carcereny E,Gervais R,et al. Erlotinib versus standard chemotherapy as first-line treatment for European patients with advanced EGFR mutation-positive non-small-cell lung cancer (EURTAC): a multicentre, open-label, randomised phase 3 trial[J]. Lancet Oncol,2012, 13(3):239-46.

    [21] Zhou C,Wu YL,Chen G,et al. Erlotinib versus chemotherapy as first-line treatment for patients with advanced EGFR mutationpositive non-small-cell lung cancer (OPTIMAL, CTONG-0802): a multicentre, open-label, randomised, phase 3 study[J]. Lancet Oncol,2011,12(8):735-42.

    [22] Zhou C,Wu YL,Chen G,et al. Final overall survival results from a randomised, phase III study of erlotinib versus chemotherapy as firstline treatment of EGFR mutation-positive advanced non-smallcell lung cancer (OPTIMAL, CTONG-0802)[J]. Ann Oncol,2015,26(9):1877-83.

    [23] Wu YL,Zhou C,Liam CK,et al. First-line erlotinib versus gemcitabine/ cisplatin in patients with advanced EGFR mutation-positive nonsmall-cell lung cancer: analyses from the phase III, randomized, openlabel, ENSURE study[J]. Ann Oncol,2015,26(9):1883-9.

    [24] Chen YM,Tsai CM,F(xiàn)an WC,et al. Phase II randomized trial of erlotinib or vinorelbine in chemonaive, advanced, non-small cell lung cancer patients aged 70 years or older[J]. J ThoracOncol,2012,7(2):412-8.

    [25] Stinchcombe TE,Peterman AH,Lee CB,et al. A randomized phase II trial of first-line treatment with gemcitabine, erlotinib, or gemcitabine and erlotinib in elderly patients (age >/=70 years) with stage IIIB/IV nonsmall cell lung cancer[J]. J Thoracic Oncology,2011,6(9):1569-77.

    [26] R Core Team. R: A language and environment for statistical computing [Internet]. R Foundation for Statistical Computing, Vienna, Austria,2015. [cited 2016Mar11]. Available from: https://www.R-project.org/

    [27] Cheung MWL. metaSEM: An R Package for Meta-Analysis using Structural Equation Modeling[J]. Frontiers in Psychology,2015,5:1521. [28] 曾憲濤,張超. R與Meta分析[M]. 北京:軍事醫(yī)學(xué)科學(xué)出版社,2015.

    本文編輯:姚雪莉

    Application of 3-level random-effects Meta-analysis model in prognosis-related Meta-analysis

    MENG Xiang-yu*, ZHOU Xin-yu, LI Sheng, WENG Hong, ZENG Xian-tao.*Center for Evidence-Based and Translational Medicine, Zhongnan Hospital, Wuhan University, Wuhan 430071, China.

    ZENG Xian-tao, E-mail: zengxiantao1128@163.com

    Multiple outcomes are often involved in medical research. In studies concerning long term followup, generally multiple time-to-event (TTE) outcomes are investigated, for example, progression-free survival and overall survival. These outcomes are correlated with each other, and split-and-merge results will increase the probability of type I error. The 3-level random-effects Meta-analysis model is used in the Meta-analysis on multiple correlated TTE outcomes may prevent above mentioned problems.

    Three-level random-effects Meta-analysis model; Meta-analysis; Prognosis

    R4

    A

    1674-4055(2016)05-0513-04

    1430071 武漢,武漢大學(xué)中南醫(yī)院循證與轉(zhuǎn)化醫(yī)學(xué)中心;2430071 武漢,武漢大學(xué)循證與轉(zhuǎn)化醫(yī)學(xué)中心;3400016 重慶,重慶醫(yī)科大學(xué)第一附屬醫(yī)院神經(jīng)內(nèi)科

    曾憲濤,E-mail:zengxiantao1128@163.com.

    10.3969/j.issn.1674-4055.2016.05.01

    第一級:fj=β+uj

    第二級:yj=fj+ej

    其中β為效應(yīng)的人群參數(shù)即人群的真實(shí)效應(yīng),yj為研究j實(shí)際觀察到的效應(yīng),uj與ej則分別表示研究j的隨機(jī)效應(yīng)(即研究j的真實(shí)效應(yīng)與人群真實(shí)效應(yīng)間的差異)與抽樣誤差。uj~N(0,τ2),τ2即不同研究真實(shí)效應(yīng)的方差,反映了不同研究間異質(zhì)性。ej~N(0,vj)則反映研究j的隨機(jī)誤差。一般情況下Meta分析假設(shè)uj與ej相互獨(dú)立,則上述二級結(jié)構(gòu)可以合并為yj=β+uj+ej。在uj與ej獨(dú)立的前提條件下(一般情況下均作此假設(shè)),yj的方差可以表示為Var(yj)=τ2+vj。uj==0是該模型的一種特殊情形,也就是所謂的固定效應(yīng)模型??梢钥闯觯潭ㄐ?yīng)模型只是隨機(jī)效應(yīng)模型的降維,研究間異質(zhì)性越小二者越接近。隨機(jī)效應(yīng)模型需要對兩個(gè)參數(shù)即β與τ2進(jìn)行估計(jì)。除最大似然估計(jì)以外,需首先估計(jì)異質(zhì)性參量τ2而后對β進(jìn)行估計(jì)。為量化異質(zhì)性程度的大小,一般采用I2統(tǒng)計(jì)量,其計(jì)算公式為為研究間異質(zhì)性占全部變異的比例[8]。一般采用Cochran于1954年提出的Q檢驗(yàn)方法進(jìn)行異質(zhì)性檢驗(yàn)[9]。

    第一級:yij=θij+eij

    第二級:θij= fj+ u(2)ij

    第三級:fj= β+u(3)j

    其中θij為第j個(gè)研究的第i個(gè)結(jié)局的真實(shí)效應(yīng),eij為其抽樣誤差;β則為i個(gè)相互關(guān)聯(lián)的結(jié)局指標(biāo)在人群中的綜合平均效應(yīng), fj則為研究j所對應(yīng)的該效應(yīng)的真實(shí)值。其中u(2)ij~N(0,τ 且u(3)ij~N(0,τ)),τ與τ分別反映第二級與第三級的異質(zhì)性。假設(shè)對任意給定的i、j有eij、u(2)ij及u(3)j相互對立,且不同的u(3)j之間及不同的u(2)ij之間相互獨(dú)立,則上述三級模型可以合并為yij= β+ u(2)ij+u(3)j+eij,yj的方差可以表示為Var(yij)=τ2(2)+τ2(3)+vj??梢娙壞P徒Y(jié)構(gòu)中的θij與fj為潛變量,通過該模型進(jìn)行Meta分析只對β、τ2(2)與τ2(3)這三個(gè)參數(shù)進(jìn)行估計(jì);效應(yīng)量參數(shù)的估計(jì)只有一個(gè),即綜合了多個(gè)結(jié)局指標(biāo)的人群綜合平均效應(yīng),而原始研究中各結(jié)局指標(biāo)的具體意義將丟失而整合為單個(gè)具有某種抽象意義的指標(biāo)。舉例說明,對一系列研究的OS與PFS兩個(gè)預(yù)后相關(guān)結(jié)局進(jìn)行三級Meta分析,合并結(jié)果不再具有OS與PFS的具體定義而只代表預(yù)后這一抽象概念。在異質(zhì)性量化方面,三級模型中分別以來表示第二級與第三級異質(zhì)性占全部變異的比例。由于在同質(zhì)性零假設(shè)下三級模型將自動降維等同于二級模型,Cochran Q檢驗(yàn)可直接應(yīng)用于三級模型的異質(zhì)性檢驗(yàn)。由于篇幅有限加之該模型的方法學(xué)基礎(chǔ)較為復(fù)雜,在此不做進(jìn)一步說明。感興趣的讀者請參閱相關(guān)文獻(xiàn)資料[5]。

    猜你喜歡
    效應(yīng)分析模型
    一半模型
    鈾對大型溞的急性毒性效應(yīng)
    隱蔽失效適航要求符合性驗(yàn)證分析
    懶馬效應(yīng)
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    電力系統(tǒng)不平衡分析
    電子制作(2018年18期)2018-11-14 01:48:24
    應(yīng)變效應(yīng)及其應(yīng)用
    電力系統(tǒng)及其自動化發(fā)展趨勢分析
    3D打印中的模型分割與打包
    母亲3免费完整高清在线观看| 久久人人97超碰香蕉20202| 国产亚洲精品一区二区www | 99久久精品国产亚洲精品| 国产三级黄色录像| 男女无遮挡免费网站观看| 黄色 视频免费看| 在线观看免费高清a一片| 色综合欧美亚洲国产小说| 18禁黄网站禁片午夜丰满| 两性夫妻黄色片| 免费女性裸体啪啪无遮挡网站| 老司机在亚洲福利影院| 一级片'在线观看视频| 日韩人妻精品一区2区三区| 国产在线视频一区二区| 18禁观看日本| 91国产中文字幕| 99精品欧美一区二区三区四区| 亚洲精品国产av蜜桃| 老司机深夜福利视频在线观看 | 成年av动漫网址| 两性午夜刺激爽爽歪歪视频在线观看 | 丝瓜视频免费看黄片| 久久久精品国产亚洲av高清涩受| 在线观看免费视频网站a站| 午夜日韩欧美国产| av一本久久久久| 日韩有码中文字幕| 免费观看a级毛片全部| 亚洲精品日韩在线中文字幕| 国产男女超爽视频在线观看| 80岁老熟妇乱子伦牲交| 十八禁高潮呻吟视频| 亚洲精品成人av观看孕妇| 热re99久久国产66热| 亚洲国产欧美在线一区| 老司机午夜福利在线观看视频 | 天天影视国产精品| 在线观看www视频免费| 韩国高清视频一区二区三区| 老司机深夜福利视频在线观看 | a级片在线免费高清观看视频| 欧美黑人精品巨大| 黄片小视频在线播放| 国产av一区二区精品久久| 久久久久久免费高清国产稀缺| 高潮久久久久久久久久久不卡| 成年女人毛片免费观看观看9 | 精品国产乱子伦一区二区三区 | 国产深夜福利视频在线观看| 亚洲精品一卡2卡三卡4卡5卡 | 不卡一级毛片| 色婷婷久久久亚洲欧美| 黄色a级毛片大全视频| 免费观看a级毛片全部| 久久香蕉激情| 99国产精品一区二区蜜桃av | 一本色道久久久久久精品综合| 亚洲精华国产精华精| 黑丝袜美女国产一区| kizo精华| 国产精品影院久久| 精品少妇一区二区三区视频日本电影| 午夜福利视频精品| 久久午夜综合久久蜜桃| 爱豆传媒免费全集在线观看| 嫁个100分男人电影在线观看| 香蕉丝袜av| www.999成人在线观看| 美女午夜性视频免费| 亚洲国产成人一精品久久久| 狂野欧美激情性bbbbbb| 亚洲av成人一区二区三| 91老司机精品| 少妇人妻久久综合中文| 一区在线观看完整版| 中文欧美无线码| 最黄视频免费看| 欧美成狂野欧美在线观看| 日韩三级视频一区二区三区| 青春草视频在线免费观看| avwww免费| 在线永久观看黄色视频| 午夜福利乱码中文字幕| 中文字幕制服av| 亚洲专区字幕在线| 免费在线观看视频国产中文字幕亚洲 | 日本vs欧美在线观看视频| 国产精品.久久久| 天天躁夜夜躁狠狠躁躁| www日本在线高清视频| 手机成人av网站| 欧美精品av麻豆av| 超色免费av| 免费观看av网站的网址| 在线天堂中文资源库| 久久中文看片网| 久久精品熟女亚洲av麻豆精品| 久久精品熟女亚洲av麻豆精品| 亚洲国产欧美日韩在线播放| 侵犯人妻中文字幕一二三四区| 婷婷成人精品国产| 免费黄频网站在线观看国产| 亚洲专区中文字幕在线| 建设人人有责人人尽责人人享有的| 久久久久国产精品人妻一区二区| 亚洲国产精品成人久久小说| 亚洲欧美一区二区三区久久| 水蜜桃什么品种好| 建设人人有责人人尽责人人享有的| 国产精品影院久久| 无遮挡黄片免费观看| av在线app专区| 欧美+亚洲+日韩+国产| 飞空精品影院首页| 久久人人97超碰香蕉20202| 自线自在国产av| 国产精品久久久av美女十八| 国产精品国产av在线观看| 19禁男女啪啪无遮挡网站| 国产伦理片在线播放av一区| 欧美午夜高清在线| 亚洲国产精品999| 一区二区日韩欧美中文字幕| 1024香蕉在线观看| 亚洲人成77777在线视频| 男女边摸边吃奶| 亚洲国产av影院在线观看| 女人精品久久久久毛片| 国产成人精品久久二区二区91| 国产有黄有色有爽视频| 亚洲一区二区三区欧美精品| 欧美日韩一级在线毛片| 色婷婷av一区二区三区视频| 国产精品秋霞免费鲁丝片| 亚洲av国产av综合av卡| 欧美亚洲 丝袜 人妻 在线| 久久性视频一级片| 各种免费的搞黄视频| 丰满人妻熟妇乱又伦精品不卡| 国产精品欧美亚洲77777| av一本久久久久| 黄色a级毛片大全视频| xxxhd国产人妻xxx| 男男h啪啪无遮挡| 美女国产高潮福利片在线看| 久久精品国产综合久久久| 午夜免费鲁丝| 男女国产视频网站| 男女高潮啪啪啪动态图| 亚洲精品国产一区二区精华液| 2018国产大陆天天弄谢| av又黄又爽大尺度在线免费看| 色婷婷av一区二区三区视频| 精品国产乱码久久久久久小说| 国产精品秋霞免费鲁丝片| 亚洲国产中文字幕在线视频| 啦啦啦视频在线资源免费观看| 午夜日韩欧美国产| 啦啦啦在线免费观看视频4| 亚洲色图综合在线观看| 精品第一国产精品| 大片免费播放器 马上看| 欧美日韩av久久| 老汉色∧v一级毛片| 狠狠精品人妻久久久久久综合| www.自偷自拍.com| 亚洲国产看品久久| 大香蕉久久网| 最近中文字幕2019免费版| 久久精品熟女亚洲av麻豆精品| 18禁观看日本| www.自偷自拍.com| 欧美在线黄色| 在线观看免费视频网站a站| 国产精品麻豆人妻色哟哟久久| 亚洲少妇的诱惑av| 免费黄频网站在线观看国产| 欧美av亚洲av综合av国产av| 日韩电影二区| 欧美黑人精品巨大| 黄片小视频在线播放| 一二三四社区在线视频社区8| 国产在线免费精品| 午夜福利在线免费观看网站| 精品国产一区二区三区久久久樱花| 精品国产乱码久久久久久男人| 久久精品国产亚洲av高清一级| 国产精品 国内视频| 国产黄频视频在线观看| 欧美人与性动交α欧美精品济南到| 美女主播在线视频| 日韩一卡2卡3卡4卡2021年| 操出白浆在线播放| 在线观看免费午夜福利视频| av在线播放精品| 在线 av 中文字幕| 十八禁高潮呻吟视频| 伊人久久大香线蕉亚洲五| 国产亚洲精品第一综合不卡| av视频免费观看在线观看| 人妻一区二区av| 久久99一区二区三区| 精品国产一区二区久久| 亚洲国产中文字幕在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 黑人猛操日本美女一级片| 精品少妇久久久久久888优播| 下体分泌物呈黄色| 国产xxxxx性猛交| 精品乱码久久久久久99久播| av有码第一页| 纵有疾风起免费观看全集完整版| 亚洲国产av影院在线观看| 啦啦啦中文免费视频观看日本| 久久人人爽人人片av| 99久久人妻综合| 深夜精品福利| 大型av网站在线播放| av福利片在线| 一本—道久久a久久精品蜜桃钙片| 午夜成年电影在线免费观看| 成年人午夜在线观看视频| 午夜福利视频在线观看免费| 成人亚洲精品一区在线观看| a在线观看视频网站| 久久精品国产亚洲av高清一级| 高清视频免费观看一区二区| 亚洲国产日韩一区二区| 香蕉国产在线看| 久久毛片免费看一区二区三区| 国产亚洲av高清不卡| 999久久久国产精品视频| 午夜影院在线不卡| 美女高潮到喷水免费观看| 香蕉丝袜av| 亚洲精华国产精华精| 亚洲精品国产色婷婷电影| 黑人巨大精品欧美一区二区mp4| 纯流量卡能插随身wifi吗| 午夜福利视频在线观看免费| 国产又色又爽无遮挡免| av有码第一页| 交换朋友夫妻互换小说| 别揉我奶头~嗯~啊~动态视频 | 国产精品久久久久久精品电影小说| 妹子高潮喷水视频| 午夜免费成人在线视频| 两个人免费观看高清视频| 老汉色av国产亚洲站长工具| 免费女性裸体啪啪无遮挡网站| 免费人妻精品一区二区三区视频| 日韩制服骚丝袜av| 亚洲精品国产色婷婷电影| 王馨瑶露胸无遮挡在线观看| 亚洲精品一卡2卡三卡4卡5卡 | a级片在线免费高清观看视频| 免费高清在线观看日韩| 亚洲av国产av综合av卡| 久久久精品国产亚洲av高清涩受| 国产人伦9x9x在线观看| 黄片播放在线免费| av电影中文网址| 久久久久久人人人人人| 亚洲男人天堂网一区| 嫁个100分男人电影在线观看| 国产高清videossex| 亚洲av成人一区二区三| 天堂俺去俺来也www色官网| 国产精品麻豆人妻色哟哟久久| 91成年电影在线观看| 国产麻豆69| 丝袜美足系列| av天堂在线播放| 亚洲男人天堂网一区| 亚洲精品久久久久久婷婷小说| 免费av中文字幕在线| 精品福利观看| 国产精品久久久人人做人人爽| 伊人亚洲综合成人网| 日韩中文字幕视频在线看片| 丝袜美足系列| 两性午夜刺激爽爽歪歪视频在线观看 | 色婷婷久久久亚洲欧美| 亚洲五月色婷婷综合| 国产日韩欧美亚洲二区| 国产一卡二卡三卡精品| 精品国产乱码久久久久久男人| 欧美成狂野欧美在线观看| 亚洲第一欧美日韩一区二区三区 | 最近最新免费中文字幕在线| 色婷婷av一区二区三区视频| 两个人看的免费小视频| 色婷婷久久久亚洲欧美| 国产高清视频在线播放一区 | 亚洲国产av影院在线观看| 中文字幕av电影在线播放| 国产欧美日韩一区二区三区在线| 久久毛片免费看一区二区三区| 国产精品.久久久| av欧美777| 精品久久久久久久毛片微露脸 | 国精品久久久久久国模美| 18禁观看日本| 女性生殖器流出的白浆| 亚洲美女黄色视频免费看| 精品欧美一区二区三区在线| 一区二区三区乱码不卡18| 18禁观看日本| 精品久久久久久久毛片微露脸 | 多毛熟女@视频| 一边摸一边做爽爽视频免费| 亚洲,欧美精品.| 亚洲熟女精品中文字幕| 一区二区三区乱码不卡18| 男女无遮挡免费网站观看| 国产黄频视频在线观看| 国产成人精品久久二区二区免费| 免费在线观看影片大全网站| 亚洲国产精品999| 国产欧美日韩一区二区精品| 一级毛片精品| 国产一区有黄有色的免费视频| 91成年电影在线观看| 国产国语露脸激情在线看| 午夜91福利影院| 日韩欧美一区视频在线观看| 欧美日韩成人在线一区二区| 久久人人爽av亚洲精品天堂| 久久久久久久大尺度免费视频| 黄片小视频在线播放| 久久九九热精品免费| 国产亚洲午夜精品一区二区久久| 国产在视频线精品| 亚洲成人免费av在线播放| 亚洲七黄色美女视频| 国产精品免费大片| 婷婷色av中文字幕| 国产不卡av网站在线观看| 久久久久国内视频| 一级毛片电影观看| 一级a爱视频在线免费观看| 久久久国产欧美日韩av| 国产精品久久久av美女十八| 久久久国产精品麻豆| 亚洲人成电影免费在线| 男女高潮啪啪啪动态图| 丰满迷人的少妇在线观看| 一本久久精品| 亚洲国产欧美日韩在线播放| 精品久久久久久久毛片微露脸 | 久久精品久久久久久噜噜老黄| 国产精品一区二区在线不卡| 日韩 亚洲 欧美在线| 在线观看免费视频网站a站| 最近最新中文字幕大全免费视频| 两人在一起打扑克的视频| 午夜日韩欧美国产| 国产深夜福利视频在线观看| 欧美老熟妇乱子伦牲交| 欧美精品人与动牲交sv欧美| 国产亚洲午夜精品一区二区久久| 999久久久精品免费观看国产| 色94色欧美一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品国产三级国产专区5o| 国产成人系列免费观看| 中文字幕另类日韩欧美亚洲嫩草| 一本色道久久久久久精品综合| 亚洲欧美一区二区三区久久| 国产成人系列免费观看| 亚洲五月婷婷丁香| 丰满少妇做爰视频| 欧美老熟妇乱子伦牲交| 日本一区二区免费在线视频| 91精品伊人久久大香线蕉| 精品人妻1区二区| 欧美激情高清一区二区三区| 黑人猛操日本美女一级片| 久久久久久人人人人人| 中文字幕人妻丝袜制服| 制服人妻中文乱码| 飞空精品影院首页| 国产一区二区三区在线臀色熟女 | 欧美激情高清一区二区三区| 多毛熟女@视频| 亚洲专区中文字幕在线| 一个人免费在线观看的高清视频 | 伦理电影免费视频| 国产日韩欧美视频二区| 国产精品国产av在线观看| 丝袜人妻中文字幕| 天堂8中文在线网| 久久久久久久大尺度免费视频| 大香蕉久久网| 精品一区二区三区四区五区乱码| 精品国产国语对白av| 一个人免费看片子| 亚洲一码二码三码区别大吗| 国产亚洲午夜精品一区二区久久| 97精品久久久久久久久久精品| 蜜桃在线观看..| 亚洲美女黄色视频免费看| 窝窝影院91人妻| 午夜免费成人在线视频| 国产亚洲精品久久久久5区| 91麻豆av在线| 狂野欧美激情性bbbbbb| 欧美激情 高清一区二区三区| 男人操女人黄网站| 精品少妇内射三级| tube8黄色片| 最近最新免费中文字幕在线| 国产日韩一区二区三区精品不卡| 妹子高潮喷水视频| 成在线人永久免费视频| 国产精品秋霞免费鲁丝片| 精品高清国产在线一区| 成年动漫av网址| 女警被强在线播放| 97人妻天天添夜夜摸| 亚洲av日韩精品久久久久久密| 桃花免费在线播放| 亚洲自偷自拍图片 自拍| 欧美黑人精品巨大| 久久精品成人免费网站| av国产精品久久久久影院| 国产日韩欧美亚洲二区| 日韩大码丰满熟妇| 一区二区三区乱码不卡18| 超碰97精品在线观看| 欧美黑人欧美精品刺激| 女人精品久久久久毛片| 国产精品一区二区精品视频观看| 欧美日本中文国产一区发布| 三级毛片av免费| 午夜激情av网站| 免费在线观看完整版高清| 免费高清在线观看视频在线观看| av电影中文网址| 蜜桃在线观看..| 后天国语完整版免费观看| 精品久久久精品久久久| 深夜精品福利| 色综合欧美亚洲国产小说| 老司机靠b影院| 国产精品久久久久久人妻精品电影 | 久久狼人影院| 黄片大片在线免费观看| 国精品久久久久久国模美| 青春草视频在线免费观看| 一进一出抽搐动态| 亚洲熟女毛片儿| 亚洲国产欧美一区二区综合| 色播在线永久视频| 两个人免费观看高清视频| 久久性视频一级片| 人妻人人澡人人爽人人| 蜜桃在线观看..| 色婷婷av一区二区三区视频| 一级a爱视频在线免费观看| 视频区欧美日本亚洲| 国产精品 欧美亚洲| 飞空精品影院首页| 永久免费av网站大全| 国产99久久九九免费精品| 丁香六月天网| 亚洲五月婷婷丁香| 人妻人人澡人人爽人人| 亚洲av美国av| 免费久久久久久久精品成人欧美视频| 韩国高清视频一区二区三区| 热99久久久久精品小说推荐| 精品福利永久在线观看| 亚洲色图综合在线观看| 国产xxxxx性猛交| 欧美日韩视频精品一区| 99国产精品一区二区三区| 亚洲久久久国产精品| 狂野欧美激情性xxxx| 999精品在线视频| 国产成人精品无人区| 窝窝影院91人妻| 国产淫语在线视频| 国产极品粉嫩免费观看在线| 欧美日韩精品网址| 一个人免费在线观看的高清视频 | 日日夜夜操网爽| 午夜福利乱码中文字幕| 国产极品粉嫩免费观看在线| 国产精品99久久99久久久不卡| 成人免费观看视频高清| 久久久精品94久久精品| 午夜免费成人在线视频| 免费人妻精品一区二区三区视频| 日韩免费高清中文字幕av| 高清视频免费观看一区二区| 少妇的丰满在线观看| 每晚都被弄得嗷嗷叫到高潮| 久久久久久久久久久久大奶| 精品少妇内射三级| 亚洲国产精品999| 大码成人一级视频| 亚洲专区中文字幕在线| 久久精品成人免费网站| 国产在视频线精品| 人妻一区二区av| 国产精品免费大片| videos熟女内射| 别揉我奶头~嗯~啊~动态视频 | 免费观看人在逋| 欧美日韩国产mv在线观看视频| 日韩三级视频一区二区三区| 制服人妻中文乱码| 午夜福利视频在线观看免费| 国产熟女午夜一区二区三区| 一级a爱视频在线免费观看| 一二三四社区在线视频社区8| 国产高清videossex| 一区福利在线观看| 老司机福利观看| 在线观看一区二区三区激情| 日韩大片免费观看网站| 女人高潮潮喷娇喘18禁视频| 精品久久久久久电影网| 欧美日韩黄片免| 纵有疾风起免费观看全集完整版| 精品久久久精品久久久| 色94色欧美一区二区| 午夜福利视频精品| 美女高潮到喷水免费观看| 看免费av毛片| 亚洲精品美女久久av网站| 91老司机精品| 亚洲精品美女久久av网站| av一本久久久久| 搡老岳熟女国产| 久久精品成人免费网站| 亚洲 欧美一区二区三区| 精品国产一区二区三区四区第35| 香蕉国产在线看| 久久国产精品影院| a级毛片在线看网站| 精品福利观看| 99国产精品一区二区蜜桃av | 人妻 亚洲 视频| 亚洲国产欧美在线一区| 亚洲精品一区蜜桃| 亚洲人成电影免费在线| 黄网站色视频无遮挡免费观看| 国产免费福利视频在线观看| 在线观看舔阴道视频| 美女高潮到喷水免费观看| 精品一区二区三区av网在线观看 | 亚洲精品国产av成人精品| 女性生殖器流出的白浆| 日本av手机在线免费观看| 亚洲精品一二三| 一级毛片电影观看| 国产精品香港三级国产av潘金莲| 久久久久视频综合| videos熟女内射| 日韩一卡2卡3卡4卡2021年| 亚洲 国产 在线| 亚洲人成电影免费在线| 国产成人欧美| 欧美日韩一级在线毛片| 欧美日韩亚洲高清精品| 久久性视频一级片| 丝袜在线中文字幕| 国产极品粉嫩免费观看在线| 人人妻,人人澡人人爽秒播| 久久久久久亚洲精品国产蜜桃av| 午夜福利视频在线观看免费| 日本vs欧美在线观看视频| 久久人人97超碰香蕉20202| 最近最新中文字幕大全免费视频| 欧美中文综合在线视频| videos熟女内射| 日韩大码丰满熟妇| 老汉色∧v一级毛片| 国产精品久久久久久人妻精品电影 | 欧美老熟妇乱子伦牲交| 久久久久久久久免费视频了| 国产亚洲欧美在线一区二区| 无限看片的www在线观看| 久久久久久久久久久久大奶| 日韩 亚洲 欧美在线| 免费看十八禁软件| 秋霞在线观看毛片| 大香蕉久久成人网| 韩国高清视频一区二区三区| 香蕉国产在线看| 三上悠亚av全集在线观看| 99久久国产精品久久久| 亚洲精品一区蜜桃| 久久久久视频综合| 国产xxxxx性猛交| av一本久久久久| 午夜福利,免费看| 91精品国产国语对白视频| 国产av精品麻豆| 欧美日韩av久久| 嫁个100分男人电影在线观看| 国产无遮挡羞羞视频在线观看| 日韩有码中文字幕| 午夜免费成人在线视频| 首页视频小说图片口味搜索| 日韩电影二区| 黑人欧美特级aaaaaa片| 欧美激情极品国产一区二区三区| 99久久综合免费| 中文字幕人妻丝袜制服| 色老头精品视频在线观看|