Hayatta kalma analizinde Tehlike Oranını hesaplamak için logrank ve Mantel-Haenszel yöntemini kullanmanın artıları ve eksileri nelerdir?


17

İki sağkalım eğrisinin karşılaştırmasını özetlemenin bir yolu, tehlike oranını (HR) hesaplamaktır. Bu değeri hesaplamak için en az iki yöntem vardır.

  • Logrank yöntemi. Kaplan-Meier hesaplamalarının bir parçası olarak, her bir gruptaki ( ve ) gözlemlenen olayların (genellikle ölümler) sayısını ve hayatta kalma farkı olmayan bir sıfır hipotezi ( ve ) varsayarak beklenen olayların sayısını hesaplayın. . Tehlike oranı şu : O b E a E b H R = ( O a / E a )ÖbirÖbEbirEb
    'HR,=(Öbir/Ebir)(Öb/Eb)
  • Mantel-Haenszel yöntemi. İlk olarak her bir zaman noktasındaki hipergeometrik varyansların toplamı olan V hesaplayın. Sonra tehlike oranını şu şekilde hesaplayın: Bu denklemleri Machin, Cheung ve Parmar, Hayatta Kalma Analizi bölüm 3'ten aldım . Bu kitap, iki yöntemin genellikle çok benzer yöntemler verdiğini ve aslında kitaptaki örnekte olduğu gibi olduğunu belirtir.
    'HR,=tecrübe((Öbir-Ebir)V)

Birisi bana iki yöntemin üç faktöre göre değiştiği bir örnek gönderdi. Bu özel örnekte, logrank tahmininin mantıklı olduğu ve Mantel-Haenszel tahmininin çok uzakta olduğu açıktır. Sorum, tehlike oranının logrank tahmininin en iyi olduğu zaman ve Mantel-Haenszel tahmininin ne zaman seçilmesinin en iyi olduğu konusunda herhangi bir genel tavsiyenin olup olmadığıdır. Örneklem büyüklüğü ile ilgili mi? Bağ sayısı? Numune boyutlarının oranı?


Bu tahminler Cox regresyonunun verdiği tahminle nasıl ilişkilidir? İK tahmininde altın standart bu olmalı.
Aniko

Cox modeli ortak değişkenler içerir. Kaplan-Meier, Nelson-Aalen, Mantel-Haenszel yöntemleri sadece yaşın bir fonksiyonu olarak tehlikeyi modellemektedir.
shabbychef

@ shabbychef: Cox PH ile tek bir ikili ortak değişken kullanın, yani referans / karşılaştırma grupları için 0/1 kodlu, sonra exp (beta) = HR kullanın.
ars

Log-rank, oransal Tehlikeler varsayımı karşılandığında Cox PH'dan daha güçlü bir testtir. Bu nedenle, tek bir 2 seviyeli eş değişkenle, bir log-rank veya Mantel-Haenszel testi tercih edilir.
Thylacoleo

cevap için aşağıya bakın ...
Thylacoleo

Yanıtlar:


11

Sanırım cevabı buldum (kendi soruma). Orantılı tehlike varsayımı doğruysa, iki yöntem tehlike oranı için benzer tahminler verir. Sanırım şimdi belirli bir örnekte bulduğum tutarsızlık, bu varsayımın şüpheli olmasından kaynaklanıyor.

