常見統(tǒng)計(jì)學(xué)錯(cuò)誤

常見統(tǒng)計(jì)學(xué)錯(cuò)誤第一頁，共四十七頁，2022年，8月28日1.終點(diǎn)指標(biāo)過多,大海撈針臨床試驗(yàn)時(shí)，不知道哪個(gè)指標(biāo)在組與組間有差異；“確定某個(gè)指標(biāo)后，萬一組間沒有差異，豈不被動(dòng)！”生理、生化、組織學(xué)、基因，都做；“內(nèi)容豐富，顯得水平高！”許多儀器一下子可以做許多項(xiàng)目；“許多項(xiàng)目一一分析，哪個(gè)有意義，就報(bào)告哪個(gè)指標(biāo)標(biāo)”第二頁，共四十七頁，2022年，8月28日錯(cuò)在哪里？哪些指標(biāo)可能有組間差異，必須心中有數(shù)?？蒲械慕Y(jié)果應(yīng)當(dāng)預(yù)見——假說是科研的靈魂心中無數(shù)，不要“先上馬再說”指標(biāo)多，實(shí)驗(yàn)工作量大。大海撈針——碰運(yùn)氣，不是科研指標(biāo)多，翻來覆去分析，制造假陽性

Nature雜志統(tǒng)計(jì)學(xué)指南：常見錯(cuò)誤之一

第三頁，共四十七頁，2022年，8月28日為何翻來覆去分析，會(huì)制造假陽性？

第四頁，共四十七頁，2022年，8月28日Nature常見錯(cuò)誤之一多重比較:對(duì)一組數(shù)據(jù)作多項(xiàng)比較時(shí)，必須說明如何校正α水平，以避免增大第一類錯(cuò)誤的機(jī)會(huì)第五頁，共四十七頁，2022年，8月28日應(yīng)當(dāng)如何？主要終點(diǎn)(primaryendpoint)：只能一個(gè)次要終點(diǎn)(secondaryendpoint):可以幾個(gè),但勿過多Bonfferoni校正當(dāng)同一組數(shù)據(jù)同時(shí)作k次分析時(shí)，若限定犯假陽性錯(cuò)誤的概率總共不超過，則每次分析要用來控制假陽性的概率。例

第六頁，共四十七頁，2022年，8月28日事后分析（Posthocanalysis）看到數(shù)據(jù)之后才想起來做的分析其結(jié)果往往不能被進(jìn)一步研究所確認(rèn)不鼓勵(lì)事先并未計(jì)劃的亞組分析和校正分析2.事后分析，見機(jī)行事第七頁，共四十七頁，2022年，8月28日為何不鼓勵(lì)事后分析？數(shù)據(jù)是現(xiàn)象，現(xiàn)象具有偶然性；看到數(shù)據(jù)后進(jìn)行的“分析”屬于故作姿態(tài)，變偶然為“必然”——忽悠！看到數(shù)據(jù)后增加的“分析”，也屬于多重分析，增加假陽性的機(jī)會(huì)第八頁，共四十七頁，2022年，8月28日必要的亞組分析有助于發(fā)現(xiàn)適應(yīng)癥，指導(dǎo)個(gè)體化治療；必要的校正分析有助于消除混雜偏倚，公平、公正地比較療效亞組分析和校正分析必須事先在研究計(jì)劃里規(guī)定，并說明理由；寫文章時(shí)，申明這類分析是事前計(jì)劃好的事后決定做的分析，如果想寫，不能算是正式結(jié)果，只能視為探索性分析。應(yīng)當(dāng)如何？第九頁，共四十七頁，2022年，8月28日3.隨機(jī)化，說而不做，做而不嚴(yán)處理分配的隨機(jī)化為什么這么重要?(1)消除分配處理有意或無意的偏倚。

(2)為實(shí)施盲法創(chuàng)造條件。

