Dear R-users:
I can not reproduce figure 4 of *Statistical Tables and Plots using S and
LATEX* by Prof. Frank Harrell with the following code:
rm(list=ls())
library(Hmisc)
getHdata(pbc)
attach(pbc)
age.groups <- cut2(age, c(45,60))
g <- function(y) apply(y, 2, quantile, c(.25,.5,.75))
y <- with(pbc, cbind(Chol=chol,Bili=bili))
# You can give new column names that are not legal S names
# by enclosing them in quotes, e.g. 'Chol (mg/dl)'=chol
vars <- with(pbc, c(label(chol), label(bili)))
label(y) <- paste(vars, collapse=' and ') # Will make nice caption in table
s3 <- summary(y ~ age.groups + ascites, fun=g, data=pbc)
s3
windows(width=10.67, height=6.60)
par(mfrow=c(1,2), oma=c(3,0,3,0))
for(ivar in 1:2)
{
isub <- (1:3)+(ivar-1)*3
print(isub)
plot(s3, which=isub, main='', xlab=vars[ivar], pch=c(91,16,93))
}
mtitle(paste('Quartiles of', label(y)), cex.m=1.5)
Could any one help me?
Thanks
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
