Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

Multilevel and Quasi Monte Carlo methods for the calculation of the Expected Value of Partial Perfect Information

Wei Fang, Zhenru Wang, Michael B. Giles, View ORCID ProfileChris H. Jackson, Nicky J. Welton, Christophe Andrieu, Howard Thom
doi: https://doi.org/10.1101/2021.03.30.21254626
Wei Fang
1Mathematical Institute, University of Oxford
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zhenru Wang
1Mathematical Institute, University of Oxford
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michael B. Giles
1Mathematical Institute, University of Oxford
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Chris H. Jackson
4MRC Biostatistics Unit, University of Cambridge
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Chris H. Jackson
Nicky J. Welton
2Population Health Science, Bristol Medical School, University of Bristol
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Christophe Andrieu
3School of Mathematics, University of Bristol
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Howard Thom
2Population Health Science, Bristol Medical School, University of Bristol
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: howard.thom@bristol.ac.uk
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Data/Code
  • Preview PDF
Loading

Abstract

The expected value of partial perfect information (EVPPI) provides an upper bound on the value of collecting further evidence on a set of inputs to a cost-effectiveness decision model. Standard Monte Carlo (MC) estimation of EVPPI is computationally expensive as it requires nested simulation. Alternatives based on regression approximations to the model have been developed, but are not practicable when the number of uncertain parameters of interest is large and when parameter estimates are highly correlated. The error associated with the regression approximation is difficult to determine, while MC allows the bias and precision to be controlled. In this paper, we explore the potential of Quasi Monte-Carlo (QMC) and Multilevel Monte-Carlo (MLMC) estimation to reduce computational cost of estimating EVPPI by reducing the variance compared with MC, while preserving accuracy. In this paper, we develop methods to apply QMC and MLMC to EVPPI, addressing particular challenges that arise where Markov Chain Monte Carlo (MCMC) has been used to estimate input parameter distributions. We illustrate the methods using a two examples: a simplified decision tree model for treatments for depression, and a complex Markov model for treatments to prevent stroke in atrial fibrillation, both of which use MCMC inputs. We compare the performance of QMC and MLMC with MC and the approximation techniques of Generalised Additive Model regression (GAM), Gaussian process regression (GP), and Integrated Nested Laplace Approximations (INLA-GP). We found QMC and MLMC to offer substantial computational savings when parameter sets are large and correlated, and when the EVPPI is large. We also find GP and INLA-GP to be biased in those situations, while GAM cannot estimate EVPPI for large parameter sets.

Competing Interest Statement

The authors have declared no competing interest.

Funding Statement

WF, ZW, and HT were supported by the Hubs for Trials Methodology Research (HTMR) network grant N79 for this work. HT and NJW were supported by the HTMR Collaboration and innovation in Difficult and Complex randomised controlled Trials In Invasive procedures (ConDuCT-II). HT and NJW were also supported by the National Institute for Health Research (NIHR) Bristol Biomedical Research Centre (BRC) for part of this work. HT was furthermore supported by MRC grant MR/S036709/1. CA would like to thank support of EPSRC EP/R018561/1 Bayes4Health. CJ was funded by the UK Medical Research Council programme MC_UU_00002/11. The directly acting oral anticoagulants for prevention of stroke in atrial fibrillation model was funded by NIHR Health Technology Assessment programme project number 11/92/17 and NIHR Senior Investigator award NF-SI-0611-10168.

We are also grateful to Mark Strong at University of Sheffield for providing his R code to estimate EVPPI, plus its SE and upward bias, using GAM and GP.

Author Declarations

I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.

Yes

The details of the IRB/oversight body that provided approval or exemption for the research described are given below:

Not necessary as modelling study synthesising published, and not individual patient, data.

All necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived.

Yes

I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).

Yes

I have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.

Yes

Data Availability

All code and data necessary to the produce the results are available in this GitHub repository: https://github.com/Bogdasayen/Example_MLMC_and_QMC_for_EVPPI The simple demo of MLMC can be found at this repository: https://github.com/Bogdasayen/mlmc_evppi_demo

https://github.com/Bogdasayen/mlmc_evppi_demo

Copyright 
The copyright holder for this preprint is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. It is made available under a CC-BY-NC-ND 4.0 International license.
Back to top
PreviousNext
Posted April 04, 2021.
Download PDF
Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
Multilevel and Quasi Monte Carlo methods for the calculation of the Expected Value of Partial Perfect Information
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Multilevel and Quasi Monte Carlo methods for the calculation of the Expected Value of Partial Perfect Information
Wei Fang, Zhenru Wang, Michael B. Giles, Chris H. Jackson, Nicky J. Welton, Christophe Andrieu, Howard Thom
medRxiv 2021.03.30.21254626; doi: https://doi.org/10.1101/2021.03.30.21254626
Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
Citation Tools
Multilevel and Quasi Monte Carlo methods for the calculation of the Expected Value of Partial Perfect Information
Wei Fang, Zhenru Wang, Michael B. Giles, Chris H. Jackson, Nicky J. Welton, Christophe Andrieu, Howard Thom
medRxiv 2021.03.30.21254626; doi: https://doi.org/10.1101/2021.03.30.21254626

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Health Economics
Subject Areas
All Articles
  • Addiction Medicine (76)
  • Allergy and Immunology (194)
  • Anesthesia (54)
  • Cardiovascular Medicine (488)
  • Dentistry and Oral Medicine (89)
  • Dermatology (56)
  • Emergency Medicine (168)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (211)
  • Epidemiology (5661)
  • Forensic Medicine (3)
  • Gastroenterology (215)
  • Genetic and Genomic Medicine (856)
  • Geriatric Medicine (88)
  • Health Economics (229)
  • Health Informatics (759)
  • Health Policy (388)
  • Health Systems and Quality Improvement (250)
  • Hematology (105)
  • HIV/AIDS (181)
  • Infectious Diseases (except HIV/AIDS) (6455)
  • Intensive Care and Critical Care Medicine (388)
  • Medical Education (116)
  • Medical Ethics (28)
  • Nephrology (90)
  • Neurology (845)
  • Nursing (44)
  • Nutrition (141)
  • Obstetrics and Gynecology (161)
  • Occupational and Environmental Health (258)
  • Oncology (514)
  • Ophthalmology (162)
  • Orthopedics (44)
  • Otolaryngology (105)
  • Pain Medicine (47)
  • Palliative Medicine (21)
  • Pathology (149)
  • Pediatrics (248)
  • Pharmacology and Therapeutics (146)
  • Primary Care Research (113)
  • Psychiatry and Clinical Psychology (959)
  • Public and Global Health (2222)
  • Radiology and Imaging (375)
  • Rehabilitation Medicine and Physical Therapy (174)
  • Respiratory Medicine (311)
  • Rheumatology (109)
  • Sexual and Reproductive Health (80)
  • Sports Medicine (82)
  • Surgery (118)
  • Toxicology (25)
  • Transplantation (34)
  • Urology (42)