Displaying 2 results from an estimated 2 matches for "lower95".
Did you mean:
lower59
2011 Nov 20
1
Cox proportional hazards confidence intervals
...ypes of endovascular interventions. I can successfully
obtain the LR, Wald, and Score test p-values from the coxph.object, as
well 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 confu...
2011 Jul 20
0
Cleveland Dot plots: tick labels and error bars
...for each point. The most direct
option I've observed is from:
http://www.unc.edu/courses/2010fall/ecol/563/001/docs/solutions/assign2.htm
It seems to use the following panel function to create 95% conf. intervals:
panel=function(x,y) {
panel.grid(v=0, h=-6, lty=3)
panel.segments(new.data$lower95, as.numeric(y), new.data$upper95,
as.numeric(y), lty=1, col=1)
panel.segments(new.data$lower50, as.numeric(y), new.data$upper50,
as.numeric(y), lty=1, lwd=4, col='grey60')
panel.xyplot(x, y, pch=16, cex=1.2, col='white')
panel.xyplot(x, y, pch=1, cex=1.1, col='black')
pa...