θ = θ Π Π Parametric counting process models θ θ θ Log-likelihood: Consider counting processes: Score functions:

Similar documents
The Poisson Process Properties of the Poisson Process

Survival Prediction Based on Compound Covariate under Cox Proportional Hazard Models

Fault Tolerant Computing. Fault Tolerant Computing CS 530 Probabilistic methods: overview

International Journal Of Engineering And Computer Science ISSN: Volume 5 Issue 12 Dec. 2016, Page No.

Continuous Time Markov Chains

Midterm Exam. Tuesday, September hour, 15 minutes

( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model

Determination of Antoine Equation Parameters. December 4, 2012 PreFEED Corporation Yoshio Kumagae. Introduction

As evident from the full-sample-model, we continue to assume that individual errors are identically and

The Linear Regression Of Weighted Segments

Assessing Normality. Assessing Normality. Assessing Normality. Assessing Normality. Normal Probability Plot for Normal Distribution.

(1) Cov(, ) E[( E( ))( E( ))]

RATIO ESTIMATORS USING CHARACTERISTICS OF POISSON DISTRIBUTION WITH APPLICATION TO EARTHQUAKE DATA

Solution set Stat 471/Spring 06. Homework 2

Final Exam Applied Econometrics

Lecture 3 Topic 2: Distributions, hypothesis testing, and sample size determination

Fundamentals of Speech Recognition Suggested Project The Hidden Markov Model

Supplement Material for Inverse Probability Weighted Estimation of Local Average Treatment Effects: A Higher Order MSE Expansion

Chapter 3: Maximum-Likelihood & Bayesian Parameter Estimation (part 1)

Density estimation III.

AML710 CAD LECTURE 12 CUBIC SPLINE CURVES. Cubic Splines Matrix formulation Normalised cubic splines Alternate end conditions Parabolic blending

Real-time Classification of Large Data Sets using Binary Knapsack

Density estimation. Density estimations. CS 2750 Machine Learning. Lecture 5. Milos Hauskrecht 5329 Sennott Square

Least squares and motion. Nuno Vasconcelos ECE Department, UCSD

14. Poisson Processes

Chapter 8. Simple Linear Regression

Cyclone. Anti-cyclone

IMPROVED PORTFOLIO OPTIMIZATION MODEL WITH TRANSACTION COST AND MINIMAL TRANSACTION LOTS

To Estimate or to Predict

Density estimation III.

Solution. The straightforward approach is surprisingly difficult because one has to be careful about the limits.

Regression Approach to Parameter Estimation of an Exponential Software Reliability Model

EE 6885 Statistical Pattern Recognition

FORCED VIBRATION of MDOF SYSTEMS

Density estimation III. Linear regression.

EE 6885 Statistical Pattern Recognition

Real-Time Systems. Example: scheduling using EDF. Feasibility analysis for EDF. Example: scheduling using EDF

Probability Bracket Notation and Probability Modeling. Xing M. Wang Sherman Visual Lab, Sunnyvale, CA 94087, USA. Abstract

Linear Regression Linear Regression with Shrinkage

The ray paths and travel times for multiple layers can be computed using ray-tracing, as demonstrated in Lab 3.

FALL HOMEWORK NO. 6 - SOLUTION Problem 1.: Use the Storage-Indication Method to route the Input hydrograph tabulated below.

QR factorization. Let P 1, P 2, P n-1, be matrices such that Pn 1Pn 2... PPA

-distributed random variables consisting of n samples each. Determine the asymptotic confidence intervals for

Least Squares Fitting (LSQF) with a complicated function Theexampleswehavelookedatsofarhavebeenlinearintheparameters

4. THE DENSITY MATRIX

Partial Molar Properties of solutions

Pricing of CDO s Based on the Multivariate Wang Transform*

Comparison of the Bayesian and Maximum Likelihood Estimation for Weibull Distribution

Mixed Integral Equation of Contact Problem in Position and Time

COMPARISON OF ESTIMATORS OF PARAMETERS FOR THE RAYLEIGH DISTRIBUTION

An Efficient Dual to Ratio and Product Estimator of Population Variance in Sample Surveys

VARIATIONAL ITERATION METHOD FOR DELAY DIFFERENTIAL-ALGEBRAIC EQUATIONS. Hunan , China,

Optimal Eye Movement Strategies in Visual Search (Supplement)

Moments of Order Statistics from Nonidentically Distributed Three Parameters Beta typei and Erlang Truncated Exponential Variables

Quantum Mechanics II Lecture 11 Time-dependent perturbation theory. Time-dependent perturbation theory (degenerate or non-degenerate starting state)

Other Topics in Kernel Method Statistical Inference with Reproducing Kernel Hilbert Space

Nature and Science, 5(1), 2007, Han and Xu, Multi-variable Grey Model based on Genetic Algorithm and its Application in Urban Water Consumption

Fully Fuzzy Linear Systems Solving Using MOLP

Queuing Theory: Memory Buffer Limits on Superscalar Processing

Use of Non-Conventional Measures of Dispersion for Improved Estimation of Population Mean

Imputation Based on Local Linear Regression for Nonmonotone Nonrespondents in Longitudinal Surveys

EDUCATION COMMITTEE OF THE SOCIETY OF ACTUARIES ADVANCED TOPICS IN GENERAL INSURANCE STUDY NOTE CREDIBILITY WITH SHIFTING RISK PARAMETERS

ENGI 3423 Simple Linear Regression Page 12-01

The Mean Residual Lifetime of (n k + 1)-out-of-n Systems in Discrete Setting

Some Probability Inequalities for Quadratic Forms of Negatively Dependent Subgaussian Random Variables

The Bernstein Operational Matrix of Integration

Suppose we have observed values t 1, t 2, t n of a random variable T.

Survival Analysis for Randomized Clinical Trials II Cox Regression. Ziad Taib Biostatistics AstraZeneca February 26, 2008

Quantitative Portfolio Theory & Performance Analysis

Multiple Choice Test. Chapter Adequacy of Models for Regression

Photon Statistics. photons/s. Light beam stream of photons Difference with classical waves relies on different statistical properties

Lecture 8: Linear Regression

STK3100 and STK4100 Autumn 2017

THE WEIBULL LENGTH BIASED DISTRIBUTION -PROPERTIES AND ESTIMATION-

ESTIMATION AND TESTING

Multiphase Flow Simulation Based on Unstructured Grid

Department of Economics University of Toronto

Asymptotic Behavior of Solutions of Nonlinear Delay Differential Equations With Impulse

Spike-and-Slab Dirichlet Process Mixture Models

STK3100 and STK4100 Autumn 2018

PubH 7440 Spring 2010 Midterm 2 April

Linear Minimum Variance Unbiased Estimation of Individual and Population slopes in the presence of Informative Right Censoring

Decompression diagram sampler_src (source files and makefiles) bin (binary files) --- sh (sample shells) --- input (sample input files)

NUMERICAL SOLUTIONS TO ORDINARY DIFFERENTIAL EQUATIONS

Chapter 5 Properties of a Random Sample

Common MidPoint (CMP) Records and Stacking

Fresnel Equations cont.

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

For the plane motion of a rigid body, an additional equation is needed to specify the state of rotation of the body.

Analysis of Variance with Weibull Data

8. Queueing systems lect08.ppt S Introduction to Teletraffic Theory - Fall

Interval Estimation. Consider a random variable X with a mean of X. Let X be distributed as X X

The t copula with Multiple Parameters of Degrees of Freedom: Bivariate Characteristics and Application to Risk Management

THE ROYAL STATISTICAL SOCIETY HIGHER CERTIFICATE

Research on portfolio model based on information entropy theory

NOTE ON SIMPLE AND LOGARITHMIC RETURN

Midterm Exam 1, section 2 (Solution) Thursday, February hour, 15 minutes

3.1 Introduction to Multinomial Logit and Probit

Analysis of a Stochastic Lotka-Volterra Competitive System with Distributed Delays

Orbital Euclidean stability of the solutions of impulsive equations on the impulsive moments

Transcription:

Paramerc coug process models Cosder coug processes: N,,..., ha cou he occurreces of a eve of eres for dvduals Iesy processes: Lelhood λ ( ;,,..., N { } λ < Log-lelhood: l( log L( Score fucos: U ( l( log λ ( ; dn λ ( ; d λ λ log ( ; dn ( ; d he MLE ˆ ( ˆ may be foud by maxmzg,..., ˆ q he log-lelhood or by solvg he lelhood equaos (umercally, f eeded U (,,..., q ˆ ( ˆ ˆ,..., q s approxmaely mulvarae ormally dsrbued aroud s rue value wh a covarace marx ha may be esmaed as Ι ( ˆ where Ι ( s he observed formao marx wh elemes Example.: Pecewse cosa hazard ad occurrece/exposure raes - Assume ha λ ( ; Y for < h ( U ( l( h h he lelhood rao, score ad Wald ess apply as usual Iroduce: O N ( N ( "occurrece" Y ( u du "exposure"

he lelhood may be wre N { } λ < N { Y } < Π ( ( exp Y O exp Log-lelhood: l( O log Score fucos U ( l( he MLEs become he occurrece/exposure raes: ˆ O he observed formao marx has elemes h ( U ( h O { h } Ι ( ( O for h for h 6 hus O O Ι ( dag,..., Now ˆ ( ˆ ˆ,..., q s approxmaely mulvarae ormally dsrbued aroud s rue value wh a covarace marx ha may be esmaed as ˆ O O Ι ( dag,..., ˆ ˆ se( ˆ ˆ O ˆ ˆ dag,..., O O I follows ha occurrece/exposure raes are approxmaely depede ad ormally dsrbued wh sadard errors ˆ 7 Posso regresso Assume proporoal hazards: α ( x α( ; exp( β x wh pecewse cosa basele hazard: α ( ; for < Iroduce: O N ( N ( Y ( u du 8

Lelhood: he lelhood s proporoal o he lelhood of "depede Posso varables" O wh "parameers" e N { } λ < β x N { } < β x Π ( e Y exp e Y β x O {( exp β x β x e ( e } O {( exp β x β x e ( e } 9 Esmao procedure: F model by GLM-sofware reag he O as "depede Posso varables" wh "parameers" wh { ψ β x log } β x e exp + + ψ log Use logarhmc l ad log as offses me erval s reaed as a caegorcal covarae A dvdual corbues oe record o he daa fle for each me erval whe he/she s a rs Posso regresso caegorcal covaraes Assume he covarae vecors x ca oly aa L dsc values: x (,, x (L ( l ( l O ( l β x ( ( β x l ( L( β, e exp e l Suffce sascs (sums are over wh x x (l l O O ad ( l ( l ( l ( l O ( l β x ( l ( ( β x l e exp( e Moraly hree Norwega coues Basele Grea compuaoal savgs for large daa ses (Number of records aggregaed daa se: 6**6

Shapes of hazard raes.8.6....... x x Moraly of paes wh myocardal farco. Dvorce raes. Moraly of cacer paes. Icdece of chldhood leuema. Geeral moraly. Icdece of mos cacers. How ca he dffere paers be erpreed? We eed o dsgush bewee he dvdual hazards ad he populao hazard Heerogeey survval aalyss Dffereces bewee dvduals are of wo ds: hose observed hrough covaraes hose ha are uobserved he laer oes are usually dsregarded survval aalyss hs may lead o dsoros as explaed he fraly heory here s a large leraure o fraly heory, bu we wll revew us a few basc ssues Populao hazard ad dvdual hazard We wll dsgush bewee he populao hazard rae ad he dvdual hazard raes (deals below he populao hazard rae s flueced by seleco: hose wh hghes rs experece he eve early he shape of he populao hazard may be erely dffere from ha of he dvdual hazards Hece he populao hazard ca o be erpreed as gvg formao o dvdual developme rs A smple example Assume a populao s composed of wo groups, a "low rs group" ad a "hgh rs group" A he ouse, he proporos he wo groups are p ad p - p Low rs group has hazard α { α u du} S exp ( Hgh rs group has hazard α { α } S exp ( u du ad survval fuco ad survval fuco 6

Illusrao: Populao survval fuco S p S + p S..... Hazard raes 6 8 me ed: low rs Blac: hgh rs....6.8. Survval fucos 6 8 me 7 Populao hazard rae ( µ S S wh ps + ps p S + p S ps S ps S + ps + ps ( ( + ( ( S p S p S S { } w α + w α w ps p S + p S 8 Illusrao: Hazard..... 6 8 me ed: low rs Blac: hgh rs Blue: populao 9