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

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

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

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

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

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

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

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

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

（強烈推薦）合格對象82例拒絕參與7例隨機分組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)費不足不能保證數(shù)據(jù)質(zhì)量需要事先估算最小樣本量論文必須報告估算最小樣本量的根據(jù)

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

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

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

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

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

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

為什麼不對？有負初衷——探討反應(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ù)、標準差；將各時間點測量值的均值用線連接，標出“誤差”線ii)在各時間點做兩組比較(t

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

先做“球形檢驗：不同時間點上數(shù)據(jù)之間沒有相關(guān)性？任何兩個時間點之間的相關(guān)性都一樣？得陰性結(jié)果才進一步作方差分析。例兔子喂養(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)疾病的患病率較高，有人懷疑與當?shù)厥覂?nèi)點香的習慣有關(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ù)的影響僅當其他變量的影響與某變量無關(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有療效嗎？由此，足以認為中西醫(yī)結(jié)合好嗎？第三十九頁，共四十七頁，2022年，8月28日

應(yī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才說：中