Thank you, Ray! That worked.
Thank you also, Yihui. I tried your subplot() suggestion, but couldn't
get it to work with a map() on a map()
Jean
Ray Brownrigg <ray.brownr...@ecs.vuw.ac.nz> wrote on 09/27/2011 07:26:34
PM:
>
> On Wed, 28 Sep 2011, Jean V Adams wrote:
> > I want to overlay a small inset map on top of another map, but I can't
> > figure out how to do it.
> > For example, here are two different maps:
> >
> > # map 1 - Ohio
> > map("state", region= "ohio")
> >
> > # map 2 - US with Ohio darkened
> > map("state")
> > map("state", region="ohio", fill=T, add=T)
> >
> > I would like to add map 2 as a small inset in the corner of map 1.
> > I have tried:
> >
> > map("state", region= "ohio")
> > par(new=TRUE, mar=c(3, 3, 15, 15))
> > map("state")
> > map("state", region="ohio", fill=T, add=T)
> >
> > but this seems to erase map 1 and replace it with a full size version
of
> > map 2.
> >
> > I can successfully overlay an unrelated plot using similar code:
> >
> > map("state", region= "ohio")
> > par(new=TRUE, mar=c(3, 3, 15, 15))
> > plot(1:10, 1:10)
> >
> > So, there must be something about the maps() function that I'm
tripping
> > over.
>
> I think the problem is that map() 'insists' (in some sense) on a
> clean frame so it can get
> the aspect ratio right.
>
> >
> > Any suggestions?
> >
> How about something like:
> map("state", region= "ohio", xlim=c(-85, -80), ylim=c(38, 42))
> par(usr=c(-216, -66, 24, 144)) # you should be able to
> 'automate' this calculation
> map("state", add=T)
> map("state", region="ohio", fill=T, add=T)
>
> HTH
> Ray Brownrigg
>
> > I am using R for Windows 2.13.0
> > and the maps package version 2.1-5.
> >
> > Jean
> >
> >
> > `�.,, ><(((�> `�.,, ><(((�> `�.,, ><(((�>
> >
> > Jean V. Adams
> > Statistician
> > U.S. Geological Survey
> > Great Lakes Science Center
> > 223 East Steinfest Road
> > Antigo, WI 54409 USA
>
>
[[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.
