John Smith wrote: > 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
John, I could reproduce the problem and this indicates a significant bug in plot.Design in some past version. I am almost certain the current version is correct because it now produces output that is consistent with the odds ratio chart. On the bottom row of the figure on p. 140 the voter turnout variable was missing, and as you mentioned the plot for % white is off. The plot for farm population shows a line that is too high. I have fixed the units error and will send you a corrected plot.Design in a separate e-mail. Thanks for reporting these problems. Frank > > [[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. > -- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University ______________________________________________ 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.
