[R] Hmisc, Design, summary.Design plot- changing confidence intervals, adding color or decreasing font size

Elizabeth Stanny estanny at me.com
Sat Jul 18 02:31:29 CEST 2009


Hi,

1. I want 95% not 99% confidence intervals in my summary.Design plot using the Design package. Putting conf.int=.95 as an argument in plot does not work.  The  default appears  to be .99 not .95 as stated in the package Design manual (p. 164). 

2.  My sweave chuck is below and my output is attached as well as linked here:
http://www.sonoma.edu/users/s/stanny/330A/project/ciplot.pdf

3.  Alternatively, how do I make the confidence interval text smaller for the second plot (the interval text is cut-off) or  color the intervals? 

4.  I have 14 plots and would like to have 2 per page.

Thanks in advance.  

<<ghg07,fig=T>>=
par(mfrow=c(2,1))
s.fghg.new.06 <- summary(fghg.new.06)
s.fghg.06 <- summary(fghg.06)
dimnames(s.fghg.new.06)[[1]] <- newoddsl
dimnames(s.fghg.06)[[1]] <- oldoddsl

plot(s.fghg.06, log=T,
     main="Disclosed Accounting Method",
     at=c(.5,1,2,4,8,20),
     cex=.65,
     cex.c=.65,
     cex.t=.65)

plot(s.fghg.new.06, log=T,
     main="Did not not Disclose CDP06",
     at=c(.5,1,2,4,8,20),
     cex= .65,
     cex.c= .65,
     cex.t= .65)
@

where:
fghg.06 <- lrm(ghg.2006 ~ghg.2005+intsalep.2006 + EPA.2006 + toper.2006 + vol.2006 + lnass.2006,
               data=dw,  x=TRUE, y=TRUE)

fghg.new.06 <- lrm(ghg ~intsalep + EPA + toper + vol + lnass,
                data=newghg,subset=year==0, x=TRUE, y=TRUE)


Elizabeth Stanny
Professor of Business
Sonoma State University
stanny at sonoma.edu


More information about the R-help mailing list