Dynamic Disease Screening
|
|
- Neil Pope
- 5 years ago
- Views:
Transcription
1 Dynamic Disease Screening Peihua Qiu Department of Biostatistics University of Florida December 10, 2014, NCTS Lecture, Taiwan p.1/25
2 Motivating Example SHARe Framingham Heart Study of NHLBI. Many residents at Framingham MA were involved. Major risk factors of cardiovascular diseases: blood pressure, total cholesterol level (TCL), smoking, obesity,... Identify patients with irregular longitudinal patterns of the disease risk factors as early as possible. Disease early detection and prevention December 10, 2014, NCTS Lecture, Taiwan p.1/25 Dynamic screening (DS) problem
3 DS Problem DS problem is popular Many products (e.g., airplanes, cars) are checked regularly or occasionally about certain variables related to their quality and/or performance. If the observed values of a product are significantly worse than the values of a typical well-functioning product of the same age, then some adjustments or interventions should be made to avoid unpleasant consequences. December 10, 2014, NCTS Lecture, Taiwan p.2/25
4 Possible Statistical Methods Confidence interval of the mean responses by longitudinal data analysis. This method uses the cross-sectional comparison approach. It does not make use of all history data of a subject. It cannot detect a shift sequentially. December 10, 2014, NCTS Lecture, Taiwan p.3/25
5 Possible Statistical Methods Statistical process control (SPC) methods Monitor each subject sequentially Use all history data of the subject They monitor subjects separately and cannot compare different subjects The process mean and variance may not be constants even when the subject is IC December 10, 2014, NCTS Lecture, Taiwan p.4/25
6 Dynamic Screening System (DySS) Estimate regular longitudinal pattern from an IC dataset Standardize observations of a new subject to monitor Monitor the standardized observations by a control chart December 10, 2014, NCTS Lecture, Taiwan p.5/25
7 References Qiu, P., and Xiang, D. (2014), Dynamic screening system: an approach for dynamically identifying irregular individuals, Technometrics, 56, Qiu, P., and Xiang, D. (2014), Surveillance of cardiovascular diseases using a multivariate dynamic screening system, revised for Statistics in Medicine. Qiu, P., Zi, X., and Zou, C. (2014), Dynamic nonparametric curve monitoring, submitted. Li, J., and Qiu, P. (2014), Nonparametric dynamic screening system for monitoring correlated longitudinal data, submitted. Xiang, D., Qiu, P., and Pu, X. (2013), Nonparametric regression analysis of multivariate longitudinal data, Statistica Sinica, 23, December 10, 2014, NCTS Lecture, Taiwan p.6/25
8 MDySS Qiu, P., and Xiang, D. (2014) Multivariate dynamic screening system (MDySS) December 10, 2014, NCTS Lecture, Taiwan p.7/25
9 Estimate regular longitudinal pattern IC data: observations of m well-functioning subjects For i = 1,2,...,m,j = 1,2,...,J i,t ij [0,1], y(t ij ) = µ(t ij )+ε(t ij ) y(t ij ) = (y 1 (t ij ),...,y q (t ij )) µ(t ij ) = (µ 1 (t ij ),...,µ q (t ij )) Regular pattern: µ(t) and Σ(s, t) = Cov(y(s), y(t)) Xiang, Qiu, and Pu (2013): estimation of µ(t) December 10, 2014, NCTS Lecture, Taiwan p.8/25 and Σ(s,t).
10 Standardize Observations New subject s y values are observed at t 1,t 2,... over [0,1]. When s/he is IC, y(t j) = µ(t j)+σ 1 2 (t j,t j)ǫ(t j) Standardized observations: ǫ(t j) = Σ ( ) 1 2 (t j,t j) y(t j) µ(t j; Σ) December 10, 2014, NCTS Lecture, Taiwan p.9/25
11 A Note By using the standardized observations of the new subject, we have actually compared its longitudinal pattern cross-sectionally with the estimated regular longitudinal pattern at the time points t 1,t 2,... December 10, 2014, NCTS Lecture, Taiwan p.10/25
12 Sequential Monitoring Zou and Qiu (2009): LASSO-based MEWMA chart MEWMA statistic U j = λ L ǫ(t j)+(1 λ L )U j 1 min α R q(u j α) (U j α)+ γ k q l=1 α l U jl Q j = max k=1,...,q W j, γk E(W j, γk ) Var(Wj, γk ) > h L December 10, 2014, NCTS Lecture, Taiwan p.11/25
13 Performance Evaluation Performance measures: IC average run length ARL 0 OC average run length ARL 1 December 10, 2014, NCTS Lecture, Taiwan p.12/25
14 Performance Evaluation (Con d) If {t j,j = 1,2,...} are unequally spaced, ARL 0 and ARL 1 may not be appropriate Basic time unit ω: largest time unit that all observation times are integer multiples of ω Define n j = t j /ω, for j = 0,1,2,..., where n 0 = t 0 = 0. t j = n jω, for all j. December 10, 2014, NCTS Lecture, Taiwan p.13/25
15 Performance Evaluation (Con d) IC: If a signal is given at the sth observation time, then E(n s) measures the IC average time to signal (ATS), denoted as ATS 0. OC: If a shift occurs at the τth observation time and a signal is given at the sth observation time with s τ, then E(n s n τ) is the OC ATS, denoted as ATS 1. December 10, 2014, NCTS Lecture, Taiwan p.14/25
16 SHARe Framingham Heart Study m = 945 non-stroke patients (IC data) 27 stroke patients (new subjects) each patient was followed 7 times (i.e., J = 7) Four medical indices: systolic blood pressure (mmhg), diastolic blood pressure (mmhg), total cholesterol level (mg/100ml), and glucose level (mg/100ml) December 10, 2014, NCTS Lecture, Taiwan p.15/25
17 SHARe Framingham Heart Study (con d) Qj Patient 1 Patient 2 Patient 3 Patient 4 Qj Patient 5 Patient 6 Patient 7 Patient 8 Qj Patient 9 Patient 10 Patient 11 Patient 12 Qj Patient 13 Patient 14 Patient 15 Patient 16 Qj Patient 17 Patient 18 Patient 19 Patient 20 Qj Patient 21 Patient 22 Patient 23 Patient j Qj Patient 25 Patient 26 Patient 27 December 10, 2014, NCTS Lecture, Taiwan p.16/ j j j
18 SHARe Framingham Heart Study (con d) DySS approach: 26 out of 27 stroke patients got signals; 131 out of 945 non-stroke patients got signals. The average signal time is years. December 10, 2014, NCTS Lecture, Taiwan p.17/25
19 Dynamic Curve Monitoring Qiu, Zi, and Zou (2014) Model µ(t i )+σ(t i )ε(t i ), for t i [0,τ], y(t i )= µ(t i )+σ(t i )g(t i )+σ(t i )ε(t i ), for t i (τ,t], After the transformation {y(t i ) µ(t i )}/σ(t i ), y(t i )= { ε(ti ), for t i [0,τ], g(t i )+ε(t i ), for t i (τ,t]. H 0 : τ > T versus H 1 : τ [0,T] December 10, 2014, NCTS Lecture, Taiwan p.18/25
20 Test and Estimation of g(t i ) Loss function: Q(t m ;λ) argmin a R m i=1 m i=1 {y(t i ) g(t i )} 2 (1 λ) t m t i {y(t i ) a} 2 (1 λ) t m t i ĝ λ (t m ) = m i=1 w i(t m )y(t i )/ m i=1 w i(t m ) Q H1 (t m ;λ) = m i=1 {y(t i) ĝ(t i )} 2 w i (t m ) Q H0 (t m ;λ) = m i=1 {y(t i)} 2 w i (t m ) December 10, 2014, NCTS Lecture, Taiwan p.19/25
21 Test and Estimation of g(t i ) Weighted GLR test statistic (WGLR): W λ (t m ) =Q H0 (t m ;λ) Q H1 (t m ;λ) m = w i (t m ){2y(t i ) ĝ(t i )}ĝ(t i ). i=1 Recursive formulas: W λ (t m ) = w m 1 (t m )W λ (t m 1 )+{2y(t m ) ĝ(t m )}ĝ(t m ), ĝ(t m ) = {α m 1 ĝ(t m 1 )+y(t m )}/α m where α m = m i=1 w i(t m ) = w m 1 (t m )α m December 10, 2014, NCTS Lecture, Taiwan p.20/25
22 Test of g(t i ) (con d) Dynamic EWMA (DEWMA) chart: W λ(t m ) = {W λ (t m ) E λ (t m )}/ V λ (t m ) > L When the observation times are equally spaced, DEWMA is the conventional EWMA chart. Benefits: Accommodate the unequally spaced observation times by using the weights (1 λ) t m t i Wλ (t m) is robust when g(t) values change December 10, 2014, NCTS Lecture, Taiwan p.21/25 much over time
23 Simulation Results t [0,1000], d = 1 IC model: µ(t) = 1+0.3t 1/2, σ 2 (t) = µ 2 (t) OC models: (I) Step Shift: g(t) = δ, for t > τ (II) Quadratic Drift: g(t) = (t τ) 2 δ (III) Sine Drift: g(t) = sin(0.003π(t τ))δ December 10, 2014, NCTS Lecture, Taiwan p.22/25
24 Model (I) Model (II) Model (III) Y δ = 0 δ = 0.05 δ = 0.1 Y δ = 0 δ = 0.5 δ = 1 Y δ = 0 δ = 0.05 δ = t (a) t (b) t (c) Y True function local linear EWMA Time (d) Y Time (e) Y December 10, 2014, NCTS Lecture, Taiwan p.23/25 Time (f)
25 Model I (equally) Model II (equally) Model III (equally) log(ats) DEWMA(λ = 0.05) DEWMA(λ = 0.2) EWMA(λ = 0.05) EWMA(λ = 0.2) log(ats) log(ats) δ (a) Model I (random) δ (b) Model II (random) δ (c) Model III (random) log(ats) log(ats) log(ats) δ (d) δ (e) December 10, 2014, NCTS Lecture, δ Taiwan p.24/25 (f)
26 Future Research Autocorrelation Nonparametric charts Accommodation of covariates December 10, 2014, NCTS Lecture, Taiwan p.25/25
Modern Statistical Process Control Charts and Their Applications in Analyzing Big Data
Modern Statistical Process Control Charts and Their Applications in Analyzing Big Data (Part III: Profile Monitoring and Dynamic Screening) Peihua Qiu pqiu@ufl.edu Department of Biostatistics University
More informationConstruction of An Efficient Multivariate Dynamic Screening System. Jun Li a and Peihua Qiu b. Abstract
Construction of An Efficient Multivariate Dynamic Screening System Jun Li a and Peihua Qiu b a Department of Statistics, University of California at Riverside b Department of Biostatistics, University
More informationUnivariate Dynamic Screening System: An Approach For Identifying Individuals With Irregular Longitudinal Behavior. Abstract
Univariate Dynamic Screening System: An Approach For Identifying Individuals With Irregular Longitudinal Behavior Peihua Qiu 1 and Dongdong Xiang 2 1 School of Statistics, University of Minnesota 2 School
More informationRejoinder. Peihua Qiu Department of Biostatistics, University of Florida 2004 Mowry Road, Gainesville, FL 32610
Rejoinder Peihua Qiu Department of Biostatistics, University of Florida 2004 Mowry Road, Gainesville, FL 32610 I was invited to give a plenary speech at the 2017 Stu Hunter Research Conference in March
More informationCONTROL CHARTS FOR MULTIVARIATE NONLINEAR TIME SERIES
REVSTAT Statistical Journal Volume 13, Number, June 015, 131 144 CONTROL CHARTS FOR MULTIVARIATE NONLINEAR TIME SERIES Authors: Robert Garthoff Department of Statistics, European University, Große Scharrnstr.
More informationModule B1: Multivariate Process Control
Module B1: Multivariate Process Control Prof. Fugee Tsung Hong Kong University of Science and Technology Quality Lab: http://qlab.ielm.ust.hk I. Multivariate Shewhart chart WHY MULTIVARIATE PROCESS CONTROL
More informationRejoinder. 1 Phase I and Phase II Profile Monitoring. Peihua Qiu 1, Changliang Zou 2 and Zhaojun Wang 2
Rejoinder Peihua Qiu 1, Changliang Zou 2 and Zhaojun Wang 2 1 School of Statistics, University of Minnesota 2 LPMC and Department of Statistics, Nankai University, China We thank the editor Professor David
More informationTECHNICAL APPENDIX WITH ADDITIONAL INFORMATION ON METHODS AND APPENDIX EXHIBITS. Ten health risks in this and the previous study were
Goetzel RZ, Pei X, Tabrizi MJ, Henke RM, Kowlessar N, Nelson CF, Metz RD. Ten modifiable health risk factors are linked to more than one-fifth of employer-employee health care spending. Health Aff (Millwood).
More informationQuantile Regression Methods for Reference Growth Charts
Quantile Regression Methods for Reference Growth Charts 1 Roger Koenker University of Illinois at Urbana-Champaign ASA Workshop on Nonparametric Statistics Texas A&M, January 15, 2005 Based on joint work
More informationStat 579: Generalized Linear Models and Extensions
Stat 579: Generalized Linear Models and Extensions Linear Mixed Models for Longitudinal Data Yan Lu April, 2018, week 12 1 / 34 Correlated data multivariate observations clustered data repeated measurement
More informationWeighted Likelihood Ratio Chart for Statistical Monitoring of Queueing Systems
Weighted Likelihood Ratio Chart for Statistical Monitoring of Queueing Systems Dequan Qi 1, Zhonghua Li 2, Xuemin Zi 3, Zhaojun Wang 2 1 LPMC and School of Mathematical Sciences, Nankai University, Tianjin
More informationStatistical Process Control for Multivariate Categorical Processes
Statistical Process Control for Multivariate Categorical Processes Fugee Tsung The Hong Kong University of Science and Technology Fugee Tsung 1/27 Introduction Typical Control Charts Univariate continuous
More informationInference for correlated effect sizes using multiple univariate meta-analyses
Europe PMC Funders Group Author Manuscript Published in final edited form as: Stat Med. 2016 April 30; 35(9): 1405 1422. doi:10.1002/sim.6789. Inference for correlated effect sizes using multiple univariate
More informationMisclassification Rates in Hypertension Diagnosis due to Measurement Errors
Misclassification Rates in Hypertension Diagnosis due to Measurement Errors Camila Friedman-Gerlicz Claremont McKenna College cgerlicz0@cmc.edu and Isaiah Lilly California State University at Sacramento
More informationCorrection for classical covariate measurement error and extensions to life-course studies
Correction for classical covariate measurement error and extensions to life-course studies Jonathan William Bartlett A thesis submitted to the University of London for the degree of Doctor of Philosophy
More informationEstimating Optimal Dynamic Treatment Regimes from Clustered Data
Estimating Optimal Dynamic Treatment Regimes from Clustered Data Bibhas Chakraborty Department of Biostatistics, Columbia University bc2425@columbia.edu Society for Clinical Trials Annual Meetings Boston,
More informationDistribution-Free Monitoring of Univariate Processes. Peihua Qiu 1 and Zhonghua Li 1,2. Abstract
Distribution-Free Monitoring of Univariate Processes Peihua Qiu 1 and Zhonghua Li 1,2 1 School of Statistics, University of Minnesota, USA 2 LPMC and Department of Statistics, Nankai University, China
More informationTime Series. Anthony Davison. c
Series Anthony Davison c 2008 http://stat.epfl.ch Periodogram 76 Motivation............................................................ 77 Lutenizing hormone data..................................................
More informationECAS Summer Course. Quantile Regression for Longitudinal Data. Roger Koenker University of Illinois at Urbana-Champaign
ECAS Summer Course 1 Quantile Regression for Longitudinal Data Roger Koenker University of Illinois at Urbana-Champaign La Roche-en-Ardennes: September 2005 Part I: Penalty Methods for Random Effects Part
More informationThe Robustness of the Multivariate EWMA Control Chart
The Robustness of the Multivariate EWMA Control Chart Zachary G. Stoumbos, Rutgers University, and Joe H. Sullivan, Mississippi State University Joe H. Sullivan, MSU, MS 39762 Key Words: Elliptically symmetric,
More informationCorrelation and Simple Linear Regression
Correlation and Simple Linear Regression Sasivimol Rattanasiri, Ph.D Section for Clinical Epidemiology and Biostatistics Ramathibodi Hospital, Mahidol University E-mail: sasivimol.rat@mahidol.ac.th 1 Outline
More informationPower and Sample Size Calculations with the Additive Hazards Model
Journal of Data Science 10(2012), 143-155 Power and Sample Size Calculations with the Additive Hazards Model Ling Chen, Chengjie Xiong, J. Philip Miller and Feng Gao Washington University School of Medicine
More informationLecture 4 Multiple linear regression
Lecture 4 Multiple linear regression BIOST 515 January 15, 2004 Outline 1 Motivation for the multiple regression model Multiple regression in matrix notation Least squares estimation of model parameters
More informationBINF 702 SPRING Chapter 8 Hypothesis Testing: Two-Sample Inference. BINF702 SPRING 2014 Chapter 8 Hypothesis Testing: Two- Sample Inference 1
BINF 702 SPRING 2014 Chapter 8 Hypothesis Testing: Two-Sample Inference Two- Sample Inference 1 A Poster Child for two-sample hypothesis testing Ex 8.1 Obstetrics In the birthweight data in Example 7.2,
More informationData Mining Stat 588
Data Mining Stat 588 Lecture 9: Basis Expansions Department of Statistics & Biostatistics Rutgers University Nov 01, 2011 Regression and Classification Linear Regression. E(Y X) = f(x) We want to learn
More informationMarginal Structural Cox Model for Survival Data with Treatment-Confounder Feedback
University of South Carolina Scholar Commons Theses and Dissertations 2017 Marginal Structural Cox Model for Survival Data with Treatment-Confounder Feedback Yanan Zhang University of South Carolina Follow
More informationDesigning Information Devices and Systems I Fall 2016 Babak Ayazifar, Vladimir Stojanovic Homework 12
EECS 6A Designing Information Devices and Systems I Fall 206 Babak Ayazifar, Vladimir Stojanovic Homework 2 This homework is due November 22, 206, at P.M. Recommended Reading: Gilbert Strang, Introduction
More informationDesigning Information Devices and Systems I Spring 2016 Elad Alon, Babak Ayazifar Homework 11
EECS 6A Designing Information Devices and Systems I Spring 206 Elad Alon, Babak Ayazifar Homework This homework is due April 9, 206, at Noon.. Homework process and study group Who else did you work with
More informationNonparametric Regression Analysis of Multivariate Longitudinal Data
Nonparametric Regression Analysis of Multivariate Longitudinal Data Dongdong Xiang 1, Peihua Qiu 2 and Xiaolong Pu 1 1 School of Finance and Statistics, East China Normal University 2 School of Statistics,
More informationEstimation of Conditional Kendall s Tau for Bivariate Interval Censored Data
Communications for Statistical Applications and Methods 2015, Vol. 22, No. 6, 599 604 DOI: http://dx.doi.org/10.5351/csam.2015.22.6.599 Print ISSN 2287-7843 / Online ISSN 2383-4757 Estimation of Conditional
More informationLecture 32: Infinite-dimensional/Functionvalued. Functions and Random Regressions. Bruce Walsh lecture notes Synbreed course version 11 July 2013
Lecture 32: Infinite-dimensional/Functionvalued Traits: Covariance Functions and Random Regressions Bruce Walsh lecture notes Synbreed course version 11 July 2013 1 Longitudinal traits Many classic quantitative
More informationPoint Formulation and Adaptive Smoothing
Statistica Sinica (2008): Preprint 1 Nonparametric Control Chart for Monitoring Profiles Using Change Point Formulation and Adaptive Smoothing Changliang Zou 1, Peihua Qiu 2 and Douglas Hawkins 2 1 Nankai
More informationNonparametric Monitoring of Multiple Count Data
Nonparametric Monitoring of Multiple Count Data Peihua Qiu 1, Zhen He 2 and Zhiqiong Wang 3 1 Department of Biostatistics University of Florida, Gainesville, United States 2 College of Management and Economics
More informationMS&E 226: Small Data
MS&E 226: Small Data Lecture 12: Logistic regression (v1) Ramesh Johari ramesh.johari@stanford.edu Fall 2015 1 / 30 Regression methods for binary outcomes 2 / 30 Binary outcomes For the duration of this
More informationA Generalized Global Rank Test for Multiple, Possibly Censored, Outcomes
A Generalized Global Rank Test for Multiple, Possibly Censored, Outcomes Ritesh Ramchandani Harvard School of Public Health August 5, 2014 Ritesh Ramchandani (HSPH) Global Rank Test for Multiple Outcomes
More informationSingle Equation Linear GMM with Serially Correlated Moment Conditions
Single Equation Linear GMM with Serially Correlated Moment Conditions Eric Zivot October 28, 2009 Univariate Time Series Let {y t } be an ergodic-stationary time series with E[y t ]=μ and var(y t )
More informationTwo sample hypothesis testing
Statistics February 26, 2014 Debdeep Pati Two sample hypothesis testing 1. Suppose we want to study the relationship between use of oral contraceptives (OC) and level of blood pressure (BP) in women. 2.
More informationMachine Learning Linear Classification. Prof. Matteo Matteucci
Machine Learning Linear Classification Prof. Matteo Matteucci Recall from the first lecture 2 X R p Regression Y R Continuous Output X R p Y {Ω 0, Ω 1,, Ω K } Classification Discrete Output X R p Y (X)
More informationDirectionally Sensitive Multivariate Statistical Process Control Methods
Directionally Sensitive Multivariate Statistical Process Control Methods Ronald D. Fricker, Jr. Naval Postgraduate School October 5, 2005 Abstract In this paper we develop two directionally sensitive statistical
More informationResiduals in the Analysis of Longitudinal Data
Residuals in the Analysis of Longitudinal Data Jemila Hamid, PhD (Joint work with WeiLiang Huang) Clinical Epidemiology and Biostatistics & Pathology and Molecular Medicine McMaster University Outline
More informationMixed Effects Multivariate Adaptive Splines Model for the Analysis of Longitudinal and Growth Curve Data. Heping Zhang.
Mixed Effects Multivariate Adaptive Splines Model for the Analysis of Longitudinal and Growth Curve Data Heping Zhang Yale University NIA Workshop on Statistical Methods for Longitudinal Data on Aging
More informationWeb Appendix for Effect Estimation using Structural Nested Models and G-estimation
Web Appendix for Effect Estimation using Structural Nested Models and G-estimation Introductory concepts and notation Anonymized authors. First, we provide some additional details on the general data framework
More informationExponentially Weighted Moving Average Control Charts for Monitoring Increases in Poisson Rate
Exponentially Weighted Moving Average Control Charts for Monitoring Increases in Poisson Rate Lianjie SHU 1,, Wei JIANG 2, and Zhang WU 3 EndAName 1 Faculty of Business Administration University of Macau
More informationSingle Equation Linear GMM with Serially Correlated Moment Conditions
Single Equation Linear GMM with Serially Correlated Moment Conditions Eric Zivot November 2, 2011 Univariate Time Series Let {y t } be an ergodic-stationary time series with E[y t ]=μ and var(y t )
More informationAn exponentially weighted moving average scheme with variable sampling intervals for monitoring linear profiles
An exponentially weighted moving average scheme with variable sampling intervals for monitoring linear profiles Zhonghua Li, Zhaojun Wang LPMC and Department of Statistics, School of Mathematical Sciences,
More informationSocial connectedness is associated with fibrinogen level in a human social network
SUPPLEMENTARY FIGURES AND TABLES FOR: Social connectedness is associated with fibrinogen level in a human social network David A. Kim, M.D., Ph.D. Emelia J. Benjamin, M.D., Sc.M. James H. Fowler, Ph.D.
More informationLecture 2: Constant Treatment Strategies. Donglin Zeng, Department of Biostatistics, University of North Carolina
Lecture 2: Constant Treatment Strategies Introduction Motivation We will focus on evaluating constant treatment strategies in this lecture. We will discuss using randomized or observational study for these
More informationBayesian causal forests: dealing with regularization induced confounding and shrinking towards homogeneous effects
Bayesian causal forests: dealing with regularization induced confounding and shrinking towards homogeneous effects P. Richard Hahn, Jared Murray, and Carlos Carvalho July 29, 2018 Regularization induced
More informationarxiv: v1 [stat.me] 14 Jan 2019
arxiv:1901.04443v1 [stat.me] 14 Jan 2019 An Approach to Statistical Process Control that is New, Nonparametric, Simple, and Powerful W.J. Conover, Texas Tech University, Lubbock, Texas V. G. Tercero-Gómez,Tecnológico
More informationDistribution-free ROC Analysis Using Binary Regression Techniques
Distribution-free Analysis Using Binary Techniques Todd A. Alonzo and Margaret S. Pepe As interpreted by: Andrew J. Spieker University of Washington Dept. of Biostatistics Introductory Talk No, not that!
More informationBIOL 51A - Biostatistics 1 1. Lecture 1: Intro to Biostatistics. Smoking: hazardous? FEV (l) Smoke
BIOL 51A - Biostatistics 1 1 Lecture 1: Intro to Biostatistics Smoking: hazardous? FEV (l) 1 2 3 4 5 No Yes Smoke BIOL 51A - Biostatistics 1 2 Box Plot a.k.a box-and-whisker diagram or candlestick chart
More informationDirectional Control Schemes for Multivariate Categorical Processes
Directional Control Schemes for Multivariate Categorical Processes Nankai University Email: chlzou@yahoo.com.cn Homepage: math.nankai.edu.cn/ chlzou (Joint work with Mr. Jian Li and Prof. Fugee Tsung)
More informationProbabilistic classification CE-717: Machine Learning Sharif University of Technology. M. Soleymani Fall 2016
Probabilistic classification CE-717: Machine Learning Sharif University of Technology M. Soleymani Fall 2016 Topics Probabilistic approach Bayes decision theory Generative models Gaussian Bayes classifier
More informationarxiv: v1 [stat.ap] 6 Apr 2018
Individualized Dynamic Prediction of Survival under Time-Varying Treatment Strategies Grigorios Papageorgiou 1, 2, Mostafa M. Mokhles 2, Johanna J. M. Takkenberg 2, arxiv:1804.02334v1 [stat.ap] 6 Apr 2018
More informationMonitoring Wafer Geometric Quality using Additive Gaussian Process
Monitoring Wafer Geometric Quality using Additive Gaussian Process Linmiao Zhang 1 Kaibo Wang 2 Nan Chen 1 1 Department of Industrial and Systems Engineering, National University of Singapore 2 Department
More informationUniversity of California, Berkeley
University of California, Berkeley U.C. Berkeley Division of Biostatistics Working Paper Series Year 2009 Paper 251 Nonparametric population average models: deriving the form of approximate population
More informationOptimal Treatment Regimes for Survival Endpoints from a Classification Perspective. Anastasios (Butch) Tsiatis and Xiaofei Bai
Optimal Treatment Regimes for Survival Endpoints from a Classification Perspective Anastasios (Butch) Tsiatis and Xiaofei Bai Department of Statistics North Carolina State University 1/35 Optimal Treatment
More informationMonitoring General Linear Profiles Using Multivariate EWMA schemes
Monitoring General Linear Profiles Using Multivariate EWMA schemes Changliang Zou Department of Statistics School of Mathematical Sciences Nankai University Tianjian, PR China Fugee Tsung Department of
More informationMS&E 226: Small Data
MS&E 226: Small Data Lecture 9: Logistic regression (v2) Ramesh Johari ramesh.johari@stanford.edu 1 / 28 Regression methods for binary outcomes 2 / 28 Binary outcomes For the duration of this lecture suppose
More informationThe concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.
The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes
More informationCOMPARING GROUPS PART 1CONTINUOUS DATA
COMPARING GROUPS PART 1CONTINUOUS DATA Min Chen, Ph.D. Assistant Professor Quantitative Biomedical Research Center Department of Clinical Sciences Bioinformatics Shared Resource Simmons Comprehensive Cancer
More information2.830J / 6.780J / ESD.63J Control of Manufacturing Processes (SMA 6303) Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 2.830J / 6.780J / ESD.63J Control of Processes (SMA 6303) Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/term
More informationA NOTE ON A NONPARAMETRIC REGRESSION TEST THROUGH PENALIZED SPLINES
Statistica Sinica 24 (2014), 1143-1160 doi:http://dx.doi.org/10.5705/ss.2012.230 A NOTE ON A NONPARAMETRIC REGRESSION TEST THROUGH PENALIZED SPLINES Huaihou Chen 1, Yuanjia Wang 2, Runze Li 3 and Katherine
More informationA DYNAMIC QUANTILE REGRESSION TRANSFORMATION MODEL FOR LONGITUDINAL DATA
Statistica Sinica 19 (2009), 1137-1153 A DYNAMIC QUANTILE REGRESSION TRANSFORMATION MODEL FOR LONGITUDINAL DATA Yunming Mu and Ying Wei Portland State University and Columbia University Abstract: This
More informationSystem Monitoring with Real-Time Contrasts
System Monitoring with Real- Contrasts HOUTAO DENG Intuit, Mountain View, CA 94043, USA GEORGE RUNGER Arizona State University, Tempe, AZ 85287, USA EUGENE TUV Intel Corporation, Chandler, AZ 85226, USA
More informationConstrained Maximum Likelihood Estimation for Model Calibration Using Summary-level Information from External Big Data Sources
Constrained Maximum Likelihood Estimation for Model Calibration Using Summary-level Information from External Big Data Sources Yi-Hau Chen Institute of Statistical Science, Academia Sinica Joint with Nilanjan
More informationEstimation in Covariate Adjusted Regression
Estimation in Covariate Adjusted Regression Damla Şentürk 1 and Danh V. Nguyen 2 1 Department of Statistics, Pennsylvania State University University Park, PA 16802, U.S.A. email: dsenturk@stat.psu.edu
More informationEfficient Control Chart Calibration by Simulated Stochastic Approximation
Efficient Control Chart Calibration by Simulated Stochastic Approximation 1/23 Efficient Control Chart Calibration by Simulated Stochastic Approximation GIOVANNA CAPIZZI and GUIDO MASAROTTO Department
More informationSample Size and Power Considerations for Longitudinal Studies
Sample Size and Power Considerations for Longitudinal Studies Outline Quantities required to determine the sample size in longitudinal studies Review of type I error, type II error, and power For continuous
More informationTA: Sheng Zhgang (Th 1:20) / 342 (W 1:20) / 343 (W 2:25) / 344 (W 12:05) Haoyang Fan (W 1:20) / 346 (Th 12:05) FINAL EXAM
STAT 301, Fall 2011 Name Lec 4: Ismor Fischer Discussion Section: Please circle one! TA: Sheng Zhgang... 341 (Th 1:20) / 342 (W 1:20) / 343 (W 2:25) / 344 (W 12:05) Haoyang Fan... 345 (W 1:20) / 346 (Th
More informationImplementation of MEWMA Control Chart in Equipment Condition Monitoring
Implementation of MEWMA Control Chart in Equipment Condition Monitoring S. S. Lampreia 1,, J. G. Requeijo 2, J. M. Dias 2 and V. Vairinhos 1,3 1 Naval Academy Mechanical Engineer Department, CINAV Centro
More informationA Nonparametric Control Chart Based On The Mann Whitney
A Nonparametric Control Chart Based On The Mann Whitney We have made it easy for you to find a PDF Ebooks without any digging. And by having access to our ebooks online or by storing it on your computer,
More informationDifferential Equations Practice: 2nd Order Linear: Nonhomogeneous Equations: Undetermined Coefficients Page 1
Differential Equations Practice: 2nd Order Linear: Nonhomogeneous Equations: Undetermined Coefficients Page 1 Questions Example (3.5.3) Find a general solution of the differential equation y 2y 3y = 3te
More informationLikelihood-Based EWMA Charts for Monitoring Poisson Count Data with Time-Varying Sample Sizes
Likelihood-Based EWMA Charts for Monitoring Poisson Count Data with Time-Varying Sample Sizes Qin Zhou 1,3, Changliang Zou 1, Zhaojun Wang 1 and Wei Jiang 2 1 LPMC and Department of Statistics, School
More informationCovariate Adjusted Varying Coefficient Models
Biostatistics Advance Access published October 26, 2005 Covariate Adjusted Varying Coefficient Models Damla Şentürk Department of Statistics, Pennsylvania State University University Park, PA 16802, U.S.A.
More informationAn Adaptive Exponentially Weighted Moving Average Control Chart for Monitoring Process Variances
An Adaptive Exponentially Weighted Moving Average Control Chart for Monitoring Process Variances Lianjie Shu Faculty of Business Administration University of Macau Taipa, Macau (ljshu@umac.mo) Abstract
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 2 Laplace Transform I 1/52
1/52 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 2 Laplace Transform I Linear Time Invariant Systems A general LTI system may be described by the linear constant coefficient differential equation: a n d n
More information171:162 Design and Analysis of Biomedical Studies, Summer 2011 Exam #3, July 16th
Name 171:162 Design and Analysis of Biomedical Studies, Summer 2011 Exam #3, July 16th Use the selected SAS output to help you answer the questions. The SAS output is all at the back of the exam on pages
More informationLecture 4: Generalized Linear Mixed Models
Dankmar Böhning Southampton Statistical Sciences Research Institute University of Southampton, UK S 3 RI, 11-12 December 2014 An example with one random effect An example with two nested random effects
More informationHarvard University. Harvard University Biostatistics Working Paper Series. Sheng-Hsuan Lin Jessica G. Young Roger Logan
Harvard University Harvard University Biostatistics Working Paper Series Year 01 Paper 0 Mediation analysis for a survival outcome with time-varying exposures, mediators, and confounders Sheng-Hsuan Lin
More informationMarginal Screening and Post-Selection Inference
Marginal Screening and Post-Selection Inference Ian McKeague August 13, 2017 Ian McKeague (Columbia University) Marginal Screening August 13, 2017 1 / 29 Outline 1 Background on Marginal Screening 2 2
More informationStudy Ch. 10.3, 67 70all, (no CI), 81, 83
GOALS: 1. Compare 2 sample means when the population standard deviations are not equal. 2. Use the distribution of the difference between the means to evaluate the samples. 3. Arrive at a conclusion: are
More informationDe-biasing the Lasso: Optimal Sample Size for Gaussian Designs
De-biasing the Lasso: Optimal Sample Size for Gaussian Designs Adel Javanmard USC Marshall School of Business Data Science and Operations department Based on joint work with Andrea Montanari Oct 2015 Adel
More informationIntroduction to Statistical Analysis
Introduction to Statistical Analysis Changyu Shen Richard A. and Susan F. Smith Center for Outcomes Research in Cardiology Beth Israel Deaconess Medical Center Harvard Medical School Objectives Descriptive
More informationCross-Validation with Confidence
Cross-Validation with Confidence Jing Lei Department of Statistics, Carnegie Mellon University WHOA-PSI Workshop, St Louis, 2017 Quotes from Day 1 and Day 2 Good model or pure model? Occam s razor We really
More informationStatistical Inference with Monotone Incomplete Multivariate Normal Data
Statistical Inference with Monotone Incomplete Multivariate Normal Data p. 1/4 Statistical Inference with Monotone Incomplete Multivariate Normal Data This talk is based on joint work with my wonderful
More informationLocal Box Cox transformation on time-varying parametric models for smoothing estimation of conditional CDF with longitudinal data
Journal of Statistical Computation and Simulation ISSN: 0094-9655 (Print) 1563-5163 (Online) Journal homepage: http://www.tandfonline.com/loi/gscs20 Local Box Cox transformation on time-varying parametric
More informationLecture 7 Time-dependent Covariates in Cox Regression
Lecture 7 Time-dependent Covariates in Cox Regression So far, we ve been considering the following Cox PH model: λ(t Z) = λ 0 (t) exp(β Z) = λ 0 (t) exp( β j Z j ) where β j is the parameter for the the
More informationMulti-model Markov Decision Processes
Multi-model Markov Decision Processes Lauren N. Steimle Department of Industrial and Operations Engineering, University of Michigan, Ann Arbor, MI 48109, steimle@umich.edu David L. Kaufman Management Studies,
More informationeappendix: Description of mgformula SAS macro for parametric mediational g-formula
eappendix: Description of mgformula SAS macro for parametric mediational g-formula The implementation of causal mediation analysis with time-varying exposures, mediators, and confounders Introduction The
More informationAnalysis of Longitudinal Data. Patrick J. Heagerty PhD Department of Biostatistics University of Washington
Analsis of Longitudinal Data Patrick J. Heagert PhD Department of Biostatistics Universit of Washington 1 Auckland 2008 Session Three Outline Role of correlation Impact proper standard errors Used to weight
More informationCross-Validation with Confidence
Cross-Validation with Confidence Jing Lei Department of Statistics, Carnegie Mellon University UMN Statistics Seminar, Mar 30, 2017 Overview Parameter est. Model selection Point est. MLE, M-est.,... Cross-validation
More informationStatistical Inference with Monotone Incomplete Multivariate Normal Data
Statistical Inference with Monotone Incomplete Multivariate Normal Data p. 1/4 Statistical Inference with Monotone Incomplete Multivariate Normal Data This talk is based on joint work with my wonderful
More informationSubset selection with sparse matrices
Subset selection with sparse matrices Alberto Del Pia, University of Wisconsin-Madison Santanu S. Dey, Georgia Tech Robert Weismantel, ETH Zürich February 1, 018 Schloss Dagstuhl Subset selection for regression
More informationWell-developed and understood properties
1 INTRODUCTION TO LINEAR MODELS 1 THE CLASSICAL LINEAR MODEL Most commonly used statistical models Flexible models Well-developed and understood properties Ease of interpretation Building block for more
More informationNonparametric meta-analysis for diagnostic accuracy studies Antonia Zapf
Nonparametric meta-analysis for diagnostic accuracy studies Antonia Zapf joint work with A. Hoyer, K. Kramer, and O. Kuss Table of contents Motivation Nonparametric approach Simulation study Application
More informationwith the usual assumptions about the error term. The two values of X 1 X 2 0 1
Sample questions 1. A researcher is investigating the effects of two factors, X 1 and X 2, each at 2 levels, on a response variable Y. A balanced two-factor factorial design is used with 1 replicate. The
More informationMAS3301 / MAS8311 Biostatistics Part II: Survival
MAS3301 / MAS8311 Biostatistics Part II: Survival M. Farrow School of Mathematics and Statistics Newcastle University Semester 2, 2009-10 1 13 The Cox proportional hazards model 13.1 Introduction In the
More informationMultivariate Process Control Chart for Controlling the False Discovery Rate
Industrial Engineering & Management Systems Vol, No 4, December 0, pp.385-389 ISSN 598-748 EISSN 34-6473 http://dx.doi.org/0.73/iems.0..4.385 0 KIIE Multivariate Process Control Chart for Controlling e
More informationMonitoring Paired Binary Surgical Outcomes Using Cumulative Sum Charts
Monitoring Paired Binary Surgical Outcomes Using Cumulative Sum Charts Stefan H. Steiner and Richard J. Cook Department of Statistics and Actuarial Sciences University of Waterloo Waterloo, Ontario Canada
More information