
一、競爭風險分析簡介
競爭風險分析是生存分析的一種延伸方法,用於處理有多種失敗事件(或風險)且事件之間互相競爭的時間-事件資料。例如,在骨髓移植的研究中,病人可能因復發或在緩解期間死亡而結束觀察,兩者為競爭事件,因為某一事件發生後,另一事件便無法再發生。如:對於復發,死亡是一個競爭風險,因為復發事件在患者死亡後不再發生,或者說患者還來不及復發就先死亡了。
在競爭風險分析中,主要有兩個重要概念:

競爭風險分析是生存分析的一種延伸方法,用於處理有多種失敗事件(或風險)且事件之間互相競爭的時間-事件資料。例如,在骨髓移植的研究中,病人可能因復發或在緩解期間死亡而結束觀察,兩者為競爭事件,因為某一事件發生後,另一事件便無法再發生。如:對於復發,死亡是一個競爭風險,因為復發事件在患者死亡後不再發生,或者說患者還來不及復發就先死亡了。
在競爭風險分析中,主要有兩個重要概念:
在進行多變量邏輯回歸(Multivariable Logistic Regression)時,大家通常知道 EPV (Events Per Variable) 一般來說要大於 10 才穩定。但最常被問到的問題是:「這裡的 Variable,到底是指變項個數,還是參數個數?」
舉例說明: 「如果你有一個『教育程度』變項,分為五個組別,它在模型中會產生4 個 Dummy variables。在計算 EPV 時,它貢獻的不是 1,而是 4。如果再加上其他變項,分母會迅速增加,這也是為什麼小樣本研究容易出現模型不穩定的原因。」
根據主流統計學文獻(如 Peduzzi et al. 1996 與 Harrell et al. 1985)以及後續的模擬研究(如 Wynants et al. 2015),關於 EPV 的計算方式,正確答案是:
EPV 的分母應為「模型中估計的參數數量(Parameters)」,意即必須計入所有的 Dummy variables。
以下為詳細的文獻理據與計算規範:
最近看到有研究者在做回溯性研究(retrospective research) 時,也會擷取未校正的Odds Ratio (OR) or Hazard Ratio (HR) 值進行統合分析 (Meta-analysis),讓我原本認為應該要用adjusted OR/HR的認知似乎不同,因此再次查閱文獻,看一下當進行統合分析中有retrospective research時應如何處理? 整理如下:
在回溯性研究的 Meta-Analysis 中,應該使用 adjusted OR 或 adjusted HR,而不是 crude 值。這是因為 adjusted 值已經考慮了混雜變數的影響,能更準確地反映變數之間的真實關係。以下是詳細的解釋與注意事項:
在社會科學、心理學及行為科學研究中,中介變數(mediator) 與 調節變數(moderator) 常常被用來探討變數之間的關係。而在醫學研究中,這兩種變數的使用可以幫助研究者更深入地了解健康結果與風險因素之間的關係,從而改善臨床實踐和公共衛生政策。因此了解中介與調節變數的差異並對於進行研究設計、資料分析與結果解釋至關重要。
中介變數 是指能夠解釋自變數(Independent variable, IV)如何影響依變數(Dependent variable, DV)的中間過程。它描述了一個「因果鏈」,即自變數先影響中介變數,而中介變數進一步影響依變數。因此,透過中介變數,我們可以更清楚地理解自變數與依變數之間的作用機制。 在醫學研究中,中介變數有助於解釋某一治療、風險因素或環境因素是如何影響患者健康結果的。例如,透過探討治療干預措施如何通過中介變數產生效果,研究者可以更清楚地理解治療的作用機制。
調節變數 是指會影響自變數與依變數之間關係強度或方向的變數。它本身不會直接影響自變數或依變數,但能夠「調整」兩者的關係,使得該關係在不同情境或條件下產生不同的效果。 在醫學研究中, 調節變數常被用來研究不同的患者特徵如何影響治療效果或風險因素與健康結果之間的關係。這種分析有助於識別哪些亞群的患者可能從特定治療中獲益更多,或哪些群體的風險因素更顯著。

通常做研究時,當會拿到原始數據,再導入統計軟件中直接分析,但偏偏遇到有些資料丟失,只有當初已統計好的平均值、標準差以及個案數,那該如何求得p 值?
若是問正在學統計學的學生,絕對沒有問題,套入公式在查表便是,但統計軟體用的行雲流水的職場統計師,可有點傷腦筋了~要再回顧學生時的場景了....這邊分幾個步驟講解,最後提供驗證過的SAS語法,之後用起來便可以信心十足拉~也可以審閱他人文章或是檢查學生結果算得對不對啦~
步驟:

當研究之連續變項沒有臨床上的切點時,大多數會將其均分成3或4組,此處用3個連續變項均分成3組來示範:
方法一
利用SAS macro可以快速找出多個變項的百分位點,此處以找最小值,第33百分位,第66百分位及最大值為例:
data a;
do i=1 to 100;
今天審稿時,看到一篇文章: 3組之間連續值比較直接做兩兩比較,所以有3個P-value: p of A vs B, P of A vs C, P of B vs C,這時問題來了~
1. p of A vs B, P of A vs C, P of B vs C 是否有經過Bonferroni or Dunn-Sidak adjustment?
(看統計方法及Table note 均無說明,應該是沒做修正了~)
2.應該要再提供一個完整的比較: Overall p-value,通常 Overall p-value 顯著了,才會再進一步做事後檢定
3.此篇探討的是Post hoc comparisons for the Kruskal-Wallis test,是針對3組以上連續資料無母數分析之案例,若是針對3組以上連續資料母數分析,應採用ANOVA and 相應的事後檢定分析

雖然我是SAS愛用者,使用SAS已有十幾年的經驗,但一直沒有去考證照,近幾年SAS公司積極推出各項課程及 SAS® 全球認證計劃在校園中推廣,發現越來越多青年學子的履歷中,有SAS證照的比例越來越高了!不禁讓我也開始排文查找相關資訊,整理如下,可以讓大家參考~

此文為兩組配對資料檢定-Wilcoxon Signed Rank Test Paired t test using SAS 之進階版
分兩個語法介紹,語法2會用到SAS macro,當重複變項很多時,可善用macro,提升工作效率
語法1
在SAS Blog看到這樣的詢問,一方面自己也使用過這兩種方法
問題: 為何 ranuni(123) 和 rand('uniform') 不同?

當兩組樣本資料屬性為連續變項,相互配對或重覆且
(1) 其差距服從常態分佈時,則採Paired t test檢定兩組差距之平均值是否相等;
(2) 其差距未服從常態分佈時,則採Wilcoxon Signed Rank Test檢定兩組差距之中位數是否相等。
目標: 檢測before & after 是否有顯著差異
Step 1: 產生10個1~50隨機配對資料;