You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
This is the code I ran: library(rworldxtra) library(rworldmap) library(rgdal) cmcpath='/eems/runeems_snps/run1' plotpath='/eems/runeems_snps/run1' coord__long_lat <- read.table(paste0('inputeems', ".coord")) projection_none<-"+proj=longlat +datum=WGS84" projection_mercator<-"+proj=merc +datum=WGS84"`
This is the error I get: Error in h(simpleError(msg, call)) : error in evaluating the argument 'x' in selecting a method for function 'nlayers': cannot do this transformation In addition: Warning message: In rgdal::rawTransform(projfrom, projto, nrow(xy), xy[, 1], xy[, : 701 projected point(s) not finite
How do I get around it? Thank you!
The text was updated successfully, but these errors were encountered:
If I read this correctly, the coordinates are in longitude/latitude but the contour plot of the migration rates is in Mercator. So that might be the issue.
To diagnose the error, I suggest that you start with the basic EEMS plot and then add components one by one. For example, first add the map, then change the projection, then add the original sampling locations. At the very least it will help you figure out what causes the error.
Hello,
I am trying to add a map to my eems plot.
This is the code I ran:
library(rworldxtra)
library(rworldmap)
library(rgdal)
cmcpath='/eems/runeems_snps/run1'
plotpath='/eems/runeems_snps/run1'
coord__long_lat <- read.table(paste0('inputeems', ".coord"))
projection_none<-"+proj=longlat +datum=WGS84"
projection_mercator<-"+proj=merc
+datum=WGS84"`eems.plots(mcmcpath='/Users/nitinr/eems/runeems_snps/run1/', plotpath= '/Users/nitinr/eems/runeems_snps/run1/eems', longlat=TRUE, m.plot.xy = { points(coord__long_lat, col = "purple", pch=18) }, q.plot.xy = { points(coord__long_lat, col = "purple", pch=18) },out.png=FALSE),
projection.in = projection_none, projection.out = projection_mercator,
add.map =TRUE, col.map ="black", lwd.map =5)
This is the error I get:
Error in h(simpleError(msg, call)) : error in evaluating the argument 'x' in selecting a method for function 'nlayers': cannot do this transformation In addition: Warning message: In rgdal::rawTransform(projfrom, projto, nrow(xy), xy[, 1], xy[, : 701 projected point(s) not finite
How do I get around it? Thank you!
The text was updated successfully, but these errors were encountered: