similar to: R-beta: contour labels

> From: Bill Simpson <wsimpson at uwinnipeg.ca> > Subject: Re: contour labelling [was "Re: [R] image ()"] > MIME-Version: 1.0 > > > We may be missing the crux here: was it to have a means of > > displaying such surfaces? > > Yes, bang on. I knew of the > existence of the Delaunay code in the R libraries. Yes I want > to make a picture that

I include the argument pty="s" to image, but still when I make an image by the method > postscript("rstuff/test.ps") > tauseq<-seq(0,1,.5) > cif2d.image(x,,y,tauseq) > dev.off() the image plot created is NOT square. I thought pty="s" would make it square. Generally it would be nice if the image saved to disk were like the one we see while in R

I have a contour plot where the surface can take on positive and negative values. I would like to use solid contour lines for positive values and dotted contour lines for negative values. I saw a plot like this produced by S-Plus before. Is there a way to do this with R? If not, can anyone suggest a way to make interpretation of the contour plot easier? (Right now, it is not clear if a set of

From: Martin Maechler <maechler at stat.math.ethz.ch> Subject: Re: R-beta: Contour package There's still no "persp" (perspective plots), even though Ross once told about a ``prototype'' hidden somewhere in his private code. Interesting coincidence... Yesterday I pulled out some very old pascal code I wrote to do perspective plots of surfaces. I hope to have a

> From: Martin Maechler <maechler at stat.math.ethz.ch> > To: royle at bearmtn-e0.cgd.ucar.edu > CC: lgygax at access.unizh.ch, r-help at stat.math.ethz.ch > Subject: contour labelling [was "Re: [R] image ()"] > Mime-Version: 1.0 (generated by tm-edit 7.106) > > >>>>> "Jeffrey" == Jeffrey A Royle <royle at bearmtn-e0.cgd.ucar.edu>

Questions on image: 1) How can I put labels on the x and y axes? 2) How to tell it to use e.g. 32 grey levels (not some colour map)? 3) How to know what the legend is (i.e. each grey level = what z value)? Thanks very much for any help. BTW I was wondering if persp was on the To Do list. That would be great! Bill Simpson

I have a question on the use of return(). (Nothing on it in the docs I have) The test code below gives the error: Error: Object "x" not found when I do: thingy2(). How should it be fixed? Thanks very much for any help! (My original solution to this sort of problem was to use global variables x<<-... y<<-...) Bill Simpson ----------------------------- thingy<-function(k)

Can someone please tell me to make a density histogram? hist makes one with count or relative frequency on the y-axis. I want the density, which is (rel freq)/(bin width) In the help I see: intensities: values f^(x[i]), as estimated density values. If `all(diff(breaks) == 1)', they are the relative frequencies `counts/n' and in general satisfy

Thanks very much Douglas for the pointer to nlm. Maybe the "Notes on R" maintainer can add at least a mention of nlm in the section on nonlinear fitting? I never did nonlinear fitting in S-Plus before, so I have nothing to unlearn, but I was hoping someone could show me how to do a least squares fit with nlm. example: x<-c(1,2,3,4,5,6) y<-.3*x^-.6 +.2 y<-y+rnorm(6,0,.01)

Thanks very much Douglas for the pointer to nlm. Maybe the "Notes on R" maintainer can add at least a mention of nlm in the section on nonlinear fitting? I never did nonlinear fitting in S-Plus before, so I have nothing to unlearn, but I was hoping someone could show me how to do a least squares fit with nlm. example: x<-c(1,2,3,4,5,6) y<-.3*x^-.6 +.2 y<-y+rnorm(6,0,.01)

I have data saved in files on disk in following format x y z ------------------ 1 1 1 1 2 2 2 1 3 2 2 4 Normally I read the file into a data frame, but suppose I do it this way: x<-seq(1,2) y<-seq(1,2) z<-seq(1,4) data<-data.frame(x,y,z) image(unique(data$x),unique(data$y),data$z, col=gray(0:19/19), xlab="x", ylab="y") I was helped at some point in the past on

Greetings, All! I want to draw contour lines in red, using contour(), but also have the contour labels (for the level-values) in black so that they will stand out against a coloured background already generated using filled.contour() (the background shades from green at low levels of "risk" to red at high levels). In any case, contour labels in red are already somewhat inconspicuous

Hi, I want to display some contour labels. It works well, but the thickness of the label is to low. I tried the "labcex" command contour(x,y,z, labcex=2) But only the size is rising and not the thickness of labels. Only if i set the value to "labcex=10" the thickness is good, but the size is to big, so i can't see anything ;) Is there any command to rise the thickness

After some paper clipping problems I checked options() and saw that the default paper size was a4; I need US Letter. So as per instructions I uncommented the R_PAPERSIZE line in config.site, R_PAPERSIZE=letter Restarting R, everything is still the same though. I guess I have to reinstall? That seems awkward. Is there a way a prompt ("a4 or letter paper?") can be inserted in the

Brian D Ripley writes: Interactive rotation of perspective plots would be a great boon, as would shading options for perspective plots, including contouring on such plots. I am envisioning a division into static and dynamic plots. The static version will concentrate on rendering as nice a picture as possible - anti-aliasing (if supported) arbitrary text rotation etc, whereas the dynamic

R-0.99.0a now plots contour labels from contour(), but the values in levels seem to be coerced to integer. Reproduce by: > set.seed(2) > contour(matrix(runif(36),nrow=6), labcex=1.2) and compare with: > set.seed(2) > contour(matrix(10*runif(36),nrow=6), labcex=1.2) Ray Brownrigg -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list --

version .62: --------------------------------------------- > ?qpois The Poisson Distribution dpois(x, lambda) ppois(q, lambda) qpois(p, lambda) rpois(n, lambda) Arguments: x: vector of (positive) quantiles. p: vector of probabilities. n: number of random values to return. lambda: vector of positive means.

I am running R-0.61.1 on linux under X. If I kill the R Graphics window (click on X box in upper right corner), then subsequently do x11(), R crashes. I find that things only work the right way if I close the graphics window using dev.off(). I was wondering if maybe it could be arranged that closing the graphics window via click would seem to R to be equivalent to typing dev.off(). BTW thanks

This surprised me. Is it the way it is supposed to be? > x<-c(1,2,3,4,5) > y<-c(1.2,1.3,1.4,1.5,NA) > x-y [1] -0.2 0.7 1.6 2.5 NA > (x-y)^2 [1] 0.04 0.49 2.56 6.25 NaN >>>so NA^2 = NaN? Why not still NA? > sum((x-y)^2) [1] NaN >>>yes that is reasonable So if you ever have a data set with missing observations (NAs), you can't do any nlm() least

If I do fit<-lm(y~x) Is it possible to access the SE of the slope? (Analogous to getting the slope like this: fit$coef[2]) If not, it would be handy. (I want SE of 1/slope, and an approx way is fit$se[2]/(fit$coef[2]^2)) Thanks for any help. Bill Simpson -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read