search for: logtest

...a parametric survival model (with aftreg {eha}) and the output tells me the overall p-value of my model is < 0.001. My simple question is: Does the result mean my model fits the data well OR does it mean my model DOES NOT fit the data well? Some side information how the p-value is calculated: logtest <- -2 * (x$loglik[1] - x$loglik[2]) pvalue <- 1-pchisq(logtest,df) with x$loglik[1] = -274 x$loglik[2] = -235 df = 25 I know the answer would probably be "read the manual" but I found different opinions on the web and want to make sure I am interpreting it correctly. Thanks Thoma...

Hi, I have a question about how to extract paramters from a fitted model. I can extract coefficients and std, but from some other statistics, I dont know how to extract. Can anyone help? Here it is an example: > coxout<-coxph(Surv(t,t.censor)~x) > coxout Call: coxph(formula = Surv(t, t.censor) ~ x) coef exp(coef) se(coef) z p x 0.349 1.42 0.257 1.36 0.17 Likelihood

...as the hazard ratio as follows: formula.obj = Surv(days, status) ~ type coxph.model = coxph(formula.obj, df) fit = summary(coxph.model) hazard.ratio = fit$conf.int[1] lower95 = fit$conf.int[3] upper95 = fit$conf.int[4] logrank.p.value = fit$sctest[3] wald.p.value = fit$waldtest[3] lr.p.value = fit$logtest[3] I had intended to report logrank P values with the hazard ratio and CI obtained from this function. In one case the P was 0.04 yet the CI crossed one, which confused me, and certainly will raise questions by reviewers. In retrospect I can see that the CI calculated by coxph is intimately rela...

Dear List, How do I extract the approximate Wald test for the frailty (in the following example 17.89 value)? What about the P-values, other Chisq, DF, se(coef) and se2? How can they be extracted? ######################################################> kfitm1 Call: coxph(formula = Surv(time, status) ~ age + sex + disease + frailty(id, dist = "gauss"), data = kidney)