similar to: sppolot: fill below minimum legend value

In the code to follow, I'm trying to label points with their corresponding values but have been unsuccessful after many variations to the following code. The code below creates the plot I want, I simply cannot get the black points ("+") to display the actual value. I'm guessing the problem is somewhere in the second to last line of code (starts with "pts<-..."). I

Hello- After looking through ?spplot, I would expect that I could specify the values of the cuts: "...?cuts? number of cuts or the actual cuts to use..." So in the following command, spplot(lzm.krige.dir["var1.pred"], scales=list(draw=TRUE), xlab="Easting",ylab="Northing",

To Whom It May Concern: I have an R package and want to put this package be part of R-project and available to anyone who is interested in. The R package is created for my paper, titled "Acceptance Sampling Plans from Truncated Life Tests Based on the Birnbaum-Saunders Distribution for Percentiles". The paper has been accepted by Communications in Statistics: Simulation and

Hi, Is there a simple way to suppress the plotting of polygon boundaries with spplot() ? # simple list of 12 colors cols <- brewer.pal(12, "Paired") # plot pile of polygons, with 12 classes: spplot(x, zcol='class2', col.regions=cols, scales=list(draw=T), xlab="Easting (m)", ylab="Northing (m)") ... seems to work well. However the polygon boundaries

Hello I'm working with a large hydrological data set stored in a netCDF format. The file stores x and y coordinates in the UTM projected coordinate system, yet when I use image to graphically display the z variable, the image is distorted in the sense that it does not plot the map in the correct spatial organization. I'm wondering if I need to define the projection of the netCDF file

Dear R useRs, I have a problem with nls.lm function of minpackl.lm package. I need to fit the Van Genuchten Model to a set of data of Theta and hydraulic conductivity with nls.lm function of minpack.lm package. For the first fit, the parameter estimates keep changing even after 1000 iterations (Th) and I have a following error message for fit of hydraulic conductivity (k); Reason for

hello mailing list! i still consider myself an R beginner, so please bear with me if my questions seems strange. i'm in the field of biology, and have done consecutive hydraulic conductivity measurements in three parallels ("Sample"), resulting in three sets of conductivity values ("PLC" for percent loss of conductivity, relative to 100%) at multiple pressures

Greetings, I have a groupedData (nmGroupedData) object created with the following syntax: Soil <- groupedData( ksat ~ conc | soil_id/sar/rep, data=soil.data, labels=list(x='Solution Concentration', y='Saturated Hydraulic Conductivity'), units=list(x='(cmol_c)', y='(cm/s)') ) the original data represents longitudinal observations in the form of:

I am trying to fit a nonlinear regression to infiltration data in order to determine saturated hydraulic conductivity and matric pressure. The original equation can be found in Bagarello et al. 2004 SSSAJ (green-ampt equation for falling head including gravity). I am also VERY new to R and to nonlinear regressions. I have searched the posts, but am still unable to determine why my data come up

Hello fellow R-users,   I’m stuck with something i think is pretty stupid, but i can’t find out where i’m wrong, and it’s turning me crazy!   I am doing a very simple linear regression with Northing/Easting data, then I plot the data as well as the regression line :   > plot(x=Dataset$EASTING, y=Dataset$NORTHING) > fit <- lm(formula = NORTHING ~ EASTING, data = Dataset) > abline(fit)

I am trying to write a loop to forecast realized volatility over successive days for the purpose of VaR prediction using the HAR-RV-CJ model which is as follows: log(RV_t+1) = ?_0 + ?_CD log(CV_t) + ?_CW log(CV_t-5) + ?_CM log(CV_t-22) + ?_JD log(J_t + 1) + ?_JW log(J_t-5 + 1) + ?_JM log(J_t-22 + 1) + e_t where RV is realized volatility, CV is continuous volatility and J is the jump which is

Full_Name: Drew Tyre Version: 1.4.0 OS: Win 2000 Submission from: (NULL) (140.253.74.225) When creating a rather complex plot, the resulting figures 'shrink' when either copied to the clipboard or saved as a metafile using the GUI menu. Also saving as postscript or pdf causes the same problem, as does printing the figure directly from the graphics window. Only the figure areas shrink,

DeaR comRades: I have a 2D spatial binomial process as shown in the data and code below. I am plotting the number of trials and the number of successes in the spatial binomial experiments and would like to draw the spatial cells were the trials and successes were counted, i.e. a partial grid in the plot only for those cells where there is a number. The cells are 2x2 km cells. The count of Trials

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, 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 have some geospatial data where two layers are thematic and the third is a percentage, so the maps need to have different themes. thumbDf <- as( stack( thumb), "SpatialGridDataFrame") names(thumbDf at data) <- c("pri", "sec", "pct") thumbDf at data$pri <- factor(thumbDf at data$pri, levels=c(0:8), labels= names( peelClasses)) thumbDf at

Hello, I want to create side-by-side maps of similar attribute data in two different cities using a single legend. To simply display side-by-side census block group boundary (non-thematic) maps for Minneapolis & Cleveland I do the following: library(rgdal) library(sp) Minneapolis=readOGR("../Minneapolis/Census/2010/Census_BlockGroup_GEO/","tl_2010_27053_bg10")

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

Hi! How does one change the background color of the map-panel in spplot()? Example: library(sp) data(meuse.grid) gridded(meuse.grid) = ~x+y spplot(meuse.grid, "part.a") How would I get another background-color for the map-panel (but not for the whole plot) here? Thank you! Marcel