Survival Data Analysis Parametric Models
|
|
- Vanessa Robinson
- 6 years ago
- Views:
Transcription
1 1 Survival Data Analysis Parametric Models January 21, 2015 Sandra Gardner, PhD Dalla Lana School of Public Health University of Toronto
2 2 January 21, 2015 Agenda Basic Parametric Models Review: hazard & cumulative hazard functions; likelihood function Proportional hazards versus accelerated failure Exponential model Weibull model Log-Normal model Log-Logistic model Checking assumptions Gamma model Goodness of fit and residuals Other Models Changepoint model (piecewise exponential model ) Reference: Matthews & Farewell 1982 Gamel-Boag (cure fraction) model Reference: Frankel & Longmate 2002 Bayesian analysis
3 3 January 21, 2015 Probability density function Random survivaltime T 0 f ( t) h( t) S ( t)
4 January 21, Hazard function Specifies the instantaneous rate of failure at T=t t t T t t T t P t h t ) ( ) ( lim 0 ) ( ) ( ) ( t S t f t h See K&M Section 2.3
5 January 21, Cumulative hazard function ) ( log ) (. ) ( ) (, ] [ ) ( 0 ) ( t S t H Note du u h t H where e t P T t S t u t H
6 Likelihood Full likelihood for parametric models Assuming censoring is independent of failure and noninformative 1 1 ) ( ) ( ), Pr( ), min( ), Pr( ) ( ( t S t f t and C X wheret t L C S x f L r n i i i R i r n D i i January 21, K&M 3.5.1
7 Likelihood ] [ ] [ ] [ ] [ ) ( exp )] [ ) exp )] [ ) exp[ ) )exp[ ) ( ( L i n i i t n i i t n i t i i n i i t H t h ds s h t h ds s h ds s h t h t S t f i i i January 21, K&M 3.5.3
8 8 January 21, 2015 Parametric Survival models Fully specified model with hazard rate a function of covariates (including intercept) Proportional Hazards (PH) constant hazard ratios across time Exponential, Weibull Accelerated Failure Models (AFT) constant time ratios across survival percentiles Exponential, Weibull, Log Normal, Log Logistic
9 PH versus AFT January 21, e x t x t TR e t x h t x h HR binary is X g e ) 0, ( ) 1, ( ) 0, ( ) 1, ( PH AFT Exponential Model
10 PH versus AFT January 21, ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ' ' ' ' ' t e S X t S e t e h X t h t S X t S e t h X t h X X X e X X PH AFT Be careful of parameterization of models in texts and software.
11 11 January 21, 2015 Sample Weibull hazard plots - HR=1.5 h(t) years x 0 1
12 12 January 21, 2015 Sample Weibull survival plots - TR=.67 (or AF=1.5) S(t) years x 0 1
13 13 January 21, 2015 Error distributions f ( ) exp( exp( )) f ( ) exp( ) Y logt X f ( ) e 1 e 2 Be careful of parameterization of models in texts and software.
14 January 21,
15 15 January 21, 2015 Exponential Model constant hazard functions both PH and AFT model underlying error function has an extreme value function with σ=1 S( t) e t h( t) ln(.5) Median 1 Mean.69
16 n i i n i i i n i i i n i t t l t mle ˆ ] [ ) ] log[ ( ) ( exp ] [ ) L( Exponential Model January 21,
17 January 21,
18 January 21,
19 January 21,
20 20 January 21, 2015 Weibull monotone increasing or decreasing hazard functions both PH and AFT model Exponential model is special case (γ=1) S( t) e t h( t) t Median 1 ln(.5) 1
21 January 21,
22 January 21,
23 January 21,
24 January 21,
25 January 21,
26 January 21,
27 Log Normal hazard functions rise to a maximum then slowly decline, AFT model only January 21, ) ln( 2 1 ) ( ) ( ) ( 2 1 ) ( ) ln( 1 ) ( e Mean e e Median t S t f t h e t t f t t S t
28 January 21,
29 January 21,
30 January 21,
31 January 21,
32 January 21,
33 January 21,
34 Log Logistic hazard functions rise to a maximum then slowly decline or are monotone decreasing, AFT model only January 21, ) ( 2 1) ( 1 1 ) ( ) ( ) ( ) (1 ) ( 1 1 ) ( Median t t t S t f t h t t t f t t S
35 January 21,
36 January 21,
37 January 21,
38 January 21,
39 January 21,
40 January 21,
41 41 January 21, 2015 Example Data Set Patients diagnosed with brain cancer are randomized to a treatment group versus placebo. N=222, with only 15 censored cases Mean age around 48 years and 64% male. Other covariates are available in data set.
42 January 21,
43 43 January 21, 2015 Overall Survival Estimated median=27.4 and mean=44.5
44 44 January 21, logs(t) Plot Plot versus t If a straight line then exponential model (H(t)=λt)
45 January 21,
46 46 January 21, 2015 log-logs(t) Plot Plot versus log(t) If a straight line then Weibull model H(t)=λt γ logh(t)=log(λ)+γlog(t)
47 January 21,
48 48 January 21, 2015 Probit Plot Plot Ф -1 (1-S(t)) versus log(t) If a straight line then Log Normal model S(t)=1-Ф((log(t)-u)/σ)
49 January 21,
50 50 January 21, 2015 Logit Plot Plot log((1-s(t)) /S(t)) versus log(t) Plot of odds of having the event by time t If a straight line then Log Logistic model S(t)=1/(1+αt γ )
51 January 21,
52 52 January 21, 2015 Other options Non-parametric smoothing of hazard function Probability plots Likelihood ratio tests of nested models (Gamma) Check distribution of t or log(t) for the noncensored cases
53 53 January 21, 2015 Smoothed Hazard Function
54 54 January 21, 2015 Smoothed Hazard Function (reduced range)
55 55 January 21, 2015 SAS Code (SAS 9.3 using ODS graphics) proc lifereg data=sda.brain; model weeks*event(0)=/d=exponential; probplot; inset; title 'LifeReg: Overall Survival - Probability Plot (Exponential)'; run; proc lifetest data=sda.brain plot=(survival(atrisk outside) hazard logsurv loglogs) notable; time weeks*event(0); title 'LifeTest: Overall Survival: hazard'; run;
56 56 January 21, 2015 Exponential Probability Plot
57 57 January 21, 2015 Weibull Probability Plot
58 58 January 21, 2015 Log Normal Probability Plot
59 59 January 21, 2015 Log Logistic Probability Plot
60 60 January 21, 2015 Gamma Model SAS fits the generalized 3-parameter model it can fit a Weibull (exponential) and log-normal model (test using likelihood ratio test) it can also fit a model with a U-shaped hazard function Survivor and hazard functions involve incomplete gamma functions
61 January 21,
62 January 21,
63 January 21,
64 64 January 21, 2015 Unadjusted model
65 65 January 21, 2015 Basic data summary Variable Sum event 207 weeks 9426 Estimated rate: 207/9426= ln( )= overall median % CI (23.14, 31.43) mean SE= 3.285
66 66 January 21, 2015 Exponential (Intercept only) -2 Log Likelihood= , AIC= Standard 95% Confidence Chi- Parameter DF Estimate Error Limits Square Pr > ChiSq Intercept <.0001 Scale Weibull Scale Weibull Shape Lagrange Multiplier Statistics Parameter Chi-Square Pr > ChiSq Scale λ = exp( ) = S(t) = exp(-λ*t) h(t) = λ Median = -ln(.5)/λ = 31.4 Mean = 1/λ = 45.5
67 67 January 21, 2015 Weibull(Intercept only) -2 Log Likelihood= AIC= Standard 95% Confidence Chi- Parameter DF Estimate Error Limits Square Pr > ChiSq Intercept <.0001 Scale * 1/shape Weibull Scale Weibull Shape * gamma λ = exp( * ) = γ = S(t) = exp(-λ*t**γ) h(t) = γ*λ*(t**(γ-1)) Median = (-ln(.5)/λ)**(1/γ) = 32.3
68 68 January 21, 2015 Log Normal(Intercept only) -2 Log Likelihood= , AIC= Standard 95% Confidence Chi- Parameter DF Estimate Error Limits Square Pr > ChiSq Intercept <.0001 Scale u = σ = S(t) = 1-Φ((ln(t)-u)/σ) f(t) = 1/(sqrt(2*π)*t*σ)*exp(-1/2*((ln(t)-u)/σ)**2) h(t) = f(t)/s(t) Median = exp(u) = 28.9 Mean = exp(u+0.5σ**2) = 45.7
69 69 January 21, 2015 Log Logistic (Intercept only) -2 Log Likelihood= , AIC= Standard 95% Confidence Chi- Parameter DF Estimate Error Limits Square Pr > ChiSq Intercept <.0001 Scale α = exp( / ) = γ = 1/ = S(t) = 1/(1+α*t**γ) f(t) = (α*γ*t**(γ-1))/(1+α*t**γ)**2 h(t) = f(t)/s(t) Median = (1/α)**(1/γ) = 27.8
70 January 21,
71 71 January 21, 2015 Gamma model (Intercept only) -2 Log Likelihood= , AIC= Standard 95% Confidence Chi- Parameter DF Estimate Error Limits Square Pr > ChiSq Intercept <.0001 Scale Shape If shape parameter is 0 then log-normal model If shape parameter is 1 then Weibull model If shape =1 and scale=1 then exponential model If shape and scale are equal, then standard gamma distribution Likelihood ratio test: Gamma vs log normal chi-square = = 7.667, p=0.006
72 72 January 21, 2015 Gamma hazard (Allison LIFEHAZ macro) Gamma-unadjusted
73 73 January 21, 2015 Comparison of Gamma model and Kaplan-Meier curve S(t) weeks PLOT Gamma KM
74 74 January 21, 2015 Model Comparison Model -2logL AIC AICC BIC Exponential Weibull LogNormal LogLogistic Gamma
75 75 January 21, 2015 Akaike Information Criteria AIC=-2log(Likelihood)+2(p+k) K&M k=1 (exponential) k=2 for Weibull, log logistic and log normal k=3 for generalized gamma In our example, AIC for gamma (606.3) is close to AIC for log-logistic (608.3). AICC AIC 2 p( n p 1) p 1 BIC 2log L p log( n)
76 76 January 21, 2015 Adjusted model Use preferred model building strategy to add covariates into the model (to be discussed further next month) Choosing two binary covariates for illustration Treated (treat=1); not treated (treat=0) Age <50 (age50=0) and age 50 (age50=1)
77 77 January 21, 2015 Covariates: median survival median treat=no : (20.57, 28.00) median treat=yes : (26.29, 37.00) median age<50 : (27.14, 39.71) median age>=50 : (19.00, 27.29)
78 January 21,
79 January 21,
80 January 21,
81 81 January 21, 2015 Exponential (treatment) Standard 95% Confidence Chi- Parameter DF Estimate Error Limits Square Pr > ChiSq Intercept <.0001 treat Scale Weibull Shape HR=exp(-beta)=exp( )=0.83 TR=exp(beta)=exp(0.1853)=1.20
82 82 January 21, 2015 Exponential (Hazard Ratio) Note closed form solution for hazard ratio can be calculated from the summary data below (unadjusted for other covariates): No treatment: 104/4300= (note log(0.0242)=-3.722) and Yes, treated: 103/5126= HR: / = Log(HR): log(0.0201/0.0242)= TR: exp(0.1856)=1.204 (output from proc means) treat Obs Variable Sum No 112 event 104 weeks 4300 Yes 110 event 103 weeks 5126
83 January 21,
84 January 21,
85 January 21,
86 86 January 21, 2015 Exponential/Weibull (age grouped) Standard 95% Confidence Chi- Parameter DF Estimate Error Limits Square Pr > ChiSq Intercept <.0001 age λ = exp(-( *age50)) Standard 95% Confidence Chi- Parameter DF Estimate Error Limits Square Pr > ChiSq Intercept <.0001 age Scale Weibull Shape λ = exp( *( *age50)) HR = exp(-beta*1.0688) = 1.69 TR = exp(beta) = 0.61 AF = 1.64
87 87 January 21, 2015 Goodness of fit Sample plots How well does model match Kaplan-Meier curves? Cox-Snell residuals Log-log(SDF) or cumulative hazard of residuals is a straight line? ˆ Other residuals: e.g. normal deviate residuals, see Nardi & Schemper r j and H ( T r Z where Hˆ is estimated j j ) distributed j exp(1) from data logti SAS output : log( S( ' xib ))
88 88 January 21, 2015 Goodness of fit Martingale residuals Klein & Moeschberger: estimate of the excess number of deaths seen in the data, but not predicted by model δ j -H(T j Z j ) i.e. δ j -r j Deviance residuals Klein & Moeschberger: more symmetric about 0 Transformed martingale residuals
89 January 21,
90 January 21,
91 January 21,
92 92 January 21, 2015 Simulated Exponential Data To show what plots look like using randomly generated data from an exponential distribution
93 January 21,
94 January 21,
95 January 21,
96 96 January 21, 2015 Simulated Exponential Data
97 97 January 21, 2015 Exponential (treatment and age) -2 Log Likelihood = Standard 95% Confidence Chi- Parameter DF Estimate Error Limits Square Pr > ChiSq Intercept <.0001 treat age Scale Weibull Shape
98 January 21,
99 99 January 21, 2015 Weibull(treatment and age) -2 Log Likelihood = Standard 95% Confidence Chi- Parameter DF Estimate Error Limits Square Pr > ChiSq Intercept <.0001 treat age Scale Weibull Shape
100 January 21,
101 101 January 21, 2015 Log Normal(treatment and age) -2 Log Likelihood = Standard 95% Confidence Chi- Parameter DF Estimate Error Limits Square Pr > ChiSq Intercept <.0001 treat age Scale
102 January 21,
103 103 January 21, 2015 Log Logistic(treatment and age) -2 Log Likelihood = Standard 95% Confidence Chi- Parameter DF Estimate Error Limits Square Pr > ChiSq Intercept <.0001 treat age Scale
104 January 21,
105 105 January 21, 2015 Residuals (SAS Code) /* Cox-Snell */ proc lifereg data=wsda.brain; model weeks*event(0)=treat age50/d=weibull; output out=wout cres=cres sres=sres p=predm std=stdm; title 'LifeReg: Treatment & Age groups - Weibull'; run; proc lifetest data=wout plots=(ls) notable; * looking for evidence that cres is exponential using the -log(s(t)) plot; * note that censoring value is maintained from original data set; time cres*event(0); title1 'Cox-Snell Residuals - Weibull'; run; /* martingale and deviance*/ lambda=exp(-( *treat *age50)); sexp=exp(-lambda*weeks); xbexp= *treat *age50; chexp=-log(sexp); martexp=event-chexp; devexp=sign(martexp)*(-2*(martexp+event*log(event-martexp)))**1/2;
106 106 January 21, 2015 Martingale Residual Plots - Weibull Model MR Weeks event 0 1
107 107 January 21, 2015 Martingale Residual Plots - Log Logistic Model MR Weeks event 0 1
108 108 January 21, 2015 Deviance Residual Plots - Weibull Model DR Weeks event 0 1
109 109 January 21, 2015 Deviance Residual Plots - Log Logistic Model DR Weeks event 0 1
110 110 January 21, 2015 Model summaries
111 January 21,
112 January 21,
113 January 21,
114 January 21,
115 January 21,
116 116 January 21, 2015 Exponential Summary Exponential Exponential Exponential Exponential Exponential t S(t), h(t), Median, S(t), h(t), (weeks) Age>=50 Treatment=1 Treatment=1 Treatment=1 Treatment=0 Treatment=0 26 Yes Yes Yes Exponential Exponential Exponential t Median, Hazard Exponential Exponential Median (weeks) Treatment=0 Ratio(t) beta(treatment); Time Ratio Ratio
117 117 January 21, 2015 Weibull Summary Weibull Weibull Weibull Weibull S(t), h(t), Weibull S(t), h(t), t Treatment= Treatment= Median, Treatment= Treatment= (weeks) Age>= Treatment= Yes Yes Yes Weibull Weibull Weibull Weibull t Median, Hazard Weibull Time Median (weeks) Treatment=0 Ratio(t) beta(treatment); Ratio Ratio
118 118 January 21, 2015 Log Normal Summary Log Normal Log Normal Log Normal Log Normal S(t), h(t), Log Normal S(t), h(t), t Treatment= Treatment= Median, Treatment= Treatment= (weeks) Age>= Treatment= Yes Yes Yes Log Normal Log Normal Log Normal t Median, Hazard Log Normal Log Normal Median (weeks) Treatment=0 Ratio(t) beta(treatment); Time Ratio Ratio beta=σ*(probit(1-slnorm0) - probit(1-slnorm1)); /* log normal scale and S(t) for each group */ TR=exp(beta)
119 119 January 21, 2015 Log Logistic Summary Log Logistic Log Logistic Log Logistic Log Logistic Log Logistic t S(t), h(t), Median, S(t), h(t), (weeks) Age>=50 Treatment=1 Treatment=1 Treatment=1 Treatment=0 Treatment=0 26 Yes Yes Yes Log Logistic Log Log Odds Log Logistic Logistic Log Logistic S(t)/(1-S(t) t Median, Hazard Log Logistic Logistic Median, (weeks) Treatment=0 Ratio(t) beta(treatment); Time Ratio Ratio Treatment= Log Logistic Odds Log S(t)/(1-S(t) Log Log Logistic Log Logistic Logisitic t, Logisitic Alpha, Alpha, Alpha (weeks) Treatment=0 Odds Ratio Treatment=1 Treatment=0 Ratio
120 120 January 21, 2015 Log logistic Odds Ratio (SAS Code) hrllog=hllog1/hllog0; if sllog1^=1 then odds1=sllog1/(1-sllog1); if sllog0^=1 then odds0=sllog0/(1-sllog0); oddsratio=odds1/odds0; alpharatio=alpha0/alpha1; beta_llog=log(oddsratio)*σ; /* log logistic scale; oddsratio=exp(beta_llog/σ) */ trllog=exp(beta_llog); mrllog=mllog1/mllog0;
121 121 January 21, 2015 Odds of survival=prob(alive)/prob(died) Odds ratio (treated to not treated)= 1.47 S(t) PLOT Treated Not Treated Prob(died, treated) Prob(alive, treated) Prob(died, not treated) Prob(alive, not treated) t
122 122 January 21, 2015 Discussion Not limited to parametric models discussed today. For example (next week): Changepoint model (piecewise exponential model) Gamel-Boag model (allows for a proportion of subjects to be long term survivors) Bayesian analysis
123 123 January 21, 2015 Changepoint model When the hazard rate is constant within in time periods and changes at known timepoint For example, brain cancer hazard rate is constant for the first year of follow up but hazard rate is reduced if patient survives at least one year. 1t S( t) e t e 1 e 2 ( t ) t
124 124 January 21, 2015 SAS code to restructure data data brain2(keep=id weeks event weeks2 event2 year1); set sda.brain; id=_n_; if weeks<=52 then do; event2=event; weeks2=weeks; year1=1; output; end; else do; event2=0; weeks2=52; year1=1; output; event2=event; weeks2=weeks-52; year1=0; output; end; run;
125 125 January 21, 2015 SAS code to fit the model proc lifereg data=brain2; model weeks2*event2(0)=year1/d=exponential; title 'Piecewise Exponential'; run; data brain3; do weeks=0 to 200 by 1; /* time frame */ lambda1=exp(-( )); * = ; lambda2=exp(-(4.533)); * = ; if weeks<=52 then sexp=exp(-lambda1*weeks); else sexp=exp(-lambda1*52)*exp(-lambda2*(weeks-52)); output; end; run;
126 126 January 21, 2015 Changepoint model Standard 95% Confidence Chi- Parameter DF Estimate Error Limits Square Pr > ChiSq Intercept <.0001 year <.0001 Scale Weibull Shape
127 January 21,
128 Gamel-Boag Model Allows for a proportion of subjects to be long term survivors. Events are modeled using log-normal model. January 21, ) (0, ~ ) ln( 1 ) ( ) ( )) ( (1 ) ( ) ( ' ' ' N e e x t e e x p x t S x p x p x t S i i i i x x f
129 129 January 21, 2015 SAS code to fit log-normal model /* log normal model using proc lifereg */ proc lifereg data=sda.brain; model weeks*event(0)=age50/d=lnormal; title 'LifeReg: Survival by age group (Log normal)'; run; /* log normal model using proc nlp (SAS/OR) - maximize loglikelihood */ proc nlp data=sda.brain tech=tr cov=2 stderr; parms int gamma sig; pi= ; u=int+gamma*age50; if event=1 then logl=log(1/(sqrt(2*pi)*weeks*sig)* exp(-1/2*((lweeks-u)/sig)**2)); if event=0 then logl=log(1-probnorm((lweeks-u)/sig)); max logl; run;
130 130 January 21, 2015 SAS code to fit Gamel-Boag model /* Gamel-Boag cure model using proc nlp (SAS/OR), maximize modified loglikelihood, reference Frankel & Longmate */ proc nlp data=sda.brain tech=tr cov=2 stderr; parms intg gamma intb beta sig; pi= ; p=exp(intb+beta*age50)/(1+exp(intb+beta*age50)); /* model proportion cured by age group */ u=intg+gamma*age50; if event=1 then logl=log((1-p)*(1/(sqrt(2*pi)*weeks*sig)* exp(-1/2*((lweeks-u)/sig)**2))); if event=0 then logl=log(p+(1-p)*(1-probnorm((lweeks-u)/sig))); max logl; run;
131 131 January 21, 2015 SAS output from Proc Lifereg Analysis of Maximum Likelihood Parameter Estimates Parameter DF Estimate Standard Error 95% Confidence Limits Chi-Square Pr > ChiSq Intercept <.0001 age Scale
132 132 January 21, 2015 SAS output from Proc NLP (1) Optimization Results Parameter Estimates N Parameter Estimate Approx Std Err t Value Approx Pr > t Gradient Objective Function 1 int E gamma sig E Value of Objective Function =
133 133 January 21, 2015 SAS output from Proc NLP (2) Optimization Results Parameter Estimates N Parameter Estimate Approx Std Err t Value Approx Pr > t Gradient Objective Function 1 intg E gamma intb E beta sig E Value of Objective Function = p1 p0 or lor
134 134 Survival: Age>=50 January 21, 2015 S(t) weeks PLOT LogNormal Gamel-Boag KM
135 135 Survival: Age<50 January 21, 2015 S(t) weeks PLOT LogNormal Gamel-Boag KM
136 136 January 21, 2015 Bayesian analysis Gibbs sampling used for the location-scale models Can add priors for model parameters Can output posterior samples proc lifereg data=sda.brain; model weeks*event(0)=age50/d=weibull; bayes WeibullShapePrior=gamma seed=1254 outpost=postweibull; run;
137 137 January 21, 2015 Analysis of Maximum Likelihood Parameter Estimates Standard 95% Confidence Parameter DF Estimate Error Limits Intercept age Scale Weibull Shape Posterior Summaries Standard Percentiles Parameter N Mean Deviation 25% 50% 75% Intercept age WeibShape Posterior Intervals Parameter Alpha Equal-Tail Interval HPD Interval Intercept age WeibShape
138 January 21,
139 139 January 21, 2015 Discussion: Why fit parametric models? Able to describe the hazard rate AF model alternative when hazard rates are non-proportional Easier and more convenient to predict outcome for a particular outcome (see Reid (1994) conversation with D.R. Cox) If underlying hazard function is correctly specified, then parametric models give more precise estimates (K & M, p.373). Applications where parametric models are compared to Cox proportional hazard models: Chapman et al (2006). Application of log-normal model which authors conclude has a major advantage over the Cox model Nardi and Schemper (2003). Authors compare Cox and parametric models in clinical settings. Carroll (2003). Author illustrates the practical benefits of a Weibullbased analysis.
140 140 January 21, 2015 References Applied Survival Analysis, D.W. Hosmer, S. Lemeshow, S. May, Wiley 2008 Survival Analysis: Techniques for Censored and Truncated Data, J.P.Klein, M.L. Moeschberger, Springer 1997 Nardi, A. and Schemper, M. (2003), Comparing Cox and parametric models in clinical studies, Statistics in Medicine, 22, J-A. W. Chapman, H.L.A Lickley, et al. Ascertaining Prognosis for Breast Cancer in Node- Negative Patients with Innovative Survival Analysis. The Breast Journal 2006, 12(1): Moran J. L., Bersten, A.D. et al. (2008). Modelling survival in acute severe illness: Cox versus accelerated failure time models. Journal of Evaluation in Clinical Practice Carroll, K. J. (2003). On the use and utility of the Weibull model in the analysis of survival data. Controlled Clinical Trials Matthews, D. E. and Farewell, V. T. (1982). On testing for constant hazard against a changepoint alternative. Biometrics 38, Frankel, P. and Longmate, J. (2002). Parametric models for accelerated and long-term survival: a comment on proportional hazards. Statistics in Medicine Reid, N. (1994). A Conversation with Sir David Cox. Statistical Science 9(3) Cox, D.R. and Oakes, D. (1984). Analysis of Survival Data. Chapman & Hall.
Estimation Procedure for Parametric Survival Distribution Without Covariates
Estimation Procedure for Parametric Survival Distribution Without Covariates The maximum likelihood estimates of the parameters of commonly used survival distribution can be found by SAS. The following
More informationSurvival Analysis APTS 2016/17 Preliminary material
Survival Analysis APTS 2016/17 Preliminary material Ingrid Van Keilegom KU Leuven (ingrid.vankeilegom@kuleuven.be) August 2017 1 Introduction 2 Common functions in survival analysis 3 Parametric survival
More informationDuration Models: Parametric Models
Duration Models: Parametric Models Brad 1 1 Department of Political Science University of California, Davis January 28, 2011 Parametric Models Some Motivation for Parametrics Consider the hazard rate:
More informationDuration Models: Modeling Strategies
Bradford S., UC-Davis, Dept. of Political Science Duration Models: Modeling Strategies Brad 1 1 Department of Political Science University of California, Davis February 28, 2007 Bradford S., UC-Davis,
More informationChapter 2 ( ) Fall 2012
Bios 323: Applied Survival Analysis Qingxia (Cindy) Chen Chapter 2 (2.1-2.6) Fall 2012 Definitions and Notation There are several equivalent ways to characterize the probability distribution of a survival
More information**BEGINNING OF EXAMINATION** A random sample of five observations from a population is:
**BEGINNING OF EXAMINATION** 1. You are given: (i) A random sample of five observations from a population is: 0.2 0.7 0.9 1.1 1.3 (ii) You use the Kolmogorov-Smirnov test for testing the null hypothesis,
More informationBayesian Multinomial Model for Ordinal Data
Bayesian Multinomial Model for Ordinal Data Overview This example illustrates how to fit a Bayesian multinomial model by using the built-in mutinomial density function (MULTINOM) in the MCMC procedure
More informationPractice Exam 1. Loss Amount Number of Losses
Practice Exam 1 1. You are given the following data on loss sizes: An ogive is used as a model for loss sizes. Determine the fitted median. Loss Amount Number of Losses 0 1000 5 1000 5000 4 5000 10000
More informationGamma Distribution Fitting
Chapter 552 Gamma Distribution Fitting Introduction This module fits the gamma probability distributions to a complete or censored set of individual or grouped data values. It outputs various statistics
More information1. You are given the following information about a stationary AR(2) model:
Fall 2003 Society of Actuaries **BEGINNING OF EXAMINATION** 1. You are given the following information about a stationary AR(2) model: (i) ρ 1 = 05. (ii) ρ 2 = 01. Determine φ 2. (A) 0.2 (B) 0.1 (C) 0.4
More informationStatistical Analysis of Life Insurance Policy Termination and Survivorship
Statistical Analysis of Life Insurance Policy Termination and Survivorship Emiliano A. Valdez, PhD, FSA Michigan State University joint work with J. Vadiveloo and U. Dias Sunway University, Malaysia Kuala
More informationConfidence Intervals for an Exponential Lifetime Percentile
Chapter 407 Confidence Intervals for an Exponential Lifetime Percentile Introduction This routine calculates the number of events needed to obtain a specified width of a confidence interval for a percentile
More informationThe Weibull in R is actually parameterized a fair bit differently from the book. In R, the density for x > 0 is
Weibull in R The Weibull in R is actually parameterized a fair bit differently from the book. In R, the density for x > 0 is f (x) = a b ( x b ) a 1 e (x/b) a This means that a = α in the book s parameterization
More informationModelling component reliability using warranty data
ANZIAM J. 53 (EMAC2011) pp.c437 C450, 2012 C437 Modelling component reliability using warranty data Raymond Summit 1 (Received 10 January 2012; revised 10 July 2012) Abstract Accelerated testing is often
More informationIntro to GLM Day 2: GLM and Maximum Likelihood
Intro to GLM Day 2: GLM and Maximum Likelihood Federico Vegetti Central European University ECPR Summer School in Methods and Techniques 1 / 32 Generalized Linear Modeling 3 steps of GLM 1. Specify the
More informationHomework Problems Stat 479
Chapter 10 91. * A random sample, X1, X2,, Xn, is drawn from a distribution with a mean of 2/3 and a variance of 1/18. ˆ = (X1 + X2 + + Xn)/(n-1) is the estimator of the distribution mean θ. Find MSE(
More informationThe comparison of proportional hazards and accelerated failure time models in analyzing the first birth interval survival data
Journal of Physics: Conference Series PAPER OPEN ACCESS The comparison of proportional hazards and accelerated failure time models in analyzing the first birth interval survival data To cite this article:
More informationThe Cox Hazard Model for Claims Data: a Bayesian Non-Parametric Approach
The Cox Hazard Model for Claims Data: a Bayesian Non-Parametric Approach Samuel Berestizhevsky, InProfix Inc, Boca Raton, FL Tanya Kolosova, InProfix Inc, Boca Raton, FL ABSTRACT General insurance protects
More informationLecture 21: Logit Models for Multinomial Responses Continued
Lecture 21: Logit Models for Multinomial Responses Continued Dipankar Bandyopadhyay, Ph.D. BMTRY 711: Analysis of Categorical Data Spring 2011 Division of Biostatistics and Epidemiology Medical University
More informationBayesian Hierarchical Modeling for Meta- Analysis
Bayesian Hierarchical Modeling for Meta- Analysis Overview Meta-analysis is an important technique that combines information from different studies. When you have no prior information for thinking any
More informationbook 2014/5/6 15:21 page 261 #285
book 2014/5/6 15:21 page 261 #285 Chapter 10 Simulation Simulations provide a powerful way to answer questions and explore properties of statistical estimators and procedures. In this chapter, we will
More informationBuilding and Checking Survival Models
Building and Checking Survival Models David M. Rocke May 23, 2017 David M. Rocke Building and Checking Survival Models May 23, 2017 1 / 53 hodg Lymphoma Data Set from KMsurv This data set consists of information
More informationPASS Sample Size Software
Chapter 850 Introduction Cox proportional hazards regression models the relationship between the hazard function λ( t X ) time and k covariates using the following formula λ log λ ( t X ) ( t) 0 = β1 X1
More informationAn Introduction to Event History Analysis
An Introduction to Event History Analysis Oxford Spring School June 18-20, 2007 Day Three: Diagnostics, Extensions, and Other Miscellanea Data Redux: Supreme Court Vacancies, 1789-1992. stset service,
More informationCalculating the Probabilities of Member Engagement
Calculating the Probabilities of Member Engagement by Larry J. Seibert, Ph.D. Binary logistic regression is a regression technique that is used to calculate the probability of an outcome when there are
More informationA Comprehensive, Non-Aggregated, Stochastic Approach to. Loss Development
A Comprehensive, Non-Aggregated, Stochastic Approach to Loss Development By Uri Korn Abstract In this paper, we present a stochastic loss development approach that models all the core components of the
More informationMaximum Likelihood Estimation
Maximum Likelihood Estimation EPSY 905: Fundamentals of Multivariate Modeling Online Lecture #6 EPSY 905: Maximum Likelihood In This Lecture The basics of maximum likelihood estimation Ø The engine that
More informationSOCIETY OF ACTUARIES/CASUALTY ACTUARIAL SOCIETY EXAM C CONSTRUCTION AND EVALUATION OF ACTUARIAL MODELS EXAM C SAMPLE QUESTIONS
SOCIETY OF ACTUARIES/CASUALTY ACTUARIAL SOCIETY EXAM C CONSTRUCTION AND EVALUATION OF ACTUARIAL MODELS EXAM C SAMPLE QUESTIONS Copyright 2008 by the Society of Actuaries and the Casualty Actuarial Society
More informationSubject CS1 Actuarial Statistics 1 Core Principles. Syllabus. for the 2019 exams. 1 June 2018
` Subject CS1 Actuarial Statistics 1 Core Principles Syllabus for the 2019 exams 1 June 2018 Copyright in this Core Reading is the property of the Institute and Faculty of Actuaries who are the sole distributors.
More informationINTRODUCTION TO SURVIVAL ANALYSIS IN BUSINESS
INTRODUCTION TO SURVIVAL ANALYSIS IN BUSINESS By Jeff Morrison Survival model provides not only the probability of a certain event to occur but also when it will occur... survival probability can alert
More informationboxcox() returns the values of α and their loglikelihoods,
Solutions to Selected Computer Lab Problems and Exercises in Chapter 11 of Statistics and Data Analysis for Financial Engineering, 2nd ed. by David Ruppert and David S. Matteson c 2016 David Ruppert and
More informationA Comparison of Univariate Probit and Logit. Models Using Simulation
Applied Mathematical Sciences, Vol. 12, 2018, no. 4, 185-204 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.818 A Comparison of Univariate Probit and Logit Models Using Simulation Abeer
More informationTo be two or not be two, that is a LOGISTIC question
MWSUG 2016 - Paper AA18 To be two or not be two, that is a LOGISTIC question Robert G. Downer, Grand Valley State University, Allendale, MI ABSTRACT A binary response is very common in logistic regression
More informationก ก ก ก ก ก ก. ก (Food Safety Risk Assessment Workshop) 1 : Fundamental ( ก ( NAC 2010)) 2 3 : Excel and Statistics Simulation Software\
ก ก ก ก (Food Safety Risk Assessment Workshop) ก ก ก ก ก ก ก ก 5 1 : Fundamental ( ก 29-30.. 53 ( NAC 2010)) 2 3 : Excel and Statistics Simulation Software\ 1 4 2553 4 5 : Quantitative Risk Modeling Microbial
More informationUsing New SAS 9.4 Features for Cumulative Logit Models with Partial Proportional Odds Paul J. Hilliard, Educational Testing Service (ETS)
Using New SAS 9.4 Features for Cumulative Logit Models with Partial Proportional Odds Using New SAS 9.4 Features for Cumulative Logit Models with Partial Proportional Odds INTRODUCTION Multicategory Logit
More informationMultivariate Cox PH model with log-skew-normal frailties
Multivariate Cox PH model with log-skew-normal frailties Department of Statistical Sciences, University of Padua, 35121 Padua (IT) Multivariate Cox PH model A standard statistical approach to model clustered
More informationA Comprehensive, Non-Aggregated, Stochastic Approach to Loss Development
A Comprehensive, Non-Aggregated, Stochastic Approach to Loss Development by Uri Korn ABSTRACT In this paper, we present a stochastic loss development approach that models all the core components of the
More informationDescribing Uncertain Variables
Describing Uncertain Variables L7 Uncertainty in Variables Uncertainty in concepts and models Uncertainty in variables Lack of precision Lack of knowledge Variability in space/time Describing Uncertainty
More informationCREDIT SCORING & CREDIT CONTROL XIV August 2015 Edinburgh. Aneta Ptak-Chmielewska Warsaw School of Ecoomics
CREDIT SCORING & CREDIT CONTROL XIV 26-28 August 2015 Edinburgh Aneta Ptak-Chmielewska Warsaw School of Ecoomics aptak@sgh.waw.pl 1 Background literature Hypothesis Data and methods Empirical example Conclusions
More informationChapter 4: Commonly Used Distributions. Statistics for Engineers and Scientists Fourth Edition William Navidi
Chapter 4: Commonly Used Distributions Statistics for Engineers and Scientists Fourth Edition William Navidi 2014 by Education. This is proprietary material solely for authorized instructor use. Not authorized
More informationModeling. joint work with Jed Frees, U of Wisconsin - Madison. Travelers PASG (Predictive Analytics Study Group) Seminar Tuesday, 12 April 2016
joint work with Jed Frees, U of Wisconsin - Madison Travelers PASG (Predictive Analytics Study Group) Seminar Tuesday, 12 April 2016 claim Department of Mathematics University of Connecticut Storrs, Connecticut
More informationSTA 4504/5503 Sample questions for exam True-False questions.
STA 4504/5503 Sample questions for exam 2 1. True-False questions. (a) For General Social Survey data on Y = political ideology (categories liberal, moderate, conservative), X 1 = gender (1 = female, 0
More informationQuantile Regression in Survival Analysis
Quantile Regression in Survival Analysis Andrea Bellavia Unit of Biostatistics, Institute of Environmental Medicine Karolinska Institutet, Stockholm http://www.imm.ki.se/biostatistics andrea.bellavia@ki.se
More informationproc genmod; model malform/total = alcohol / dist=bin link=identity obstats; title 'Table 2.7'; title2 'Identity Link';
BIOS 6244 Analysis of Categorical Data Assignment 5 s 1. Consider Exercise 4.4, p. 98. (i) Write the SAS code, including the DATA step, to fit the linear probability model and the logit model to the data
More informationPhd Program in Transportation. Transport Demand Modeling. Session 11
Phd Program in Transportation Transport Demand Modeling João de Abreu e Silva Session 11 Binary and Ordered Choice Models Phd in Transportation / Transport Demand Modelling 1/26 Heterocedasticity Homoscedasticity
More informationME3620. Theory of Engineering Experimentation. Spring Chapter III. Random Variables and Probability Distributions.
ME3620 Theory of Engineering Experimentation Chapter III. Random Variables and Probability Distributions Chapter III 1 3.2 Random Variables In an experiment, a measurement is usually denoted by a variable
More informationEconometric Methods for Valuation Analysis
Econometric Methods for Valuation Analysis Margarita Genius Dept of Economics M. Genius (Univ. of Crete) Econometric Methods for Valuation Analysis Cagliari, 2017 1 / 25 Outline We will consider econometric
More informationSome Characteristics of Data
Some Characteristics of Data Not all data is the same, and depending on some characteristics of a particular dataset, there are some limitations as to what can and cannot be done with that data. Some key
More informationCase Study: Applying Generalized Linear Models
Case Study: Applying Generalized Linear Models Dr. Kempthorne May 12, 2016 Contents 1 Generalized Linear Models of Semi-Quantal Biological Assay Data 2 1.1 Coal miners Pneumoconiosis Data.................
More informationProcess capability estimation for non normal quality characteristics: A comparison of Clements, Burr and Box Cox Methods
ANZIAM J. 49 (EMAC2007) pp.c642 C665, 2008 C642 Process capability estimation for non normal quality characteristics: A comparison of Clements, Burr and Box Cox Methods S. Ahmad 1 M. Abdollahian 2 P. Zeephongsekul
More informationSYSM 6304 Risk and Decision Analysis Lecture 2: Fitting Distributions to Data
SYSM 6304 Risk and Decision Analysis Lecture 2: Fitting Distributions to Data M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas Email: M.Vidyasagar@utdallas.edu September 5, 2015
More informationEstimation Appendix to Dynamics of Fiscal Financing in the United States
Estimation Appendix to Dynamics of Fiscal Financing in the United States Eric M. Leeper, Michael Plante, and Nora Traum July 9, 9. Indiana University. This appendix includes tables and graphs of additional
More informationEX-POST VERIFICATION OF PREDICTION MODELS OF WAGE DISTRIBUTIONS
EX-POST VERIFICATION OF PREDICTION MODELS OF WAGE DISTRIBUTIONS LUBOŠ MAREK, MICHAL VRABEC University of Economics, Prague, Faculty of Informatics and Statistics, Department of Statistics and Probability,
More informationClark. Outside of a few technical sections, this is a very process-oriented paper. Practice problems are key!
Opening Thoughts Outside of a few technical sections, this is a very process-oriented paper. Practice problems are key! Outline I. Introduction Objectives in creating a formal model of loss reserving:
More informationBackground. opportunities. the transformation. probability. at the lower. data come
The T Chart in Minitab Statisti cal Software Background The T chart is a control chart used to monitor the amount of time between adverse events, where time is measured on a continuous scale. The T chart
More informationCS 361: Probability & Statistics
March 12, 2018 CS 361: Probability & Statistics Inference Binomial likelihood: Example Suppose we have a coin with an unknown probability of heads. We flip the coin 10 times and observe 2 heads. What can
More informationINSTITUTE AND FACULTY OF ACTUARIES. Curriculum 2019 SPECIMEN EXAMINATION
INSTITUTE AND FACULTY OF ACTUARIES Curriculum 2019 SPECIMEN EXAMINATION Subject CS1A Actuarial Statistics Time allowed: Three hours and fifteen minutes INSTRUCTIONS TO THE CANDIDATE 1. Enter all the candidate
More informationjoint work with K. Antonio 1 and E.W. Frees 2 44th Actuarial Research Conference Madison, Wisconsin 30 Jul - 1 Aug 2009
joint work with K. Antonio 1 and E.W. Frees 2 44th Actuarial Research Conference Madison, Wisconsin 30 Jul - 1 Aug 2009 University of Connecticut Storrs, Connecticut 1 U. of Amsterdam 2 U. of Wisconsin
More informationFrequency Distribution Models 1- Probability Density Function (PDF)
Models 1- Probability Density Function (PDF) What is a PDF model? A mathematical equation that describes the frequency curve or probability distribution of a data set. Why modeling? It represents and summarizes
More informationCommonly Used Distributions
Chapter 4: Commonly Used Distributions 1 Introduction Statistical inference involves drawing a sample from a population and analyzing the sample data to learn about the population. We often have some knowledge
More informationLoss Simulation Model Testing and Enhancement
Loss Simulation Model Testing and Enhancement Casualty Loss Reserve Seminar By Kailan Shang Sept. 2011 Agenda Research Overview Model Testing Real Data Model Enhancement Further Development Enterprise
More informationis the bandwidth and controls the level of smoothing of the estimator, n is the sample size and
Paper PH100 Relationship between Total charges and Reimbursements in Outpatient Visits Using SAS GLIMMIX Chakib Battioui, University of Louisville, Louisville, KY ABSTRACT The purpose of this paper is
More informationThe Cox Hazard Model for Claims Data
The Cox Hazard Model for Claims Data By Samuel Berestizhevsky and Tanya Kolosova ABSTRACT Claim management requires applying statistical techniques in the analysis and interpretation of the claims data.
More informationPersonalized screening intervals for biomarkers using joint models for longitudinal and survival data
Personalized screening intervals for biomarkers using joint models for longitudinal and survival data Dimitris Rizopoulos, Jeremy Taylor, Joost van Rosmalen, Ewout Steyerberg, Hanneke Takkenberg Department
More informationQuantile Regression. By Luyang Fu, Ph. D., FCAS, State Auto Insurance Company Cheng-sheng Peter Wu, FCAS, ASA, MAAA, Deloitte Consulting
Quantile Regression By Luyang Fu, Ph. D., FCAS, State Auto Insurance Company Cheng-sheng Peter Wu, FCAS, ASA, MAAA, Deloitte Consulting Agenda Overview of Predictive Modeling for P&C Applications Quantile
More informationLogistic Regression. Logistic Regression Theory
Logistic Regression Dr. J. Kyle Roberts Southern Methodist University Simmons School of Education and Human Development Department of Teaching and Learning Logistic Regression The linear probability model.
More informationAppendix A. Selecting and Using Probability Distributions. In this appendix
Appendix A Selecting and Using Probability Distributions In this appendix Understanding probability distributions Selecting a probability distribution Using basic distributions Using continuous distributions
More informationMonitoring Accrual and Events in a Time-to-Event Endpoint Trial. BASS November 2, 2015 Jeff Palmer
Monitoring Accrual and Events in a Time-to-Event Endpoint Trial BASS November 2, 2015 Jeff Palmer Introduction A number of things can go wrong in a survival study, especially if you have a fixed end of
More information9. Logit and Probit Models For Dichotomous Data
Sociology 740 John Fox Lecture Notes 9. Logit and Probit Models For Dichotomous Data Copyright 2014 by John Fox Logit and Probit Models for Dichotomous Responses 1 1. Goals: I To show how models similar
More informationA Survival Analysis of the Approval of U.S. Patent Applications
Econometrics Working Paper EWP0707 ISSN 1485-6441 Department of Economics A Survival Analysis of the Approval of U.S. Patent Applications Ying Xie & David E. Giles* Department of Economics, University
More informationMaximum Likelihood Estimates for Alpha and Beta With Zero SAIDI Days
Maximum Likelihood Estimates for Alpha and Beta With Zero SAIDI Days 1. Introduction Richard D. Christie Department of Electrical Engineering Box 35500 University of Washington Seattle, WA 98195-500 christie@ee.washington.edu
More informationChapter 3 Statistical Quality Control, 7th Edition by Douglas C. Montgomery. Copyright (c) 2013 John Wiley & Sons, Inc.
1 3.1 Describing Variation Stem-and-Leaf Display Easy to find percentiles of the data; see page 69 2 Plot of Data in Time Order Marginal plot produced by MINITAB Also called a run chart 3 Histograms Useful
More informationIntroduction to General and Generalized Linear Models
Introduction to General and Generalized Linear Models Generalized Linear Models - IIIb Henrik Madsen March 18, 2012 Henrik Madsen () Chapman & Hall March 18, 2012 1 / 32 Examples Overdispersion and Offset!
More informationList of Examples. Chapter 1
REFERENCES 485 List of Examples Chapter 1 1.1 : 1.1: Bayes theorem in Case Control studies. DATA: imaginary. Page: 4. 1.2 : 1.2: Goals scored by the national football team of Greece in Euro 2004 (Poisson
More informationIntroduction to the Maximum Likelihood Estimation Technique. September 24, 2015
Introduction to the Maximum Likelihood Estimation Technique September 24, 2015 So far our Dependent Variable is Continuous That is, our outcome variable Y is assumed to follow a normal distribution having
More informationOne-Sample Cure Model Tests
Chapter 713 One-Sample Cure Model Tests Introduction This module computes the sample size and power of the one-sample parametric cure model proposed by Wu (2015). This technique is useful when working
More informationWC-5 Just How Credible Is That Employer? Exploring GLMs and Multilevel Modeling for NCCI s Excess Loss Factor Methodology
Antitrust Notice The Casualty Actuarial Society is committed to adhering strictly to the letter and spirit of the antitrust laws. Seminars conducted under the auspices of the CAS are designed solely to
More informationREINSURANCE RATE-MAKING WITH PARAMETRIC AND NON-PARAMETRIC MODELS
REINSURANCE RATE-MAKING WITH PARAMETRIC AND NON-PARAMETRIC MODELS By Siqi Chen, Madeleine Min Jing Leong, Yuan Yuan University of Illinois at Urbana-Champaign 1. Introduction Reinsurance contract is an
More informationAn Improved Skewness Measure
An Improved Skewness Measure Richard A. Groeneveld Professor Emeritus, Department of Statistics Iowa State University ragroeneveld@valley.net Glen Meeden School of Statistics University of Minnesota Minneapolis,
More informationCHAPTER 12 EXAMPLES: MONTE CARLO SIMULATION STUDIES
Examples: Monte Carlo Simulation Studies CHAPTER 12 EXAMPLES: MONTE CARLO SIMULATION STUDIES Monte Carlo simulation studies are often used for methodological investigations of the performance of statistical
More informationHierarchical Generalized Linear Models. Measurement Incorporated Hierarchical Linear Models Workshop
Hierarchical Generalized Linear Models Measurement Incorporated Hierarchical Linear Models Workshop Hierarchical Generalized Linear Models So now we are moving on to the more advanced type topics. To begin
More informationFitting parametric distributions using R: the fitdistrplus package
Fitting parametric distributions using R: the fitdistrplus package M. L. Delignette-Muller - CNRS UMR 5558 R. Pouillot J.-B. Denis - INRA MIAJ user! 2009,10/07/2009 Background Specifying the probability
More informationLAST SECTION!!! 1 / 36
LAST SECTION!!! 1 / 36 Some Topics Probability Plotting Normal Distributions Lognormal Distributions Statistics and Parameters Approaches to Censor Data Deletion (BAD!) Substitution (BAD!) Parametric Methods
More informationBy-Peril Deductible Factors
By-Peril Deductible Factors Luyang Fu, Ph.D., FCAS Jerry Han, Ph.D., ASA March 17 th 2010 State Auto is one of only 13 companies to earn an A+ Rating by AM Best every year since 1954! Agenda Introduction
More informationLecture Note: Analysis of Financial Time Series Spring 2017, Ruey S. Tsay
Lecture Note: Analysis of Financial Time Series Spring 2017, Ruey S. Tsay Seasonal Time Series: TS with periodic patterns and useful in predicting quarterly earnings pricing weather-related derivatives
More informationHigh-Frequency Data Analysis and Market Microstructure [Tsay (2005), chapter 5]
1 High-Frequency Data Analysis and Market Microstructure [Tsay (2005), chapter 5] High-frequency data have some unique characteristics that do not appear in lower frequencies. At this class we have: Nonsynchronous
More informationModeling Credit Risk of Portfolio of Consumer Loans
ing Credit Risk of Portfolio of Consumer Loans Madhur Malik * and Lyn Thomas School of Management, University of Southampton, United Kingdom, SO17 1BJ One of the issues that the Basel Accord highlighted
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -5 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -5 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Moments of a distribubon Measures of
More informationBasic notions of probability theory: continuous probability distributions. Piero Baraldi
Basic notions of probability theory: continuous probability distributions Piero Baraldi Probability distributions for reliability, safety and risk analysis: discrete probability distributions continuous
More informationAustralian Journal of Basic and Applied Sciences. Conditional Maximum Likelihood Estimation For Survival Function Using Cox Model
AENSI Journals Australian Journal of Basic and Applied Sciences Journal home page: wwwajbaswebcom Conditional Maximum Likelihood Estimation For Survival Function Using Cox Model Khawla Mustafa Sadiq University
More informationMAINTAINABILITY DATA DECISION METHODOLOGY (MDDM)
TECHNICAL REPORT NO. TR-2011-19 MAINTAINABILITY DATA DECISION METHODOLOGY (MDDM) June 2011 APPROVED FOR PUBLIC RELEASE; DISTRIBUTION IS UNLIMITED U.S. ARMY MATERIEL SYSTEMS ANALYSIS ACTIVITY ABERDEEN PROVING
More informationAlastair Hall ECG 790F: Microeconometrics Spring Computer Handout # 2. Estimation of binary response models : part II
Alastair Hall ECG 790F: Microeconometrics Spring 2006 Computer Handout # 2 Estimation of binary response models : part II In this handout, we discuss the estimation of binary response models with and without
More informationStatistics & Flood Frequency Chapter 3. Dr. Philip B. Bedient
Statistics & Flood Frequency Chapter 3 Dr. Philip B. Bedient Predicting FLOODS Flood Frequency Analysis n Statistical Methods to evaluate probability exceeding a particular outcome - P (X >20,000 cfs)
More informationInternational Journal of Scientific and Research Publications, Volume 6, Issue 12, December ISSN
International Journal of Scientific and Research Publications, Volume 6, Issue 12, December 2016 61 Frequentist Comparison of the Bayesian Credible and Maximum Likelihood Confidence for the Median of the
More informationStatistics and Finance
David Ruppert Statistics and Finance An Introduction Springer Notation... xxi 1 Introduction... 1 1.1 References... 5 2 Probability and Statistical Models... 7 2.1 Introduction... 7 2.2 Axioms of Probability...
More informationModel fit assessment via marginal model plots
The Stata Journal (2010) 10, Number 2, pp. 215 225 Model fit assessment via marginal model plots Charles Lindsey Texas A & M University Department of Statistics College Station, TX lindseyc@stat.tamu.edu
More informationMultinomial Logit Models for Variable Response Categories Ordered
www.ijcsi.org 219 Multinomial Logit Models for Variable Response Categories Ordered Malika CHIKHI 1*, Thierry MOREAU 2 and Michel CHAVANCE 2 1 Mathematics Department, University of Constantine 1, Ain El
More informationFiring Costs, Employment and Misallocation
Firing Costs, Employment and Misallocation Evidence from Randomly Assigned Judges Omar Bamieh University of Vienna November 13th 2018 1 / 27 Why should we care about firing costs? Firing costs make it
More informationMultiple Regression and Logistic Regression II. Dajiang 525 Apr
Multiple Regression and Logistic Regression II Dajiang Liu @PHS 525 Apr-19-2016 Materials from Last Time Multiple regression model: Include multiple predictors in the model = + + + + How to interpret the
More informationWeb Appendix. Are the effects of monetary policy shocks big or small? Olivier Coibion
Web Appendix Are the effects of monetary policy shocks big or small? Olivier Coibion Appendix 1: Description of the Model-Averaging Procedure This section describes the model-averaging procedure used in
More information