end-stage renal disease

Similar documents
Section 3.2 Least-Squares Regression

IAPT: Regression. Regression analyses

Eating and Sleeping Habits of Different Countries

MODELING DISEASE FINAL REPORT 5/21/2010 SARAH DEL CIELLO, JAKE CLEMENTI, AND NAILAH HART

Choosing Life: empowerment, Action, Results! CLEAR Menu Sessions. Adherence 1: Understanding My Medications and Adherence

Chapter 1 Introduction

Chapter 3: Examining Relationships

Following is a list of topics in this paper:

Lab 4 (M13) Objective: This lab will give you more practice exploring the shape of data, and in particular in breaking the data into two groups.

About Health TV with Jeanne Blake Prostate cancer update

Making Your Treatment Work Long-Term

PESTICIDE TOXICOLOGY. Introduction. Materials and Methods. Warning about handling pesticides. Experimental design

Unit 8 Day 1 Correlation Coefficients.notebook January 02, 2018

Clincial Biostatistics. Regression

Choosing Life: Empowerment, Action, Results! CLEAR Menu Sessions. Substance Use Risk 5: Drugs, Alcohol, and HIV

1.4 - Linear Regression and MS Excel

Political Science 15, Winter 2014 Final Review

COPING WITH A CANCER DIAGNOSIS. Tips for Dealing with What Comes Next

MBios 478: Systems Biology and Bayesian Networks, 27 [Dr. Wyrick] Slide #1. Lecture 27: Systems Biology and Bayesian Networks

Chapter 23. Inference About Means. Copyright 2010 Pearson Education, Inc.

Sheila Barron Statistics Outreach Center 2/8/2011

Quadratic Functions I

3.2 Least- Squares Regression

Stepwise method Modern Model Selection Methods Quantile-Quantile plot and tests for normality

Standard Deviation and Standard Error Tutorial. This is significantly important. Get your AP Equations and Formulas sheet

Testing for. Prostate Cancer

Chapter 12. The One- Sample

appstats26.notebook April 17, 2015

Regressions usually happen at the following ages 4 months, 8 months, 12 months, 18 months, 2 years.

Graphical Exploration of Statistical Interactions. Nick Jackson University of Southern California Department of Psychology 10/25/2013

(2) In each graph above, calculate the velocity in feet per second that is represented.

Analysis Plans in Economics

Two-Way Independent ANOVA

1. Scatter or bar? 2. Scatter or bar? 3. Scatter or bar? 4. Scatter or bar? Response rate, ms. Average Ht, ft. Temperature, C. Sex

Difficult conversations. Dr Amy Waters MBBS, FRACP Staff Specialist in Palliative Medicine, St George Hospital Conjoint Lecturer, UNSW

Chunyang Ding. A Sweet Conclusion: Comparing the Osmotic Potentials of Yams and Sweet Potatoes. 2 October Mr. Allen. AP/IB Biology P.

t-test for r Copyright 2000 Tom Malloy. All rights reserved

Choosing Life: Empowerment, Action, Results! CLEAR Menu Sessions. Substance Use Risk 2: What Are My External Drug and Alcohol Triggers?

Mathematical Modeling of Infectious Diseases

NEXTGEN ICD10 TIPS DEMONSTRATION

Overdose. Everything you need to know. Includes free DVD!

3.2A Least-Squares Regression

Building Friendships: Avoid Discounting

1. The figure below shows the lengths in centimetres of fish found in the net of a small trawler.

MiSP Solubility Lab L3

I don t want to be here anymore. I m really worried about Clare. She s been acting different and something s not right

Camel Milk and Cancer: A Possible Cure?

Take new look emotions we see as negative may be our best friends (opposite to the script!)

Diets Cause Weight Re-GAIN Here s How

PROTEIN LAB BASED ON THE RESEARCH OF DR. RICHARD LONDRAVILLE

Keys to a Successful Legislative Visit

STP 231 Example FINAL

Selection at one locus with many alleles, fertility selection, and sexual selection

Safer injecting Reducing your risk. Includes free DVD for injectors!

AP Human Geography Kuby: Tracking the AIDS Epidemic. Mapping the Diffusion of AIDS

PSY 205 Module 3 Supplement. Comparing Correlation, Ex post facto, and Experimental Approaches to Research

Correlation and regression

What You Need to Know About. TB Infection TUBERCULOSIS

Sample Exam Paper Answer Guide

Reverse Engineering a Regression Table. Achim Kemmerling Willy Brandt School Erfurt University

