John Smith
2007-Nov-24 15:41 UTC
[R] Hmisc: can not reproduce figure 4 of Statistical Tables and Plots using S and LATEX
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]]
Frank E Harrell Jr
2007-Nov-25 15:56 UTC
[R] Hmisc: can not reproduce figure 4 of Statistical Tables and Plots using S and LATEX
John Smith wrote:> 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))Add ,reset.par=FALSE before the final ) above. This will be correct in the next version of the document. Next time please follow the posting guide with regard to specifying information about the version of R, platform, and the description of the problem you encountered. Note that you did reproduce the code in the document exactly as it produced a pdf file, not an on-screen graphic. Frank> } > mtitle(paste('Quartiles of', label(y)), cex.m=1.5) > > > Could any one help me? > > Thanks > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at 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