More changes to
survHE. Earlier this week, I gave a talk at the Virtual Ispor US conference (our panel was W11). In preparation for the panel (and other things I was doing anyway), I got to talk to colleagues and thought a bit more of how survival modelling in HTA is kind of weird and “special”…
In particular, I realised how
survHE doesn’t give by default any real information about the hazard and the cumulative hazard functions — this isn’t entirely true, because once the survival curve has been estimated (which is the whole point of
survHE!) then \(h(t)\) and \(H(t)\) can be derived analytically in many cases, or at least numerically.
But I had not planned for a graphical representation of these two functions, which in fact is not really good practice, because they can be very helpful in understanding the model, particularly in the extrapolation portion (more on this here and in the NICE guidelines, among other places).
So, I’ve changed the code on the
devel branch of the
survHE GitHub repository so that the
plot function now has an extra option
what, which allows the user to choose whether the graph should depict \(S(t)\), \(h(t)\) or \(H(t)\).
library(survHE) # Fits a simple model to the built-in dataset, using a Weibull distribution m=fit.models(Surv(time,censored)~as.factor(arm),data,"wei") # Plots the survival curve plot(m)
# Plots the hazard curve plot(m,what="hazard")
# Plots the cumulative hazard curve plot(m,what="cumhazard")
The hazard curve is approximated by computing numerically the derivative of the cumulative hazard function.