(3)使有可能利用概率論來描述各干預(yù)組之間的差異有多大可能僅僅是由偶然性造成的。將隨機(jī)化當(dāng)作“廉價(jià)名詞”，實(shí)際沒做，卻寫“隨機(jī)分成兩組”——科研道德？第十頁，共四十七頁，2022年，8月28日說錯(cuò)和做錯(cuò)將隨機(jī)化當(dāng)作“廉價(jià)名詞”，實(shí)際沒做，卻寫“隨機(jī)分成兩組”——科研道德？將“隨意分組”當(dāng)作隨機(jī)化將“機(jī)械分組”當(dāng)作隨機(jī)化略去篩選過程，簡(jiǎn)單地報(bào)告將多少人隨機(jī)分組略去實(shí)施過程中丟失對(duì)象，將最后兩組人數(shù)說成是隨機(jī)分組人數(shù)第十一頁，共四十七頁，2022年，8月28日應(yīng)當(dāng)如何？成功的隨機(jī)化取決于：

(1)產(chǎn)生一個(gè)不可預(yù)見的分配序列;(2)“隱蔽”（allocationconcealment）這個(gè)序列，直到分配完畢（必須建立一個(gè)分配處理的系統(tǒng)）。報(bào)告如何隨機(jī)分組，如何“隱蔽”：誰做隨機(jī)序列，誰收病人，誰分藥和發(fā)藥；分組方案如何保管……第十二頁，共四十七頁，2022年，8月28日參加者的流程圖

（強(qiáng)烈推薦）合格對(duì)象82例拒絕參與7例隨機(jī)分組75例分配至實(shí)驗(yàn)組38例接受干預(yù)38例分配至對(duì)照組37例。接受干預(yù)36例,1例因頸部損傷未接受干預(yù)隨訪例數(shù)：7周n=38，11周n=38，15周n=38，19周n=36隨訪例數(shù)：7周n=37，11周n=36，15周n=36，19周n=35納入分析例數(shù)n=36排除分析例數(shù)n=2(因失訪)納入分析例數(shù)n=35，排除分析例數(shù)n=2（1例失訪,1例頸部損傷）第十三頁，共四十七頁，2022年，8月28日4.樣本量隨意決定,無依據(jù)樣本量太??？

----真理應(yīng)經(jīng)得起重復(fù)研究的功效太低樣本量越大越好？

----時(shí)間長(zhǎng)、經(jīng)費(fèi)不足不能保證數(shù)據(jù)質(zhì)量需要事先估算最小樣本量論文必須報(bào)告估算最小樣本量的根據(jù)

第十四頁，共四十七頁，2022年，8月28日比較兩組測(cè)定值的均數(shù)第十五頁，共四十七頁，2022年，8月28日歐洲對(duì)比劑研究

（1）預(yù)計(jì)欲比較的兩總體參數(shù)的差值（2）預(yù)計(jì)總體標(biāo)準(zhǔn)差（3）允許出現(xiàn)假陽性結(jié)果的機(jī)會(huì)（4）允許出現(xiàn)假陰性結(jié)果的機(jī)會(huì)第十六頁，共四十七頁，2022年，8月28日歐洲研究的樣本量估算決定每組含61名患者。第十七頁，共四十七頁，2022年，8月28日比較兩組發(fā)生某結(jié)局的百分比（1）預(yù)計(jì)兩組發(fā)生某結(jié)局的百分比約為（2）允許犯假陽性錯(cuò)誤的機(jī)會(huì)（3）允許犯假陰性錯(cuò)誤的機(jī)會(huì)

第十八頁，共四十七頁，2022年，8月28日南韓研究南韓對(duì)比劑研究第十九頁，共四十七頁，2022年，8月28日南韓研究的樣本量估算第二十頁，共四十七頁，2022年，8月28日5．劑量-反應(yīng)關(guān)系

不能作均數(shù)比較或回歸例有人分析蛇毒因子（CVF）的劑量對(duì)血液白細(xì)胞噬菌率的影響，得如下數(shù)據(jù)，欲討論劑量-反應(yīng)關(guān)系。第二十一頁，共四十七頁，2022年，8月28日

做法1：?jiǎn)我蛩胤讲罘治?!F=0.701，P>0.5,

均數(shù)間差別無統(tǒng)計(jì)學(xué)意義

為什麼不對(duì)？有負(fù)初衷——探討反應(yīng)隨劑量變化的趨勢(shì)*由多個(gè)劑量組的比較只能得知均數(shù)間是否有差異*有統(tǒng)計(jì)學(xué)差異也不等于有劑量-反應(yīng)關(guān)系第二十二頁，共四十七頁，2022年，8月28日

做法2：

