Preciso estimar a função de risco de linha de base em um modelo Cox dependente do tempo
Enquanto eu fazia o curso de Sobrevivência, lembro que a derivada direta da função de risco cumulativo ( ) não seria um bom estimador, porque o estimador Breslow fornece uma função de etapa.
Então, existe alguma função no R que eu possa usar diretamente? Ou alguma referência a esse tópico?
Não tenho certeza se vale a pena abrir outra pergunta, portanto, apenas adiciono alguns antecedentes sobre por que a função de risco de linha de base é importante para mim. A fórmula abaixo estima a probabilidade de que o tempo de sobrevivência para um sujeito seja maior que o outro. Sob uma configuração do modelo Cox, a função de risco de linha de base é necessária.
Respostas:
Um modelo de Cox foi explicitamente projetado para poder estimar as taxas de risco sem ter que estimar a função de risco da linha de base. Esta é uma força e uma fraqueza. O ponto forte é que você não pode cometer erros nas funções que não estima. Essa é uma força real e é a razão pela qual as pessoas se referem a ela como "semi-paramétrica" e é em grande parte responsável por sua popularidade. No entanto, também é uma fraqueza real, pois, uma vez que você quer saber algo diferente da taxa de risco, muitas vezes precisará da função de risco de linha de base e isso derrota o próprio objetivo de um modelo de Cox.
Portanto, costumo usar os modelos Cox apenas quando estou interessado em taxas de risco e nada mais. Se eu quiser saber outras coisas, normalmente passo para outros modelos, como os discutidos aqui: http://www.stata.com/bookstore/flexible-parametric-survival-analysis-stata/
fonte
A função de risco da linha de base pode ser estimada em R usando a função "basehaz". O arquivo "help" afirma que é a função "sobrevivência prevista" que claramente não é. Se alguém inspeciona o código, é claramente a função de risco cumulativo de um
survfit
objeto. Para bobagens adicionais, a configuração padrão écentered=TRUE
a) a) não é uma função de risco de linha de base (como o nome sugere) eb) emprega previsão nos meios que são amplamente desacreditados como válidos em qualquer sentido prático.E ao seu ponto anterior: sim, esta função faz uso da função step. Você pode transformar essa saída em uma função de perigo usando suavização. A pior parte de tudo, qual é o intervalo de incerteza para essa previsão? Você pode receber uma medalha Fields se puder obtê-la. Acho que nem sabemos se o bootstrap funciona ou não.
Como um exemplo:
fonte