[R] xYplot, error bars, log scale
Jacob Etches
jacob.etches at utoronto.ca
Tue Feb 5 23:26:07 CET 2008
I'm new to Hmisc and trying to 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|malef*ftf, groups=classval0,
as.table=T,
ylim=c(.8,4), pch=1:4, type="b", method="alt bars", label.curves=F,
lty.bar=3,
data=tphreg1_parameterestimates1987, subset=ft != 2 & parameter ==
"finc_mean_qu",
scales=list(
#y=list(log="e",at=c(1,1.25,1.5,2,2.5,3)),
x=list(at=seq(3,18,by=3))
),
key=list(columns=4,
text=list(c("Q0: Poorest","Q1","Q2","Q3")),
points=list(pch = 1:4,col = trellis.par.get("superpose.symbol")
$col[1:4])
)
)
Doctoral candidate, Epidemiology Program
Department of Public Health Sciences, University of Toronto Faculty of
Medicine
Research Associate
Institute for Work & Health
800-481 University Avenue, Toronto, Ontario, Canada M5G 2E9
T: 416.927.2027 ext. 2290
F: 416.927.4167
jetches at iwh.on.ca
www.iwh.on.ca
More information about the R-help
mailing list