similar to: help color coding map in R

Hola, no me funciona: Te contesto lo mismo que le contesté a otro compañero de la lista, espero ser más claro. Les comparto la carpeta con la capa y los datos que estoy usando. https://www.dropbox.com/s/umy0evt3qm1wr4d/dissolve.zip Esto es lo que estoy haciendo: library(rgdal) library(maptools) library(sp) library(classInt) library(RColorBrewer) zm <- readOGR(?.?, ?zmdis?) data <-

Hola, Intento representar en un mapa participaciones porcentuales de los sectores económicos y no logro hacerlo con spplot. He intentado con spplot(zm, c(?part88?, ?part93?), cuts=4, col.regions=brewer.pal(4, ?Set3?)) Lo pude hacer utilizando el base graphics de R, definiendo: > plotvar88 <- zm$part88 > nclr <- 8 > plotclr <- brewer.pal(nclr, "PuOr") > plotclr

Hola, El problema con la propuesta de Olivier es que los intervalos son diferentes para cada variable. La forma sencilla y rápida es: spplot(zm["part88"], col.regions=plotclr, at=class$brks) Pero para que quede más elegante hay que dar algunos pasos más: ## Intervalos en forma character op <- options(digits=4) tab <- print(class) options(op) intChar <- names(tab) ## Indice

Hola, Lo hice de la forma en que propones y no me sale, no sé si algo esté mal con R, porque no me respeta ni el ?key.space?, me lo sigue poniendo en la parte izquierda del gráfico Saludos El 06/05/2014, a las 02:21, Olivier Nuñez <onunez en unex.es> escribió: > El segundo argumento en spplot ha de definir el número/nombre del intervalo al que > pertenece el area a pintar. >

I'm using the following code to create quartile intervals of my dataset and it is working like it is suppose too. In my final output of the data I would like currency amounts to be displayed with commas. I have not found an option either through FormatC or sprintf that will display the data with commas, they only seem to add decimals. Is there a format that will display commas without me

Hello: I have a problem when I tried to color the USA map with different colors. The following is my data (I only used the second column of data): alaska, 1, 2 Hawaii, 0, 0 USA, 5, 5 And here is my code: library("latticeExtra") library("mapproj") state<-read.table("C:\\usaclass.txt",sep=",") state

Dear R-listers, I am using R for Windows Vista 64-bit. I have no experience working with maps, but now I have to plot a map of Denmark in which each Danish county is color-coded according to its incidence rate of a particular disease. My problem is that I don't know where to start. I have read the help files to packages like 'mapplot', and their examples indicate that I am on the

I create a lot of plots with a dollar amount as the y-axis. How do I change the formatting so comma and $ are in place. Like $1000,000 or to use common dollar abbreviations like 100K? Chris Anderson Data Analyst Medical Affairs wk: 925-677-4870 cell: 707-315-8486 Fax:925-677-4670 This electronic message transmission, including any attachments, contains information which may be confidential,

I am trying to create multiple dataset by group like the following using either a loop or vectorization: vocallGp1<-subset(vocallsub, Group==1) vocallGp2<-subset(vocallsub, Group==2) vocallGp3<-subset(vocallsub, Group==3) vocallGp4<-subset(vocallsub, Group==4) vocallGp5<-subset(vocallsub, Group==5) When I create the loop I get the following error:

I'm trying to display my deciles without scientific notation, but have not found an option that will allow me to do so. According to web searches, the options(scipen=999) should remove scientific notation, but it seems not too. Does this option work with quantcut function? Is there any other option that can be used? I am using verison 2.10.0. library(gtools) library(gdata)

I notice that R basic installation does not take advantage of multi-thread or multi-core processors. Is there an option for multi-threaded processors? I currently run all of my analysis on thinkpad laptop and most processing is very quick, however, when I run stepAIC with interaction terms it can take a couple of hours to complete. I have tried running the same process and a quad core machine 64

I have a process where I am creating a effects plot similar to the cowles effect example. I would like to add the point estimates to the effects plot, can someone show me the correct syntax. I have included the "R" effects example, so you can show me the correct syntax. Thanks mod.cowles <- glm(volunteer ~ sex + neuroticism*extraversion, data=Cowles, family=binomial)

r-help lattice adepts: I have a question which is somewhat geospatial, so I posted to r-sig-geo rather than here: https://stat.ethz.ch/pipermail/r-sig-geo/2012-November/016757.html > summary: How to overlay a geographical map on each panel in a lattice > (or Trellis), e.g., of levelplot's? Note I am not inquiring about > creating choropleth maps[,] which Sarkar 2008 covers quite

Is there a way to generate maps in R. Specifically, I have calculated estimates of intra-regional inequality for US states, and would like to project that information onto a map. Thanks, Nirmala

I need to make a fairly complex animated graphic and decided to use grid for it. A very simple example of what I need: ##============================================================================== library(grid) grid.newpage() pushViewport(plotViewport()) pushViewport(viewport(xscale = extendrange(c(0, 100)), yscale = extendrange(c(0, 100)))) grid.xaxis() grid.yaxis()

Hello dear developers, I find myself often having the result of "range" oder "extendrange", which I want to create a sequence with. But "seq" needs two seperate arguments "from" and "two". Could an argument "range" be added? Otherwise I will have to create an object with the range (may come from a longer calculation), index twice from

Hi, I often need to extend the plot range to the right, but not to the left (for example: not below 0 so that log = "x" still works...). This could be a handy improvement of extendrange(): --- utils.R 2015-08-25 18:18:20.000000000 -0400 +++ utils.R 2018-04-25 17:21:25.000000000 -0400 @@ -30,6 +30,6 @@ ## Purpose: extend a range by a factor 'f' - on each side

Dear all, How can I make an index plot with lattice, that is plotting a vector simply against its particular index in the vector, i.e. something similar to y <- rnorm(10) plot(y) I don't want to specify the x's manually, as this could become cumbersome when having multiple panels. I tried something like library(lattice) mp <- function(x, y, ...) { x <- 1:length(y)

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Symptoms: curve with log x axis gets the wrong x limits (in 2.6.1, I believe introduced in this version). Credit goes to Mike McCoy for detecting the problem. Demonstration: x = 1:5 plot(x,5*exp(-x),xlim=c(0.1,1000),log="x") xvec = 10^seq(-1,3,length=100) lines(xvec,5*exp(-xvec)) curve(5*exp(-x),add=TRUE,col=2,lwd=3) I believe

Dear R-experts, I have two problems: PROBLEM (1) ----------- I want to produce a very small png file (35 x 18 px) that contains a histogram without a figure region or margins, only the pure heights. In the base graphic system this is simple: png(filename = "hist.png", res = 72, width=35, height=18) par(mar=c(0,0,0,0), oma=c(0,0,0,0)) hist(rnorm(100),