A mathematical model of tumor dynamics during stereotactic body radiation therapy for non-small cell lung cancer
|
|
- Noah Jacobs
- 5 years ago
- Views:
Transcription
1 A mathematical model of tumor dynamics during stereotactic body radiation therapy for non-small cell lung cancer Russell Injerd, Emma Turian Abstract Image guided radiotherapy allows tumor volume dynamics to be measured at intervals during treatment. This has improved our ability to model certain types of tumors, like non-small cell lung cancer, during treatment. Previous one and two population ODE models of tumor volume dynamics in response to stereotactic ablative radiotherapy (SABR) were created using exponential and logistic growth. Previous studies indicated that a two population exponential model provided the best balance between t and mathematical complexity and may serve a functional role in clinical practice. Our study reevaluates previous ndings and tests the suitability of Gompertz growth models by assessing the goodness of t versus its mathematical complexity. The exponential, Velhurst, and Gompertz models are calibrated using data from 11 patients extracted from a previous study and evaluated using statistical methods, such as the Akaike information criterion and the Wilcoxon rank sum test. Model comparison indicates that Gompertz growth model does not signicantly improve model t, however it may oer advantages as an inferential tool for clinical use. Keywords Image-guided radiotherapy, dynamic modeling, model selection, non-small cell lung cancer, stereotactic ablative radiotherapy. 1 Introduction Tumor modeling has undergone many changes that have resulted in improvements in the description and understanding of tumor dynamics. Complex models can include morphology, pressure and nutrient gradients, and genetic probability [3]. These methods have been most luminary when applied to in vitro spheroids, but are increasingly used in multi-phase and in vitro models with success. The complexity and variability of tumor growth models in vivo require specic patient information which makes regular clinical use problematic. There remains a use for simple mathematical models that make use This study was a requirement for Math 371, Spring Department of Mathematics, NEIU, Chicago, IL, 60625, math@neiu.edu. 1
2 R. Injerd 2 of regularly collected clinical data. Image guided radiotherapy, like stereotactic ablative radiotherapy (SABR), provides precise delivery of radiation using imaging during treatment. This allows the use of higher fractional doses of radiation and reduces the time patients are exposed to radiation [4]. As a result of these new therapies, tumor volume is recorded at several points during treatment. This provides researchers the data to apply mathematical models to describe tumor volume dynamics. Because tumor volume has been found to be a signicant predictor in cancer prognosis [5], simple models, like ordinary dierential equations, that incorporate tumor volume measurements collected during treatment could be a practical clinical tool. A previous study examined 11 patients undergoing SABR treatment and evaluated exponential and logistic models using the Akaike Information Criterion [5]. The study demonstrated that simple ODE models can adequately describe tumor volume of non-small cell lung cancer (NSCLC) patients receiving SABR. This study reassesses these earlier models and introduces Gompertz growth to model the tumor volume dynamics. 2 Data Patient data was extracted from a previous clinical study of SABR treatment for NSCLC [5]. Patients received targeted doses of radiation, higher than the conventional 2 Gy per fraction, aided by advanced image guidance techniques [4]. Images of the tumor are captured prior to each fractional dose. Gross tumor volume (GTV) is estimated by an experienced medical professional by a combination of automation and manual adjustment techniques. A second physicist conducts a second independent estimation of the GTV. The average variability between each clinician's GTV measurements was 3.7% [5]. Patients with issues during GTV contouring due to tissue density classication were excluded from the study. Patients receiving less than 5 fractional doses were also excluded. The study includes 11 patients with NSCLC receiving 5 or 8 fractional doses of SABR treatment with tumors that allow accurate contouring. Patient characteristics are summarized in table 1. A detailed table of GTV can be found in appendix A.
3 R. Injerd 3 Table 1: Patient characteristics Patient no. Age Gender No. of fractions Dose per fraction (Gy) Initial GTV (cc) 1 81 Male Male Female Male Female Female Female Female Female Female Male Mathematical Modeling The main objective of this study is to assess the Gompertz growth function as an alternative model to describe NSCLC tumor volume in patients receiving SABR. Previous studies [5] have assessed exponential and Velhurst models as pragmatic methods to describe tumor dynamics in NSCLC SABR patients. These non-spatialized models make use of available data and omit highly detailed microscopic processes to balance model accuracy and model complexity. Gompertz growth is also an accepted and widely used model for describing tumor growth [3]. Gompertz growth does not introduce any new parameters. Using it to describe NSCLC tumors receiving SABR may improve model accuracy without introducing additional complexity. Previous models were reassessed alongside Gompertz models using both one and two population methods. 3.1 One Population Models Figure 1: One population tumor growth models. growth Tumor Cells Fractional Dose Surviving Fraction GTV One population models assume tumors are composed of homogeneous cells that can be described by a single variable. Though it has been shown through
4 R. Injerd 4 CT imaging that tumor volume is heterogeneous in terms of density, it is practical and useful to assume homogeneity for model simplicity. The rst model is simple exponential growth given by the following equation: dn kn (1) dt where t is time, and k > 0 is the proliferation rate. The initial tumor volume will be taken from the rst CT measurement. This model assumes that all vital nutrients are available in abundance and places no constraints on growth. This assumption is reasonable for small tumors like stage 1 NSCLC with tumor diameter less than 5cm [5]. The linear-quadratic equation will be used to model the eects of radiotherapy. After each fractional dose a proportion of the cells will die and the remainder of the cells survive treatment. S expp αd βd 2 q (2) where S represents the surviving fraction, d is the dose per fraction, and α and β are the linear and quadratic coecients for cell killing. To limit the possible combinations of α and β the ratio is xed to 10Gy 1, thereby requiring estimation of α only [5]. After a dose of radiotherapy at time t, an innitesimal time later p tq the tumor volume will be: N pt tq N ptqs (3) which will continue to proliferate until the next treatment. Equations (1)-(3) complete the rst model which will be referred to as "POP1EXP". When nutrients are thought to be limited, which is more realistic for tumor cell proliferation, growth can be shown to slow as cells reach a carrying capacity. This can be modeled using the Velhurst growth model: and the Gompertz growth model: dn dt kn 1 N θ dn dt kn ln N θ where k > 0 is the proliferation rate and θ > 0 is the carrying capacity of the cell population. These models have fast initial growth that slow as the carrying capacity is approached. Equations (4), (2), and (3) complete model 2 which will be referred to as POP1LOG. Equations (5), (2), and (3) complete model 3 which will be referred to as POP1GOM. The owchart shown in Fig. 1 illustrates the change in GTV using one population models. (4) (5)
5 R. Injerd Two population models. Figure 2: Two population tumor growth model. Living Cells Fractional Dose Surviving Fraction growth GTV Dead Cells clearance Dead Cells not cleared from tumor Modeling tumors using a single population shows instantaneous shrinking of the tumor volume after a fractional dose of radiotherapy as shown in Eq.(3). Actual reduction in tumor volume after radiotherapy requires dead cells to be cleared from the tumor site, which may take appreciable time [5]. To model this process, it becomes useful to model two populations of cells, living N l and dead N d. Using this methodology exponential growth becomes: the Velhurst growth: or the Gompertz growth: dn l dt dn l dt dn l dt kn l (6) kn l 1 N l θ kn l ln Nl θ After receiving a fractional dose, a proportion of the cells survive irradiation shown by Eq.(2) and the rest become dead cells: (7) (8) N l pt tq N l ptqs (9) N d pt tq N d ptq N l ptqr1 Ss (10) Dead cells will be cleared from the tumor exponentially: dn d dt cn d (11)
6 R. Injerd 6 where c > 0 is the clearance rate. The total tumor volume is the sum of the living and dead cells: N N l N d (12) Using two populations is a useful method to incorporate a delay into dead cell clearance. However, CT images cannot distinguish living and dead cells. The only volume collected at each measurement is the gross tumor volume. Using two populations is merely a mathematical technique to incorporate delay into volume reduction caused by radiotherapy. Equations (6), (2), and (9)-(12) complete model 4 which will be referred to as POP2EXP. Equations (7), (2), and (9)-(12) complete model 5 which will be referred to as POP2LOG. Equations (8), (2), and (9)-(12) complete model 6 which will be referred to as POP2GOM. The owchart shown in Fig. 2 illustrates the change in GTV using two population models. 3.3 Model Calibration & Selection The models are calibrated to the GTV of each patient by estimating the parameters which minimize the sum of squared errors (SSE): SSE Ţ t1 N ptq ˆN ptq 2 (13) where N ptq is the measured GTV right before the t th fraction, T is the total number of fractions for an individual patient, and ˆN ptq is the corresponding modeled GTV. Constraints are placed on each of the parameters to omit impossible values, while providing enough exibility to allow the tting algorithm to nd a suitable t [5]. Sensitivity analysis will be conducted using Monte Carlo simulations. The coecient of variation CV will be calculated to determine which parameters the modeled GTV are most sensitive to. Table 2: Model parameters Parameter Description Units Constraints k proliferation rate d 1 [ ] α cell killing coecient Gy 1 [ ] θ carrying capacity cc [N p0q - 100] c clearance rate for dead cells d 1 [ ] To assess goodness of t the % root mean squared error (%RMSE) is calculated: RMSE a SSE{T (14) %RMSE RMSE{pNmax Nminq 100% (15)
7 R. Injerd 7 The Shapiro-Wilks normality test will be conducted on the residuals. Graphical analysis will also be conducted using quantile-quantile plots (QQ plots). Care must be taken when using models as inferential tools if residuals are not normal; it implies that the predictive value of the model changes across the range of data. Pearson's correlation analysis will also be conducted to determine whether the parameters are independent or dependent on each other. Dependent parameters may not be fully identiable from available data and inferential use may be problematic. To balance model t and model complexity Akaike Information Criterion (AIC) is calculated using the mean sum of dierences (MSD): MSD 1 T Ţ t1 2 ˆN ptq N ptq (16) N ptq AIC T ln pmsdq 2p (17) where p is the number of free parameters. The AIC rewards goodness of t, while penalizing the number of free parameters. The Wilcoxon paired rank sum test will be used to determine whether the %RMSE and the AIC values calculated for the models dier signicantly. 4 Results Generally, ts are similar in the one population models. Median %RM SE for POP1EXP, POP1LOG, and POP1GOM are %, %, and % respectively. The Wilcoxon paired rank sum test shows signicant dierence in the medians of POP1GOM and POP1EXP with a p-value below.05. Two population models appear to improve in t moving from POP2EXP, to POP2LOG, and to POP2GOM with median % RMSE of %, %, and % respectively. The Wilcoxon paired rank sum test did not produce p-values low enough to reject the null hypothesis of equal medians. There was signicant improvement in t moving from one population models to two population models. The t data is summarized in Fig. (3) and Wilcoxon p-values can be found in appendix B.
8 R. Injerd 8 Figure 3: Box plots of each model comparing t for all 11 patients When complexity is accounted for by penalizing free parameters, one and two population models appear not to benet from the inclusion of a carrying capacity. There is improvement in AIC scores of two population models over one population models. Median AIC values for POP1EXP, POP1LOG, and POP1GOM are , , and respectively. Median values for POP2EXP, POP2LOG, and POP2GOM are , , respectively. Few of the models had p-values low enough to reject the null hypothesis of equal medians. The AIC data is summarized in Fig. (4) and Wilcoxon p-values can be found in appendix C.
9 R. Injerd 9 Figure 4: Box plots of each model comparing t and complexity for all 11 patients Sensitivity analysis conducted using 10 7 iterations for each model indicate α produces the most variation in 4 of the 6 models. The results of the sensitivity analysis are shown in appendix D. The Pearson correlation analysis indicate that most of the parameters are dependent with p-values above.05. The results of the Pearson correlation analysis are shown in appendix F. The Shapiro-Wilks test and QQ plots indicate that the residuals of several patients are not normally distributed. The results of the Shapiro-Wilks test are found in appendix E and QQ plots are included in appendix H. 5 Concluding remarks As reported in previous studies, two population models t the patient data better than one population models [5]. Two population models also produce α values closer to reported values. The POP2EXP had the lowest AIC value indicating that the added θ parameter included in POP2LOG and POP2GOM increases information loss. Though POP2GOM did not outperform POP2EXP when accounting for complexity, it appeared to provide a better t when judging by the %RM SE alone. The Shapiro-Wilks test and the QQ plots also suggest that residuals are more consistently close to being normally distributed for POP2GOM. Depending on the patient, using a two population Gompertz model may improve the utility of inferential tests. Because the study only in-
10 R. Injerd 10 cluded 11 patients and the Wilcoxon rank sum tests failed to show signicance in median dierences, more research and data are necessary to determine which two population model is best suited for clinical application. 6 Acknowledgments This study was produced during the Cancer Modeling Seminar, Math 371, as part of its required work load.
11 R. Injerd 11 A Patient Tumor Volume Day GT V 1 GT V 2 GT V 3 GT V 4 GT V 5 GT V 6 GT V 7 GT V 8 GT V 9 GT V 10 GT V B RMSE% p-values Paired Wilcoxon Rank Sum POP1EXP POP2EXP POP1LOG POP2LOG POP1GOM POP2GOM POP1EXP POP2EXP POP1LOG POP2LOG POP1GOM POP2GOM C AIC p-values Paired Wilcoxon Rank Sum POP1EXP POP2EXP POP1LOG POP2LOG POP1GOM POP2GOM POP1EXP POP2EXP POP1LOG POP2LOG POP1GOM POP2GOM
12 R. Injerd 12 D Sensitivity Analysis µ x σ CV POP1EXP α k POP2EXP α k c POP1LOG α k θ POP2LOG α k θ c POP1GOM α k θ POP2GOM α k θ c E Residual Shapiro-Wilks Test Statistic Patient no. POP1EXP POP2EXP POP1LOG POP2LOG POP1GOM POP2GOM Omitted values indicate a p-value<.05 indicating likely non-normality
13 R. Injerd 13 F Parameter Correlation α θ c POP1EXP k POP2EXP α k POP1LOG α k θ POP2LOG α k θ POP1GOM α k θ POP2GOM α k θ Omitted values indicate a p-value<.05 indicating likely independence
14 R. Injerd 14 G Patient Plots Figure 5: Patient 1 Figure 6: Patient 2
15 R. Injerd 15 Figure 7: Patient 3 Figure 8: Patient 4
16 R. Injerd 16 Figure 9: Patient 5 Figure 10: Patient 6
17 R. Injerd 17 Figure 11: Patient 7 Figure 12: Patient 8
18 R. Injerd 18 Figure 13: Patient 9 Figure 14: Patient 10
19 R. Injerd 19 Figure 15: Patient 11
20 R. Injerd 20 H QQ Plots Figure 16: Patient 1 POP1EXP Figure 19: Patient 1 POP2LOG Figure 17: Patient 1 POP2EXP Figure 20: Patient 1 POP1GOM Figure 18: Patient 1 POP1LOG Figure 21: Patient 1 POP2GOM
21 R. Injerd 21 Figure 22: Patient 2 POP1EXP Figure 25: Patient 2 POP2LOG Figure 23: Patient 2 POP2EXP Figure 26: Patient 2 POP1GOM Figure 24: Patient 2 POP1LOG Figure 27: Patient 2 POP2GOM
22 R. Injerd 22 Figure 28: Patient 3 POP1EXP Figure 31: Patient 3 POP2LOG Figure 29: Patient 3 POP2EXP Figure 32: Patient 3 POP1GOM Figure 30: Patient 3 POP1LOG Figure 33: Patient 3 POP2GOM
23 R. Injerd 23 Figure 34: Patient 4 POP1EXP Figure 37: Patient 4 POP2LOG Figure 35: Patient 4 POP2EXP Figure 38: Patient 4 POP1GOM Figure 36: Patient 4 POP1LOG Figure 39: Patient 4 POP2GOM
24 R. Injerd 24 Figure 40: Patient 5 POP1EXP Figure 43: Patient 5 POP2LOG Figure 41: Patient 5 POP2EXP Figure 44: Patient 5 POP1GOM Figure 42: Patient 5 POP1LOG Figure 45: Patient 5 POP2GOM
25 R. Injerd 25 Figure 46: Patient 6 POP1EXP Figure 49: Patient 6 POP2LOG Figure 47: Patient 6 POP2EXP Figure 50: Patient 6 POP1GOM Figure 48: Patient 6 POP1LOG Figure 51: Patient 6 POP2GOM
26 R. Injerd 26 Figure 52: Patient 7 POP1EXP Figure 55: Patient 7 POP2LOG Figure 53: Patient 7 POP2EXP Figure 56: Patient 7 POP1GOM Figure 54: Patient 7 POP1LOG Figure 57: Patient 7 POP2GOM
27 R. Injerd 27 Figure 58: Patient 8 POP1EXP Figure 61: Patient 8 POP2LOG Figure 59: Patient 8 POP2EXP Figure 62: Patient 8 POP1GOM Figure 60: Patient 8 POP1LOG Figure 63: Patient 8 POP2GOM
28 R. Injerd 28 Figure 64: Patient 9 POP1EXP Figure 67: Patient 9 POP2LOG Figure 65: Patient 9 POP2EXP Figure 68: Patient 9 POP1GOM Figure 66: Patient 9 POP1LOG Figure 69: Patient 9 POP2GOM
29 R. Injerd 29 Figure 70: Patient 10 POP1EXP Figure 73: Patient 10 POP2LOG Figure 71: Patient 10 POP2EXP Figure 74: Patient 10 POP1GOM Figure 72: Patient 10 POP1LOG Figure 75: Patient 10 POP2GOM
30 R. Injerd 30 Figure 76: Patient 11 POP1EXP Figure 79: Patient 11 POP2LOG Figure 77: Patient 11 POP2EXP Figure 80: Patient 11 POP1GOM Figure 78: Patient 11 POP1LOG Figure 81: Patient 11 POP2GOM
31 R. Injerd 31 References [1] Araujo RP, McElwain LS. (2004). A History of the Study of Solid Tumour Growth: The Contribution of Mathematical Modelling. Bulletin of Mathematical Biology, 66, [2] Chvestov A V, Palta J J, Nagata Y (2008). Time-dependent cell disintegration kinetics in lung tumors after irradiation. Physics in Medicine & Biology, 53, [3] Enderling H, Chaplain M. (2013). Mathematical Modelling of Tumour Growth and Treatment. Current pharmaceutical design, November, [4] Jain P, Baker A, Distefano G, Scott AJ, Webster GJ, Hatton MQ. (2013).Stereotactic ablative radiotherapy in the UK: current status and developments. BR J Radiol, 86, [5] Tariq I, Humbert-Vidan L, Chen T, South C, Ezhil V, Kirkby N, Jena R, Nisbet A. (2014). Mathematical modelling of tumour volume dynamics in response to stereotactic ablative radiotherapy fo non-small cell lung cancer. Physics in Medicine & Biology, 60,
Mathematical models of tumor growth and radiation response
Mathematical models of tumor growth and radiation response Yoichi Watanabe, Ph.D. Department of Radiation Oncology University of Minnesota North Central Chapter of AAPM Meeting, Wisconsin Dells, WI April
More informationA mathematical model of tumor growth and its response to single irradiation
Watanabe et al. Theoretical Biology and Medical Modelling (2016) 13:6 DOI 10.1186/s12976-016-0032-7 RESEARCH Open Access A mathematical model of tumor growth and its response to single irradiation Yoichi
More informationChapter 11: Advanced Remedial Measures. Weighted Least Squares (WLS)
Chapter : Advanced Remedial Measures Weighted Least Squares (WLS) When the error variance appears nonconstant, a transformation (of Y and/or X) is a quick remedy. But it may not solve the problem, or it
More informationMathematics Meets Oncology
.. Mathematics Meets Oncology Mathematical Oncology Philippe B. Laval Kennesaw State University November 12, 2011 Philippe B. Laval (Kennesaw State University)Mathematics Meets Oncology November 12, 2011
More informationMathematical Modeling of PDGF-Driven Glioblastoma Reveals Optimized Radiation Dosing Schedules
Mathematical Modeling of PDGF-Driven Glioblastoma Reveals Optimized Radiation Dosing Schedules Kevin Leder, Ken Pittner, Quincey LaPlant, Dolores Hambardzumyan, Brian D. Ross, Timothy A. Chan, Eric C.
More informationSimple Linear Regression the model, estimation and testing
Simple Linear Regression the model, estimation and testing Lecture No. 05 Example 1 A production manager has compared the dexterity test scores of five assembly-line employees with their hourly productivity.
More informationEC352 Econometric Methods: Week 07
EC352 Econometric Methods: Week 07 Gordon Kemp Department of Economics, University of Essex 1 / 25 Outline Panel Data (continued) Random Eects Estimation and Clustering Dynamic Models Validity & Threats
More informationSensitivity, specicity, ROC
Sensitivity, specicity, ROC Thomas Alexander Gerds Department of Biostatistics, University of Copenhagen 1 / 53 Epilog: disease prevalence The prevalence is the proportion of cases in the population today.
More informationMultiple Bivariate Gaussian Plotting and Checking
Multiple Bivariate Gaussian Plotting and Checking Jared L. Deutsch and Clayton V. Deutsch The geostatistical modeling of continuous variables relies heavily on the multivariate Gaussian distribution. It
More informationRisk of secondary cancer induced by radiotherapy
Risk of secondary cancer induced by radiotherapy Iuliana Toma-Dasu Medical Radiation Physics Stockholm University and Karolinska Institutet Radiation - the two-edged sword Risk of secondary cancer induced
More informationResearch Article A Mathematical Model of Tumor Volume Changes during Radiotherapy
The Scientific World Journal Volume 203, Article ID 8070, 5 pages http://dx.doi.org/0.55/203/8070 Research Article A Mathematical Model of Tumor Volume Changes during Radiotherapy Ping Wang and Yuanming
More informationCHAPTER TWO MECHANISMS OF RADIATION EFFECTS
10-2 densely ionizing radiation CHAPTER TWO MECHANISMS OF RADIATION EFFECTS 2.0 INTRODUCTION Cell survival curves describe the relationship between the fractional survival, S, of a population of radiated
More information11/18/2013. Correlational Research. Correlational Designs. Why Use a Correlational Design? CORRELATIONAL RESEARCH STUDIES
Correlational Research Correlational Designs Correlational research is used to describe the relationship between two or more naturally occurring variables. Is age related to political conservativism? Are
More informationLec 02: Estimation & Hypothesis Testing in Animal Ecology
Lec 02: Estimation & Hypothesis Testing in Animal Ecology Parameter Estimation from Samples Samples We typically observe systems incompletely, i.e., we sample according to a designed protocol. We then
More informationOverview of Non-Parametric Statistics
Overview of Non-Parametric Statistics LISA Short Course Series Mark Seiss, Dept. of Statistics April 7, 2009 Presentation Outline 1. Homework 2. Review of Parametric Statistics 3. Overview Non-Parametric
More informationHeterogeneity and statistical signi"cance in meta-analysis: an empirical study of 125 meta-analyses -
STATISTICS IN MEDICINE Statist. Med. 2000; 19: 1707}1728 Heterogeneity and statistical signi"cance in meta-analysis: an empirical study of 125 meta-analyses - Eric A. Engels *, Christopher H. Schmid, Norma
More informationModelling and Bayesian adaptive prediction of individual patients tumour volume
Modelling and Bayesian adaptive prediction of individual patients tumour volume change during radiotherapy Imran Tariq 1, Tao Chen 1,*, Norman F. Kirkby 2, Rajesh Jena 3 1 Department of Chemical and Process
More informationRanking radiotherapy treatment plans: physical or biological objectives?
Ranking radiotherapy treatment plans: physical or biological objectives? Martin Ebert Department of Radiation Oncology, Sir Charles Gairdner Hospital, Western Australia, Australia Background. The ranking
More informationChapters from Clinical Oncology
Chapters from Clinical Oncology Lecture notes University of Szeged Faculty of Medicine Department of Oncotherapy 2012. 1 RADIOTHERAPY Technical aspects Dr. Elemér Szil Introduction There are three possibilities
More informationProject for Math. 224 DETECTION OF DIABETES
Project for Math. 224 DETECTION OF DIABETES Diabetes is a disease of metabolism which is characterized by too much sugar in the blood and urine. Because of the lack of insulin (a hormone), the patient
More informationPitfalls in Linear Regression Analysis
Pitfalls in Linear Regression Analysis Due to the widespread availability of spreadsheet and statistical software for disposal, many of us do not really have a good understanding of how to use regression
More informationClinical Implementation of patient-specific dosimetry, comparison with absorbed fraction-based method
Clinical Implementation of patient-specific dosimetry, comparison with absorbed fraction-based method George Sgouros, Ph.D. Russell H. Morgan Dept of Radiology & Radiological Science Johns Hopkins University,
More informationDisclosure SBRT. SBRT for Spinal Metastases 5/2/2010. No conflicts of interest. Overview
Stereotactic Body Radiotherapy (SBRT) for Recurrent Spine Tumors Arjun Sahgal M.D., F.R.C.P.C. Assistant Professor Princess Margaret Hospital Sunnybrook Health Sciences Center University of Toronto Department
More informationBayesian Confidence Intervals for Means and Variances of Lognormal and Bivariate Lognormal Distributions
Bayesian Confidence Intervals for Means and Variances of Lognormal and Bivariate Lognormal Distributions J. Harvey a,b, & A.J. van der Merwe b a Centre for Statistical Consultation Department of Statistics
More informationMathematical Modeling of the Role of Survivin on Dedifferentiation and Radioresistance in Cancer
Mathematical Modeling of the Role of Survivin on Dedifferentiation and Radioresistance in Cancer Adam Rhodes, Thomas Hillen April 30, 2016 Centre for Mathematical Biology, University of Alberta, Canada
More informationA Brief (very brief) Overview of Biostatistics. Jody Kreiman, PhD Bureau of Glottal Affairs
A Brief (very brief) Overview of Biostatistics Jody Kreiman, PhD Bureau of Glottal Affairs What We ll Cover Fundamentals of measurement Parametric versus nonparametric tests Descriptive versus inferential
More informationDr. Kelly Bradley Final Exam Summer {2 points} Name
{2 points} Name You MUST work alone no tutors; no help from classmates. Email me or see me with questions. You will receive a score of 0 if this rule is violated. This exam is being scored out of 00 points.
More information11/24/2017. Do not imply a cause-and-effect relationship
Correlational research is used to describe the relationship between two or more naturally occurring variables. Is age related to political conservativism? Are highly extraverted people less afraid of rejection
More informationA kinetic model of tumor growth and its radiation response with an application to Gamma Knife stereotactic radiosurgery
A kinetic model of tumor growth and its radiation response with an application to Gamma Knife stereotactic radiosurgery Yoichi Watanabe 1, Erik L. Dahlman 1, Kevin Z. Leder 2, and Susanta K. Hui 1 1 Department
More informationInvestigating mathematical models of immunointeractions with early-stage cancer under an agent-based modelling perspective
Figueredo et al. BMC Bioinformatics, (Suppl ):S http://www.biomedcentral.com/7-5//s/s PROCEEDINGS Open Access Investigating mathematical models of immunointeractions with early-stage cancer under an agent-based
More informationCentre for Education Research and Policy
THE EFFECT OF SAMPLE SIZE ON ITEM PARAMETER ESTIMATION FOR THE PARTIAL CREDIT MODEL ABSTRACT Item Response Theory (IRT) models have been widely used to analyse test data and develop IRT-based tests. An
More informationMostly Harmless Simulations? On the Internal Validity of Empirical Monte Carlo Studies
Mostly Harmless Simulations? On the Internal Validity of Empirical Monte Carlo Studies Arun Advani and Tymon Sªoczy«ski 13 November 2013 Background When interested in small-sample properties of estimators,
More informationWorkplace smoking ban eects in an heterogenous smoking population
Workplace smoking ban eects in an heterogenous smoking population Workshop IRDES June 2010 1 Introduction 2 Summary of ndings 3 General population analysis The French ban had no impact on overall smoking
More informationModel-free detection of synchrony in neuronal spike trains, with an application to primate somatosensory cortex
Neurocomputing 32}33 (2000) 1103}1108 Model-free detection of synchrony in neuronal spike trains, with an application to primate somatosensory cortex A. Roy, P.N. Steinmetz, K.O. Johnson, E. Niebur* Krieger
More informationPerformance of Median and Least Squares Regression for Slightly Skewed Data
World Academy of Science, Engineering and Technology 9 Performance of Median and Least Squares Regression for Slightly Skewed Data Carolina Bancayrin - Baguio Abstract This paper presents the concept of
More informationMathematics for Infectious Diseases; Deterministic Models: A Key
Manindra Kumar Srivastava *1 and Purnima Srivastava 2 ABSTRACT The occurrence of infectious diseases was the principle reason for the demise of the ancient India. The main infectious diseases were smallpox,
More informationPET in Radiation Therapy. Outline. Tumor Segmentation in PET and in Multimodality Images for Radiation Therapy. 1. Tumor segmentation in PET
Tumor Segmentation in PET and in Multimodality Images for Radiation Therapy Wei Lu, Ph.D. Department of Radiation Oncology Mallinckrodt Institute of Radiology Washington University in St. Louis Outline
More informationANOVA in SPSS (Practical)
ANOVA in SPSS (Practical) Analysis of Variance practical In this practical we will investigate how we model the influence of a categorical predictor on a continuous response. Centre for Multilevel Modelling
More informationDesigning and Implementing a Model for Ethical Decision Making
Designing and Implementing a Model for Ethical Decision Making 1 Introductory Remarks One of us [GR] has been involved in developing a tool for computing the "relative evil" of pairs of military actions
More informationEcological Statistics
A Primer of Ecological Statistics Second Edition Nicholas J. Gotelli University of Vermont Aaron M. Ellison Harvard Forest Sinauer Associates, Inc. Publishers Sunderland, Massachusetts U.S.A. Brief Contents
More informationMathematical Microbiologists: Why we have to return to our square roots to uncover uncertainty (of measurement) in Quantitative PCR (qpcr)
Mathematical Microbiologists: Why we have to return to our square roots to uncover uncertainty (of measurement) in Quantitative PCR (qpcr) J. Ian Stuart George Zahariadis Marina Salvadori Objectives 1.
More informationAppendix III Individual-level analysis
Appendix III Individual-level analysis Our user-friendly experimental interface makes it possible to present each subject with many choices in the course of a single experiment, yielding a rich individual-level
More informationRadiation Damage Comparison between Intensity Modulated Radiotherapy (IMRT) and Field-in-field (FIF) Technique In Breast Cancer Treatments
Radiation Damage Comparison between Intensity Modulated Radiotherapy () and Field-in-field (FIF) Technique In Breast Cancer Treatments Huisi Ai 1 and Hualin Zhang 2 1. Department of Radiation Oncology,
More informationMonte Carlo Analysis of Univariate Statistical Outlier Techniques Mark W. Lukens
Monte Carlo Analysis of Univariate Statistical Outlier Techniques Mark W. Lukens This paper examines three techniques for univariate outlier identification: Extreme Studentized Deviate ESD), the Hampel
More informationCopyright 2008 Society of Photo Optical Instrumentation Engineers. This paper was published in Proceedings of SPIE, vol. 6915, Medical Imaging 2008:
Copyright 28 Society of Photo Optical Instrumentation Engineers. This paper was published in Proceedings of SPIE, vol. 695, Medical Imaging 28: Computer Aided Diagnosis and is made available as an electronic
More informationKnowledge-Based IMRT Treatment Planning for Prostate Cancer: Experience with 101. Cases from Duke Clinic. Deon Martina Dick
Knowledge-Based IMRT Treatment Planning for Prostate Cancer: Experience with 101 Cases from Duke Clinic by Deon Martina Dick Department of Medical Physics Duke University Date: Approved: Joseph Lo, Chair
More informationBayesian and Frequentist Approaches
Bayesian and Frequentist Approaches G. Jogesh Babu Penn State University http://sites.stat.psu.edu/ babu http://astrostatistics.psu.edu All models are wrong But some are useful George E. P. Box (son-in-law
More informationThe impact of pre-selected variance inflation factor thresholds on the stability and predictive power of logistic regression models in credit scoring
Volume 31 (1), pp. 17 37 http://orion.journals.ac.za ORiON ISSN 0529-191-X 2015 The impact of pre-selected variance inflation factor thresholds on the stability and predictive power of logistic regression
More informationSTATISTICS AND RESEARCH DESIGN
Statistics 1 STATISTICS AND RESEARCH DESIGN These are subjects that are frequently confused. Both subjects often evoke student anxiety and avoidance. To further complicate matters, both areas appear have
More informationPart I. Boolean modelling exercises
Part I. Boolean modelling exercises. Glucose repression of Icl in yeast In yeast Saccharomyces cerevisiae, expression of enzyme Icl (isocitrate lyase-, involved in the gluconeogenesis pathway) is important
More informationResearch Article Impact of Dose and Sensitivity Heterogeneity on TCP
Computational and Mathematical Methods in Medicine, Article ID 182935, 7 pages http://dx.doi.org/10.1155/2014/182935 Research Article Impact of Dose and Sensitivity Heterogeneity on TCP Kristin Wiklund,
More informationPsychology Research Process
Psychology Research Process Logical Processes Induction Observation/Association/Using Correlation Trying to assess, through observation of a large group/sample, what is associated with what? Examples:
More informationHeroin Epidemic Models
Heroin Epidemic Models icholas A. Battista Intro to Math Biology Project School of Mathematical Sciences, Rochester Institute of Technology, 85 Lomb Memorial Drive, Rochester, Y 14623-5603, USA May 21,
More informationImproving ecological inference using individual-level data
STATISTICS IN MEDICINE Statist. Med. 2006; 25:2136 2159 Published online 11 October 2005 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/sim.2370 Improving ecological inference using individual-level
More informationEvaluation of Monaco treatment planning system for hypofractionated stereotactic volumetric arc radiotherapy of multiple brain metastases
Evaluation of Monaco treatment planning system for hypofractionated stereotactic volumetric arc radiotherapy of multiple brain metastases CASE STUDY Institution: Odette Cancer Centre Location: Sunnybrook
More informationPairwise- and mixed treatment comparison models in multicriteria benefit-risk analysis
Master Thesis Pairwise- and mixed treatment comparison models in multicriteria benefit-risk analysis Author: Jing Zhao jingzhao@kth.se Supervisors : Henrik Boström, Department of Computer and Systems Sciences,
More informationAdaptive modelling of tumour volume dynamics under. radiotherapy
Adaptive modelling of tumour volume dynamics under radiotherapy Thesis by Imran Tariq Submitted for the Degree of Doctor of Philosophy Department of Chemical and Process Engineering Faculty of Engineering
More informationSelection of Linking Items
Selection of Linking Items Subset of items that maximally reflect the scale information function Denote the scale information as Linear programming solver (in R, lp_solve 5.5) min(y) Subject to θ, θs,
More informationJSM Survey Research Methods Section
Methods and Issues in Trimming Extreme Weights in Sample Surveys Frank Potter and Yuhong Zheng Mathematica Policy Research, P.O. Box 393, Princeton, NJ 08543 Abstract In survey sampling practice, unequal
More informationDosimetric Analysis of 3DCRT or IMRT with Vaginal-cuff Brachytherapy (VCB) for Gynaecological Cancer
Dosimetric Analysis of 3DCRT or IMRT with Vaginal-cuff Brachytherapy (VCB) for Gynaecological Cancer Tan Chek Wee 15 06 2016 National University Cancer Institute, Singapore Clinical Care Education Research
More informationOut-of-field doses of CyberKnife in stereotactic radiotherapy of prostate cancer patients
EURAMED,MELODI and EURADOS Workshop, Helsinki 19.-20.4.2017 Out-of-field doses of CyberKnife in stereotactic radiotherapy of prostate cancer patients Jan Seppälä, Chief Physicist Tuomas Virén, Medical
More informationList of Figures. List of Tables. Preface to the Second Edition. Preface to the First Edition
List of Figures List of Tables Preface to the Second Edition Preface to the First Edition xv xxv xxix xxxi 1 What Is R? 1 1.1 Introduction to R................................ 1 1.2 Downloading and Installing
More informationResearch Methods in Forest Sciences: Learning Diary. Yoko Lu December Research process
Research Methods in Forest Sciences: Learning Diary Yoko Lu 285122 9 December 2016 1. Research process It is important to pursue and apply knowledge and understand the world under both natural and social
More informationAssignment #6. Chapter 10: 14, 15 Chapter 11: 14, 18. Due tomorrow Nov. 6 th by 2pm in your TA s homework box
Assignment #6 Chapter 10: 14, 15 Chapter 11: 14, 18 Due tomorrow Nov. 6 th by 2pm in your TA s homework box Assignment #7 Chapter 12: 18, 24 Chapter 13: 28 Due next Friday Nov. 13 th by 2pm in your TA
More informationEpidemiological Model of HIV/AIDS with Demographic Consequences
Advances in Applied Mathematical Biosciences. ISSN 2248-9983 Volume 5, Number 1 (2014), pp. 65-74 International Research Publication House http://www.irphouse.com Epidemiological Model of HIV/AIDS with
More informationDOSE DISTRIBUTION ANALYZE OF THE BODY STI USED MONTE CARLO METHOD
Proceedings of the Tenth EGS4 Users' Meeting in Japan, KEK Proceedings 2002-18, p.65-73 DOSE DISTRIBUTION ANALYZE OF THE BODY STI USED MONTE CARLO METHOD N. Tohyama, H. Saitoh, T. Fujisaki 1, S. Abe 1
More informationReview of TG-186 recommendations
Review of TG-186 recommendations Implementation of advanced brachytherapy dose calculation algorithms beyond TG-43 Rowan M. Thomson Carleton Laboratory for Radiotherapy Physics Carleton University Ottawa
More informationCOMPARING SEVERAL DIAGNOSTIC PROCEDURES USING THE INTRINSIC MEASURES OF ROC CURVE
DOI: 105281/zenodo47521 Impact Factor (PIF): 2672 COMPARING SEVERAL DIAGNOSTIC PROCEDURES USING THE INTRINSIC MEASURES OF ROC CURVE Vishnu Vardhan R* and Balaswamy S * Department of Statistics, Pondicherry
More informationLINEAR REGRESSION FOR BIVARIATE CENSORED DATA VIA MULTIPLE IMPUTATION
STATISTICS IN MEDICINE Statist. Med. 18, 3111} 3121 (1999) LINEAR REGRESSION FOR BIVARIATE CENSORED DATA VIA MULTIPLE IMPUTATION WEI PAN * AND CHARLES KOOPERBERG Division of Biostatistics, School of Public
More informationMotivation Empirical models Data and methodology Results Discussion. University of York. University of York
Healthcare Cost Regressions: Going Beyond the Mean to Estimate the Full Distribution A. M. Jones 1 J. Lomas 2 N. Rice 1,2 1 Department of Economics and Related Studies University of York 2 Centre for Health
More informationInvestigating the robustness of the nonparametric Levene test with more than two groups
Psicológica (2014), 35, 361-383. Investigating the robustness of the nonparametric Levene test with more than two groups David W. Nordstokke * and S. Mitchell Colp University of Calgary, Canada Testing
More informationA Learning Method of Directly Optimizing Classifier Performance at Local Operating Range
A Learning Method of Directly Optimizing Classifier Performance at Local Operating Range Lae-Jeong Park and Jung-Ho Moon Department of Electrical Engineering, Kangnung National University Kangnung, Gangwon-Do,
More informationBLADDER RADIOTHERAPY PLANNING DOCUMENT
A 2X2 FACTORIAL RANDOMISED PHASE III STUDY COMPARING STANDARD VERSUS REDUCED VOLUME RADIOTHERAPY WITH AND WITHOUT SYNCHRONOUS CHEMOTHERAPY IN MUSCLE INVASIVE BLADDER CANCER (ISRCTN 68324339) BLADDER RADIOTHERAPY
More informationMathematical Model of Vaccine Noncompliance
Valparaiso University ValpoScholar Mathematics and Statistics Faculty Publications Department of Mathematics and Statistics 8-2016 Mathematical Model of Vaccine Noncompliance Alex Capaldi Valparaiso University
More informationAutomatic Definition of Planning Target Volume in Computer-Assisted Radiotherapy
Automatic Definition of Planning Target Volume in Computer-Assisted Radiotherapy Angelo Zizzari Department of Cybernetics, School of Systems Engineering The University of Reading, Whiteknights, PO Box
More informationA Comparison of Robust and Nonparametric Estimators Under the Simple Linear Regression Model
Nevitt & Tam A Comparison of Robust and Nonparametric Estimators Under the Simple Linear Regression Model Jonathan Nevitt, University of Maryland, College Park Hak P. Tam, National Taiwan Normal University
More informationImaging of Scattered Radiation for Real Time Tracking of Tumor Motion During Lung SBRT
Imaging of Scattered Radiation for Real Time Tracking of Tumor Motion During Lung SBRT April 25 nd, 2015 Lung Cancer Lung cancer is the most lethal cancer: Over 224,000 new diagnoses in the U.S. predicted
More informationRESPONSE SURFACE MODELING AND OPTIMIZATION TO ELUCIDATE THE DIFFERENTIAL EFFECTS OF DEMOGRAPHIC CHARACTERISTICS ON HIV PREVALENCE IN SOUTH AFRICA
RESPONSE SURFACE MODELING AND OPTIMIZATION TO ELUCIDATE THE DIFFERENTIAL EFFECTS OF DEMOGRAPHIC CHARACTERISTICS ON HIV PREVALENCE IN SOUTH AFRICA W. Sibanda 1* and P. Pretorius 2 1 DST/NWU Pre-clinical
More informationEssentials in Bioassay Design and Relative Potency Determination
BioAssay SCIENCES A Division of Thomas A. Little Consulting Essentials in Bioassay Design and Relative Potency Determination Thomas A. Little Ph.D. 2/29/2016 President/CEO BioAssay Sciences 12401 N Wildflower
More informationPossible Consequences of Inhomogeneous Suborgan Distribution of Dose and the Linear No-Threshold Dose-Effect Relationship
Possible Consequences of Inhomogeneous Suborgan Distribution of Dose and the Linear No-Threshold Dose-Effect Relationship Balázs G. Madas, Imre Balásházy Centre for Energy Research, Hungarian Academy of
More informationImage Guided Stereotactic Radiotherapy of the Lung
Image Guided Stereotactic Radiotherapy of the Lung Jamie Marie Harris, MS DABR Avera McKennan Radiation Oncology September 25, 2015 Stereotactic Body Radiotherapy - Clinical Dose/Fractionation - Normal
More informationAppendix (unedited, as supplied by the authors)
Appendix (unedited, as supplied by the authors) Details of the mathematical model The model used for the simulations in this manuscript was adapted from a previously described model of HCV transmission
More informationThe Fairest of Them All: Using Variations of Beta-Binomial Distributions to Investigate Robust Scoring Methods
The Fairest of Them All: Using Variations of Beta-Binomial Distributions to Investigate Robust Scoring Methods Mary Good - Bluffton University Christopher Kinson - Albany State University Karen Nielsen
More informationUnderstandable Statistics
Understandable Statistics correlated to the Advanced Placement Program Course Description for Statistics Prepared for Alabama CC2 6/2003 2003 Understandable Statistics 2003 correlated to the Advanced Placement
More informationT-Statistic-based Up&Down Design for Dose-Finding Competes Favorably with Bayesian 4-parameter Logistic Design
T-Statistic-based Up&Down Design for Dose-Finding Competes Favorably with Bayesian 4-parameter Logistic Design James A. Bolognese, Cytel Nitin Patel, Cytel Yevgen Tymofyeyef, Merck Inna Perevozskaya, Wyeth
More informationNon-classical radiobiology relevant to high-doses per fraction
Non-classical radiobiology relevant to high-doses per fraction Michael Joiner Wayne State University Radiation Oncology Detroit, Michigan joinerm@wayne.edu Why reconsider high dose fractions? Because we
More informationarxiv: v2 [math.oc] 15 Nov 2016
Optimizing chemoradiotherapy to target multi-site metastatic disease and tumor growth arxiv:1603.00349v2 [math.oc] 15 Nov 2016 Hamidreza Badri, 1 Ehsan Salari, 2 Yoichi Watanabe, 3 Kevin Leder 1 1 Department
More informationThe Impact of Cobalt-60 Source Age on Biologically Effective Dose in Gamma Knife Thalamotomy
The Impact of Cobalt-60 Source Age on Biologically Effective Dose in Gamma Knife Thalamotomy BH Kann, JB Yu, J Bond, C Loiselle, VL Chiang, RS Bindra, JL Gerrard, DJ Carlson Leksell Gamma Knife Society
More informationStereotactic MR-guided adaptive radiation therapy (SMART) for locally advanced pancreatic tumors
Stereotactic MR-guided adaptive radiation therapy (SMART) for locally advanced pancreatic tumors Anna Bruynzeel, Radiation Oncologist VU University Medical Center, Amsterdam, The Netherlands Current standard
More informationUPPER MIDWEST MARKETING AREA ANALYSIS OF COMPONENT LEVELS AND SOMATIC CELL COUNT IN INDIVIDUAL HERD MILK AT THE FARM LEVEL 2007
UPPER MIDWEST MARKETING AREA ANALYSIS OF COMPONENT LEVELS AND SOMATIC CELL COUNT IN INDIVIDUAL HERD MILK AT THE FARM LEVEL 2007 Staff Paper 08-01 Prepared by: Corey Freije 2008 Federal Milk Market Administrator
More informationVariable Dose Rate Dynamic Conformal Arc Therapy (DCAT) for SABR Lung: From static fields to dynamic arcs using Monaco 5.10
Variable Dose Rate Dynamic Conformal Arc Therapy (DCAT) for SABR Lung: From static fields to dynamic arcs using Monaco 5.10 Simon Goodall Radiation Oncology Physicist Genesis Care Western Australia Introduction
More informationLecture Outline. Biost 590: Statistical Consulting. Stages of Scientific Studies. Scientific Method
Biost 590: Statistical Consulting Statistical Classification of Scientific Studies; Approach to Consulting Lecture Outline Statistical Classification of Scientific Studies Statistical Tasks Approach to
More informationSummary Report for HIV Random Clinical Trial Conducted in
Summary Report for HIV Random Clinical Trial Conducted in 9-2014 H.T. Banks and Shuhua Hu Center for Research in Scientific Computation North Carolina State University Raleigh, NC 27695-8212 USA Eric Rosenberg
More information1 Introduction Synchronous ring of action potentials amongst multiple neurons is a phenomenon that has been observed in a wide range of neural systems
Model-free detection of synchrony in neuronal spike trains, with an application to primate somatosensory cortex A. Roy a, P. N. Steinmetz a 1, K. O. Johnson a, E. Niebur a 2 a Krieger Mind/Brain Institute,
More informationA Survey of Techniques for Optimizing Multiresponse Experiments
A Survey of Techniques for Optimizing Multiresponse Experiments Flavio S. Fogliatto, Ph.D. PPGEP / UFRGS Praça Argentina, 9/ Sala 402 Porto Alegre, RS 90040-020 Abstract Most industrial processes and products
More informationResults & Statistics: Description and Correlation. I. Scales of Measurement A Review
Results & Statistics: Description and Correlation The description and presentation of results involves a number of topics. These include scales of measurement, descriptive statistics used to summarize
More informationApplication of the concept of biologically effective dose (BED) to patients with Vestibular Schwannomas treated by radiosurgery
Jour. of Radiosurgery and SBRT, Vol. 2, pp. 257 271 Reprints available directly from the publisher Photocopying permitted by license only 2013 Old City Publishing, Inc. Published by license under the OCP
More informationKidane Tesfu Habtemariam, MASTAT, Principle of Stat Data Analysis Project work
1 1. INTRODUCTION Food label tells the extent of calories contained in the food package. The number tells you the amount of energy in the food. People pay attention to calories because if you eat more
More informationEmpirical Formula for Creating Error Bars for the Method of Paired Comparison
Empirical Formula for Creating Error Bars for the Method of Paired Comparison Ethan D. Montag Rochester Institute of Technology Munsell Color Science Laboratory Chester F. Carlson Center for Imaging Science
More informationChanging Paradigms in Radiotherapy
Changing Paradigms in Radiotherapy Marco van Vulpen, MD, PhD Mouldroomdag-2015 Towards the elimination of invasion 1 NIH opinion on the future of oncology Twenty-five years from now,i hope that we won
More information