David Phillippo
@dmphillippo.bsky.social
44 followers
110 following
67 posts
Statistician at University of Bristol | Bayesian, meta-analysis and evidence synthesis, #rstats
Posts
Media
Videos
Starter Packs
David Phillippo
@dmphillippo.bsky.social
· Jun 30
R for Health Technology Assessment
R for Health Technology Assessment discusses the use of proper statistical software, specifically R, to perform the whole pipeline of analytic modelling in health technology assessment (HTA). It has b...
www.routledge.com
David Phillippo
@dmphillippo.bsky.social
· Jun 17
David Phillippo
@dmphillippo.bsky.social
· Jun 17
David Phillippo
@dmphillippo.bsky.social
· Jun 17
PSA: multinma is temporarily unavailable from CRAN
A small memory allocation bug in Stan tripped some additional CRAN checks, which needs to be patched by rstan.
multinma is still fully functional and passes all tests. In the meantime:
A small memory allocation bug in Stan tripped some additional CRAN checks, which needs to be patched by rstan.
multinma is still fully functional and passes all tests. In the meantime:
David Phillippo
@dmphillippo.bsky.social
· Jun 17
Misaligned address sanitizer errors with `csr_matrix_times_vector()` - leading to CRAN package failing additional tests · Issue #1111 · stan-dev/rstan
Summary: Sparse matrix arithmetic using csr_matrix_times_vector() seems to trigger sanitizer "misaligned address" errors. Description: My package multinma that fits models using rstan has been flag...
github.com
David Phillippo
@dmphillippo.bsky.social
· Jun 17
David Phillippo
@dmphillippo.bsky.social
· Jun 17
David Phillippo
@dmphillippo.bsky.social
· Jun 17
Multilevel network meta-regression for general likelihoods: synthesis of individual and aggregate data with applications to survival analysis
Network meta-analysis combines aggregate data (AgD) from multiple randomised controlled trials, assuming that any effect modifiers are balanced across populations. Individual patient data (IPD) meta-regression is the "gold standard" method to relax this assumption, however IPD are frequently only available in a subset of studies. Multilevel network meta-regression (ML-NMR) extends IPD meta-regression to incorporate AgD studies whilst avoiding aggregation bias, but currently requires the aggregate-level likelihood to have a known closed form. Notably, this prevents application to time-to-event outcomes.
We extend ML-NMR to individual-level likelihoods of any form, by integrating the individual-level likelihood function over the AgD covariate distributions to obtain the respective marginal likelihood contributions. We illustrate with two examples of time-to-event outcomes, showing the performance of ML-NMR in a simulated comparison with little loss of precision from a full IPD analysis, and demonstrating flexible modelling of baseline hazards using cubic M-splines with synthetic data on newly diagnosed multiple myeloma.
ML-NMR is a general method for synthesising individual and aggregate level data in networks of all sizes. Extension to general likelihoods, including for survival outcomes, greatly increases the applicability of the method. R and Stan code is provided, and the methods are implemented in the multinma R package.
arxiv.org
David Phillippo
@dmphillippo.bsky.social
· Jun 17
David Phillippo
@dmphillippo.bsky.social
· Jun 17
David Phillippo
@dmphillippo.bsky.social
· Jun 17
David Phillippo
@dmphillippo.bsky.social
· Jun 17
David Phillippo
@dmphillippo.bsky.social
· Jun 17