similar to: Contour plot for arbitrary (x,y,z)

Hello. I have an R program that calls gawk (GNU Awk 3.06 for Windows) from within pipe() to preprocess a large file before it is read into a data frame with read.table(). I've recently upgraded from Win98SE to WinXP, and have also upgraded from R1.5.0 to R1.6.1 over the past month or so. This program worked before the upgrade(s), but now fails. I observe the following sort of behavior with

(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

Full_Name: Ronnen Levinson Version: 1.5.1 OS: Mac OS 10.1 Submission from: (NULL) (12.232.201.92) The value returned by strptime behaves badly after arithmetic operation and/or combination. > a=strptime("2002-06-01 12:15:01","%Y-%m-%d %H:%M:%S") > a [1] "2002-06-01 12:15:01" > a+0 [1] "1932-04-25 21:46:45" > a+3600 [1] "1932-04-25

Hi. The following statement works fine in R-win 1.8.0, but yields a syntax error in R-win 2.0.1 (and possibly in other versions after 1.8.0): plot(c(1,2),main=expression(a==b==c)) I note that the following workaround executes successfully in both versions of R... plot(c(1,2),main=expression(a*"="*b*"="*c)) ...but I don't really understand

Hi. Can anyone suggest a simple way to obtain in R a list of vector coordinates of the following form? The code below is Mathematica. In[5]:= Flatten[Table[{i,j,k},{i,3},{j,4},{k,5}], 2] Out[5]= {{1,1,1},{1,1,2},{1,1,3},{1,1,4},{1,1,5},{1,2,1},{1,2,2},{1,2,3},{1 ,2,4},{1,2, 5},{1,3,1},{1,3,2},{1,3,3},{1,3,4},{1,3,5},{1,4,1},{1,4,2},{1,4,3}, {1,4,

Hello. Is it possible to set specify the font used by legend()? I would like to specify a fixed-width font so that I can line up parts of vertically stacked curve labels. For example, it would be nice if I could align the names, ages, and weights in the following three curve labels: Bob age=7 weight=100 Alexander age=13 weight=150 Susan age=20 weight=130 Is there perhaps a clever

Hello. I have an XYY series that I would like to graph with matplot() or some other single function that will do the trick. The X in question is a vector of POSIXt values obtained from strptime(). Is it possible to tell matplot() how to handle POSIXt x values? I have examined the examples at http://lark.cc.ukans.edu/~pauljohn/R/statsRus.html#5.22 , but would prefer not have to overlay the

Hi. As shown in the simplified example below, I'm having trouble allocating memory for large vectors, even though it would appear that there is more than enough memory available. That is, even with a memory limit of 1500 MB, R 2.6.1 (Win) will allocate memory for a first vector of 285 MB, but not for a second vector of the same size. Forcing garbage collection does not seem

Hi. I often wind up with many help windows cluttering my RGui screen when running Windows R 2.0.1. Is there an R instruction to close one or more help windows, or an RGui command to close all help windows? Yours truly, /Ronnen. /P.S. E-mailed CC:s of posted replies appreciated. [[alternative HTML version deleted]]

Hi. Is there an analog to abline() that can be used to plot a polynomial fit? For example, I can draw the straight-line fit fit <- lm(y ~ x) via abline(coef=fit$coef) but I'm not sure how to draw the polynomial fit fit <- lm(y ~ poly(x,2)) I do see the function curve(), but not how to prepare an expr for curve() based on the coefficients returned by the polynomial

(1) With MSVCRT.DLL (version 4.20.6164) from ftp://ftp.microsoft.com/softlib/mslfiles/msvcrt.exe installed in rw1051\bin directory, as per RWinFAQ 2.14: RGUI caused an invalid page fault in module MSVCRT.DLL at 017f:78014b90. Registers: EAX=0093007a CS=017f EIP=78014b90 EFLGS=00010246 EBX=017a2360 SS=0187 ESP=0091eee4 EBP=0091ef10 ECX=ffffffb5 DS=0187 ESI=00000000 FS=5bbf EDX=81d7d520 ES=0187

Hi. I'm having a little difficulty with R1.5.1. To wit: (1) Can anyone explain the following behavior? > a=strptime("2002-06-01 12:15:01","%Y-%m-%d %H:%M:%S") > a [1] "2002-06-01 12:15:01" > a+0 [1] "1932-04-25 21:46:45" > a+3600 [1] "1932-04-25 22:46:45" (2) Can anyone recommend an introductory R text other than that

Hello, I have read recent posts on this topic (Dr. Ronnen Levinson's Monday 02:39:55 pm), but before I install a 64 bit system, and purchase more RAM, I want to make sure I understand this interesting issue. I was attempting to plot a data frame containing Dow Jones stock information: > length(Date) [1] 19957 > length(Close) [1] 19957 > head(DowJones) Date Open High

Folks, A few days ago, I had asked a question on this mailing list about making a contour plot where a function z(x,y) is evaluated on a grid of (x,y) points, and the data structure at hand is a simple table of (x,y,z) points. As usual, R has wonderful resources (and subtle complexity) in doing this, and the gurus of the list showed me the way. Here's a complete working example. One might

Hi, one of the nice things about contour is that it labels all contour lines. contourplot only labels each particular elevation a single time. i.e., if there are two contour lines corresponding to z = 45, it will only label one of them. is there a way to get contourplot to automatically label all the contour lines, even those that are repeated? i ask because i want to plot contour plots *w/

Phil, You write: >> Using x,y, and z variables in the enclosed test.RData, the command >> >> contour(x,y,z) >> >> produces a lovely 2D plot, whereas >> >> contourplot(z~x*y) >> >> produces the error.. >> >> Error in if (lim[1] > lim[2]) stop("Improper value of limit") : >> missing value where

Hi, I have done the interpolation for my data and I was able to create the contours in multipanel with the help of Pascal. Now, I want to clip the contour with the shapefile. I want only the portion of contour to be displayed which falls inside the boundary of the shapefile. The data mydata.csv can be found on https://www.dropbox.com/s/khi7nv0160hi68p/mydata.csv The data for shapefile can be

Dear Users, I helped to install R Commander and FactoMineR to one of my collegaues. He wanted to do an MCA. Selecting three variables and using the default settings results in only one graph, the variables representation, where he gets three points for the three variables (which is totally fine). Running the code (output of the point-and-click method) EuTop100.MCA<-EuTop100[,

Hi All, I'm having difficulty making a contour plot a would like some help. A standard contourplot can be made by having an x,y, and some matrix (shape x*y) with contents: x = 1:10 y = 1:10 cont <- matrix(runif(100,min=1,max=2),nrow=10,ncol=10) filled.contour(x,y,cont) Looks very nice. Sometimes data comes in a different format however. Basically an x and an y, and a vector z with all

I have four large two-dimensional matrices of which I want to create contour plots. Something like filled.contour(<matrix>) contourplot(<matrix>) works but only gives me one plot at a time. If I combine the four matrices into one three-dimensional matrix, which I'll name "seven", there should be a way of doing something like this contourplot(seven[,,k] for k in 1 to 4)