Dear David,
Thank you very much for your replies! I didn't know about
par('usr').
I get different coordinates though:
[1] -19.75966 54.75966 33.60000 71.40000
But the arrow is not at the bottom of the map. I will keep playing with
this. Thanks again!
Sincerely,
Milu
On Mon, Apr 11, 2016 at 12:00 AM, David Winsemius <dwinsemius at
comcast.net>
wrote:
>
>
> > On Apr 10, 2016, at 1:45 PM, Miluji Sb <milujisb at gmail.com>
wrote:
> >
> > Dear David,
> >
> > The device was the issue. The quartz() device works fine but pdf()
does
> not. Now I just need to figure out the limits for map for Europe. Thanks
> for all your help and patience.
>
>
> After plotting a map of Europe with base graphics the coordinates of the
> lower-left and upper-right corners are obtained by par('usr')
>
> > par('usr')
> [1] -12.20000 47.20000 25.89375 79.10625
>
> --
> David.
>
>
> >
> > Sincerely,
> >
> > Milu
> >
> > On Sun, Apr 10, 2016 at 7:10 PM, David Winsemius <dwinsemius at
comcast.net>
> wrote:
> >
> > > On Apr 10, 2016, at 4:12 AM, Miluji Sb <milujisb at
gmail.com> wrote:
> > >
> > > Hello David,
> > >
> > > This is exactly what I want but I still can't get the arrows.
R and R
> studio is updated. Thanks again!
> >
> > I didn't try it in Rstudio until just now (and I don't
remember that you
> ever mentioned RStudio as a possible issue.) The plotting I see in the
> default graphics Rstudio window is rather different than what I see in the
> default plotting window for the R.app GUI. I'm guessing you are
reporting
> results from viewing plots in that IDE's viewing window.
> >
> > I would try plotting with the png() or pdf() devices and see if the
> results are more predictable. I just tried with pdf() from RStudio and the
> results were much closer to what I was seeing with saving from R.app. The
> `quartz()` device seems to deliver consistent results for me. The RStudio
> device is something they call RStudioGD, and I don't have sufficient
> experience to explain its quirks.
> >
> > --
> > David.
> >
> > >
> > > Sincerely,
> > >
> > > Milu
> > >
> > > On Sat, Apr 9, 2016 at 10:29 PM, David Winsemius <
> dwinsemius at comcast.net> wrote:
> > >
> > > > On Apr 9, 2016, at 1:27 PM, David Winsemius <dwinsemius
at comcast.net>
> wrote:
> > > >
> > > >
> > > >> On Apr 9, 2016, at 11:18 AM, David Winsemius <
> dwinsemius at comcast.net> wrote:
> > > >>
> > > >>
> > > >>> On Apr 9, 2016, at 10:46 AM, Miluji Sb <milujisb
at gmail.com> wrote:
> > > >>>
> > > >>> Dear David,
> > > >>>
> > > >>> Thank you for your answer. Sorry for the
embarrassing mistake.
> However, even with when I generate a map for the whole world using:
> > > >>>
> > > >>> eps <- mapCountryData(n,
nameColumnToPlot="eps_score",
> mapTitle="EPS Score",colourPalette=colourPalette,
> > > >>>
catMethod="fixedWidth", missingCountryCol > "white",
addLegend=FALSE)
> > > >>>
> > > >>> And then use:
> > > >>>
> > > >>> do.call(addMapLegend, c(eps,
legendLabels="none",
> arrows(-100,-140,100,-140,code=3)))
> > > >>
> > > >> I do get an arrow using same version of R and OSX. See
attached. (I
> think that png images will be accepted by the mailserver.)
> > > >
> > > > Nope I was wrong, but the copy to Milugi did arrive.
> > > >
> > > > Here's a (somewhat larger) pdf:
> > >
> > >
> > >
> > >
> > > >
> > > >
> > > > --
> > > > David
> > > >>
> > > >>
> > > >>
> > > >> --
> > > >> David.
> > > >>>
> > > >>> Only a legend with the colours is generated, no
arrows. My session
> info is below. Thanks again!
> > > >>>
> > > >>> R version 3.2.4 (2016-03-10)
> > > >>> Platform: x86_64-apple-darwin13.4.0 (64-bit)
> > > >>> Running under: OS X 10.11.2 (El Capitan)
> > > >>>
> > > >>> locale:
> > > >>> [1]
en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
> > > >>>
> > > >>> attached base packages:
> > > >>> [1] stats graphics grDevices utils datasets
methods
> base
> > > >>>
> > > >>> other attached packages:
> > > >>> [1] countrycode_0.18 ggplot2_2.1.0
RColorBrewer_1.1-2
> foreign_0.8-66 maptools_0.8-39 rworldmap_1.3-6
> > > >>> [7] sp_1.2-0
> > > >>>
> > > >>> loaded via a namespace (and not attached):
> > > >>> [1] Rcpp_0.12.4 lattice_0.20-33 grid_3.2.4
plyr_1.8.3
> gtable_0.2.0 spam_1.3-0
> > > >>> [7] scales_0.4.0 tools_3.2.4 munsell_0.4.3
maps_3.1.0
> fields_8.3-6 colorspace_1.2-6
> > > >>>
> > > >>> On Sat, Apr 9, 2016 at 7:34 PM, David Winsemius <
> dwinsemius at comcast.net> wrote:
> > > >>>
> > > >>>> On Apr 9, 2016, at 8:13 AM, Miluji Sb
<milujisb at gmail.com> wrote:
> > > >>>>
> > > >>>> Forgot to copy the list
> > > >>>>
> > > >>>> Dear Jim,
> > > >>>>
> > > >>>> Thank you for your reply. I must be doing
something wrong, If
> this is my
> > > >>>> command to plot a map of Europe:
> > > >>>>
> > > >>>> eps_europe <- mapCountryData(n,
nameColumnToPlot="eps_score",
> mapTitle="EPS
> > > >>>> Score -
Europe",colourPalette=colourPalette,
> > > >>>>
catMethod="fixedWidth",
> missingCountryCol > > > >>>> "white",
mapRegion="Europe", addLegend=FALSE)
> > > >>>>
> > > >>>> The following command does not seem to add the
arrow. What am I
> doing wrong?
> > > >>>>
> > > >>>> do.call(addMapLegend, c(eps_europe,
legendLabels="none",
> > > >>>> arrows(-100,-140,100,-140,code=3)))
> > > >>>>
> > > >>>
> > > >>> Your earlier question had a full world map. That was
the context
> for Jim's reply, which did plot a two headed arrow above the legend in
your
> earlier question. Now you have restricted the plot region to Europe so the
> coordinates of -100,-140,100,-140 no longer are on the visible plot area.
> You need to decide where you want the arrows using sensible coordinates.
> > > >>>
> > > >>> --
> > > >>> David.
> > > >>>
> > > >>>
> > > >>>> Thank you again. I really appreciate it.
> > > >>>>
> > > >>>> Sincerely,
> > > >>>>
> > > >>>> Milu
> > > >>>>
> > > >>>> On Sat, Apr 9, 2016 at 12:20 PM, Jim Lemon
<drjimlemon at gmail.com>
> wrote:
> > > >>>>
> > > >>>>> Hi Miluji,
> > > >>>>> Try this:
> > > >>>>>
> > > >>>>> arrows(-100,-140,100,-140,code=3)
> > > >>>>>
> > > >>>>> Jim
> > > >>>>>
> > > >>>>>
> > > >>>>> On Fri, Apr 8, 2016 at 10:24 PM, Miluji Sb
<milujisb at gmail.com>
> wrote:
> > > >>>>>> I am trying to draw maps for the world
using:
> > > >>>>>>
> > > >>>>>> library(rworldmap)
> > > >>>>>> library(maptools)
> > > >>>>>> library(RColorBrewer)
> > > >>>>>>
> > > >>>>>>
> > > >>>>>> tmp2<- dput(head(pece,10))
> > > >>>>>> structure(list(iso3 = c("AUS",
"AUT", "BEL", "CAN", "CHE",
> "CHL",
> > > >>>>>> "CZE", "DEU",
"DNK", "ESP"), eps_score = c(0.877343773841858,
> > > >>>>>> 2.68984365463257, 1.31406247615814,
1.98046875,
> 2.61666655540466,
> > > >>>>>> NA, 1.44414067268372, 2.34257817268372,
2.89687490463257,
> > > >>>>> 2.15937495231628
> > > >>>>>> ), gov_eff = c(1.76499999562899,
1.85666667421659,
> 1.74500000476837,
> > > >>>>>> 1.88416666785876, 1.99181815710935,
1.21499997377396,
> 0.865833342075348,
> > > >>>>>> 1.64999999602636, 2.15416664878527,
1.36833332975705),
> sh_va_enint > > > >>>>>> c(13.4375638961792,
> > > >>>>>> 8.90904521942139, 10.368335723877,
14.0469560623169, NA, NA,
> > > >>>>>> 13.5679216384888, 9.67090892791748,
10.5978908538818,
> 8.34146690368652
> > > >>>>>> ), rd_in_va = c(2.17547988891602,
2.47147130966187,
> 2.53955459594727,
> > > >>>>>> 2.01138758659363, NA, NA,
1.49587619304657, 2.72330951690674,
> > > >>>>>> 2.5316367149353, 1.48551619052887)),
datalabel = "", time.stamp
> = " 9 Mar
> > > >>>>>> 2016 17:43", .Names =
c("iso3",
> > > >>>>>> "eps_score",
"gov_eff", "sh_va_enint", "rd_in_va"), formats
> c("%9s",
> > > >>>>>> "%8.0g", "%10.0g",
"%9.0g", "%9.0g"), types = c(6L, 254L, 255L,
> > > >>>>>> 254L, 254L), val.labels =
c("", "", "", "", ""),
var.labels > c("",
> > > >>>>>> "(mean) eps_score",
"(mean) gov_eff", "(mean) sh_va_enint",
> "(mean)
> > > >>>>>> rd_in_va"
> > > >>>>>> ), expansion.fields =
list(c("_dta", "ReS_i", "countrycode"),
> > > >>>>>> c("_dta",
"ReS_ver", "v.2"), c("_dta", "ReS_j",
"year"),
> > > >>>>>> c("_dta",
"ReS_str", "0"), c("_dta", "ReS_Xij",
"a_"),
> c("_dta",
> > > >>>>>> "__JVarLab",
"ACT"), c("_dta", "__XijVarLabrdd_", "(sum)
rdd"
> > > >>>>>> ), c("_dta",
"__XijVarLabp", "Value"), c("_dta",
> "__XijVarLabpop",
> > > >>>>>> "Population"),
c("_dta", "__XijVarLabest_lu_f", "Source of lu"
> > > >>>>>> ), c("_dta",
"__XijVarLablu", "Percentage of No Schooling"
> > > >>>>>> ), c("_dta",
"__XijVarLabest_lp_f", "Source of lp"), c("_dta",
> > > >>>>>> "__XijVarLablp",
"Percentage of Primary"), c("_dta",
> "__XijVarLablh",
> > > >>>>>> "Percentage of Tertiary"),
c("_dta", "__XijVarLabest_lh_f",
> > > >>>>>> "Source of lh"),
c("_dta", "__XijVarLabls", "Percentage of
> Secondary"
> > > >>>>>> ), c("_dta",
"__XijVarLabest_ls_f", "Source of ls"), c("_dta",
> > > >>>>>> "__XijVarLabvalue",
"Value"), c("_dta", "_TStvar", "year"
> > > >>>>>> ), c("_dta",
"_TSpanel", "id2"), c("_dta",
"_TSdelta",
> > > >>>>>> "+1.0000000000000X+000"
> > > >>>>>> ), c("_dta",
"_TSitrvl", "1"), c("_dta", "tis",
"year"),
> > > >>>>>> c("_dta", "iis",
"id2")), version = 12L, row.names = c("1",
> > > >>>>>> "2", "3",
"4", "5", "6", "7", "8",
"9", "10"), class > "data.frame")
> > > >>>>>> n <- joinCountryData2Map(pece,
joinCode="ISO3",
> nameJoinColumn="iso3")
> > > >>>>>> n <-
n[-which(row.names(n)=='Antarctica'),]
> > > >>>>>>
> > > >>>>>> # EPS
> > > >>>>>> colourPalette <- rev(brewer.pal(7,
"RdYlGn"))
> > > >>>>>>
> > > >>>>>> eps <- mapCountryData(n,
nameColumnToPlot="eps_score",
> mapTitle="EPS
> > > >>>>>> Score",colourPalette=colourPalette,
> > > >>>>>>
catMethod="fixedWidth", missingCountryCol > > >
>>>>> "white",
> > > >>>>>> addLegend=FALSE)
> > > >>>>>> do.call(addMapLegend, c(eps,
legendLabels="all",
> legendWidth=0.5))
> > > >>>>>>
> > > >>>>>> Instead of adding numeric based legend,
I would like to add a
> two-headed
> > > >>>>>> arrow with some text. I would be
grateful for any help. Thank
> you!
> > > >>>>>>
> > > >>>>>> Sincerely,
> > > >>>>>>
> > > >>>>>> Milu
> > > >>>>>>
> > >
> > >
> > >
> >
> > David Winsemius
> > Alameda, CA, USA
> >
> >
>
> David Winsemius
> Alameda, CA, USA
>
>
[[alternative HTML version deleted]]