Dear R-users:

I am trying to use the following code to reproduce the figures on page 140
of Prof. Frank Harrell's book 'Regression Modeling Strategies':

rm(list=ls())
options(width=128)
library(Hmisc)
library(Design)

getHdata(counties)
counties$older <- counties$age6574 + counties$age75
label(counties$older) <- '% age >= 65, 1990'
counties$pdensity <- log10(counties$pop.density+1)
label(counties$pdensity) <- 'log 10 of 1992 pop per 1990 miles^2'

dd <- datadist(counties)
options(datadist='dd')

f <- ols(democrat ~ rcs(pdensity,4) + rcs(pop.change,3) + rcs(older,3) +
crime + rcs(college,5)
         + rcs(income,4) + rcs(college,5) %ia% rcs(income,4) + rcs(farm,3) +
rcs(white,5) +
         rcs(turnout,3), data=counties)
f

r <- resid(f)

windows(width=10.67, height=6.60)
par(mfrow=c(3,3))
par(mar=c(4, 4, 2, 2))
plot(f, ylim=c(20,70))



But I have following problems:
1): the last two subfigures are different with those from the book;
2): I have a error reads:
Error in if (units == "") g(label) else if (label == "") g(units) else if
(plotmath &&  :
  missing value where TRUE/FALSE needed
3): I think there should be another subfigure on variable "turnout", but is
missing in book also

I am using most current version of R, Design and Hmisc.

Thanks

	[[alternative HTML version deleted]]

