Displaying 4 results from an estimated 4 matches for "logtest".
Did you mean:
logdest
2010 Oct 01
2
Small p-value good or bad?
...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...
2008 Apr 29
2
Help on extract paramters from fitted models
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
2011 Nov 20
1
Cox proportional hazards confidence intervals
...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...
2007 Apr 17
3
Extracting approximate Wald test (Chisq) from coxph(..frailty)
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)