Here's an example. The problem is that hexbin uses grid graphics, not
the base graphics that plot.svysmooth() uses, so you need to do the
plotting yourself with grid. The necessary viewport is returned by
svyplot()
hexstuff<-svyplot(api00~api99, design=dstrat, style="hex", xlab="1999
API",ylab="2000 API")
smoothh<-svysmooth(api00~api99, design=dstrat, method="quantreg",quantile=0.75)
smoothl<-svysmooth(api00~api99, design=dstrat, method="quantreg",quantile=0.25)
pushHexport(hexstuff$plot.vp)
grid.lines(smoothh$api99$x,smoothh$api99$y,default.units="native",gp=gpar(col="purple",lwd=2))
grid.lines(smoothl$api99$x,smoothl$api99$y,default.units="native",gp=gpar(col="purple",lwd=2))
popViewport()
-thomas
On Fri, Oct 19, 2012 at 6:49 AM, Durant, James T. (ATSDR/DCHI/SSB)
<h...@cdc.gov> wrote:
> Hi all-
>
> So sorry to bother you all with something pretty basic.
>
> I am trying to add the lines method output from svysmooth to a svyplot with
> style="grayhex". However, the line either appears in the wrong place or if I
> am running in R Studio it causes the system to crash.
>
> I know this is something to do with Lattice graphics, but for the life of me
> I can not figure out how. Dr. Lumley in his excellent book on page 118
> mentions that there is code on his website to do this, but I can not find it.
>
>
>
> So for example:
>
> library(survey)
>
> data(api)
> dclus2<-svydesign(id=~dnum+snum, fpc=~fpc1+fpc2, data=apiclus2)
> svyplot(api00~api99, dclus2)
>
> s1 <-svysmooth(api00~api99, dclus2)
>
> lines(s1)
>
> #works
>
> svyplot(api00~api99, dclus2, style="grayhex")
> lines(s1)
>
> #does not work (line either appears in the wrong position in RGui or crashes
> RStudio).
>
>
> VR
>
> James
>
> James T. Durant, MSPH CIH
> Environmental Health Scientist
> US Agency for Toxic Substances and Disease Registry
> Atlanta, GA 30341
> 770-488-0668
>
>
>
>
> [[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.
--
Thomas Lumley
Professor of Biostatistics
University of Auckland
______________________________________________
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.
