Dear List, Must be a silly question, but I was wondering whether there is a direct way of calculating "standard error of a HR or exp(coef)" from coxph objects x <- coxph(Surv(time, status) ~ age + inst, lung)> x coef exp(coef) se(coef) z page 0.0190 1.02 0.00925 2.06 0.04inst -0.0104 0.99 0.01028 -1.01 0.31 cheers, Ehsan [[alternative HTML version deleted]]
On Jul 23, 2011, at 1:41 AM, Ehsan Karim wrote:> > Dear List, > Must be a silly question, but I was wondering whether there is a > direct way of calculating "standard error of a HR or exp(coef)" from > coxph objects > x <- coxph(Surv(time, status) ~ age + inst, lung)> x coef > exp(coef) se(coef) z page 0.0190 1.02 0.00925 2.06 > 0.04inst -0.0104 0.99 0.01028 -1.01 0.31 >?predict.coxph (The third example is precisely what you are requesting, at least under the assumption that you are considering a discrete predictor or want the HR for a one unit change in a continuous one.)> cheers, > EhsanDavid Winsemius, MD West Hartford, CT
Dear List: Wondering how to get around the following problem: any suggestions are welcome. Cheers, Ehsan> install.packages("geepack")Warning in install.packages("geepack") :? argument 'lib' is missing: using '/home/grad/student/Library/R/2.11/library'--- Please select a CRAN mirror for use in this session ---Loading Tcl/Tk interface ... doneError in structure(.External("dotTclObjv", objv, PACKAGE = "tcltk"), class = "tclObj") :?? [tcl] invalid command name "toplevel".In addition: Warning message:In fun(...) : no display name and no $DISPLAY environment variable> R.Version()$platform[1] "x86_64-apple-darwin9.8.0"$arch[1] "x86_64"$os[1] "darwin9.8.0"$system[1] "x86_64, darwin9.8.0"$status[1] "Under development (unstable)"$major[1] "2"$minor[1] "11.0"$year[1] "2010$month[1] "03"$day[1] "14"$`svn rev`[1] "51276"$language[1] "R"$version.string[1] "R version 2.11.0 Under development (unstable) (2010-03-14 r51276)"