similar to: map issues

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"

Dear all, I would like to produce a map with information about the patenting activity in German districts, by coloring districts with different degrees of patenting activity in different colors. I work with the packages maptools, maps and spdep. The map data is read from an external .shp file (+ the corresponding .shx and .dbf files). Plotting a map with the IDs or the patenting indicator itself

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

Hi, I was plotting or creating a map for Texas school districts using the shape file of Texas. I could not find any other helpful mail in the mailing list. txshp<-read.shape(system.file("S:\\Districts_10_11.shp", package="maptools")) Error- read.shape no found. But read.shape is there in maptools. If anyone can help me out it will be great. Thanks in advance. Shant

I am working on a project to unite several local school districts. We will have 14 different districts, every district would have their own asterisk box on location. We all have fiber lines running to a central location at our isd. This provides connectivity to all the districts. 1. would it be wiser to install a asterisk box at the central location of the ISD and have all asterisk boxes

Dear all, I want to build a model in R based on animal collection data, that look like the following Nr Village District Site Survey Species Count 1 AX A F Dry B 0 2 AY A V Wet A 5 3 BX B F Wet B 1 4 BY B V Dry B 0 Each data point shows one collection unit in a certain Village, District, Site, and Survey for a certain Species. 'Count' is the number of animals collected in that

Hi is there a way in R to create a new column vector with the numbers of 2 others columns simply appended? Example : I have a column with provinces codes (1 to 19) I have a second column with districts codes (1 to ##, depending on the province) I want to create a third column with a code that would be unique for each district: Province  District   District_Unique 1             1            11

Hello! I am doing an analysis on a questionnaire of hunters taken in 4 different districts of some mysterious foreign country. The aim of the study was to gather info on the factors that determine the hunting success of a peculiarly beautiful bird in that area. All variables are factors, i.e. they are variables such as "Use of Guns - yes / no", "Use of Dogs - yes / no" and

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

I am trying to use the value of an ID variable in an if statement and not getting the results I expected. # ID values for two school districts > with(rf, tapply(DistrictID, DistrictName, min) ) Aberdeen School Dist. # 58 Buhl Joint School District 59340 53409 This creates DNAME as I expected ...

I have just learned how to play with map, but something weird (or not) is happening. Suppose I want to draw a map of two countries (that have disconnected components), like Argentina and Brazil. If I command: library(maps) library(mapdata) map("worldHires", c("Argentina", "Brazil")) It works fine. However, if I want to _colour_ the interior:

No, you're quite right, I was wrong. I must have been very confused when I put up the first 'bug' report, #886. I noticed my error a couple of days after and put up a comment to that effect on R-bugs, but it seems to have vanished (or alternatively I did that wrong as well). I agree with you on the criterion for the choice of districts as well. There remains a minuscule discrepancy

Hello R Folks, I am new to R. I have been struggling to overlay a barchart with a xyplot together on one plot but did not get this worked out. Any help and idea are greatly appreciated. I attached R scripts for barchart and xyplot below and also data I used. What I am trying to do is just to put the barchart and xyplot together on one plot. Huapeng <<dispersal_infestation.csv>>

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

Selling analog GSM gateways at 119.00 USD. Polarity reverse for "call answer" and "call end". Compatible with all Digium/Asterisk analog FXO cards. Suitable for low-traffic call termination, inbound calls from cellular network or for phone bill optimization applications. Main Functions - Provide Accurate Anti-Polarity Signal - Incoming Calls Display - Add Prefixal

We have had another instance of this..... since I am forwarding to other lists, "this" involves a lost file, due to accidental deletion. in this case, we had a backup, but from the backup time, till deletion time, a lot of data had been lost. So, we have not enough disk space to do hourly backups, novell allowed recovery of a lost file like this, so: is there a filesystem that we

Hi, Is there any way to plot a Map with different shading patterns, instead of color intensity, in R2.0.1 over Linux ? Thanks. D ----------------------------------------- Stay ahead of the information curve. Receive GIS news and jobs on your desktop daily. Subscribe today to the GIS CafeNews newsletter. [ http://www10.giscafe.com/nl/newsletter_subscribe.php ] It's informative and

I have setup SAMBA to authenticate to a 2003 (probably SP1) domain and I am having two very distinct problems. I have searched the archives, but each time I see this question is ask I never see any answers and I need an answer or I am going to have to shut this off. I work in the Kentucky school systems and we have one AD forest containing one domain for each county in the state. That totals up

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

For a sanitation project in Bangladesh I need to design a three stage sample survey to be representative of around 40 million people. I find myself suddenly with several challenges with which I struggle and would be gratefully for any help. As the questions are linked I kept them together rather than creating multiple posts 1) SURVEY DESIGN So far I designed mainly two stage cluster surveys but