similar to: Overlay line on a bar plot - multiple axis

Dear all, I am trying to make a similar plot - https://peltiertech.com/images/2013-09/BarLineSampleChart4.png. I have data for two variables; count and z by city and week. I would like to have a horizontal bar plot of *count* by city and a line plot of weekly average of the variable *z*. I have tried the following: ggplot() + geom_bar(data=dat, aes(x=city, y=count),

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. 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: > >

> 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

Hello David, This is exactly what I want but I still can't get the arrows. R and R studio is updated. Thanks again! 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,

> 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

Hi Milu, My fault here. As I don't have the data to make the map and try out my suggestions I mixed up the x and y coordinates. Try this: par(xpd=TRUE) arrows(-19.75966,53,33.60000,53,code=3) par(xpd=FALSE) Jim On Tue, Apr 12, 2016 at 10:11 PM, Miluji Sb <milujisb at gmail.com> wrote: > Hello Jim, > > Thanks again. I am getting the two-headed arrow but I cannot seem to get

> 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 <-

Hello Jim, Thanks again. I am getting the two-headed arrow but I cannot seem to get the coordinates right for the arrow to appear beneath the map. These coordinates puts the arrow on the left hand side. Thanks again! Sincerely, Milu On Tue, Apr 12, 2016 at 1:15 PM, Jim Lemon <drjimlemon at gmail.com> wrote: > Hi Milu, > There is a two-headed arrow on the image you sent, and it

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:

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",

> 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.

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 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, >

> 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 -

xtabs does this automatically if your cross classifying variables are factors with levels all the cities (sorted, if you like): > x <- sample(letters[1:5],8, rep=TRUE) > y <- sample(letters[1:5],8,rep=TRUE) > xtabs(~ x + y) y x c d e a 1 0 0 b 0 0 1 c 1 0 0 d 1 1 1 e 1 1 0 > lvls <- sort(union(x,y)) > x <- factor(x, levels = lvls) > y <-

or in base R : ?xtabs ?? as in: xtabs(~previous_location + current_location,data=x) (You can convert the 0s to NA's if you like) Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, May 8, 2018 at 9:21 AM, Huzefa

Dear Bert and Huzefa, Apologies for the late reply, my account got hacked and I have just managed to recover it. Thank you very much for your replies and the solutions. Both work well. I was wondering if there was any way to ensure (force) that all possible combinations show up in the output. The full dataset has 25 cities but of course people have not moved from Boston to all the other 24

Ista, et. al: efficiency? (Note: I needed to correct my previous post: do.call() is required for pmax() over the data frame) > x <- data.frame(matrix(runif(12e6), ncol=12)) > system.time(r1 <- do.call(pmax,x)) user system elapsed 0.049 0.000 0.049 > identical(r1,r2) [1] FALSE > system.time(r2 <- apply(x,1,max)) user system elapsed 2.162 0.045 2.207 ##

Dear William and Ben, Thank you for your replies and elegant solutions. I am having trouble with the fact that two of the previous locations do not appear in current locations (that is no one moved to OKC and Dallas from other cities), so these two cities are not being included in the output. I have provided a better sample of the data and the ideal output (wide form - a 10x10 bilateral matrix)

Hi Milu, There is a two-headed arrow on the image you sent, and it seems to be where you specified. Did you want it beneath the map, as: par(xpd=TRUE) arrows(-22,54.75,-22,74,code=3) par(xpd=FALSE) Jim On Tue, Apr 12, 2016 at 7:58 PM, Miluji Sb <milujisb at gmail.com> wrote: > Dear Jim, > > Thanks again! I do want the arrows at the bottom (beneath the map). This is > what I am