An Introduction to Huntington s Disease

MULTIPLE LINEAR REGRESSION 24.1 INTRODUCTION AND OBJECTIVES OBJECTIVES

The Pretest! Pretest! Pretest! Assignment (Example 2)

Lose Weight. without dieting.

Chapter 1. Dysfunctional Behavioral Cycles

Overcoming Perfectionism

Mathematics Meets Oncology

K to 12 BASIC EDUCATION CURRICULUM SENIOR HIGH SCHOOL CORE SUBJECT

Living Life with Persistent Pain. A guide to improving your quality of life, in spite of pain

Unraveling Recent Cervical Cancer Screening Updates and the Impact on Your Practice

STAT 201 Chapter 3. Association and Regression

Subliminal Messages: How Do They Work?

Estimation. Preliminary: the Normal distribution

Lecture 12: Normal Probability Distribution or Normal Curve

Ordinary Least Squares Regression

Business Statistics Probability

Chapter 7: Descriptive Statistics

An Interview with a Chiropractor

CBT+ Measures Cheat Sheet

5 To Invest or not to Invest? That is the Question.

Designing and Analyzing RCTs. David L. Streiner, Ph.D.

9 INSTRUCTOR GUIDELINES

HIV/AIDS Patient Involvement. in Antiretroviral Treatment Decisions

Gage R&R. Variation. Allow us to explain with a simple diagram.

BOOKLET ONE. Introduction to Behavioural Activation for Depression

Pancreatic Cancer: Associated Signs, Symptoms, Risk Factors and Treatment Approaches

Anxiety. Top ten fears. Glossophobia fear of speaking in public or of trying to speak

Total English Placement Test

Teens. Self-Talk. Be positive. Practice. Try it and see. Using your thoughts, feelings and actions PATIENT EDUCATION

Constructing a Bivariate Table:

Economics 345 Applied Econometrics

Day One: After you ve tested positive

Step Five. Admitted to ourselves and another human being the exact nature of our wrongs.

Regression CHAPTER SIXTEEN NOTE TO INSTRUCTORS OUTLINE OF RESOURCES

DOCTOR: The last time I saw you and your 6-year old son Julio was about 2 months ago?

Teachers Students Total

Complex Regression Models with Coded, Centered & Quadratic Terms

Michael Norman International &

Still important ideas

Presentation outline. Issues affecting African Communities in New Zealand. Key findings Survey. Findings cont... Findings cont..

Transcription:

Case study: AIDS and end-stage renal disease Robert Smith? Department of Mathematics and Faculty of Medicine The University of Ottawa

AIDS and end-stage renal disease ODEs Curve fitting AIDS End-stage renal disease Mortality Prevalence Mortality Prevalence Curve fitting Curve fitting Determining impact Determining impact Model fitting Fit lines to data Estimating parameters Appendix G Lab Determine impact Revise model Predicting future outcomes Explore drug effects Realistic outcomes Effects of different treatment efficacies

End-stage renal disease Many patients with AIDS develop end-stage renal disease One of the opportunistic infections that kills you Basically, kidney failure In the US, this is particularly prevalent in African Americans Thus, we ll focus on this subset of the population.

Antiretroviral drugs HAART has drastically changed the face of HIV Reduced the number of AIDS deaths Made HIV a disease it s possible to live with Not clear what effect HAART has had on the prevalence of AIDS or end-stage renal disease......so we ll investigate it ourselves.

Our questions 1. Has HAART had an impact on the prevalence of AIDS? 2. Has HAART had an impact on the prevalence of end-stage renal disease? 3. If aggressive treatment is initiated now, with different effects, what will the long-term outcome be?

Our approach We ll need to formulate a model fit parameters to data draw conclusions predict the future This combines various strands of modelling while incorporating real-world data.

Impact of HAART on AIDS mortality? Mortality data from the CDC: These are the number of deaths due to AIDS for African Americans in the US To see it a bit more clearly, let s plot it. 1991 10475 1992 11946 1993 15460 1994 17844 1995 18971 1996 15909 1997 10333 1998 8744 1999 9097 2000 8723 2001 9085 2002 8927 2003 9077 2004 9302 2005 8562

HAART has clearly had an effect 20,000 18,000 16,000 14,000 AIDS deaths 12,000 10,000 8,000 HAART 6,000 4,000 2,000 0 1990 1995 2000 2005 time

