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]]

