1、關(guān)于常見統(tǒng)計(jì)學(xué)錯(cuò)誤第1頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三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)”第2頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三錯(cuò)在哪里？哪些指標(biāo)可能有組間差異，必須心中有數(shù)。 科研的結(jié)果應(yīng)當(dāng)預(yù)見 假說是科研的靈魂 心中無數(shù)，不要“先上馬再說”指標(biāo)多，實(shí)驗(yàn)工作量大。 大海撈針 碰運(yùn)氣，不是科研指標(biāo)多，翻來覆去分析，制造

2、假陽性 Nature雜志統(tǒng)計(jì)學(xué)指南：常見錯(cuò)誤之一 第3頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三為何翻來覆去分析，會(huì)制造假陽性？第4頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三Nature常見錯(cuò)誤之一多重比較: 對(duì)一組數(shù)據(jù)作多項(xiàng)比較時(shí)，必須說明如何校正 水平，以避免增大第一類錯(cuò)誤的機(jī)會(huì)第5頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三應(yīng)當(dāng)如何？主要終點(diǎn)(primary end point) ：只能一個(gè) 次要終點(diǎn)(secondary end point) : 可以幾個(gè), 但勿過多Bonfferoni 校正 當(dāng)同一組數(shù)據(jù)同時(shí)作k次分析時(shí)，若限定犯假陽性錯(cuò)誤

3、的概率總共不超過 ，則每次分析要用 來控制假陽性的概率。例 第6頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三事后分析（Post hoc analysis） 看到數(shù)據(jù)之后才想起來做的分析 其結(jié)果往往不能被進(jìn)一步研究所確認(rèn)不鼓勵(lì)事先并未計(jì)劃的亞組分析和校正分析2. 事后分析，見機(jī)行事第7頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三為何不鼓勵(lì)事后分析？數(shù)據(jù)是現(xiàn)象，現(xiàn)象具有偶然性； 看到數(shù)據(jù)后進(jìn)行的“分析”屬于故作姿態(tài)，變偶然為“必然” 忽悠！看到數(shù)據(jù)后增加的“分析”，也屬于多重分析，增加假陽性的機(jī)會(huì)第8頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三必要的亞組分

4、析有助于發(fā)現(xiàn)適應(yīng)癥，指導(dǎo) 個(gè)體化治療；必要的校正分析有助于消除 混雜偏倚，公平、公正地比較療效亞組分析和校正分析必須事先在研究計(jì)劃 里規(guī)定，并說明理由；寫文章時(shí)，申明這 類分析是事前計(jì)劃好的事后決定做的分析，如果想寫，不能算是 正式結(jié)果，只能視為探索性分析。應(yīng)當(dāng)如何？第9頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三3. 隨機(jī)化，說而不做，做而不嚴(yán)處理分配的隨機(jī)化為什么這么重要? (1) 消除分配處理有意或無意的偏倚。 (2) 為實(shí)施盲法創(chuàng)造條件。 (3) 使有可能利用概率論來描述各干預(yù)組之間的差異有多大可能僅僅是由偶然性造成的。將隨機(jī)化當(dāng)作“廉價(jià)名詞”，實(shí)際沒做，卻寫“隨機(jī)分成兩

5、組” 科研道德？第10頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三說錯(cuò)和做錯(cuò)將隨機(jī)化當(dāng)作“廉價(jià)名詞”，實(shí)際沒做，卻寫“隨機(jī)分成兩組” 科研道德？將“隨意分組”當(dāng)作隨機(jī)化將“機(jī)械分組”當(dāng)作隨機(jī)化略去篩選過程，簡單地報(bào)告將多少人隨機(jī)分組略去實(shí)施過程中丟失對(duì)象，將最后兩組人數(shù)說成是隨機(jī)分組人數(shù)第11頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三應(yīng)當(dāng)如何？成功的隨機(jī)化取決于： (1) 產(chǎn)生一個(gè)不可預(yù)見的分配序列; (2) “隱蔽” （allocation concealment ）這個(gè)序列，直到分配完畢（必須建立一個(gè)分配處理的系統(tǒng)） 。報(bào)告如何隨機(jī)分組，如何“隱蔽” ：誰做隨

