search for: geosphere

Hey, I'm trying to build something like this http://flowingdata.com/2011/05/11/how-to-map-connections-with-great-circles/ but with my own data in csv files. The code runs well if I use the same csv files as the author, but with mine , this is what I get *Code* library(maps) library(geosphere) map("world") xlim <- c(-180.00, 180.00) ylim <- c(-90.00, 90.00) map("world", col = "#f2f2f2", fill = TRUE, bg = "white", lwd = 0.05,xlim = xlim, ylim = ylim) airports <- read.csv("/Users/shabnam/Desktop/airports.csv", as.is=TRUE, h...

...ot;euclidean" in (for example) corSpher would be incorrect. I would be grateful for help on how to write a new distance metric for the corSpatial function. I believe there are several ways that distances on a sphere can be calculated in R, for example the "distMeeus" function in the geosphere library. However, I have no idea how to write this into a corSpatial function. The aim is to end up with a metric = "sphere" option that calculates great circle distances between points using latitude and longitude. Many thanks, Dan

...ta on a map. Everything is working well except the size of the map is really too small to show the data effectively and I can't seem to figure out how to make the output map larger. Do I need to change the device characteristics or is it a map.??? call. Here is the code: library(maps) library(geosphere) airports <- read.csv("airports.csv", header=TRUE) flights <- read.csv("FSRflights.csv", header=TRUE, as.is=TRUE) checkDateLine <- function(l){ n<-0 k<-length(l) k<-k-1 for (j in 1:k){ n[j] <- l[j+1] - l[j] } n <- abs(n) m<-max(n, r...

Hi all, I am a very recent user of R. Mine is probably a basic question. I would like to know how do we create a certain set of points with geographic coordinates, between two endpoints. Please explain me how to do this using R. Which are the packages to be used for doing this. Thanks in advance, Ravichandra Mondreti. Ravichandra Mondreti PhD student/ Doctorant Spatial

Dear R-listers, Here, I would like to hearing helps from you. I have GPS data (multiple points in the geographic scale) in longitude/latitude. I intend to calculate distance (in kilometer) among such points and output the distance matrix in dist class. I have gotten some progress, but I still can not get final goal. Could please give me any directions/advice? This email cc. to Mr. Pierre, the

I'm in the process of converting some S3 methods to S4 methods. I have this function : setGeneric("enrichmentCalc", function(rs, organism, seqLen, ...){standardGeneric("enrichmentCalc")}) setMethod("enrichmentCalc", c("GenomeDataList", "BSgenome"), function(rs, organism, seqLen, ...) { ... ... ... })

...(2.9), catspec (0.96), cda (1.1.3), clue (0.3-44), coin (1.0-21), crs (0.15-13), cvTools (0.3.0), Deducer (0.6-2), directlabels (2.4), doRNG (1.2.1), dplR (1.5.3), DSL (0.1-1), eqtl (1.1-7), ergm (3.0-1), ergm.userterms (3.0-1), fda.usc (0.9.6), frailtypack (2.2-22), games (1.0-5), gdsfmt (0.9.6), geosphere (1.2-27), GeoXp (1.5.6), GLDEX (1.0.4.3), gMCP (0.7-8), GSM (1.3), hive (0.1-14), HLMdiag (0.1.5), HumMeth27QCReport (1.2.13), ibr (1.4.1), LaplacesDemon (12.03.05), lattice (0.20-5), lattice (0.20-6), lifecontingencies (0.9.1), lifecontingencies (0.9), lubridate (1.1.0), Mangrove (1.1), matrixStat...

...8), ftsa (1.9), FunNet (1.00-9), gamlss (4.0-0), gamlss.add (4.0-0), gamlss.cens (4.0.0), gamlss.data (4.0-0), gamlss.dist (4.0-0), gamlss.mx (4.0-0), gamlss.nl (4.0-0), gamlss.tr (4.0-0), gamlss.util (4.0-0), gap (1.0-23), gbm (1.6-3.1), gclus (1.3), gdata (2.8.0), geoR (1.6-28), geoRglm (0.8-28), geosphere (1.2-4), geozoo (0.4.1), GExMap (1.1), ggplot2 (0.8.8), ghyp (1.5.3), glmnet (1.4), gmm (1.3-2), gplots (2.8.0), gputools (0.21), gRain (0.8.5), grImport (0.6-1), gss (1.1-5), gsubfn (0.5-3), gtools (2.6.2), gWidgets (0.0-41), gWidgetsRGtk2 (0.0-67), gWidgetsrJava (0.0-19), gWidgetstcltk (0.0-38),...

Perdon por no se lo suficientemente claro :( Tu codigo produce `validPairs` que tiene 7 variables y 360 observaciones. Donde > validPairs[1,] V1 V2 V3 V4 V5 V6 valid 60 e1 g1 c1 e1 g1 c2 Valid indica que un maestro tiene asignado c1 y c2 en la escuela e1 y el grado g1. Correcto? Si es asi, esto es casi lo que queira producir y creo que puedo llegar a donde quiero usando tu codigo de base.