similar to: package (in development)

There is an update for the `foreign' package now on the CRAN master and soon on a mirror near you. This adds support for DBF files (by Nicholas Lewin-Koh, Roger Bivand and myself) Systat files written on MS-DOS/Windows (by Roger Bivand) The reasons for announcing this update are 1) People quite often ask for DBF support. 2) There is a `maptools' update that depends on getting this

There is an update for the `foreign' package now on the CRAN master and soon on a mirror near you. This adds support for DBF files (by Nicholas Lewin-Koh, Roger Bivand and myself) Systat files written on MS-DOS/Windows (by Roger Bivand) The reasons for announcing this update are 1) People quite often ask for DBF support. 2) There is a `maptools' update that depends on getting this

I want to read data from a number of files into R. Reading individual files one by one requires writing enormous amount of code that will look something like the following. **************** maptools:::dbf.read("wb-01vc.dbf")->dist1 maptools:::dbf.read("wb-02vc.dbf")->dist2 maptools:::dbf.read("wb-03vc.dbf")->dist3

The main limitation of the shapefiles package that I put together is that it does not create shapefiles from R objects - rather it only writes shapefiles that have been read into R and manipulated within the constraints of the existing file structure. By this I mean that for example you can change the coordinates of points and write them back out. Or you can add a bunch of blank columns in the

Dear R users, I use R 2.4.0 on an iMac running Mac OS X 10.4.8, with a 2.16GHz Intel Core 2 Duo and 2GB 667 MHz DDR2 SDRAM. > sessionInfo() R version 2.4.0 (2006-10-03) i386-apple-darwin8.8.1 locale: es_ES.UTF-8/es_ES.UTF-8/es_ES.UTF-8/C/es_ES.UTF-8/es_ES.UTF-8 attached base packages: [1] "methods" "stats" "graphics" "grDevices"

Hi, I meet a strange trouble writing shapefiles with write.pointShape() from maptools. > write.pointShape(mylines[,5:9], file="Sxtplinshp", mycoordspt) Error in write.pointShape(mylines[, 5:9], file = "Sxtplinshp", mycoordspt) : shapefile names must conform to the 8.3 format It seems that the function request a special shapefile name spelling (conformed to 8.3

Hi, I am trying to compile a package i am developing on my home computer where i am running Windows NT 4.0 service pack 6. Everything works fine on my other computer at work running windows 98. I am using the cygwin 1.1 bash shell and tolls with the MingW compiler, version 2.95.2-1 installed seperately and first in the path. I have active perl 6.16 as well. Here is the error I get when i run

Hi, we all know that R is powerful in dealing with graphics. I am now trying to use R to show some geographic attributes for a specific geographic area on a graph. For example, if I would like to generate a graph showing the population densities (in terms of the intensity of colors) for Illinois State by zip codes. I did some research; and it seems that package "maptools" is

Dear R users, I would like to draw map and import it in maptools/spatstat packages. The 'raw data' I have come from a web page (<map>...</map>) and are basically a list of coordinates of a polygon. I would like to know how to import them in R; I checked the maptools packages, but all the examples use existing .dbf files. I just have a (serie of) text file(s) looking like

thanks to all the help so far I have code now that works on windows (albeit very clunky) that can read the pertinent info from the shape header and return a list with all the necessary parameters to read info out of the shape file. Now come the big design questions 1) I have considered building a map object, basically a list that would hold the info from the shape file and a data frame as part of

Hi ! I have a shapefile that I can easily read into R using library(maptools). My problem stems from some warning messages that come even though everything seems to work fine. library(maptools) districts <- read.shape(filen = "a_ds", dbf.data = TRUE) length(districts$Shapes) so far so good. when I try and plot this I get a nice plot and a warning message: > plot(districts)

Hi all, I've written some code to generate 4 maps per screen and write the output to a jpeg. The output is fairly quick at the start (about 5 jpegs per minute) and then slows down greatly (1-2 jpegs per minute). Is there some way to speed it up? One of my thoughts is to keep the base map static on the screen and just update the points that are being plotted on the map (with the exception

Hi all, I have a set of 54 files that I need to convert from ASCII grid format to .shp files to .bnd files for BayesX. I have the following R code to operate on those files: library(maptools) library(Grid2Polygons) library(BayesX) library(BayesXsrc) library(R2BayesX) readfunct <- function(x) { u <- readAsciiGrid(x) } modfilesmore <- paste0("MaxFloodDepth_", 1:54,

Hello, I would like to place region names (or abbreviations if space is too limited) on a map. I have tried the following, copying and adapting various examples I have seen but without success. The first example is self-contained, but the second one requires the .shp file, which can retrieved at http://www.gadm.org/country. ----- Example 1 ----- library(sp) con <-

My sincerest apologies, as this is a very elementary problem, but I have searched through archives and FAQs and on the web, and I am at the end of my own resources. I need to do analysis of a spatial point process occuring with four counties, using spatstat. I've had no trouble importing the shapefiles, but I can't seem to figure out what I need to do in order to use the boundary as an

Hi, The following function does work optimize.W<-function(W,k,G,Groups,cph.call,z){ n<-length(Groups) grp.wt<-rep(0,n) for(i in 1:(length(G))){ ind<-Groups == G[i] if(G[i]!=k){ grp.wt[ind]<-W[i] } elsegrp.wt[ind]<-1 } z<-data.frame(cbind(z,grp.wt=grp.wt)) #needed to make the case weights #part of the data

Hola a todos. Estoy utilizando R , más concretamente los paquetes maptools y sp, para realizar mapas de España por comunidades y provincias y de Andalucía por provincias y municipios. El caso es que para el mapa de municipios en Andalucía, una vez leído el fichero shp con readShapePoly , si luego por ejemplo , quiero poner todos los polígonos de un color utilizando library(maps)

Hi, I am having trouble with the survival library, particualrily the coxph function. the following works coxph(jtree9$cph.call,z,rep(1,dim(z)[1])) Call: coxph(formula = jtree9$cph.call, data = z, weights = rep(1, dim(z)[1])) coef exp(coef) se(coef) z p SM 0.2574 1.294 0.0786 3.274 1.1e-03 Sex -0.1283 0.880 0.1809 -0.709

I downloaded the latest version of R (1.5.0) and issued a read.table() command that used to work with (1.3.1). Unfortunately I erased 1.3.1 from my hard drive after I downloaded R 1.5.0. Here's the command and the error message I get under 1.5.0: > golub1<-read.table("C:/Microarrays/code/data_set_ALL_AML_train.txt",sep="\t",q

Hi, I have been thinking about your questions and comments. And I think I have the answers. :D > the way you draw the polygons has the consequence that if the larger one > is drawn after the smaller one, > the color of the smaller one is "overwritten" > and therefore the map is really wrong in the relevant area. > so this needs fixing. Another approach is sorting the