similar to: Drawing several pictures on a plot

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

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 from bottom to top. I can't find the function in text() to write

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 but failed. Thanks in advance, Marcin

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.1,102,112,112.5,115) xy<-cbind(x,y) library(cluster) UPGMA.orig<-agnes(x)

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

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 bad (non-updated) legend. Compare it with: hexbinplot(y ~ x | a, mixdata, colramp = BTC) What

Hello, I''m Alexandra Well, i have been working with cbq around 3 month ago without any problems, what i do is control the bandwith to 64Kbit for input and output and it works well. But now I don''t want control the Bandwith, i only need prioritize udp trafic over tcp trafic Can somebody to help me with my requirement? Thanks Alexandra

Hi, I'm currently using a Belkin UPS ( http://www.belkin.com/au/IWCatProductPage.process?Product_Id=137904) with nut, and my notify script doesn't receive a low battery warning. I receive ONBATT and ONLINE, but not any battery low message. I've run upsmon -D and I've seen it polls the ups.status property. The problem seems to be that this UPS doesn't put LB (I guess that

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 file. But when printing a pdf file, the background is white as it should be. What

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: "Relative efficiency for sample size n = n_i", where n should be in italics, and of course n_i

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 interested in your opinion on these books - which you would recommend. The quality of writing and

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 what's going on with that? Thanks in advance, Marcin

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 (from 0 to k, k being the number of objects). Is it at all possible? (Descriptions of both functions

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! So, maybe the limits? par(pty="s"); plot(1:10, 1:10, xlim = c(0,11),

I would like to draw a circle on a graph I am plotting in R and use the current plot's coordinate system. The most basic functionality I am looking for is the ability to plot a circle on an already existing graph by simply providing the xy coordinates for it's centre and the radius. I am also looking for a way to obtain x,y values of a point (of a circle) by specifying only the origin x,y

hello, i want to stack two lattice plots beneath each other using one x-axis and sharing the same text-panels, like: ##################################################################### library(lattice) y1 <- rnorm(100,100,10) y2 <- rnorm(100,10,1) facs<-expand.grid(Sites=rep(c("Site I","Site II"),25),Treatment=c("A","B")) pl1<-dotplot(y1 ~

Hello again, First I want to thank all the people who answered my question about line width in graphs. I promise I will learn the 'par' help page by heart for the end of the month ! I now want to trace some ellipses to emphasize groups of data. I found how to trace circles with 'symbols()', but no ellipse. I'm planning on writing my own function based on

Does anyone have any recommendations for producing dotplots with error bars? Are there packages available for this? I searched far and wide and cannot find a suitable option. I am trying to produce publication-quality figures for my thesis results. Dotplots (Cleveland dotplots) are a much better form of summarizing barchart-type data. It does not appear that any of the main plotting packages in r

Hi, I would like to use the curve() function to draw the predicted curve from an nls() object. for example: dd<-read.table("dd.txt",sep='\t',header=T,row.names=1) obj<-nls(y~c+(d-c)/(1+(x/e)^b),data=dd,start=list(b=-1, c=0, d=100, e=150)) coef(obj) b c d e -1.1416422 0.6987028 102.8613176 135.9373131

Hello everyone, I have the following code which keeps giving me an error. The code is: dat<-read.table(file="C:\\Documents and Settings\\Owner\\My Documents\\eisen.txt",header=T,row.names=1,blank.lines.skip=F,na..strings="NA") dimnames(dat)((1)) <-as.character(dat(,1)) dat<-dat(,-1) dat<-as.data.frame(dat) file.show(file="C:\\Documents and Settings\\Owner\\My