similar to: Plotting from different data sources on the same plot (with ggplot2)

Hi everyone, Hi Hadley, I am a heavy user of coord_equal() in ggplot2 since most of my data is spatial, on x,y coordinates. Everything works. However by enforcing an aspect ratio of 1 for the plotting region, coord_equal() usually wastes a lot of space if the region of interest is not a perfect square. For example: x=runif(10) a=data.frame(x=x*3,y=x) ggplot(data=a, aes(x=x,y=y)) +

Buenos días. Hago mis mapas con el código que os copio abajo, pero me gustaría que me dibujase las líneas del mapa después de poner los puntos para que se vean. Puedo hacer los puntos más pequeños, pero entonces quedan separados y yo quiero que rellenen el mapa. Gracias, como siempre, Manuel print(ggplot(legend=FALSE)+geom_path(data=map_data('world'), aes(x=long,

Buenos días. Tengo dos mapas que hago con ggplot (abajo) y necesito superponer el 2º al 1º. En el 1º se representa la variable PA, binomial (0, 1), en gris, y en el 2º la variable ExtTG7085, que es numérica y va de 0 a 200. La 2ª es un gradiente en rojo. Lo ideal sería que donde coinciden el rojo de la 2ª se viera más oscuro, pero supongo que eso no es fácil. Podría también ser en otro color, pero

I'm sure this is smack-head moment, but I haven't been able to find an example of this on Nabble or SO, so thought I'd ask. This works: michigan <- map_data('county', 'michigan') mich_points <- data.frame(x = rnorm(n = 200, median(michigan[,1]), 0.75), y = rnorm(n = 200, median(michigan[,2]), 0.75)) ggplot() + geom_path(aes(long, lat, group = group), data =

Gracias a los tres, Raúl, Marcelino y Carlos. Lo del "get" de Marcelino me da la respuesta a lo que yo exactamente preguntaba, y funciona, pero ahora tengo problemas con el for, por lo que probablemente recurra al eval parse de Raúl o Carlos, que ya tienen el for. Aún así, lo intento 1º con el get. Con subset(df, subset=get(GT[i])>0) el problema es que en el for hago un

Gracias Marcelino. Si, (i in 1:length(GT)), lo he utilizado mil veces, pero se me sigue olvidando de una vez a otra. Lo iba a mirar, pero me centré primero en que me hiciera bien el mapa. He probado el for y me da este error: Error in aes(x = lon, y = lat, color = get(GT[i]), size = 2) + scale_colour_gradient(low = ("white"), : non-numeric argument to binary operator

Dear all I want to display 4 dimensional space by some suitable way. I searched CRAN and found miscellaneous 3 dim graphics packages which I maybe can modify but anyway I am open to any hint how to efficiently display data like: longitude, latitude, height, value Thank you Petr Pikal petr.pikal at precheza.cz

Muy buenas. Quiero hacer un mapa que me pinte en azul los valores negativos y en rojo los positivos. Los negativos llegan hasta -400 y los positivos hasta 200. Si pongo limits=c(-200, 200), me colorea bien los positivos (rojo), pero los negativos por debajo de -200 me los pone gris. Si pongo limits=c(-400, 200), me colorea bien los negativos (azul), pero los positivos apenas se ven. Necesito que

Hi, QUESTION TOPIC #1 I have some data I want to plot on a map. But what I have are home addresses: street, City, State, complete postal code--i.e 95377-1234. Is there a way to plot this data or do I need latitudinal and longitude coordinates? If so how do I convert them? Is there a package that will do the conversion in R? QUESTION TOPIC #2 I was trying to experiment with this code that I found

Hi R-help I am trying to mask the ocean from an image plot I have made. Here is some example code: library(mapdata) image(x=110:155, y =-40:-10, z = outer(1:45, 1:30, "+"), xlab = "lon", ylab = "lat") outline <- map("worldHires", plot=FALSE) # returns a list of x/y coords xrange <- range(outline$x, na.rm=TRUE) # get bounding box yrange

Dear R users, I need to produce rainfall maps using R. I know that this is possible, I looked though the web, I found the example below reported (the author is Andrew Tredennick). I would ask you if this is the most performing way to make rainfall maps; if yes would someone be able to give me an example of how file.asc and pointfile.csv should be? If no would somebody please show me another way

json_dir is a list of JSON lists mapping lat/long route points between locations using CloudMade's API. post_url is the URL of the HTTP request for (n in json_dir) { i = i + 1 if (typeof(json_dir[[i]]) != "NULL") { if (i == 1) { dat_add <- ldply(json_dir[[i]], function(x) t(data.frame(x)), .progress = "text")

Manuel: I changed your variable names from x to 'long' and y to 'lat' on the riqueza_out.csv file. The code below should do what you want. Also, since the legend title is kind of long, I broke it down into three lines so you can see more plot area. I am cc'ing the other groups so more people use it if needed. library(rgdal) library(ggplot2) library(sp) library(maptools)

To all, I'm working with code below to produce a map with station data plotted in points, but right now I'm having trouble with the mapping portion of this code states <- data.frame(map("state", plot=FALSE,xlim= c(-85,-75),ylim=c(33,37))[c("x","y")]) usamap<- ggplot(states)+geom_path(aes(x,y)) usamap When I plot this the problem is that the bounds of

Gracias Emilio. Si, son asimétricos porque se mueven en rangos distintos. De hecho, miden cosas distintas, pero una muestra no puede ser positiva y negativa a la vez. Me planteo multiplicar las positivas por 2 y al editar el mapa corregirlo para que indique los valores reales. Eso funcionaría, pero son muchos mapas y quizás haya una función en ggplot que evite el problema directamente. Manuel El

Hello everyone, I am (happily) using ggplot2 for all my plotting now and I wondered is there is an easy way to specify xlim and ylim somewhere when using the ggplot syntax, as opposed to the qplot syntax. Eg. qplot(data=mtcars,y=wt, x=qsec,xlim=c(0,30)) <-> ggplot(mtcars, aes(y=wt, x=qsec)) + geom_point() + ??? Indeed the ggplot syntax is in general more flexible and powerful and

Hi everyone! I'm digging into ggplot for some while now, I must say it's great! But - as some others have posted before, and hadley knows very well himself - the documentation is lacking some bits... So I have to pose that question directly: I'd like to produce a series of maps with different data on, but exactly the same extent in each plot. Is there a way of switching the

No sé muy bien, Carlos, como hacerlo, pero no parece más sencillo que lo de multiplicar por 2 los positivos y después editar el mapa. De lo que decías, Emilio, de cambiar el color del punto central creo que haría más confusa la interpretación. El vie., 20 sept. 2019 a las 10:33, Carlos Ortega (<cof en qualityexcellence.es>) escribió: > Hola, > > Una alternativa que suele

Muy buenas. Quiero hacer un loop en el que en cada iteración se hace un subset con el que se queda con las muestras para la que cierta variable es positiva. Si hago esto, sale bien: df2<-subset(df, subset = var1>0) Pero he probado así (y de no sé cuantas formas más), antes de hacer el for, y no sale: GT<- c("var1","var2", ? ) df2<-subset(df,

Dear R users, #---------------------------------------------------------------------------------------------------------------------------- #There is table containing 1000 (lat, lon, time) GPS positions, wchich should be recognized and labeled in every row of that #recognized position by label from defined POI vector GPS$Lat =c(9,2.2,2,3,4,5,6,7,2,9,2.1,10,3.3,0 , 6, 8, 9) GPS$Lon =c(0,1