2013-01-15 27 views
5

Verilerim için en iyi Cox Orantısal Tehlikeler modellerini seçerken, R paketindeki sıçrama paketindeki regsubsets'e benzer bir işlevi kullanmaya çalışıyorum. Mümkün mü? ve eğer böyle bir işlev zaten mevcutsa? Ben zaten aşina tahmin ediyorumTüm alt kümeler değişken seçimini R? 'Daki ​​Cox Oransal Tehlikeler Modeli için yapabilir miyim?

+2

Önce bunun iyi bir fikir olup olmadığını ciddi olarak sorgulamalısınız. –

+0

Hangi model seçim prosedürünü önerirsiniz? Bir modeli "en iyi" olarak kabul etmeden önce rakip modellerin ne kadar yakın olduğunu görmek uygun değil mi? – marcellt

+0

Öncelikle yatırım alanında temel bilim düşünmelisiniz. Hedefler açık olmadığı sürece "en iyi prosedür" önermek mümkün değildir. –

cevap

0

aşağıdaki ... sonra bu makul bir başlangıç ​​noktası olacağını 'top model' için kriter olarak AIC kullanıyorsanız:

library(survival) 
data(colon) 
c1 <- coxph(Surv(time=time, event=status) ~ 
    as.factor(extent) + age + sex, data=colon) 
step(c1) 

Dikkat Eğer eksik değerleriniz varsa (NA). Elbette bu yöntemle bulunmayan daha iyi bir model olabilir, ancak az sayıda olası tahmin ediciyle bunu kaçırmanız mümkün değildir. Bilinçli bir görüşün daha güvenilir olabileceği sayısal yöntemler kullanılarak kullanım hakkında yukarıdaki gibi (teşekkürler @DWin).