Outline. Recall that Aalen additive hazards model and the semiparametric version
|
|
- Eleanore Benson
- 6 years ago
- Views:
Transcription
1 Outlne Clustered survval data (addtve models) Addtve-Multplcatve hazards model Advanced Survval Analyss 21 Copenhagen Gulana Cortese Clustered survval data Margnal addtve models Addtve-Multplcatve hazards model Cox-Aalen model Proportonal excess hazards models Aalen addtve model Clustered survval data Recall that Aalen addtve hazards model and the semparametrc verson p λ(t) = Y (t)x T (t)β(t) = Y (t) X j(t)β j(t), (1) and the semparametrc verson j=1 λ (t) = Y [ (t) X T (t)β(t) + Z T (t)γ ]. (2) These are models for the ntensty of the observed countng process, e.. λ(t)dt = E [dn(t) σ(n(s), X (s), Z(s), Y (s), s [, t) ] In clustered data, falure tmes wthn each cluster may be correlated. Many examples: to blndness n both eyes Twn studes Recurrent events Here, choce of survval model depends on the ams of the study, such as 1) correlaton of falure tmes wthn clusters (condtonal models) 2) comparson of falure tmes across clusters (margnal models)
2 Clustered survval data Clustered survval data Two basc approaches : 1) Fralty models: gven random eects "Z" the survval tmes are ndependent wth hazards λ(t)z exp(x T β). 2) Margnal models: Gven covarates the margnal rate (hazard) s Y (t)λ(t) exp(x T β) = E(dN(t) X, Y (t) = 1)/dt For k = 1,..., K, = 1,..., n, let T and C be the falure and censorng tmes for the th ndvdual n the kth cluster, and let X (t) be a p-vector of covarates. Put T k = ( T1k,..., Tnk), C k = (C1k,..., C nk), X k(t) = (X1k(t),..., X nk(t)). We assume that ( Tk, C k, X k( )), k = 1,..., K are ndependent and dentcally dstrbuted varables and these varables follow the model descrbed n the followng. The rght-censored falure tme s denoted T = T C and as usual we let Y (t) = 1(T t). Margnal models Rght-censored falure tmes T = T C, Y (t) = 1(T t), N (t) = 1(T t, T = T) Margnal (ntensty) model s a model Cox model: F t = σ{n (s), Y (s), X (s) : s t}, (3) λ F (t) = Y (t)λ(t) exp (X T (t)β). (4) It s mportant to note that (4) s not the ntensty wth respect to the observed ltraton F t = F k t, (5) k F k t = σ{n (s), Y (s), X (s) : = 1, n, s t} s the nformaton generated by observng all the ndvduals n the kth cluster. Margnal addtve models Instead of the Cox model we assume the Aalen addtve model for the margnal ntenstes, λ F (t) = Y (t)x T (t)β(t). (6) Possble nteracton terms are ncluded n the covarate vector X T (t). Let Xk(t) be the n p-matrx wth th row and let (Y (t)x 1(t),..., Y (t)x p(t)). N k(t) = (N1k(t),..., N nk(t)) T M k(t) = N k(t) X k(s)db(s). The -th component of M k(t), M (t) s a martngale wth respect to, but M k(t) s not a martngale wth respect to the observed ltraton. F t
3 Margnal addtve models Margnal addtve models One may also show that (7) has the same lmt dstrbuton as The (unweghted) workng ndependence estmator of B(t) s We have that K K ˆB(t) = [ X T k (s) Xk(s)] 1 X T k (s)dn k(s) K K 1/2 (ˆB(t) B(t)) = K 1/2 K [K 1 X T k (s) Xk(s)] 1 X T k (s)dm k(s) K = K 1/2 ɛ k(t) (7) whch s essentally a sum of..d. components (replace K 1 K X T (t) Xk(t) by ts lmt n probablty). k ˆɛ k(t) = K 1/2 K K 1/2 K ˆɛ k(t)g k, (8) K [K 1 X T k (s) Xk(s)] 1 T X k (s)d ˆMk(s) and ˆM k(t) = N k(t) X k(s)d ˆB(s). G1,..., GK are ndependent standard normal. The asymptotc covarance matrx of K 1/2 ( ˆB(t) B(t)) s estmated consstently by K K 1 ˆɛ 2 (t). k Resamplng from (8) may be used to construct condence bands and obtan p-values n tests about B(t). Example: Dabetc retnopathy data Example: Dabetc retnopathy data Subset of 197 patents wth dbetc retnopathy. Purpose of study was to assess the ecacy of laser treatment n delayng onset of blndness. Treatment assgned randomly to one eye of each patent, the other eye was wthout treatment. d tme status trteye treat adult > adult.treat<-(dabetes$adult==2)*(dabetes$treat) > ft<-aalen(surv(tme,status) ~adult+treat+adult.treat, dabetes,cluster=dabetes$d) > plot(ft) > summary(ft) Addtve Aalen Model Test for nonparametrc terms Test for non-sgnfcant effects sup hat B(t)/SD(t) p-value H_: B(t)= (Intercept) adult treat adult.treat Test for tme nvarant effects sup B(t) - (t/tau)b(tau) p-value H_: B(t)=b t (Intercept) adult treat adult.treat nt (B(t)-(t/tau)B(tau))^2dt p-value H_: B(t)=b t (Intercept) adult treat adult.treat
4 Example: Dabetc retnopathy data Example: Dabetc retnopathy data (Intercept) adult The margnal addtve model may be smpled by assumng a tme-constant eect for the nteracton term. Cumulatve coeffcents Cumulatve coeffcents > ft.s<-aalen(surv(tme,status) ~adult+treat+const(adult.treat), dabetes,cluster=dabetes$d) > par(mfrow=c(2,2));plot(ft,sm.c=2) Test for non-sgnfcant effects Supremum-test of sgnfcance p-value H_: B(t)= (Intercept) adult treat adult.treat treat Cumulatve coeffcents Cumulatve coeffcents Test for tme nvarant effects Kolmogorov-Smrnov test p-value H_:constant effect (Intercept) adult treat Parametrc terms : Coef. SE Robust SE z P-val const(adult.treat) Fgure: Dabetes-data. Cumulatve regresson estmators along wth 95% condence ntervals (full lnes) and unform bands (broken lnes). Margnal semparametrc addtve models What are the conclusons about the nteracton? What does ts regresson coecent express? Multplcatve and addtve hazards models Multplcatve ntensty models One can also consder a semparametrc addtve model the margnal ntenstes are assumed to be λ F (t) = Y [ (t) X T (t)β(t) + Z T (t)γ ], (9) γ s an unknown q-dmensonal vector of tme-constant coecents. Workng ndependence estmator of B(t) and γ can be found smlarly to the nonparametrc model. Agan, an d decomposton for K 1/2 (B(t, ˆγ) B(t)) and for K 1/2 (ˆγ γ) can be obtaned and used for resamplng. Addtve models If one compares the two models: α (t) = Y (t) exp (β(t) T X (t)γ T Z (t)) α (t) = Y (t)(β(t) T X (t) + γ T Z (t)) Multplcatve model very attractve but smoothng s needed! Multplcatve model leads to a relatve-rsk type summary. Addtve model very appealng because of lack of smoothng parameter, but hazard predctons may be negatve. Addtve model leads to excess rsk nterpretaton. Some eects are addtve (competng rsks model), and some wll be multplcatve (smokng). Goodness-of-t wll seldom clearly decde whch model s the correct one.
5 Survval Estmaton for proportonal/addtve models Survval estmate s S (t) = exp( Therefore estmate λ (s)ds: Cox: exp(βt X )dλ(s). Multplcatve: exp(β(s)t X )dλ(s). Aalen: B(s) T X. Cox-Aalen: (B(s) T X ) exp(γ T Z). λ (s)ds). In multplcatve model S (t) needs the estmates of beta(s), and thus t depends on bandwdth. Asymptotcs are more complcated. In addtve model S (t) does not need smoothng, problems of negatve hazards predctons can be handled wth post-processng. Multplcatve-addtve hazards models The addtve and multplcatve hazards models may be combned as follows: Cox-Aalen model λ(t) = Y (t) [ X T (t)α(t) ] exp{z T (t)β}. Here there s an addtve model for the baselne based on covarate X. The nterpretaton of relatve rsks for Z s the same as n Cox model, whle the coecents α(t) represent the excess baselne due to presence of X. Proportonal excess hazards model λ(t) = Y (t)x T (t)α(t)+ρ(t)λ(t) exp{z T (t)β}, Y (t) and ρ(t) are at rsk ndcators, λ(t) s the baselne hazard of the excess rsk term. Estmaton n Cox-Aalen Estmaton n Cox-Aalen λ(t) = Y (t) [ X T (t)α(t) ] exp{z T (t)β}. The log-lelhood functon leads to the score equatons for β and da(t) Z(t) T { dn Y (β, t) T da(t) } =, Y (β, t) T W (t) { dn Y (β, t) T da(t) } =, W (t) = dag(w (t)) wth w (t) = Y (t) λ = Y (t) exp( Z (t) T β), (t) X (t) T α(t) for = 1,..., n. For known β ths leads to the estmator of the cumulatve ntensty Â(β, t) = Y (β, s)dn(s), Y (β, t) = (Y (β, t) T W (t)y (β, t)) 1 Y (β, t) T W (t) Insertng ths estmator nto the score equaton for β gves U(β) = wth τ U(β) = U(β, τ) = Z T (t)g(β, t)dn(t), (1) G(β, t) = I Y (β, t)y (β, t) s the projecton onto the orthogonal space spanned by the columns of Y (β, t). Easy to see that U(β, t) = s a square ntegrable martngale. Z T (s)g(β, s)dm(s)
6 Example: PBC data Example: PBC data > ft <- cox.aalen(surv(tme/365, status) ~ prop(age) + edema + + prop(log(bl)) + prop(log(alb)) + log(protme), pbc, + maxtme = 3/365) > summary(ft) (Intercept) Edema pro.m Test for tme nvarant effects sup B(t) - (t/tau)b(tau) p-value H_: B(t)=b t (Intercept) edema log(protme) Proportonal Cox terms : Coef. Std. Error Robust SE D2log(L)^ prop(age) prop(log(bl)) prop(log(alb)) Score Tests for Proportonalty sup hat U(t) p-value H_ Cumulatve regresson functon Cumulatve regresson functon Cumulatve regresson functon prop(age) prop(log(bl)) prop(log(alb)) Proportonal excess hazard models Example: Melanoma data For model 11 It can be seen as an excess rsk model (Martnussen & schee, 22): λ(t) = Y (t)x T (t)α(t)+ρ(t)λ(t) exp{z T (t)β}, (11) data(melanoma) lt<-log(melanoma$thck) excess<-(melanoma$thck>=21) # excess rsk for thck tumors ft<-prop.excess(surv(days/365,status==1)~ sex + ulc + cox(sex)+ cox(ulc) + cox(lt), melanoma, excess=excess, n.sm=1) For model 12 ρ = 1, all, gves a mx of Aalen's and Cox's models. The Cox term represents the excess rsk for an exposed subject wth covarates Z. ρ s an excess ndcator, e.g. I (d > ) wth d doses for th subject. The addtve term X T (t)α(t) can be estmated from the study. It can also be replaced by a known functon of X, α(t, X ), whch represents the background mortalty rate of a control populaton (often n cancer studes) (Sasen, 1996): λ(t) = Y (t) [ α(t, X )+λ(t) exp{z T (t)β} ]. (12) data(mela.pop) out<-pe.sasen(surv(start,stop,status==1)~ age + sex, mela.pop, d=1:25, max.tme=7, offsets=mela.pop$rate, n.sm=1) summary(out) Proportonal terms: coef se(coef) z p age sex Test for Proportonalty sup hat U(t) p-value H_ age sex ul<-out$cum[,2]+1.96*out$var.cum[,2]^.5 ll<-out$cum[,2]-1.96*out$var.cum[,2]^.5 plot(out$cum,type="s",ylm=range(ul,ll)) lnes(out$cum[,1],ul,type="s"); lnes(out$cum[,1],ll,type="s")
7 Exercse on Dabetes data Consder the Dabetes data n the tmereg package. Start from the followng model ft.s<-aalen(surv(tme,status) ~adult+treat+const(adult.treat), dabetes,cluster=dabetes$d) Try to smplfy further the above margnal semparametrc addtve model, after successve GOF tests. What s a nal possble model? Solve Exercse 9.1 of the book Martnussen & Schee (26). Exercse on TRACE data The TRACE study group conducted a study amed at assessng the prognostc mportance of varous rsk factors on mortalty for approxmately 6 patents wth myocardal nfarcton. The patents had varous rsk factors recorded such as age, sex (male=1), congestve heart falure (chf) (present=1), dabetes (present=1) and ventrcular brllaton (vf) (present=1). Consder the TRACE data n the tmereg package (a subset of patents from the orgnal data) wth the above covarates. ft.tr <- cox.aalen(surv(tme,status==9)~ chf + age + sex + prop(dabetes) + vf, TRACE) plot(ft.tr, sm.c=2,xlab=" (years)") Conclude what s/are the best model/models (Cox, addtve, Cox-Aalen) from ttng the TRACE data and arrve to a nal concluson by usng GOF procedures. Compare the nal model wth results from a Cox model wth "vf" as a stratcaton factor. Can ths be obtaned from cox.aalen? Plot survval predctons, gven some covarate values, for your nal model and the above strated Cox model. By comparng estmated curves, what may we conclude?
STK4080/9080 Survival and event history analysis
SK48/98 Survval and event hstory analyss Lecture 7: Regresson modellng Relatve rsk regresson Regresson models Assume that we have a sample of n ndvduals, and let N (t) count the observed occurrences of
More informationGoodness-Of-Fit for Cox s Regression Model. Extensions of Cox s Regression Model. Survival Analysis Fall 2004, Copenhagen
Outline Cox s proportional hazards model. Goodness-of-fit tools More flexible models R-package timereg Forthcoming book, Martinussen and Scheike. 2/38 University of Copenhagen http://www.biostat.ku.dk
More informationChapter 13: Multiple Regression
Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to
More informationChapter 20 Duration Analysis
Chapter 20 Duraton Analyss Duraton: tme elapsed untl a certan event occurs (weeks unemployed, months spent on welfare). Survval analyss: duraton of nterest s survval tme of a subject, begn n an ntal state
More informationRegression Analysis of Clustered Failure Time Data under the Additive Hazards Model
A^VÇÚO 1 33 ò 1 5 Ï 217 c 1 Chnese Journal of Appled Probablty and Statstcs Oct., 217, Vol. 33, No. 5, pp. 517-528 do: 1.3969/j.ssn.11-4268.217.5.8 Regresson Analyss of Clustered Falure Tme Data under
More informationStat 642, Lecture notes for 01/27/ d i = 1 t. n i t nj. n j
Stat 642, Lecture notes for 01/27/05 18 Rate Standardzaton Contnued: Note that f T n t where T s the cumulatve follow-up tme and n s the number of subjects at rsk at the mdpont or nterval, and d s the
More informationTime to dementia onset: competing risk analysis with Laplace regression
Tme to dementa onset: competng rsk analyss wth Laplace regresson Gola Santon, Debora Rzzuto, Laura Fratglon 4 th Nordc and Baltc STATA Users group meetng, Stockholm, November 20 Agng Research Center (ARC),
More informationNon-Mixture Cure Model for Interval Censored Data: Simulation Study ABSTRACT
Malaysan Journal of Mathematcal Scences 8(S): 37-44 (2014) Specal Issue: Internatonal Conference on Mathematcal Scences and Statstcs 2013 (ICMSS2013) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal
More informationSTAT 511 FINAL EXAM NAME Spring 2001
STAT 5 FINAL EXAM NAME Sprng Instructons: Ths s a closed book exam. No notes or books are allowed. ou may use a calculator but you are not allowed to store notes or formulas n the calculator. Please wrte
More information1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands
Content. Inference on Regresson Parameters a. Fndng Mean, s.d and covarance amongst estmates.. Confdence Intervals and Workng Hotellng Bands 3. Cochran s Theorem 4. General Lnear Testng 5. Measures of
More informationThe Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction
ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also
More informationBIO Lab 2: TWO-LEVEL NORMAL MODELS with school children popularity data
Lab : TWO-LEVEL NORMAL MODELS wth school chldren popularty data Purpose: Introduce basc two-level models for normally dstrbuted responses usng STATA. In partcular, we dscuss Random ntercept models wthout
More informationChapter 9: Statistical Inference and the Relationship between Two Variables
Chapter 9: Statstcal Inference and the Relatonshp between Two Varables Key Words The Regresson Model The Sample Regresson Equaton The Pearson Correlaton Coeffcent Learnng Outcomes After studyng ths chapter,
More informationDummy variables in multiple variable regression model
WESS Econometrcs (Handout ) Dummy varables n multple varable regresson model. Addtve dummy varables In the prevous handout we consdered the followng regresson model: y x 2x2 k xk,, 2,, n and we nterpreted
More informationx i1 =1 for all i (the constant ).
Chapter 5 The Multple Regresson Model Consder an economc model where the dependent varable s a functon of K explanatory varables. The economc model has the form: y = f ( x,x,..., ) xk Approxmate ths by
More informationStatistics for Economics & Business
Statstcs for Economcs & Busness Smple Lnear Regresson Learnng Objectves In ths chapter, you learn: How to use regresson analyss to predct the value of a dependent varable based on an ndependent varable
More informationEcon107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)
I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes
More informationUNIVERSITY OF TORONTO Faculty of Arts and Science. December 2005 Examinations STA437H1F/STA1005HF. Duration - 3 hours
UNIVERSITY OF TORONTO Faculty of Arts and Scence December 005 Examnatons STA47HF/STA005HF Duraton - hours AIDS ALLOWED: (to be suppled by the student) Non-programmable calculator One handwrtten 8.5'' x
More informationChapter 8 Indicator Variables
Chapter 8 Indcator Varables In general, e explanatory varables n any regresson analyss are assumed to be quanttatve n nature. For example, e varables lke temperature, dstance, age etc. are quanttatve n
More informationRELIABILITY ASSESSMENT
CHAPTER Rsk Analyss n Engneerng and Economcs RELIABILITY ASSESSMENT A. J. Clark School of Engneerng Department of Cvl and Envronmental Engneerng 4a CHAPMAN HALL/CRC Rsk Analyss for Engneerng Department
More informationPredictive Analytics : QM901.1x Prof U Dinesh Kumar, IIMB. All Rights Reserved, Indian Institute of Management Bangalore
Sesson Outlne Introducton to classfcaton problems and dscrete choce models. Introducton to Logstcs Regresson. Logstc functon and Logt functon. Maxmum Lkelhood Estmator (MLE) for estmaton of LR parameters.
More informationParametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010
Parametrc fractonal mputaton for mssng data analyss Jae Kwang Km Survey Workng Group Semnar March 29, 2010 1 Outlne Introducton Proposed method Fractonal mputaton Approxmaton Varance estmaton Multple mputaton
More informationComparison of Regression Lines
STATGRAPHICS Rev. 9/13/2013 Comparson of Regresson Lnes Summary... 1 Data Input... 3 Analyss Summary... 4 Plot of Ftted Model... 6 Condtonal Sums of Squares... 6 Analyss Optons... 7 Forecasts... 8 Confdence
More informationComputation of Higher Order Moments from Two Multinomial Overdispersion Likelihood Models
Computaton of Hgher Order Moments from Two Multnomal Overdsperson Lkelhood Models BY J. T. NEWCOMER, N. K. NEERCHAL Department of Mathematcs and Statstcs, Unversty of Maryland, Baltmore County, Baltmore,
More informationDepartment of Quantitative Methods & Information Systems. Time Series and Their Components QMIS 320. Chapter 6
Department of Quanttatve Methods & Informaton Systems Tme Seres and Ther Components QMIS 30 Chapter 6 Fall 00 Dr. Mohammad Zanal These sldes were modfed from ther orgnal source for educatonal purpose only.
More informationStatistics for Business and Economics
Statstcs for Busness and Economcs Chapter 11 Smple Regresson Copyrght 010 Pearson Educaton, Inc. Publshng as Prentce Hall Ch. 11-1 11.1 Overvew of Lnear Models n An equaton can be ft to show the best lnear
More informationComparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method
Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method
More informationMaximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models
ECO 452 -- OE 4: Probt and Logt Models ECO 452 -- OE 4 Maxmum Lkelhood Estmaton of Bnary Dependent Varables Models: Probt and Logt hs note demonstrates how to formulate bnary dependent varables models
More informationThe Geometry of Logit and Probit
The Geometry of Logt and Probt Ths short note s meant as a supplement to Chapters and 3 of Spatal Models of Parlamentary Votng and the notaton and reference to fgures n the text below s to those two chapters.
More informationCS-433: Simulation and Modeling Modeling and Probability Review
CS-433: Smulaton and Modelng Modelng and Probablty Revew Exercse 1. (Probablty of Smple Events) Exercse 1.1 The owner of a camera shop receves a shpment of fve cameras from a camera manufacturer. Unknown
More informationStatistics for Managers Using Microsoft Excel/SPSS Chapter 13 The Simple Linear Regression Model and Correlation
Statstcs for Managers Usng Mcrosoft Excel/SPSS Chapter 13 The Smple Lnear Regresson Model and Correlaton 1999 Prentce-Hall, Inc. Chap. 13-1 Chapter Topcs Types of Regresson Models Determnng the Smple Lnear
More informationA joint frailty-copula model between disease progression and death for meta-analysis
CSA-KSS-JSS Specal Invted Sessons 4 / / 6 A jont fralty-copula model between dsease progresson and death for meta-analyss 3/5/7 Takesh Emura Graduate Insttute of Statstcs Natonal Central Unversty TAIWAN
More informationLecture Notes on Linear Regression
Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume
More information/ n ) are compared. The logic is: if the two
STAT C141, Sprng 2005 Lecture 13 Two sample tests One sample tests: examples of goodness of ft tests, where we are testng whether our data supports predctons. Two sample tests: called as tests of ndependence
More information8/25/17. Data Modeling. Data Modeling. Data Modeling. Patrice Koehl Department of Biological Sciences National University of Singapore
8/5/17 Data Modelng Patrce Koehl Department of Bologcal Scences atonal Unversty of Sngapore http://www.cs.ucdavs.edu/~koehl/teachng/bl59 koehl@cs.ucdavs.edu Data Modelng Ø Data Modelng: least squares Ø
More informationHere is the rationale: If X and y have a strong positive relationship to one another, then ( x x) will tend to be positive when ( y y)
Secton 1.5 Correlaton In the prevous sectons, we looked at regresson and the value r was a measurement of how much of the varaton n y can be attrbuted to the lnear relatonshp between y and x. In ths secton,
More informationLecture 9: Linear regression: centering, hypothesis testing, multiple covariates, and confounding
Recall: man dea of lnear regresson Lecture 9: Lnear regresson: centerng, hypothess testng, multple covarates, and confoundng Sandy Eckel seckel@jhsph.edu 6 May 8 Lnear regresson can be used to study an
More informationLecture 9: Linear regression: centering, hypothesis testing, multiple covariates, and confounding
Lecture 9: Lnear regresson: centerng, hypothess testng, multple covarates, and confoundng Sandy Eckel seckel@jhsph.edu 6 May 008 Recall: man dea of lnear regresson Lnear regresson can be used to study
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationBasically, if you have a dummy dependent variable you will be estimating a probability.
ECON 497: Lecture Notes 13 Page 1 of 1 Metropoltan State Unversty ECON 497: Research and Forecastng Lecture Notes 13 Dummy Dependent Varable Technques Studenmund Chapter 13 Bascally, f you have a dummy
More informationStatistics for Managers Using Microsoft Excel/SPSS Chapter 14 Multiple Regression Models
Statstcs for Managers Usng Mcrosoft Excel/SPSS Chapter 14 Multple Regresson Models 1999 Prentce-Hall, Inc. Chap. 14-1 Chapter Topcs The Multple Regresson Model Contrbuton of Indvdual Independent Varables
More informationIntroduction to the R Statistical Computing Environment R Programming
Introducton to the R Statstcal Computng Envronment R Programmng John Fox McMaster Unversty ICPSR 2018 John Fox (McMaster Unversty) R Programmng ICPSR 2018 1 / 14 Programmng Bascs Topcs Functon defnton
More informationHomework Assignment 3 Due in class, Thursday October 15
Homework Assgnment 3 Due n class, Thursday October 15 SDS 383C Statstcal Modelng I 1 Rdge regresson and Lasso 1. Get the Prostrate cancer data from http://statweb.stanford.edu/~tbs/elemstatlearn/ datasets/prostate.data.
More informationDepartment of Statistics University of Toronto STA305H1S / 1004 HS Design and Analysis of Experiments Term Test - Winter Solution
Department of Statstcs Unversty of Toronto STA35HS / HS Desgn and Analyss of Experments Term Test - Wnter - Soluton February, Last Name: Frst Name: Student Number: Instructons: Tme: hours. Ads: a non-programmable
More information= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.
Chapter Matlab Exercses Chapter Matlab Exercses. Consder the lnear system of Example n Secton.. x x x y z y y z (a) Use the MATLAB command rref to solve the system. (b) Let A be the coeffcent matrx and
More informationDO NOT OPEN THE QUESTION PAPER UNTIL INSTRUCTED TO DO SO BY THE CHIEF INVIGILATOR. Introductory Econometrics 1 hour 30 minutes
25/6 Canddates Only January Examnatons 26 Student Number: Desk Number:...... DO NOT OPEN THE QUESTION PAPER UNTIL INSTRUCTED TO DO SO BY THE CHIEF INVIGILATOR Department Module Code Module Ttle Exam Duraton
More informationLinear Correlation. Many research issues are pursued with nonexperimental studies that seek to establish relationships among 2 or more variables
Lnear Correlaton Many research ssues are pursued wth nonexpermental studes that seek to establsh relatonshps among or more varables E.g., correlates of ntellgence; relaton between SAT and GPA; relaton
More informationSemiparametric Methods of Time Scale Selection
1 Semparametrc Methods of Tme Scale Selecton Therry Duchesne Unversty of Toronto, Toronto, Canada Abstract: In several relablty applcatons, there may not be a unque plausble scale n whch to analyze falure
More informationCS 2750 Machine Learning. Lecture 5. Density estimation. CS 2750 Machine Learning. Announcements
CS 750 Machne Learnng Lecture 5 Densty estmaton Mlos Hauskrecht mlos@cs.ptt.edu 539 Sennott Square CS 750 Machne Learnng Announcements Homework Due on Wednesday before the class Reports: hand n before
More informationHowever, since P is a symmetric idempotent matrix, of P are either 0 or 1 [Eigen-values
Fall 007 Soluton to Mdterm Examnaton STAT 7 Dr. Goel. [0 ponts] For the general lnear model = X + ε, wth uncorrelated errors havng mean zero and varance σ, suppose that the desgn matrx X s not necessarly
More informationStat 543 Exam 2 Spring 2016
Stat 543 Exam 2 Sprng 206 I have nether gven nor receved unauthorzed assstance on ths exam. Name Sgned Date Name Prnted Ths Exam conssts of questons. Do at least 0 of the parts of the man exam. I wll score
More informationCredit Card Pricing and Impact of Adverse Selection
Credt Card Prcng and Impact of Adverse Selecton Bo Huang and Lyn C. Thomas Unversty of Southampton Contents Background Aucton model of credt card solctaton - Errors n probablty of beng Good - Errors n
More informationLecture 6 More on Complete Randomized Block Design (RBD)
Lecture 6 More on Complete Randomzed Block Desgn (RBD) Multple test Multple test The multple comparsons or multple testng problem occurs when one consders a set of statstcal nferences smultaneously. For
More information1 Binary Response Models
Bnary and Ordered Multnomal Response Models Dscrete qualtatve response models deal wth dscrete dependent varables. bnary: yes/no, partcpaton/non-partcpaton lnear probablty model LPM, probt or logt models
More informationLecture 6: Introduction to Linear Regression
Lecture 6: Introducton to Lnear Regresson An Manchakul amancha@jhsph.edu 24 Aprl 27 Lnear regresson: man dea Lnear regresson can be used to study an outcome as a lnear functon of a predctor Example: 6
More informationDiagnostics in Poisson Regression. Models - Residual Analysis
Dagnostcs n Posson Regresson Models - Resdual Analyss 1 Outlne Dagnostcs n Posson Regresson Models - Resdual Analyss Example 3: Recall of Stressful Events contnued 2 Resdual Analyss Resduals represent
More informationAn Introduction to Censoring, Truncation and Sample Selection Problems
An Introducton to Censorng, Truncaton and Sample Selecton Problems Thomas Crossley SPIDA June 2003 1 A. Introducton A.1 Basc Ideas Most of the statstcal technques we study are for estmatng (populaton)
More informationSimulation and Random Number Generation
Smulaton and Random Number Generaton Summary Dscrete Tme vs Dscrete Event Smulaton Random number generaton Generatng a random sequence Generatng random varates from a Unform dstrbuton Testng the qualty
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationOnline Classification: Perceptron and Winnow
E0 370 Statstcal Learnng Theory Lecture 18 Nov 8, 011 Onlne Classfcaton: Perceptron and Wnnow Lecturer: Shvan Agarwal Scrbe: Shvan Agarwal 1 Introducton In ths lecture we wll start to study the onlne learnng
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationChapter 11: Simple Linear Regression and Correlation
Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests
More informationLogistic regression models 1/12
Logstc regresson models 1/12 2/12 Example 1: dogs look lke ther owners? Some people beleve that dogs look lke ther owners. Is ths true? To test the above hypothess, The New York Tmes conducted a quz onlne.
More informationBasic Business Statistics, 10/e
Chapter 13 13-1 Basc Busness Statstcs 11 th Edton Chapter 13 Smple Lnear Regresson Basc Busness Statstcs, 11e 009 Prentce-Hall, Inc. Chap 13-1 Learnng Objectves In ths chapter, you learn: How to use regresson
More informationCorrelation and Regression. Correlation 9.1. Correlation. Chapter 9
Chapter 9 Correlaton and Regresson 9. Correlaton Correlaton A correlaton s a relatonshp between two varables. The data can be represented b the ordered pars (, ) where s the ndependent (or eplanator) varable,
More informationPhase I Monitoring of Nonlinear Profiles
Phase I Montorng of Nonlnear Profles James D. Wllams Wllam H. Woodall Jeffrey B. Brch May, 003 J.D. Wllams, Bll Woodall, Jeff Brch, Vrgna Tech 003 Qualty & Productvty Research Conference, Yorktown Heghts,
More informationChapter 11: I = 2 samples independent samples paired samples Chapter 12: I 3 samples of equal size J one-way layout two-way layout
Serk Sagtov, Chalmers and GU, February 0, 018 Chapter 1. Analyss of varance Chapter 11: I = samples ndependent samples pared samples Chapter 1: I 3 samples of equal sze one-way layout two-way layout 1
More informationSTAT 3008 Applied Regression Analysis
STAT 3008 Appled Regresson Analyss Tutoral : Smple Lnear Regresson LAI Chun He Department of Statstcs, The Chnese Unversty of Hong Kong 1 Model Assumpton To quantfy the relatonshp between two factors,
More information6 Supplementary Materials
6 Supplementar Materals 61 Proof of Theorem 31 Proof Let m Xt z 1:T : l m Xt X,z 1:t Wethenhave mxt z1:t ˆm HX Xt z 1:T mxt z1:t m HX Xt z 1:T + mxt z 1:T HX We consder each of the two terms n equaton
More informationLecture 3: Probability Distributions
Lecture 3: Probablty Dstrbutons Random Varables Let us begn by defnng a sample space as a set of outcomes from an experment. We denote ths by S. A random varable s a functon whch maps outcomes nto the
More informationA Robust Method for Calculating the Correlation Coefficient
A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal
More informationStat 543 Exam 2 Spring 2016
Stat 543 Exam 2 Sprng 2016 I have nether gven nor receved unauthorzed assstance on ths exam. Name Sgned Date Name Prnted Ths Exam conssts of 11 questons. Do at least 10 of the 11 parts of the man exam.
More informationNegative Binomial Regression
STATGRAPHICS Rev. 9/16/2013 Negatve Bnomal Regresson Summary... 1 Data Input... 3 Statstcal Model... 3 Analyss Summary... 4 Analyss Optons... 7 Plot of Ftted Model... 8 Observed Versus Predcted... 10 Predctons...
More informationANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)
Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of
More informationDimension Reduction and Visualization of the Histogram Data
The 4th Workshop n Symbolc Data Analyss (SDA 214): Tutoral Dmenson Reducton and Vsualzaton of the Hstogram Data Han-Mng Wu ( 吳漢銘 ) Department of Mathematcs Tamkang Unversty Tamsu 25137, Tawan http://www.hmwu.dv.tw
More informationCS286r Assign One. Answer Key
CS286r Assgn One Answer Key 1 Game theory 1.1 1.1.1 Let off-equlbrum strateges also be that people contnue to play n Nash equlbrum. Devatng from any Nash equlbrum s a weakly domnated strategy. That s,
More informationJanuary Examinations 2015
24/5 Canddates Only January Examnatons 25 DO NOT OPEN THE QUESTION PAPER UNTIL INSTRUCTED TO DO SO BY THE CHIEF INVIGILATOR STUDENT CANDIDATE NO.. Department Module Code Module Ttle Exam Duraton (n words)
More informationFaculty of Health Sciences. Cox regression. Torben Martinussen. Department of Biostatistics University of Copenhagen. 20. september 2012 Slide 1/51
Faculty of Health Sciences Cox regression Torben Martinussen Department of Biostatistics University of Copenhagen 2. september 212 Slide 1/51 Survival analysis Standard setup for right-censored survival
More informationLecture 2: Prelude to the big shrink
Lecture 2: Prelude to the bg shrnk Last tme A slght detour wth vsualzaton tools (hey, t was the frst day... why not start out wth somethng pretty to look at?) Then, we consdered a smple 120a-style regresson
More informationDurban Watson for Testing the Lack-of-Fit of Polynomial Regression Models without Replications
Durban Watson for Testng the Lack-of-Ft of Polynomal Regresson Models wthout Replcatons Ruba A. Alyaf, Maha A. Omar, Abdullah A. Al-Shha ralyaf@ksu.edu.sa, maomar@ksu.edu.sa, aalshha@ksu.edu.sa Department
More informationLINEAR REGRESSION ANALYSIS. MODULE VIII Lecture Indicator Variables
LINEAR REGRESSION ANALYSIS MODULE VIII Lecture - 7 Indcator Varables Dr. Shalabh Department of Maematcs and Statstcs Indan Insttute of Technology Kanpur Indcator varables versus quanttatve explanatory
More information4.1. Lecture 4: Fitting distributions: goodness of fit. Goodness of fit: the underlying principle
Lecture 4: Fttng dstrbutons: goodness of ft Goodness of ft Testng goodness of ft Testng normalty An mportant note on testng normalty! L4.1 Goodness of ft measures the extent to whch some emprcal dstrbuton
More informationChapter 4: Regression With One Regressor
Chapter 4: Regresson Wth One Regressor Copyrght 2011 Pearson Addson-Wesley. All rghts reserved. 1-1 Outlne 1. Fttng a lne to data 2. The ordnary least squares (OLS) lne/regresson 3. Measures of ft 4. Populaton
More informationBasic R Programming: Exercises
Basc R Programmng: Exercses RProgrammng John Fox ICPSR, Summer 2009 1. Logstc Regresson: Iterated weghted least squares (IWLS) s a standard method of fttng generalzed lnear models to data. As descrbed
More informationSee Book Chapter 11 2 nd Edition (Chapter 10 1 st Edition)
Count Data Models See Book Chapter 11 2 nd Edton (Chapter 10 1 st Edton) Count data consst of non-negatve nteger values Examples: number of drver route changes per week, the number of trp departure changes
More informationChapter 15 - Multiple Regression
Chapter - Multple Regresson Chapter - Multple Regresson Multple Regresson Model The equaton that descrbes how the dependent varable y s related to the ndependent varables x, x,... x p and an error term
More informationScientific Question Determine whether the breastfeeding of Nepalese children varies with child age and/or sex of child.
Longtudnal Logstc Regresson: Breastfeedng of Nepalese Chldren PART II GEE models (margnal, populaton average) covered last lab Random Intercept models (subject specfc) Transton models Scentfc Queston Determne
More informationLinear Regression Analysis: Terminology and Notation
ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented
More informationIntroduction to Econometrics (3 rd Updated Edition, Global Edition) Solutions to Odd-Numbered End-of-Chapter Exercises: Chapter 13
Introducton to Econometrcs (3 rd Updated Edton, Global Edton by James H. Stock and Mark W. Watson Solutons to Odd-Numbered End-of-Chapter Exercses: Chapter 13 (Ths verson August 17, 014 Stock/Watson -
More informationSample Size Calculation Based on the Semiparametric Analysis of Short-term and Long-term Hazard Ratios. Yi Wang
Sample Sze Calculaton Based on the Semparametrc Analyss of Short-term and Long-term Hazard Ratos Y Wang Submtted n partal fulfllment of the requrements for the degree of Doctor of Phlosophy under the Executve
More informationβ0 + β1xi. You are interested in estimating the unknown parameters β
Ordnary Least Squares (OLS): Smple Lnear Regresson (SLR) Analytcs The SLR Setup Sample Statstcs Ordnary Least Squares (OLS): FOCs and SOCs Back to OLS and Sample Statstcs Predctons (and Resduals) wth OLS
More informationPsychology 282 Lecture #24 Outline Regression Diagnostics: Outliers
Psychology 282 Lecture #24 Outlne Regresson Dagnostcs: Outlers In an earler lecture we studed the statstcal assumptons underlyng the regresson model, ncludng the followng ponts: Formal statement of assumptons.
More informationA note on almost sure behavior of randomly weighted sums of φ-mixing random variables with φ-mixing weights
ACTA ET COMMENTATIONES UNIVERSITATIS TARTUENSIS DE MATHEMATICA Volume 7, Number 2, December 203 Avalable onlne at http://acutm.math.ut.ee A note on almost sure behavor of randomly weghted sums of φ-mxng
More informationFirst Year Examination Department of Statistics, University of Florida
Frst Year Examnaton Department of Statstcs, Unversty of Florda May 7, 010, 8:00 am - 1:00 noon Instructons: 1. You have four hours to answer questons n ths examnaton.. You must show your work to receve
More informationStrong Markov property: Same assertion holds for stopping times τ.
Brownan moton Let X ={X t : t R + } be a real-valued stochastc process: a famlty of real random varables all defned on the same probablty space. Defne F t = nformaton avalable by observng the process up
More informationFinancing Innovation: Evidence from R&D Grants
Fnancng Innovaton: Evdence from R&D Grants Sabrna T. Howell Onlne Appendx Fgure 1: Number of Applcants Note: Ths fgure shows the number of losng and wnnng Phase 1 grant applcants over tme by offce (Energy
More informationx = , so that calculated
Stat 4, secton Sngle Factor ANOVA notes by Tm Plachowsk n chapter 8 we conducted hypothess tests n whch we compared a sngle sample s mean or proporton to some hypotheszed value Chapter 9 expanded ths to
More informationGeneralized Linear Methods
Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set
More informationPHYS 450 Spring semester Lecture 02: Dealing with Experimental Uncertainties. Ron Reifenberger Birck Nanotechnology Center Purdue University
PHYS 45 Sprng semester 7 Lecture : Dealng wth Expermental Uncertantes Ron Refenberger Brck anotechnology Center Purdue Unversty Lecture Introductory Comments Expermental errors (really expermental uncertantes)
More informationInterval Estimation in the Classical Normal Linear Regression Model. 1. Introduction
ECONOMICS 35* -- NOTE 7 ECON 35* -- NOTE 7 Interval Estmaton n the Classcal Normal Lnear Regresson Model Ths note outlnes the basc elements of nterval estmaton n the Classcal Normal Lnear Regresson Model
More information