similar to: plot.Map Maptools

Hi, I have 142 observations off different variables and I'm trying to do a OLS. And I get this error. Any Ideas ? > f <- ols(lnmigr3~popden78+income+modern+spareha+rain) Error in La.chol2inv(x, size) : size cannot exceed nrow(x) = 1 D. ----------------------------------------- Stay ahead of the information curve. Receive GIS news and jobs on your desktop daily. Subscribe today

Hi all, I observed that the 'link' between R and BUGS (winBugs/linBugs) has totally disappeared. I'm wondering what package can be use to run bayesian models specified using the BUGS language in R (specifically under Linux). Is there any other option besides JAGS ? Thanks ----------------------------------------- Stay ahead of the information curve. Receive EDA news and jobs on

Is there any way to read and decompress inside R tar and bzip2 files ? Thanks D ----------------------------------------- Stay ahead of the information curve. Receive EDA news and jobs on your desktop daily. Subscribe today to the EDA CafeNews newsletter. [ http://www10.edacafe.com/nl/newsletter_subscribe.php ] It's informative and essential. This message was sent to you from a machine

Installed samba 3.0.4.0 and libiconv 1.9.1 (download from bullfreeware) on AIX 5.2. When I run smbclient, it says error: exec(): 0509-036 Cannot load program ./smbclient because of the following errors: 0509-150 Dependent module libreadline.a(libreadline.so) could not be loaded. 0509-022 Cannot load module libreadline.a(libreadline.so). 0509-026 System error: A file or directory in the path name

Hi I'm working with the function Kest in the package SpatStat (under LINUX with R 2.1.0). In order to evaluate the statistical significance of my point pattern I'm doing 999 Montecarlo replications. The script that use the Kest function runs OK for most of the different point patterns that I have but for a particular point pattern, which have only 17 points, it runs until the 34th

Dear Roger, I am trying to use the write.polylistShape() function of maptools for the first time and realize that it handles list of polygons of class 'polylist'. However, it seems that no as.polylist() function exist in the package. The question behind that is: in your opinion, which would be the best way to convert a list of matrix of polygon nodes coordinates into an object of

Hi everybody, I've been looking around, but can't seem to find the answer. I get a list of names (which vary, so I never know them in advance), and transform them into variables, each of which is a dataframe, using assign: polyList <- c("rs123", "rs124", "rs555", "rs000") numPoly <- length(polyList) for (k in 1:numPoly) {

Hi, I'd like to get an argument (I think it's the right term) dynamically from a list, but cannot manage to do it. Here is the code I use: #output given by database BOB <- c('A/A', 'C/C', '15/27') MARY <- c('A/A', NA, '13/12') JOHN <- c('A/A', 'C/A', '154/35') CLIFF <- c('A/C', 'C/C',

Hello, I'm fairly new to R, and despite spending quite some time exploring, testing, and looking for answers, I still have a couple of questions remaining... The first one is about creating a dataframe using a list of the variable names . I get this output file from a database: BOB <- c('A/A', 'C/C', '15/27') MARY <- c('A/A', NA, '13/12') JOHN

Dear All, In the moment i'm using the map and maptools package to read shapefiles and display the maps. I'm looking for the possibility to draw points (randomly positioned or positioned according to a grid) into the polygons instead of filling the polygons with colors. For example: a map (shapefile) with 10 countries, 15 points in the polygon of country A, 20 points in the

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

Dear listers, I am working since a while with the sp package and still wonder how the plot methods are managed with sp spatial objects. For instance, SpatialPolygonsDataFrame objects have obviously a plot method. However it cannot be found in the list provided by methods(plot) . Furthermore ?plot.SpatialPolygonsDataFrame, nor ?plot.SpatialPolygons, etc.. provide a help, though the lattice

Greetings: I would like to kindly ask for a little help. The rough code is: #________________________________________________________ dat=data.frame(read.delim(file="all.txt", header = TRUE, sep = "\t", quote="\"", dec=".",na.strings = "NA")) nc=read.shape("astae.shp", dbf.data=TRUE, verbose=TRUE) mappolys=Map2poly(nc) submap

What tools and documentation exist for producing maps and working with polygon and polylist objects in R beyond the following: * Packages maps, mapproj, maptools, and mapdata. * Becker and Wilks (1993, 1995), listed as references in the documentation on the "map" function. Thanks, Spencer Graves

Hi Guys Thanks for this! As I am a beginner, I bet I'm running into some really basic problems. Using the example from the Map2poly function in maptools: try <- read.shape("euadmll.shp") mappolys <- Map2poly(try) # this produces 14 warnings() of the type: # 1: level 2 circularity at 62 in: .mtlbuild(.mtafters(r1), rD) # etc # 14: From next release, default hole handling

Hi all, I''m thinking of useing sub-domains as account keys, something I''ve played with previously. I have accounts and users, fairly standard, except that users can belong to many accounts, and accounts can have many users. So I figure I need to have an area on the main domain, that allows users to view their account memberships and things. Sort of like a dashboard for each

Using the maptools function "map2poly" I have created a map/polylist object - a map of the Great Lakes. My goal is to passively plot additional data on this map. Unfortunately I am not able to change (distort) the relative scale of the axes for the map, and so the sample points do not line up correctly on the map. In the code below, notice that I have set the same axis ranges on each

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)

I''m trying to add a "node type" variable that can be interpolated into source values. In this, my latest attempt, $nodetype always ends up as "base": # manifests/nodes $nodetype = "base" node "base" { include postfix ...snip... } node "www01", "www02", ... , "www10" inherits "base" { $nodetype

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"