similar to: political maps world maps in R, wrld_simpl

Hi useR's, I want to overlay an image plot over a world map and I can do it, but just not the way I need to do it. Here is the code I am using (with data file attached) to create my baseline map: library(sp) load("TM_WORLD_BORDERS_SIMPL-0.2.RData") par(bty="l") plot(wrld_simpl, axes = TRUE, ylim = c(-90, 90), xlim=c(-180, 180), asp=1.5) lim <- par("usr")

Dear All, When I tried to run the following code (taken from the *spdplyr* package vignettes): library(spdplyr) library(maptools) data(wrld_simpl) worldcorner <- wrld_simpl %>% mutate(lon = coordinates(wrld_simpl)[,1], lat = coordinates(wrld_simpl)[,2]) %>% filter(lat < -20, lon > 60) %>% dplyr::select(NAME) I got the following error messages: Error in (function (cl,

Muchas gracias, Carlos. Estaba emperrado con usar OSM. Pero sevuramente los de Google presentaran alguna solucion. Intento probar pronto. Juan Abasolo. telefomotik, barjatu trebezi falta Carlos Ortega <cof en qualityexcellence.es> igorleak hau idatzi zuen (2019 mar. 19, ar. 23:13): > Hola, > > No sé si esto seguirá siendo vigente, con los cambios en las APIs de > Google...

Buenas! Estoy intentando hacer un trabajo con mapas y llegué a un punto en el que me atasco: Quiero generar el dibujo de un mapa. A ese mapa después le señalaría unos puntos y en, base a esos puntos, dividiría el mapa en secciones (usando diagramas Voronoi, que creo que sabría hacerlo, por más que todavía me quedan raro) Todavía lo que sé hacer es muy ínfimo (marcar puntos y etiquetar con

Hola a todos. Soy nuevo en R, estoy aprendiendo con tutoriales y lo hago sin un fin comercial y/o laboral. Necesito ayuda para mostrar un mapa y sobre ese mapa poner un punto dado por una coordenada. con esto consigo mostrar el mapa. library(raster) arg<-getData('GADM', country='ARG', level=0) plot(arg) En ese mapa quiero mostrar un punto que está en el sistema de

Good day, For lack of a better solution (or perhaps I am ignorant to something more elegant), I have been bootstrapping panel data by hand so to speak and I would like to know if there is a way to define multiple variables in a loop using the loop variable. I found a post (here: https://stat.ethz.ch/pipermail/r-help/2002-October/026305.html ) that discussed naming multiple variables but it

I'm having difficulty overlaying maps when writing to a file graphics device. My command sequence has the structure plot(map1) par(new = T) plot(map2) On the screen device, it works fine. When I attempt something like png(file = "map.png") plot(map1) par(new = T) plot(map2) dev.off() only the last map appears, the previous ones having been cleared. Can someone clarify? Thanks,

Allows you to safely use nbdkit-offset-filter on top of a plugin supporting extents. --- filters/offset/offset.c | 43 +++++++++++++++++++++++++++++++++++++++++ 1 file changed, 43 insertions(+) diff --git a/filters/offset/offset.c b/filters/offset/offset.c index 058571d..4e3f74d 100644 --- a/filters/offset/offset.c +++ b/filters/offset/offset.c @@ -37,6 +37,7 @@ #include <stdlib.h>

On 3/20/19 10:31 AM, Richard W.M. Jones wrote: > I think the extents map is just too complicated and is unnecessarily > so. How about instead we define the plugin interface to be: > > int can_extents (void *handle); // as before > int extents (void *handle, uint32_t count, uint64_t offset, uint32_t flags, > struct nbdkit_extents_list *list); > > and

Dear all, I need to compute tensor product of B-spline defined over equi-spaced break-points. I wrote my own program (it works in a 2-dimensional setting) library(splines) # set the break-points Knots = seq(-1,1,length=10) # number of splines M = (length(Knots)-4)^2 # short cut to splineDesign function bspline = function(x) splineDesign(Knots,x,outer.ok = T) # bivariate tensor product of

Hi all, I am running a bivariate regression with the following: p1=c(184,155,676,67,922,22,76,24,39) p2=c(1845,1483,2287,367,1693,488,435,1782,745) I1=c(1530,1505,2505,204,2285,269,1271,298,2023) I2=c(8238,6247,6150,2748,4361,5549,2657,3533,5415) R1=I1-p1 R2=I2-p2 x1=cbind(p1,R1) y1=cbind(p2,R2) fit1=lm(y1~-1+x1) summary(fit1) Response 2: Coefficients: Estimate Std. Error t value

Dear R-users, I learned today that there exists an interesting topic in numerical analysis names "best polynomial approximation" (BSA). Given a function f the BSA of degree k, say pk, is the polynomial such that pk=arginf sup(|f-pk|) Although given some regularity condition of f, pk is unique, pk IS NOT calculated with least square. A quick google tour show a rich field of research

Dear all, is there any way to transform a .Data directory created in S-plus 6.1 for windows in a .RData file? Thanks in advance, Patrizio Frederic

On Wed, Mar 20, 2019 at 11:20:15AM -0500, Eric Blake wrote: > On 3/20/19 10:31 AM, Richard W.M. Jones wrote: > > I think the extents map is just too complicated and is unnecessarily > > so. How about instead we define the plugin interface to be: > > > > int can_extents (void *handle); // as before > > int extents (void *handle, uint32_t count, uint64_t

hi R-users: I have some data related with geographical information. I wonder how to plot some data (in contour,shaded or scatter) on a map. Could anyone tell me which package or reference or examples I can learn ? thank you . -- TANG Jie Email: totangjie@gmail.com [[alternative HTML version deleted]]

Dear R-helpers, I'm having a problem in using plot.design in Design Library. Tho following example code produce the error: > n <- 1000 # define sample size > set.seed(17) # so can reproduce the results > age <- rnorm(n, 50, 10) > blood.pressure <- rnorm(n, 120, 15) > cholesterol <- rnorm(n, 200, 25) > sex <-

summary: spatial data to be input to a regional-scale environmental model must (1) be converted to netCDF and then (2) "regridded" (cropped, projected, increased resolution). In a public git repository https://github.com/TomRoche/GEIA_to_NetCDF I have R code (with bash drivers) that does the conversion step, and plots the converted output, apparently correctly. However attempts to plot

library(mvtnorm) x = seq(-4,4,length=201) xy = expand.grid(x,x) sigma = (diag(c(1,1))+1)/2 d2 = matrix(dmvnorm(xy,sigma=sigma),201) xsamp = rmvnorm(200,sigma=sigma) contour(x,x,d2) points(xsamp,col=3,pch=16) pdf("pdftry.pdf") contour(x,x,d2) points(xsamp,col=3,pch=16) dev.off() postscript("pstry.ps") contour(x,x,d2) points(xsamp,col=3,pch=16) dev.off() # I can see

hi all, I have a data frame such as: 1 blue 0.3 1 NA 0.4 1 red NA 2 blue NA 2 green NA 2 blue NA 3 red 0.5 3 blue NA 3 NA 1.1 I wish to find the last non-missing value in every 3ple: ie I want a 3 by 3 data.frame such as: 1 red 0.4 2 blue NA 3 blue 1.1 I have written a little script data = structure(list(V1 = c(1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L ), V2 = structure(c(1L, NA,

Hi, I'm trying to build some maps of the US by county that will have the following characteristics: Feature/Map Map 1 Map2 Both Broken out by county Yes Yes Yes Heatmaps of US Census Data for income by county Yes No Yes Heatmaps of US Census Data for race by county (recoded as white and &non-white, with each county color coded based on the majority) No Yes No Polygon