6、機(jī)序列，誰收病人，誰分藥和發(fā)藥；分組方案如何保管第12頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三參加者的流程圖 （強(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例頸部損傷）第13頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三4. 樣本

7、量隨意決定, 無依據(jù)樣本量太??？ - 真理應(yīng)經(jīng)得起重復(fù) 研究的功效太低樣本量越大越好？ - 時(shí)間長、經(jīng)費(fèi)不足 不能保證數(shù)據(jù)質(zhì)量需要事先估算最小樣本量 論文必須報(bào)告估算最小樣本量的根據(jù) 第14頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三比較兩組測定值的均數(shù)第15頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三歐洲對(duì)比劑研究 （1）預(yù)計(jì)欲比較的兩總體參數(shù)的差值 （2）預(yù)計(jì)總體標(biāo)準(zhǔn)差 （3）允許出現(xiàn)假陽性結(jié)果的機(jī)會(huì) （4）允許出現(xiàn)假陰性結(jié)果的機(jī)會(huì) 第16頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三歐洲研究的樣本量估算決定每組含61名患者。第17頁，共50頁，20

8、22年，5月20日，13點(diǎn)24分，星期三比較兩組發(fā)生某結(jié)局的百分比（1）預(yù)計(jì)兩組發(fā)生某結(jié)局的百分比約為（2）允許犯假陽性錯(cuò)誤的機(jī)會(huì)（3）允許犯假陰性錯(cuò)誤的機(jī)會(huì) 第18頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三南韓研究南韓對(duì)比劑研究第19頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三南韓研究的樣本量估算第20頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三5劑量-反應(yīng)關(guān)系 不能作均數(shù)比較或回歸例 有人分析蛇毒因子（CVF）的劑量對(duì)血液白細(xì)胞噬菌率的影響，得如下數(shù)據(jù)，欲討論劑量-反應(yīng)關(guān)系。第21頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三 做法

9、1：單因素方差分析?!F=0.701，P0.5, 均數(shù)間差別無統(tǒng)計(jì)學(xué)意義 為什麼不對(duì)？有負(fù)初衷 探討反應(yīng)隨劑量變化的趨勢* 由多個(gè)劑量組的比較只能得知均數(shù)間是否有差異* 有統(tǒng)計(jì)學(xué)差異也不等于有劑量-反應(yīng)關(guān)系第22頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三 做法 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.886 log(劑量)決定系數(shù)： R2=0.914。第23頁，共50頁，2022年，5月20日，13點(diǎn)24

10、分，星期三正確作法：用個(gè)體資料作回歸分析回歸方程： Y = 61.782-6.884 log(劑量)決定系數(shù)： R2=0.095回歸方程無統(tǒng)計(jì)學(xué)意義，無劑量-反應(yīng)關(guān)系！第24頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三6. 重復(fù)測量資料不能時(shí)點(diǎn)間兩兩比較例 各取7只兔子，分別以正常食物和待研究食物喂養(yǎng)，在實(shí)驗(yàn)前、喂養(yǎng)5周、10周后，各取血測量其中膽固醇濃度，自然對(duì)數(shù)轉(zhuǎn)換后, 數(shù)據(jù)見表22.1, 問血清膽固醇濃度隨時(shí)間變化的趨勢是否受該食物影響。第25頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三例 某藥物有新、舊兩種劑型。為比較兩種劑型的代謝情況，對(duì)16例某病患者服藥

11、后0、4、8、12小時(shí)的血藥濃度作了測量，問該藥新舊兩種劑型的血藥濃度-時(shí)間曲線的差別是否具有統(tǒng)計(jì)學(xué)意義。第26頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三第27頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三常見的不妥i)對(duì)每個(gè)時(shí)間點(diǎn)，分別計(jì)算均數(shù)、標(biāo)準(zhǔn)差； 將各時(shí)間點(diǎn)測量值的均值用線連接，標(biāo)出“誤差”線ii)在各時(shí)間點(diǎn)做兩組比較 ( t 檢驗(yàn)或非參數(shù)檢驗(yàn)）第28頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三為什麼不對(duì)？（1）連接各時(shí)間點(diǎn)測量值均數(shù)的方法將掩蓋個(gè)體曲線位置和形狀的特點(diǎn)： 各均數(shù)所對(duì)應(yīng)的點(diǎn)連起來形成的曲線形狀可能與諸個(gè)體的曲線形狀毫不相干；

12、各時(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ù)測量值間具有相關(guān)性。（關(guān)鍵所在）第29頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三應(yīng)當(dāng)如何？1）綜合指標(biāo)法（summary measures approach） 或 派生變量法（derived variable approach）采用少數(shù)獨(dú)立的綜合指標(biāo)來概括每個(gè)個(gè)體多個(gè)時(shí)間點(diǎn)的測量值；用單變量方法比較各組的差異第30頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三綜合指標(biāo)第31頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三2）重復(fù)測量資料的方差分析 先做

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

14、，經(jīng)抽樣調(diào)查，得數(shù)據(jù)如下第33頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三合并后作檢驗(yàn) ？！ 第34頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三為什麼不對(duì)？i) 點(diǎn)香與不點(diǎn)香者患病率不可比！ 年齡組合并掩蓋了: 點(diǎn)香（B1）與不點(diǎn)香（B2）者年齡結(jié)構(gòu)不同ii) 不同年齡組患病率不可比！ 點(diǎn)香與不點(diǎn)香者合并掩蓋了: 不同年齡組中點(diǎn)香與不點(diǎn)香者比例的不同第35頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三正確作法： 對(duì)數(shù)線性模型或 logistic 回歸i) 對(duì)數(shù)線性模型平等地討論A，B，C三因素對(duì)各格子中頻數(shù)的影響僅當(dāng)其他變量的影響與某變量無關(guān)時(shí)，才可通過

