[R] auto.key and simpleTheme

Dear all, parallel(~iris[1:4], groups = Species, iris, par.settings = simpleTheme(lwd = c(1,3,1), lty = c(1,1,2), col.line = 1), auto.key = T) Despite the use of par.settings and simpleTheme, the lines in the key and graph are not the same. Any suggestions why? Regards, Marcin [[alterna

[R] Double expressions in plot labels

Hi all, I wanted to do something of this sort: library(lattice) label <- expression(alpha) xx <- substitute(expression(lab %+-% type), list(lab = label, type = " type")) stripplot(Species ~ jitter(Petal.Length), iris, ylab = xx, horizontal=T) (The example plot itself does not make much sense, bu

Re: [R] Plot's aspect ratio and pty

gt; > > On 13.12.2010 07:30, Marcin Kozak wrote: >> >> Dear All, >> >> I've been playing with pty, and it seems it does not produce square >> plots as it is expected to (or at least as I expect it to). Consider >> this simple example: >> >> p

[R] Plot's aspect ratio and pty

Dear All, I've been playing with pty, and it seems it does not produce square plots as it is expected to (or at least as I expect it to). Consider this simple example: par(pty="s"); plot(1:10, 1:10) This should produce a square plot, right? Well, if you have a look at the graph, it is not square

[R] The width argument of stem()

Dear all, The help page of stem() says that the width argument gives the desired width of plot. However, I cannot come up with any idea what this width stays for. This becomes especially difficult if you run these couple of examples: The decimal point is 3 digit(s) to the right of the | 0 | 0

[R] Updating a hexbinplot

Dear all, Considering this simple example of hexbinplot: mixdata <- data.frame(x = c(rnorm(5000), rnorm(5000,4,1.5)), y = c(rnorm(5000), rnorm(5000,2,3)), a = gl(2, 5000)) fig <- hexbinplot(y ~ x | a, mixdata) print(fig) update(fig, colramp = BTC) produces a b

[R] Setting base level for contrasts with lme

Hi all, Note: lm(Yield ~ Block + C(Variety, base = 2), Alfalfa) equals i <- 2; lm(Yield ~ Block + C(Variety, base = i), Alfalfa) However, lme(Yield ~ C(Variety, base = 2), Alfalfa, random=~1|Block) which is fine, does not equal i <- 2; lme(Yield ~ C(Variety, base = i), Alfalfa, random=~1|Block)

[R] Comparing regression lines

Hi, How would you compare several (simple) regression lines in R? Have you heard of any implementation of the Johnson-Neyman procedure or anything of this sort? Any suggestions? Best, Marcin __ R-help@r-project.org mailing list https://stat.ethz.ch/mai

[R] xyplot: no ticks for a factor scale?

Hi, I would like to have no ticks on a scale that represents a factor. The tick.number argument from scales does not work in such a situation, as the help page as well as this simple (fairly stupid) code show: require(lattice) fact<-gl(4,1,labels=LETTERS[1:4]) y<-c(1,4,3,2) xyplot(y~fact,scales=l

[R] Truncated labels with mosaicplot

How can I deal with truncated labels in the mosaicplot()? Look at the example: mosaicplot(~ gear + carb, data = mtcars, color = TRUE) Look at "carb" (the number of carburetors): the label "8" is truncated. How this might be dealt with? Thanks in advance, Marcin __

[R] Loess fitting with bisquare

Hi, It seems there is no straightforward way to carry out in R the loess fitting with bisquare, as given by William Cleveland in his "Visualizing Data". Am I right? Thanks in advance, Marcin -- "Build up your weaknesses until they become your strong points" -- Knute Rockne _

[R] Dendrogram with the UPGMA method

Hi, I am clustering objects using the agnes() function and the UPGMA clustering method (function = "average"). Everything works well, but apparently something is wrong with the dendrogram. For example: x<-c(102,102.1,112.5,113,100.3,108.2,101.1,104,105.5,106.3) y<-c(110,111,110.2,112.1,119.5,122.

Re: [R] A ternary graph's scales

Gregory (Greg) L. Snow Ph.D. > Statistical Data Center > Intermountain Healthcare > [EMAIL PROTECTED] > (801) 408-8111 > > > >> -Original Message- >> From: [EMAIL PROTECTED] >> [mailto:[EMAIL PROTECTED] On Behalf Of Marcin Kozak >> Sent: Thursday, S

[R] A ternary graph's scales

Hi all, I am drawing a ternary graph. Everything is fine with both ternaryplot (package vcd) and triangle.plot (package ade4), but I want to present scales in neither percents nor from 0 to 1 (this is actually the only option I found in both functions). I want the scales to be in a natural scale (

Re: [R] Point size problem

Sorry. I'm using Windows XP and R 2.7.0, and the same problem occurred for various graphics devices (windows, pdf, jpeg). Marcin On Thu, Jun 26, 2008 at 11:54 AM, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > On Thu, 26 Jun 2008, Marcin Kozak wrote: > >> Hi all, >>

[R] Point size problem

Hi all, I have a simple question and couldn't find any post on this. When plotting simple scatterplots (other plots as well), e.g., x<-rnorm(30, 10, 1) y<-rnorm(30, 10, 1) plot(x, y, pch = 15, cex = 1), the points, even those close to each other, may have visibly different sizes. Do you know wha

[R] Book on R graphics

Hi all, I am looking for a book from which one could learn a wide range of graphics in R, from the very beginning topics to those advanced (though not necessarily concerned with a particular method or topic). I don't look for a list of such books since this is on R's web page, I am rather interest

[R] Italics in plot main title

Hi, I am drawing several plots and want to have italics in a main title; this is easy with expression(). However, I want also to add a value to it, say n_i, that depends on an ith plot. For this I am using paste(). An example: n_i = 10, 20, 30; I want to draw a plot for each i with the title: "Rel

[R] EM algorithm for multiple-locus haplotypes frequencies

Hi all, I've been looking in R for an EM algorithm adjusted for multiple-locus haplotypes frequencies, but failed in 100%. Has anyone heard of anything of this kind in R? Thanks in advance, Marcin __ R-help@r-project.org mailing list https://stat.ethz.

[R] Labels in plot()

Hi, Look here: xy<-matrix(rnorm(100,100,10),50) plot(xy,ann=F,labels=F) Even though the error is reported, the "labels" do what I want them to. Any reason for this? Thanks, Marcin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/list

[R] Drawing several pictures on a plot

Hi, I want to draw several pictures on one plot and do as follows: xy<-matrix(rnorm(100,100,10),50) op <- par(mfrow = c(2, 2), pty = "s") plot(xy,ann=F);plot(xy,ann=F);plot(xy,ann=F);plot(xy,ann=F) par(op) What I need now is to set the size of the pictures within this plot since there is too muc

[R] Controling width of boxes in boxplots

Hi, I want to add boxplots to a scatterplot: plot(x,y, xlim=c(80,120),ylim=c(80,120)) boxplot(y,add=TRUE,at=118) boxplot(x,add=TRUE,at=118,horizontal=TRUE) How can I control the width of the boxes (say, I'd like them to be of width 3 in the variables' scales). I've tried the "width" parameter bu

Re: [R] Drawing a loess line

rs, geom = c("point", "smooth")) > > or with a colored line: > >qplot(speed, dist, data = cars) + geom_smooth(colour = "red") > > > > On Feb 3, 2008 3:29 PM, Marcin Kozak <[EMAIL PROTECTED]> wrote: > > Dear all, > > > &

Re: [R] Drawing a loess line

lines(predict(cars.lo), col="blue") > > > > > On 03/02/2008, Marcin Kozak <[EMAIL PROTECTED]> wrote: > > Dear all, > > > > To draw a lowess line on a plot was a piece of cake; to draw a loess > > line, however, seems not that easy. Is the loes

[R] Drawing a loess line

Dear all, To draw a lowess line on a plot was a piece of cake; to draw a loess line, however, seems not that easy. Is the loess plotting implemented at all in relation to the loess function, or do I have to look in add-on packages? Thanks, Marcin __ R-

[R] Vertical text in a plot

Hi, Consider this simple plot: > plot(1:25,runif(25,0,1),ylab="First Y-axis label",xaxt="n") I want to add an additional axis as > axis(4,at=seq(0.2,1,.2), labels=1:5) I have no idea how to add now the title of the new axis as "Second Y-axis label". I want this text to be vertically directed fro

[R] Printing to jpeg (dev.print)

Hi, This may be an easy question, but let me ask it. When writing a plot to a jpeg file: > plot(runif(30)) > dev.print(file="test.jpeg", device=jpeg, width=600) the plot I receive has gray background, the result on no account I want to receive. The same situation occurs when printing to a bmp fi