常見統(tǒng)計學(xué)錯誤

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

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

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

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

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

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

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

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

（強(qiáng)烈推薦）合格對象82例拒絕參與7例隨機(jī)分組75例分配至實驗組38例接受干預(yù)38例分配至對照組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ù)研究的功效太低樣本量越大越好？

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

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

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

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

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

做法1：單因素方差分析?!F=0.701，P>0.5,

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

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

做法2：

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

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

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

常見的不妥

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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