similar to: Grid building in R

I am trying to produce a simple plot where the points in the plot are actually integers from my data. That is, I am making a length/width plot and I would like the points, rather dots or triangles, appear as different numbers corresponding to a column in my data. I have tried using this: ...,pch="column title" but this merely returns the first letter of the column title, B. Is there

Dear R-Listers, I have spatial data on an equidistant rectangular grid, similar to topographic data. I know that there are quite a few R-packages or base functions that provide nice iso-contours plot, but I don't want a plot, just the smoothed isocontour line of ONE level (e.g. 10 mm). Data sets are large, so it would be preferable if the availability of regular grid data could be exploited,

Hello, I would like to plot a wireframe of a function which is defined on the unit square using the lattice library (for trellis-like plots). The plot contains z-values of about 100 (only in the neighborhood of zero) although most of the z-values are in the range of -let's say- 0 to 10. If I evaluate this function on an equidistant grid of 25 points on the unit square the plot quality is not

Hi, I build a graph taking into account the times: 1, 2,4,6,8,10,12,15,18,21,24,28,32 and 48. Be that the scale of the X axis does not look right. It seems equidistant. (graph attached) What changes have I to do in the following commands so that the scale be correct? interaction.plot(Tempo, Trat, Valor, ylim=c(0, 2), las=1, lty=c(1,2,3,4), lwd=3, bty='l',

I am having difficulty installing packages. For example, here I have tried to install "gplots" but to no avail. I have tried multiple mirrors and different packages (errbar) but after installation I seem to be unable to access any of the commands or help pages. I am a little concerned about the 'lib' missing message but this has been present for previous successful package

Full_Name: Dr. Bernd Winkler Version: 1.2.2 OS: Linux Submission from: (NULL) (194.59.179.176) Hi all, having updated to version 1.2.2 last week I experienced some strange behaviour of the hist function. I want to plot a histogram of frequencies resp. counts, for example poisson random numbers, on a predefined and fixed x scale > hist(rpois(500, 2), breaks=0:10) But instead of the

Hi, if I dermine the break points using the hist() function and then try to re-use these in a new histogram, R fails. Here is an example of the problem: ##First, plot a histogram: data(islands) foo <- hist(islands,freq=T) ##Now, try plot it again, with the previously determined break points: hist(islands,breaks=foo$breaks,freq=T) ##... this lead to the warning message: Warning message:

Having to collect hourly electricity loads and quarter-of-an-hour electricity production data for some years I think that the tidiest way of doing it is to resort to ts but I don't know how to define such a frequency starting from a set date. Leafing through r-help mail archives I've found this *ALMOST* satisfactory message: ==========================================================

(I'm reposting this message because the original has not appeared after about 2 days. Sorry if it shows up twice.) Hello. First, thanks to those who responded to my recent inquiry about using contour() over arbitrary (x,y) by mentioning the interp() function in the akima package. That worked nicely. Now for a new question: I would like to use a pair of prediction intervals to

Suppose I have two variables, x and y. For a fixed number of knots, I want to create a spline transformation of x such that a loss function is minimized. Presumably, this loss function would be least squares, i.e. sum (f(x)-y)^2. The spline transformations would be linear, quadratic or cubic. I know I can solve this problem using some optimization function in R, but I was wondering if anyone

In an experimental setup we obtain z-data samples at equidistant grid points. The surf.gls (Kriging) algorithm produces an error under this circumstance when performing the Choleski decomposition. A workaround is to dither the grid coordinates using (x <- rnorm(length(x)) ; y<- rnowm(length(y))). Question: Is this an expected behaviour of the surf.gls function ? Regards, Berthold

Sorry, to solve your question I had tried: data(faithful) kmeans(faithful[c(1:20),1],10) Error: empty cluster: try a better set of initial centers But when I run this a second time it will be ok. It seems, that kmeans has problems to initialize good starting points, because of the random choose of these starting initial points. With kmeans(data,k,centers=c(...) the problem can be solved.

Dear R users, i try to produce the fold change versus fold change plot where i have the values for x and y ranging from 0.01 to 100. So i start with plot(x,y,xlim=c(0.01,100),ylim=c(0.01,100), axes=F). Then i would like both axes to have tick marks as c(0.01,0.1,1,10,100) but they should appear equidistant. How should i manage this? Thank you for your help, Alla.

Hello, when I tried to create 3D graph from attached file (at first: columns 1, 2, and 3 (longitude, latitude, altitude), then columns 1, 2 and 4 , then the same for column 5 etc) using Akima library for 3D interpolation (the data are neither equidistant nor regularly distributed), the console crashed by SIGABRT fault response. Do you know what is the main problem? Otherwise, can you advise me the

I think there are a couple of things in ?hist that are not quite as clear as they could be. (1) freq: logical; if 'TRUE', the histogram graphic is a representation of frequencies, the 'counts' component of the result; if 'FALSE', _relative_ frequencies ("probabilities"), component 'density', are plotted. Defaults to

Dear All, I was wondering whether someone has created a kernel density evaluator that estimates the density at given specified points. The regular density() function evaluates the kernel at equidistant points, but I am interested in doing such evaluation along a list of values existing in a pre-specified vector. (Similar to the option at() in the kdensity command in Stata). This question has

Dear All, I am a undergraduate using R for the first time. It seems like an excellent program and one that I look forward to using a lot over the next few years, but I have hit a very basic problem that I can't solve. I want to produce a standardised histogram, i.e. one where the area under the graph is equal to 1. I look at the manual for the histogram function and find this: freq:

Hello, my name is Carsten. This ist my first post to R-help mailing list. I estimate densities with the function "density" out of the package "stats". A simplified example: #generation of test data n=10 z = rnorm(n) #density estimation f=density(z,kernel="epanechnikov",n=n) #evaluation print(f$y[5]) Here I can only evaluate the estimation at given

Hello, I am trying to define a different interval for a "year". In hydrology, a "water year" is defined as the period between October 1st and September 30 of the following year. I was wondering how I might do this in R. Say I have a data.frame like the following and I want to extract a variable with the water year specs as defined above:

Hi all, I am wondering whether there is a small bug in the binom.test function of the ctest library (I'm using R 1.6.0 on windows 2000, but Splus 2000 seems to have the same behaviour). Or perhaps I've misunderstood something. the command binom.test(11,100,p=0.1) and binom.test(9,100,p=0.1) give different p-values (see below). As 9 and 11 are equidistant from 10, the mean of the