Il survival
pacchetto R
sembra focalizzarsi su modelli di sopravvivenza a tempo continuo. Sono interessato a stimare una versione temporale discreta di un modello di rischio proporzionale, il modello log-log complementare. Ho un modello di sopravvivenza abbastanza semplice, con una semplice censura giusta.
So che un modo per stimare questo modello è quello di creare un set di dati con una riga separata per ogni osservazione per ogni periodo in cui non è "morto". Quindi, è possibile utilizzare un glm
modello con il cloglog
collegamento.
Questo approccio sembra molto memoria inefficiente; in effetti, probabilmente produrrebbe un set di dati che è troppo grande per la memoria sulla mia macchina.
Un secondo approccio sarebbe codificare il MLE me stesso. Sarebbe abbastanza semplice, ma spero che esista un pacchetto con questo modello di sopravvivenza. Sarebbe più semplice collaborare ed evitare errori di codifica per usare un pacchetto.
Qualcuno sa di un tale pacchetto?
cloglog
link, comunque.
coxph(ties="exact")
, nelsurvival
pacchetto standard , renda il modello "un modello logistico condizionale, ed è appropriato quando i tempi sono un piccolo insieme di valori discreti". Questo non funzionerebbe per te? È b / c che non userebbe ilcloglog
link?