similar to: About filled.contour function

Dear all I use the library(netCDF) to read in NCEP data. The file I want to read has size 113 Mb. When i try to read it I get the following message: Error: cannot allocate vector of size 221080 Kb In addition: Warning message: Reached total allocation of 255Mb: see help(memory.size) I get a similar message when I try to read a file with 256 Mb in a PC with 2 GigaByte RAM. Is there something

Dear R-project users I would like to calculate a linear trend versus time taking into account a first order autoregressive process of a single time series (e.g. data$S80 in the following example) using th gls function. gls(S80 ~ tt,data=data,corAR1(value, form, fixed)) My question is what number to set in the position of value within corAR1? Should it be the acf at lag 1? I look forward for

Dear all I use the gls function but in contrast to the lm function in which when I type summary(lm(...))$coef I receive all the coefficients (estimate, Std. Error, t-value and pvalue), with gls when I type summary(gls(...))$coef I only receive the estimate of the reg. coefficient without std. error and t- and p-values. Dou you have any suggestion how to solve my problem? With kind regards

Hi everyone I have a matrix (let's say that it contains the values of elevation) and want to plot its values on a map using a function such as image.plot or filled.contour. The problem is that my grid is not rectilinear, it is bended. Here is an example lon<-matrix(0,20,25) lat<-matrix(0,20,25) elev<-matrix(0,20,25) for (i in 1:20) { for (j in 1:25) {

The following: library(grid) grid.newpage() example(Grid) has the yaxis label partly cut off and the x axis label does not appear at all. Also ?grid.multipanel in that example brings up documentation for grid-internal but this would not seem to be internal if its part of an example. I am using: > R.version.string # Windows XP [1] "R version 2.2.0, 2005-09-20"

Hi all, My question is not really urgent. I can write a loop and solve the problem. But I know that I'll be in a similar situation many more times so it would be useful to find out the answer Is there a fast way to perform linear fit to all the columns of a matrix? (or in the one dimension of a multi-dimensional array.) I'm talking about many single linear fits, not about a multiple fit.

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Estimado Javier, adjunto el grafico retocado con GIMP, borre las lineas que sobran en cada panel. Obviamente, no es la idea hacerlo de este modo, sino con R, pues quiero cambiar el color de la recta y hacer varios otros graficos, cuya edicion manual requeriria mucho tiempo. Espero entiendas mejor mi explicacion ahora. Cada una de esas rectas es un

I cannot achieve a multipanel display of several levelplots. Here is part of my code. At first, create the layout: push.viewport(viewport(layout=grid.layout(1, 2))) Then for the left panel: push.viewport(viewport(layout.pos.col=1,layout.pos.row=1)) push.viewport(viewport(width=0.6, height=0.6)) levelplot(z~x*y, grid) pop.viewport() Similar code for the right panel. However, there is always

Hi. I'm a big user of level plots. Because filled.contour() didn't do exactly what I wanted, I hacked at it a bit to add these options: (*) Even spacing in the key, even if the contours are not evenly spaced. (*) Horizontal key (*) options to not use layout() and to plot only the level plot _or_ only the key. This means you can use layout yourself to (e.g.) have several

Dear R-users, I am trying to create a multipanel scatter plot which displays weighted regression lines for each subgroup. I am having diffuculty figuring out how the Trellis multipanel plotting functions work, in general. Does someone know where there are some beginners examples on how to create and use panel.xxx functions in conjunction with the xyplot functions? Back to the matter at hand

Is there a simple way to specify a theme or trellis (lattice) parameters so that, in a multipanel (conditioned) plot, there is no color and in the strips there is no shading? This is the effect achieved on page 124 of Deepayan Sarkar's "Lattice" (figure 7.2). I managed to trick lattice into making a grayscale plot on my interactive display as follows: > graphics.off() >

Full_Name: Karel 'Clock' Kulhavy Version: 1.9.0 OS: GNU/Linux Submission from: (NULL) (212.71.168.94) par(mfrow=c(2,1) is ignored by filled.contour().

Dear All, I tried to draw 3 filled.contour plots on one page, but failed. It returns "Error in plot.new() : Figure margins too large". Can anyone help me to do this. BTW, does anyone has the experience to draw a isosurface with R? Dr Zhiqiang Cui School of Geography and Environmental Sciences The University of Birmingham Edgbaston, Birmingham B15 2TT

I'm trying to add a set of reference lines to a multipanel xyplot xyplot(y ~ x | Visit, panel = function(x, y, ...){ panel.xyplot(x, y, ...) abline(v = c(0.5, 1)) }) However, the reference lines are different for different visits. For example, for the first 2 visits, I'd like vertical lines at x = 0.5 and 1. For visits 3 and four, I'd like vertical lines at x = 1 and 1.5. I can

Hello, Is possible to place a label on a multipanel figure outside the boundaries of any single one of the existing graphs in a panel? Specifically, I have a function that creates a panel with quartz() and then fills it with 6 graphs consisting of exactly the same plot for 6 different groups. I would like to add one title for the entire figure at the top of the figure, but I cannot figure out

Derar Rxperts, I have a strange situation.. I see curly brackets around "strip.levels" in multipanel strips while using lattice::xyplot. .How do I get rid of the curly brackets? For some reason, I am not able to reproduce the problem using an example below... Any suggestions are highly welcome! Thanks, Santosh q <-

Hi, I'm making a multipanel lattice densityplot figure with 2 factors (3 and 20 classes in each factor) with the following statement (the type="percent" is there to prevent plotting the actual points which detract from the figure - is there another way of doing this?): densityplot(~End-Begin | Type * Chromosome, data=Mon, layout=c(5,12), xlab="Element

I want to create a simple plot containing three levelplots with one colorbar. I used the "Three levelplots" code below, but the third levelplot is drawn smaller than the first two. However, if I try the "Two levelplots" code below it works well. Can anybody tell me how could I draw three levelplots (of the same size) with one colorbar. Thanks in advance, Joaquin ### Three

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Estimada comunidad, tengo que hacer un grafico con 4 paneles (4 graficos pegados juntos) y sobre este tengo que dibujar 4 rectas ajustadas a los datos de cada panel. Intente hacerlo con lattice usando el siguiente codigo : library(lattice) library(latticeExtra) xYplot(Cbind(mf[,3],mf[,5],mf[,6])~dia|factor(sol),data=mf,xlab="Time

I would like to know if there is a way to use plot(zoo) and suppress the label text to increase the amount of the screen that if for the plot and not the labels. I tried xy.labels=FALSE, but that had no effect. Thanks [[alternative HTML version deleted]]