similar to: masking a regular lat/lon grid to extract map boundaries

Hi all Just this morning I upgraded to R 2.12.0 (for Mac OS X 10.6.4). All went well until I needed to run a help() or help.search() in my session, which I'm running within Emacs (ESS 5.3.7). Say I need help with the command 'density'. When I type help(density) or ?density the ESS help buffer opens, it is titled *help[R](density)* but it contains only a couple of lines saying,

Dear R-helpers, I'm plotting geophysical data in the form of contours using "filled.contour". The display would be much more effective if the areas with negative values could be color coded by -- say -- "cold colors" in the range of blue to green, and conversely the areas with positive values got plotted with "warm colors", from yellow to red. Right now if I use

Hello I'm trying to render in 3D what I usually plot by image(), or image.plot() from the library fields, followed by a map("world",add=TRUE) type of command. More concretely, I have a field of temperature values, for a given geographic area, and I would like to plot a 3D surface whose x and y axes consist of longitude and latitude values and whose height and color-coding correspond

Hi, I'm trying to plot a field obtained from the atmospheric model WRF-NMM which uses a "Rotated Lat-Lon? map projection. The WRF documentation mentions that: ? Rotates the earth's lat/lon grid such that the intersection of the equator and prime meridian is at the center of the model domain. ? Within the rotated framework the grid spacing is constant, but in an earth-relative

Given a set of latitude and longitude coordinates pairs (stored in variables "latitudevals" and "longitudevals"), I would like to plot them onto the image of a equirectangular world map. I would like to plot each coordinate pair with a red circle, if possible. Does anyone have any suggestions as to how I go about doing this, whether using R or using another program like Google

Dear R-listers, Here, I would like to hearing helps from you. I have GPS data (multiple points in the geographic scale) in longitude/latitude. I intend to calculate distance (in kilometer) among such points and output the distance matrix in dist class. I have gotten some progress, but I still can not get final goal. Could please give me any directions/advice? This email cc. to Mr. Pierre, the

Hi, I'm trying to read a subset of a netcdf file into R, but although I'm relatively experienced using R, I'm still new to netCDF files, so this may be a very simple/stupid question! I've included an example of the type of file I'm looking at here. www.met.reading.ac.uk/~swp06hg/ccd1983_01-dk1_20.nc (~7Mb) It's a 2D array of the variable CCD along with its lat and

Here is the answer to an old question of mine, regarding maps in R Hope this help. claudia ---------- Forwarded message ---------- Date: Mon, 25 Jun 2001 13:48:05 -0400 From: Kieran Healy <kjhealy at Princeton.EDU> To: Claudia Tebaldi <tebaldi at ucar.edu> Subject: Re: [R] maps in R Hi Claudia -- you can get the maps() library from here:

Hi All, I have a netcdf gridded file with LCC projection. I can easily use image.plot to visualize it. However, as the axises are in X,Y, not Lat and Lon, I could not add state or country maps onto it (or lat lon information). I do have a grid2d file that describes the lat and lon for each (X,Y) grid, but the lat and lon are not regularly spaced, so I could not use image.plot. Does anyone know

I'm running R-1.4.1 on a linux debian (version 2.2) box. I've tried to upgrade to R-1.6.0, and the execution of make halts with the following error messages: make[4]: Entering directory `/home/tebaldi/R/R-1.6.0/src/library/methods' dumping R code in package 'methods' Error in testRversion(descfile) : This package has not been installed properly See the Note in ?library

Hello, I am trying to match lat/lon from one dataset with the lat/lon from a second dataset and use that rows data for calculations. I am using match, but this is finding the first match and not comparing the pair, how can I determine if the lat/lon are the same? See example below. Is there a better way to determine to a matching pair of lat/lon values? Example Datasets: > data2

Two questions: 1. Is the order of the dimensions shown what is shown if you look at str(ncin) - I mean shown at the end where it describes the variable and its dimensions? 2. Is you problem how to subset the netcdf file, how to write to the .csv file, or both? -Roy > On Aug 28, 2017, at 2:21 PM, Eeusha Nafi <eshad002 at gmail.com> wrote: > > I have a series of nertCDF files

Hi all I have a dataset whose response is a categorical variable of the ordinal scale type (6 levels). I'm interested in building classification models, and I'm wondering if there is something implemented in R (or its packages) that I'm not aware of, to treat the ordinal scale measurements straightforwardly. I can think of the alternative of building a conditional (hierarchical)

Hello I have several replicates (30) of a -- short -- time series (19 obs) and I'd like to estimate the first few lags autocorrelation function. Is there an immediate way of handling multiple realizations when computing the ACF in R? I'm using R-1.4.1 on a Linux Debian box. Thanks in advance claudia ------------------------------------------------------------------------- claudia

There is probably an elegant and efficient way to do this -- perhaps already implemented, and some of you knows it. Myself, I'm getting stuck in long loops ... Here is the question: I have a long (80,000 + obs) binary time series. For each of its elements I'd like to compute and save in a vector "the number of steps 'til the first switch of state happens". For example,

I have a series of nertCDF files containing global data for a particular variable, e.g. tmin/tmax/precipiation/windspeed/relative humuidity/radiation etc. I get the following information when using *nc_open* function in R: datafile: https://www.dropbox.com/s/xpo7zklcmtm3g5r/gfdl_preci.nc?dl=0 File gfdl_preci.nc (NC_FORMAT_NETCDF4_CLASSIC): 1 variables (excluding dimension variables):

Hi I'm trying to print several plots to a multi-page postscript file. Some plots are output of basic commands, some are produced by "trellis" commands (from the package lattice). Right now I'm not able to get a straightforward black and white color scheme for the latter kind: when I open the postscript file the "trellis" plots are invisible, or parts of them are

#################################################################################### We have used maptools to construct state, county, township, census-tract, and zipcode level R maps with an Albers Conical Equal Area projection. We would like to be able to plot the location of weather stations or other point locations on the maps. The data the point locations are in latitude-longitude units

Hi, I have a 1.000 observations with 10 attributes (of different types: numeric, dicotomic, categorical ecc..) and a measure M. I need to cluster these observations in order to assign a new observation (with the same 10 attributes but not the measure) to a cluster. I want to calculate for the new observation a measure as the average of the meausures M of the observations in the cluster

Hi all I'm running lm()/leaps() kind of computations within a loop, and I submitted the entire thing by sourcing a file of R code in interactive mode. I'd like to interrupt the evaluation of this (it's taking way too long) but I am not able to. I'm running R-1.3.0 from an emacs window (ESS-5.0) and I'm on a linux box, (Debian). I tried Control-c Control-c (nothing happens),