What about prevalence? We expect prevalence to increase not as many people are dying other people are progressing from HIV to AIDS But perhaps not as sharply as it did before HAART Prevalence data from the CDC: Again, we ll plot this. 1991 14561 1992 15897 1993 60649 1994 71847 1995 81317 1996 92319 1997 105464 1998 117890 1999 112483 2000 121903 2001 181475 2002 193814 2003 204466 2004 214017 2005 225270

Prevalence of AIDS 250000 200000 AIDS prevalence 150000 100000 50000 0 1990 1995 2000 2005 time

Formulating a null hypothesis Q. How can we tell if HAART has made a difference? A. Construct a null hypthesis and test against it: N0: HAART has had no significant impact on the prevalence of AIDS (note that we are not assuming the prevalence goes up or down).

Testing the null hypothesis Q. How can we test the impact? A. Fit curves to pre- and post-haart data Compare with entire data set The data is approximately linear, so this will make life easier.

Prevalence pre- and post-haart 250000 200000 prevalence 150000 100000 HAART 50000 0 1990 1995 2000 2005 time

Initial thoughts The slope of the first line is steeper than the second Both lines are good fits: r = 0.949 and r = 0.967, respectively So it looks like HAART may have reduced the rate of increase of the prevalence (as we d hope) However, this is only half the story We still need to compare to the fit overall.

Prevalence overall 250000 200000 prevalence 150000 100000 50000 0 1992 1994 1996 1998 2000 2002 2004 time

Impact of HAART on AIDS prevalence Is this a better or worse fit? The eye can t tell, so we need to rely on the regressional coefficient In this case r = 0.981...higher than either r from before! Thus, we can t reject the null hypothesis It follows that HAART has had no significant impact on prevalence of AIDS among African Americans in the US (It has drastically reduced mortality, though).

Impact on mortality? Mortality data from the US Renal Data System: These are the number of deaths due to end-stage renal disease for African Americans in the US Again, let s plot it. 1991 88 1992 126 1993 159 1994 176 1995 255 1996 185 1997 120 1998 141 1999 149 2000 131 2001 126 2002 143 2003 135 2004 128

HAART has clearly had an effect 250 end!stage renal disease mortality 200 150 100 50 HAART 0 1990 1995 2000 2004 time

What about prevalence? We again expect prevalence to increase Prevalence data from the US Renal Data System: Again, we ll plot this. 1991 14561 1992 15897 1993 60649 1994 71847 1995 81317 1996 92319 1997 105464 1998 117890 1999 112483 2000 121903 2001 181475 2002 193814 2003 204466 2004 214017 2005 225270

Prevalence of AIDS 2500 end!stage renal disease prevalence 2000 1500 1000 500 0 1990 1995 2000 2004 time

Formulating another null hypothesis N0: HAART has had no significant impact on the prevalence of end-stage renal disease Again, we ll fit linear curves to pre- and post- HAART, as well as the entire data set (it helps that the data is approximately linear, of course).

Prevalence pre- and post-haart 3000 2500 end!stage renal disease prevalence 2000 1500 1000 500 HAART 0 1990 1995 2000 2005 time

Initial thoughts The slope of the second line is steeper than the first Both lines are even better fits: r = 0.98887 and r = 0.98683, respectively So it looks like HAART may have increased the prevalence of end-stage renal disease (not out of the question, as many more people are alive because of HAART).

Prevalence overall 3000 2500 end!stage renal disease prevalence 2000 1500 1000 500 0 1990 1995 2000 2005 time

Impact on prevalence Once again, the eye can t tell, so we need to rely on the regressional coefficient In this case r = 0.99352...higher than either r from before! Thus, we can t reject the null hypothesis It follows that HAART has had no significant impact on prevalence of either AIDS or endstage renal disease among African Americans in the US (It has drastically reduced mortality, though).

Model fitting What we ve done in each case is fit models to data True, they were simple, linear models, but we still made choices We can now use these linear fits to estimate parameters and construct a more complex differential equation model.

How to construct such a model? We have two variables of interest: AIDS prevalence and end-stage renal disease prevalence Since end-stage renal disease doesn t cause AIDS, we can consider AIDS in isolation This makes our first equation much easier.

The equation for AIDS prevalence We now have faith that it s a linear fit, so let s construct a linear differential equation Prevalence is increasing, so the derivative will be positive Thus, we could write da dt = g This is simple and we could solve it if we wanted to, but we won t We need estimates for g and A(0).

Estimating g and A(0) Using our linear fit, we can estimate g = 15133 (the slope) A(0) = 14959 (the intercept) Technical note: Time really starts at 1991, so we need to transpose the x-axis by 1991. prevalence 250000 200000 150000 100000 50000 0 1992 1994 1996 1998 2000 2002 2004 time