反應(yīng)的均數(shù)關(guān)于劑量作回歸分析?!

為什麼不對(duì)？均數(shù)做因變量造成“好”的假象!*回歸方程是否有統(tǒng)計(jì)學(xué)意義與反應(yīng)的變異狀況有關(guān)*以諸個(gè)體反應(yīng)值的均數(shù)作回歸計(jì)算,掩蓋變異性回歸方程：Y=61.786-6.886log(劑量)決定系數(shù)：R2=0.914。第二十三頁，共四十七頁，2022年，8月28日正確作法：用個(gè)體資料作回歸分析回歸方程：Y=61.782-6.884log(劑量)決定系數(shù)：R2=0.095回歸方程無統(tǒng)計(jì)學(xué)意義，無劑量-反應(yīng)關(guān)系！第二十四頁，共四十七頁，2022年，8月28日6.重復(fù)測(cè)量資料不能時(shí)點(diǎn)間兩兩比較例各取7只兔子，分別以正常食物和待研究食物喂養(yǎng)，在實(shí)驗(yàn)前、喂養(yǎng)5周、10周后，各取血測(cè)量其中膽固醇濃度，自然對(duì)數(shù)轉(zhuǎn)換后,數(shù)據(jù)見表22.1,問血清膽固醇濃度隨時(shí)間變化的趨勢(shì)是否受該食物影響。第二十五頁，共四十七頁，2022年，8月28日例某藥物有新、舊兩種劑型。為比較兩種劑型的代謝情況，對(duì)16例某病患者服藥后0、4、8、12小時(shí)的血藥濃度作了測(cè)量，問該藥新舊兩種劑型的血藥濃度-時(shí)間曲線的差別是否具有統(tǒng)計(jì)學(xué)意義。

第二十六頁，共四十七頁，2022年，8月28日第二十七頁，共四十七頁，2022年，8月28日

常見的不妥

i)對(duì)每個(gè)時(shí)間點(diǎn)，分別計(jì)算均數(shù)、標(biāo)準(zhǔn)差；將各時(shí)間點(diǎn)測(cè)量值的均值用線連接，標(biāo)出“誤差”線ii)在各時(shí)間點(diǎn)做兩組比較(t

檢驗(yàn)或非參數(shù)檢驗(yàn)）第二十八頁，共四十七頁，2022年，8月28日為什麼不對(duì)？（1）連接各時(shí)間點(diǎn)測(cè)量值均數(shù)的方法將掩蓋個(gè)體曲線位置和形狀的特點(diǎn)：各均數(shù)所對(duì)應(yīng)的點(diǎn)連起來形成的曲線形狀可能與諸個(gè)體的曲線形狀毫不相干；各時(shí)間點(diǎn)標(biāo)上相應(yīng)的標(biāo)準(zhǔn)差也會(huì)誤導(dǎo)；（2）上述方法無法體現(xiàn)不同時(shí)間點(diǎn)的數(shù)據(jù)來自同一個(gè)體，而每一個(gè)體的多次重復(fù)測(cè)量值間具有相關(guān)性。（關(guān)鍵所在）第二十九頁，共四十七頁，2022年，8月28日應(yīng)當(dāng)如何？1）綜合指標(biāo)法（summarymeasuresapproach）或派生變量法（derivedvariableapproach）采用少數(shù)獨(dú)立的綜合指標(biāo)來概括每個(gè)個(gè)體多個(gè)時(shí)間點(diǎn)的測(cè)量值；用單變量方法比較各組的差異第三十頁，共四十七頁，2022年，8月28日綜合指標(biāo)第三十一頁，共四十七頁，2022年，8月28日2）重復(fù)測(cè)量資料的方差分析

先做“球形檢驗(yàn)：不同時(shí)間點(diǎn)上數(shù)據(jù)之間沒有相關(guān)性？任何兩個(gè)時(shí)間點(diǎn)之間的相關(guān)性都一樣？得陰性結(jié)果才進(jìn)一步作方差分析。例兔子喂養(yǎng)重復(fù)測(cè)量資料分析處理因素（group）：F=31.69，P=0.0001—食物對(duì)家兔血清膽固醇有影響；2)測(cè)量時(shí)間（time）：F=11.93，P=0.0003—家兔血清膽固醇濃度隨時(shí)間變化；3)交互效應(yīng)（time×group）：F=10.57，P=0.0005—不同食物，血清膽固醇濃度變化趨勢(shì)有所不同。第三十二頁，共四十七頁，2022年，8月28日7.多維列聯(lián)表資料不可反復(fù)做檢驗(yàn)例某地區(qū)呼吸系統(tǒng)疾病的患病率較高，有人懷疑與當(dāng)?shù)厥覂?nèi)點(diǎn)香的習(xí)慣有關(guān)，經(jīng)抽樣調(diào)查，得數(shù)據(jù)如下第三十三頁，共四十七頁，2022年，8月28日合并后作檢驗(yàn)？！

