Modelo de riscos proporcionais de Cox censurado por intervalo em R

Dado o tempo de sobrevivência com intervalo de censura, como faço para executar um modelo Cox PH com intervalo de censura R? Uma pesquisa rseek exibe o pacote intcox, que não existe mais no Rrepositório. Estou quase certo de que a coxphfunção no survivalpacote não pode lidar com dados de sobrevivência censurados por intervalo.

Além disso, não quero imputar os dados e depois usar a coxphfunção Este método subestima os erros padrão dos coeficientes porque você está ignorando a incerteza da censura de intervalo.

Como mencionado acima, você pode usar a função survreg. Uma observação: este não é estritamente um modelo de Cox PH, mas sim modelos em escala de localização. Usando a transformação de log padrão, este é o modelo de ré. No caso da distribuição exponencial, os riscos proporcionais e o modelo à ré são equivalentes; portanto, se a distribuição for definida como exponencial, este será um modelo de riscos proporcionais com uma linha de base exponencial. Da mesma forma, se um modelo de linha de base de distribuição Weibull a ré for usado, as estimativas de parâmetro são apenas uma transformação linear daquelas usadas no modelo de riscos proporcionais com a distribuição de linha de base Weibull. Mas, em geral, o survreg não se encaixa no modelo de Cox PH.

Se um modelo semi-paramétrico for desejado, conforme encontrado implementado no intcox, uma palavra de cautela: há vários problemas com a versão atual do intcox (o algoritmo normalmente termina prematuramente significativamente longe do MLE, falha totalmente com observações sem censura, sem erros padrão apresentado automaticamente).

Uma nova alternativa que você pode usar é o pacote "icenReg".

Admissão de viés: esse é o autor do icenReg.

Para fazer uma análise censurada por intervalo em R, você deve criar um objeto Surv e, em seguida, usar survfit (). Se você tiver mais de uma variável, o pacote intcox resolve o problema.