Equation for end-stage renal disease This is a bit trickier A proportion s of people with AIDS develop end-stage renal disease People with end-stage renal disease die at rate, proportional to the prevalence of end-stage renal disease dn = sa δn dt Solving this is harder (and depends on A(t)) We need estimates for s, and N(0).

Estimating N(0) Using our linear fit, we can estimate N(0) = 268 (intercept) slope = 179 Slope is the derivative, so we can use this to estimate s and. end!stage renal disease prevalence 3000 2500 2000 1500 1000 500 0 1990 1995 2000 2005 time

Estimating s and Picking two points that are close to the linear fit, we have dn = sa δn = 179 dt s(20564) δ(1287) = 179 s(117890) δ(1521) = 179 [ ] [ s 20564 1287 = δ 117890 1521 [ ] 0.0048 =. 0.2563 ] 1 [ 179 179 ]

Summarising parameter estimates We thus have A(0) = 14959 N(0) = 268 g = 15133 s = 0.0048 = 0.2563 Note that we didn t solve either equation, even though one was easy and the other was doable.

Predicting future outcome Now that we have our model, we can use it to predict the future Unfortunately, by definition, we don t have actual evidence about the future (and if we wait for the evidence to arrive, it won t be the future any more) To compensate for this, we ll make a range of predictions.

Initiating aggressive HAART now Currently, treatment hasn t done much to slow the epidemic of end-stage renal disease However, treatment also hasn t been applied as aggressively as it could Especially in disadvantaged groups like African Americans Can HAART eliminated end-stage renal disease?

Effects of treatment We can represent treatment by a factor (1-h) If h=0, then treatment has no effect on progression to end-stage renal disease If h=1, treatment completely suppresses progression to end-stage renal disease Our equation thus becomes dn dt = s(1 h)a δn

A range of blocking effects Let s consider a number of values of h: h = 0.38, 0.65, 0.80, 0.95, 1 We ll run the original equation from 1991 til 2007, then each of the new equations from 2007 until 2035 We ll plot each of them on the same graph, so we can compare the effects.

Effects of aggressive HAART 4000 38% projected end!stage renal disease prevalence 3500 3000 2500 65% 2000 1500 80% 1000 500 95% 0 100% 1990 1995 2000 2005 2010 2015 2020 2025 2030 2035 time

What does this tell us? The only way to eliminate end-stage renal disease is if HAART is 100% effective at blocking progression (unlikely) All other therapies have an initial dip and then rise again Even 95% effective therapy will eventually lead to an increase in prevalence. projected end!stage renal disease prevalence 4000 38% 3500 3000 2500 65% 2000 1500 80% 1000 500 95% 0 100% 1990 1995 2000 2005 2010 2015 2020 2025 2030 2035 time

So what was the point? This doesn t mean our model doesn t tell us anything useful Even therapy that was only 38% effective would not result in an increase in prevalence for about 25 years If we re looking for eradication, then we d be disappointed But these delays give us time to come up with new strategies and hold back the disease.

What s the take-home message? Our model tells us that putting all our energy into perfecting treatment might not be the best use of our time Unless we can have 100% effective treatment, we re not going to eradicate the disease However, even fairly ineffective therapy can do a lot of good in the meantime.

Why model? In this way, modelling gives us useful information about whether to proceed or not, knowing likely outcomes 4000 And we did this with 3500 nothing more 3000 sophisticated than 2500 linear regression and 2000 simple ordinary 1500 differential equations. 1000 projected end!stage renal disease prevalence 500 38% 65% 80% 95% 0 100% 1990 1995 2000 2005 2010 2015 2020 2025 2030 2035 time

AIDS and end-stage renal disease ODEs Curve fitting AIDS End-stage renal disease Mortality Prevalence Mortality Prevalence Curve fitting Curve fitting Determining impact Determining impact Model fitting Fit lines to data Estimating parameters Appendix G Lab Determine impact Revise model Predicting future outcomes Explore drug effects Realistic outcomes Effects of different treatment efficacies

Lab work Use the mortality data on AIDS and endstage renal disease to fit lines to pre- HAART and post-haart data Compare this to a linear fit to the entire data Adjust the model to include the effect of HAART reducing the prevalence of AIDS Explore different effects: 20%, 50%, 80% Use the model to interpret biological outcomes as treatment approaches (or exceeds) 100%.