Orantılı tehlikeler varsayımı doğruysa, o zaman log (zaman) ve log (-log (St)) (St'nin t zamanında oransal sağkalım olduğu) grafiği iki paralel çizgi göstermelidir. Sorunlu veri kümesinden oluşturulan grafik aşağıdadır. Doğrusal olmaktan uzak görünüyor. Orantılı tehlike varsayımı geçerli değilse, tehlike oranı kavramı anlamsızdır ve bu nedenle tehlike oranını hesaplamak için hangi yöntemin kullanıldığı önemli değildir.

alternatif metin

Logrank ve Mantel-Haenszel'in tehlike oran tahminleri arasındaki tutarsızlığın orantılı tehlike varsayımını test etmek için bir yöntem olarak kullanılıp kullanılamayacağını merak ediyorum.


7

Yanılmıyorsam, başvuruda bulunduğunuz log-rank tahmincisi Pike tahmincisi olarak da bilinir. Genelde HR <3 için önerildiğine inanıyorum çünkü bu aralıkta daha az önyargı sergiliyor. Aşağıdaki makale ilginizi çekebilir (kağıdın O / E olarak adlandırıldığını unutmayın):

[...] O / E yöntemi önyargılıdır, ancak klinik çalışmalarda söz konusu tehlike oranlarının oranı değerleri aralığında, ortalama kare hatası açısından CML veya Mantel-Haenszel'den daha etkilidir. en büyük denemeler hariç tüm yöntem. Mantel-Haenszel yöntemi minimum düzeyde önyargılıdır, CML kullanılarak elde edilenlere çok yakın cevaplar verir ve tatmin edici yaklaşık güven aralıkları sağlamak için kullanılabilir.


Bu makaleye kısa bir göz attıktan sonra, düşündükleri tahminlerin, sorucının denklemleri ile aynı olduğundan emin değilim. Soru altındaki yorumlara katılıyorum - belki 1981'de yaklaşık yöntemler yararlıydı, ancak bu günlerde Cox regresyonunu kullanmamak için açık bir neden yok.
onestop

@onestop: hmm, yukarıda unutulan günlükle O / E == LR tanımını mı düşünüyorsunuz? Cox PH hakkında söylediklerinize katılıyorum - cevaplamaya çalıştığım soru bu değil, ancak tavsiyeniz daha geniş bağlamda daha iyi.
ars

Bernstein ve diğ. ark. iki yöntemin yanlış veya farklı olmasına neden olan bazı nedenleri (küçük n, bağlar) gösterir. Ancak gösterdikleri tüm tutarsızlıklar azdır. Bu yüzden bu makaledeki hiçbir şeyin bu soruyu doğuran üç kat farklılığını açıkladığını sanmıyorum. Geldiğim cevap için aşağıya bakın.
Harvey Motulsky

7

Aslında birkaç yöntem daha var ve seçim genellikle en erken farkları mı, daha sonra farkları mı yoksa - log-rank testi ve Mantel-Haenszel testinde olduğu gibi - tüm zaman noktalarına eşit ağırlık verip vermediğinize bağlıdır.

Eldeki soruya. Log-rank testi aslında hayatta kalma verilerine uygulanan Mantel-Haenszel testinin bir şeklidir. Mantel-Haenszel testi genellikle tabakalı durum tablolarında bağımsızlığı test etmek için kullanılır.

MH testini sağkalım verilerine uygulamaya çalışırsak, her başarısızlık zamanındaki olayların bağımsız olduğunu varsayarak başlayabiliriz. Daha sonra başarısızlık zamanına göre katmanlaşırız. MH yöntemlerini, her başarısızlık süresini bir katman haline getirerek kullanırız. Şaşırtıcı olmayan bir şekilde, genellikle aynı sonucu verirler.

Birden fazla olay aynı anda gerçekleştiğinde istisna oluşur - tam olarak aynı zaman noktasında çoklu ölümler. O zaman tedavinin nasıl farklı olduğunu hatırlayamıyorum. Bence log-rank testi, bağlı başarısızlık zamanlarının olası sıralamaları üzerinde ortalamalar.

Böylece log-rank testi hayatta kalma verileri için MH testidir ve bağlarla başa çıkabilir. Hayatta kalma verileri için MH testini hiç kullanmadım.


3

Ben bir web sitesi genelinde tökezledi düşündüm ve tam olarak bu soru ile ilgilenen referans:

http://www.graphpad.com/faq/viewfaq.cfm?faq=1226 "Karşılaştırılan iki yöntem" ile başlayın.

Site (yukarıda) bağlantılı Berstein kağıt ars referansı:

http://www.jstor.org/stable/2530564?seq=1

Site, Berstein ve arkadaşlarının sonuçlarını güzel bir şekilde özetliyor, bu yüzden alıntı yapacağım:

İkisi genellikle aynı (veya neredeyse aynı) sonuçları verir. Ancak, birkaç kişi aynı anda öldüğünde veya tehlike oranı 1.0'dan uzak olduğunda sonuçlar değişebilir.

Bernsetin ve arkadaşları simüle verileri her iki yöntemle analiz etmişlerdir (1). Tüm simülasyonlarında orantılı tehlike varsayımı doğruydu. İki yöntem çok benzer değerler verdi. Logrank yöntemi (O / E yöntemi olarak adlandırdıkları), özellikle tehlike oranı büyük veya örnek boyutu büyük olduğunda, gerçek Tehlike Oranından 1.0'a yakın değerleri bildirir.

Bağlar olduğunda, her iki yöntem de daha az doğrudur. Logrank yöntemleri, 1.0'a daha yakın olan tehlike oranlarını bildirme eğilimindedir (bu nedenle rapor edilen tehlike oranı, tehlike oranı 1,0'dan büyük olduğunda çok küçük ve tehlike oranı 1,0'dan düşük olduğunda çok büyüktür). Buna karşılık Mantel-Haenszel yöntemi, 1.0'dan daha büyük tehlike oranları bildirir (bu nedenle rapor edilen tehlike oranı, 1.0'dan büyük olduğunda çok büyük ve tehlike oranı 1.0'dan düşük olduğunda çok küçüktür).

İki yöntemi orantılı tehlike varsayımının doğru olmadığı yerlerde simüle edilmiş verilerle test etmediler. İki İK tahmininin çok farklı olduğu bir veri seti gördüm (üç faktörle) ve bu veriler için orantılı tehlike varsayımı şüpheli. Görünüşe göre Mantel-Haenszel yöntemi, geç zaman noktalarındaki tehlikedeki farklılıklara daha fazla ağırlık verirken, logrank yöntemi her yerde eşit ağırlık verir (ancak bunu ayrıntılı olarak araştırmadım). İki yöntemle çok farklı HR değerleri görüyorsanız, orantılı tehlike varsayımının makul olup olmadığını düşünün. Bu varsayım makul değilse, elbette tüm eğriyi tanımlayan tek bir tehlike oranı kavramının tamamı anlamlı değildir.

Site aynı zamanda, "İK'nın iki tahmininin çok farklı olduğu (üç faktöre göre)" veri kümesine atıfta bulunur ve PH varsayımının kilit bir husus olduğunu gösterir.

Sonra "Siteyi kim yazdı?" Diye düşündüm. Biraz arama yaptıktan sonra Harvey Motulsky olduğunu gördüm. Bu yüzden Harvey, kendi sorunuzu cevaplarken size referans vermeyi başardım. Otorite oldunuz!

"Sorun veri kümesi" herkese açık bir veri kümesi midir?


Cevabı iki gün önce anladım ve buraya yeni bir cevap olarak gönderdim. Daha sonra bulduğunuz graphpad.com adresindeki web sayfasını genişletip güncelledim. Sorunlu verileri içeren bir Excel dosyasına bağlantı eklemek için bu sayfayı tekrar düzenledim ( graphpad.com/faq/file/1226.xls ). Verileri üreten adamdan izin alana kadar bunu yapamadım (anonim olmak istiyor ve veriler belirsiz bir şekilde etiketleniyor).
Harvey Motulsky
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.