Posts

flexsurv 2

Previously, I started discussing the flexsurv package. I used it to fit a Weibull model. This is implemented as an accelerated failure …
Nov 1, 2019 5 min read

flexsurv

I’m going to write about some of my favourite R packages. I’ll start with flexsurv (https://github.com/chjackson/flexsurv-dev) by Chris …
Oct 28, 2019 4 min read

Longitudinal hurdle models 3

In the last post on longitudinal hurdle models, I had just taken samples from the marginal mean \[\begin{align} g(\theta; x) & = …
Oct 25, 2019 3 min read

Longitudinal hurdle models 2

In a previous post I fit a longitudinal hurdle model using the brms package. library(brms) summary(fit_hurdle) ## Family: …
Oct 23, 2019 5 min read

Longitudinal hurdle models

Data Recently, I have been modelling data that is longitudinal, contains excess zeros, and where the non-zero data is right-skewed and …
Oct 16, 2019 4 min read

How to calculate the log-rank statistic

Suppose we have the following data. df <- dplyr::tibble(patient_id = as.character(1:12), treatment = rep(c("C", …
Jun 3, 2019 2 min read