15、合并消除該變量ii) logistic回歸 反應(yīng)變量: 患病與否，解釋變量: 點(diǎn)香和年齡交互效應(yīng): 點(diǎn)香年齡第36頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三8. 中西醫(yī)結(jié)合療效研究，應(yīng)設(shè) 4 組不少人作兩組比較： 試驗(yàn)組：西藥A和中藥B 對(duì)照組：西藥A若 (西藥A和中藥B) 療效 西藥A 療效能說 中西醫(yī)結(jié)合好嗎? 第37頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三不妥之一 如果 (西藥A和中藥B) 療效 西藥A 療效能說：“在西藥A存在情形下，中藥B有療效”嗎？人們質(zhì)疑：“中藥B的療效可能是安慰效應(yīng)!?”于是，修改為：試驗(yàn)組：西藥A和中藥B 對(duì)照組：西藥A和中藥

16、安慰劑第38頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三不妥之二若 (西藥A和中藥B) 療效 (西藥A和中藥安慰劑) 療效則可以說： “在西藥A存在情形下，中藥B有療效！”但是，人們要問：在西藥A不存在情形下，中藥B有療效嗎？由此，足以認(rèn)為中西醫(yī)結(jié)合好嗎？第39頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三 應(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ì)照療效

17、 則說：中西醫(yī)結(jié)合療效比單獨(dú)療效之和還大！第40頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三 換一個(gè)角度看若 (西藥A和中藥B)療效 - 對(duì)照療效 (西藥A和中藥安慰劑)療效 - 對(duì)照療效 + (中藥B和西藥安慰劑)療效 - 對(duì)照療效 改寫為：(西藥A和中藥B)療效-(中藥B和西藥安慰劑)療效 (西藥A和中藥B安慰劑)療效- 對(duì)照療效中藥B存在情形下，西藥A的療效 中藥B不存在情形下，西藥A的療效則說：中藥B的存在使西藥A的療效有所增加！第41頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三 析因設(shè)計(jì)可以告訴我們什么？西藥A的作用 中藥B的作用西藥A和中藥B的協(xié)同作用若

18、 西藥A和中藥B的協(xié)同作用 0才說： 中西醫(yī)結(jié)合好！第42頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三例：某中西結(jié)合研究 22析因試驗(yàn)結(jié)果 (均數(shù)) 第43頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三兩因素間交互效應(yīng)（interaction）第44頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三小結(jié)終點(diǎn)指標(biāo)不宜過多, 大海撈針，很可能撈 上個(gè)假陽性。馬后炮分析不足取，事先計(jì)劃的分析才可信；不論干預(yù)性研究還是觀察性研究，都要事先到國際網(wǎng)站注冊(cè)，以接受同行監(jiān)督。隨機(jī)化是國際公認(rèn)的消除選擇偏倚的好辦法，不可作為廉價(jià)名詞，隨便說說；說而不做，做而不嚴(yán)，自欺欺人。第45頁，共50頁，2022年，5月20日，13點(diǎn)24分，星期三小結(jié)4. 干預(yù)研究的論文要報(bào)告樣本量估算的依據(jù)和方法；對(duì)不同的設(shè)計(jì)，樣本量估算的方法不同。5. 劑量-反應(yīng)關(guān)系既不能通過均數(shù)間比較，也不能用均數(shù)對(duì)劑量的回