search for: tphreg1_parameterestimates1987

...get the following to work, but if I un- comment the y-scale list (in order to get a log-scale for the hazard ratio), the error bars become strangely large. The dataframe is simply ODS output from TPHREG in SAS. Can someone point me towards what I'm sure is a naive error? Many thanks! tphreg1_parameterestimates1987$LCL <- exp(tphreg1_parameterestimates1987$estimate - qnorm(. 975)*tphreg1_parameterestimates1987$stderr) tphreg1_parameterestimates1987$UCL <- exp(tphreg1_parameterestimates1987$estimate + qnorm(. 975)*tphreg1_parameterestimates1987$stderr) xYplot(Cbind(exp(estimate),LCL,UCL) ~ lag|mal...