第三十四頁，共四十七頁，2022年，8月28日為什麼不對(duì)？i)點(diǎn)香與不點(diǎn)香者患病率不可比！年齡組合并掩蓋了:

點(diǎn)香（B1）與不點(diǎn)香（B2）者年齡結(jié)構(gòu)不同ii)不同年齡組患病率不可比！點(diǎn)香與不點(diǎn)香者合并掩蓋了:

不同年齡組中點(diǎn)香與不點(diǎn)香者比例的不同第三十五頁，共四十七頁，2022年，8月28日正確作法：

對(duì)數(shù)線性模型或logistic

回歸i)對(duì)數(shù)線性模型平等地討論A，B，C三因素對(duì)各格子中頻數(shù)的影響僅當(dāng)其他變量的影響與某變量無關(guān)時(shí)，才可通過合并消除該變量ii)logistic回歸

反應(yīng)變量:患病與否，解釋變量:點(diǎn)香和年齡交互效應(yīng):點(diǎn)香×年齡第三十六頁，共四十七頁，2022年，8月28日8.中西醫(yī)結(jié)合療效研究，應(yīng)設(shè)4組不少人作兩組比較：試驗(yàn)組：西藥A和中藥B

對(duì)照組：西藥A若(西藥A和中藥B)療效>西藥A療效能說中西醫(yī)結(jié)合好嗎?

第三十七頁，共四十七頁，2022年，8月28日不妥之一

如果(西藥A和中藥B)療效>西藥A療效能說：“在西藥A存在情形下，中藥B有療效”嗎？人們質(zhì)疑：“中藥B的療效可能是安慰效應(yīng)!?”于是，修改為：試驗(yàn)組：西藥A和中藥B

對(duì)照組：西藥A和中藥安慰劑第三十八頁，共四十七頁，2022年，8月28日不妥之二若(西藥A和中藥B)療效

>(西藥A和中藥安慰劑)療效則可以說：“在西藥A存在情形下，中藥B有療效！”但是，人們要問：在西藥A不存在情形下，中藥B有療效嗎？由此，足以認(rèn)為中西醫(yī)結(jié)合好嗎？第三十九頁，共四十七頁，2022年，8月28日

應(yīng)當(dāng)如何？

——析因設(shè)計(jì)！試驗(yàn)組1：西藥A和中藥B試驗(yàn)組2：西藥A和中藥安慰劑試驗(yàn)組3：中藥B和西藥安慰劑對(duì)照組：中藥安慰劑和西藥安慰劑若(西藥A和中藥B)療效-對(duì)照療效>[(西藥A和中藥安慰劑)療效-對(duì)照療效]+[(中藥B和西藥安慰劑)療效-對(duì)照療效]則說：中西醫(yī)結(jié)合療效比單獨(dú)療效之和還大！第四十頁，共四十七頁，2022年，8月28日

換一個(gè)角度看若(西藥A和中藥B)療效-對(duì)照療效>[(西藥A和中藥安慰劑)療效-對(duì)照療效]+[(中藥B和西藥安慰劑)療效-對(duì)照療效]改寫為：(西藥A和中藥B)療效-(中藥B和西藥安慰劑)療效

>(西藥A和中藥B安慰劑)療效-對(duì)照療效中藥B存在情形下，西藥A的療效

>中藥B不存在情形下，西藥A的療效則說：中藥B的存在使西藥A的療效有所增加！第四十一頁，共四十七頁，2022年，8月28日

析因設(shè)計(jì)可以告訴我們什么？西藥A的作用中藥B的作用西藥A和中藥B的協(xié)同作用若西藥A和中藥B的協(xié)同作用>0才說：中