similar to: Superposed histograms

Why does > curve(function(y) y^2, from=0,to=1) not work, whereas > myf <- function(y) y^2 > curve(myf, from=0,to=1) work? For the former, I get the error message Error in curve(function(y) y^2, from = 0, to = 1) : 'expr' must be a function or an expression containing 'x' I'm using R1.7.0 under Windows XP. Damon Wischik.

I want to plot biplots. I have seen the function biplot, but there are some extra features I would like, that I do not know how to achieve. 1. My observations, and my variables, fall into groups. Is there a way to, say, plot the observations in several different colours, according to which group the observation falls into? Similarly with the variables? I am used to the lattice idiom, in which one

When I alter the levels of a factor, why does it alter the names too? f <- factor(c(A="one",B="two",C="one",D="one",E="three"), levels=c("one","two","three")) names(f) -- gives [1] "A" "B" "C" "D" "E" levels(f) <-

I am having trouble with ESS+R. I don't know if it is an ESS problem or an R problem, so I'm posting to this mailing list in the first instance. I am using R 1.6.2, Windows XP (latest updates installed), XEmacs 21.4, and ESS 5.1.21. I have experienced the same problem under Windows 2000 and recent versions of Emacs and ESS (though my current Windows XP installation is from scratch, and

There have been a number of posts to this list by people having trouble with the HTML search engine. Often these troubles are caused by incorrect setups (user hasn't installed Java properly, or Java is disabled, or Javascript is disabled). Sometimes the trouble persists even when Java is installed properly. I have written an alternative HTML search engine, which is based on Javascript rather

I would be grateful for advice about the following problem. It's not directly R-related, but I'm hoping that R will help me analyse the following data. I have a table which indicates the progression of a certain age-related disease. At a certain point in time, a population was sampled; and I have measurements for the age of each individual, and their disease stage. (Disease stage is an

----- Original Message ----- From: "Richard A. O'Keefe" <ok at cs.otago.ac.nz> To: <paul at datavore.com> Sent: Thursday, September 04, 2003 2:56 AM Subject: Re: [R] Overlaying graphs > I do not know how to overlay the curve graphic on top of hist graphic. > > Do you know about the "add=TRUE" option for plot()? > > I am hoping to show visually

I am having trouble with optim. It claims to have converged to a minimum, yet it has in the course of the optimization visited many points which are closer to optimal. I would be grateful for any explanation of this behaviour. I'm trying to estimate the parameters in the model X ~ Binomial(1,p) * NegBin(mu,theta). So I define a log likelihood function, and invoke optim thus: o <- optim

Hi, I am having some trouble trying to plot a histogram in more than one colour. What I want to do is, plot two vectors in the same histogram, but with different colours, for instance: > x <- rnorm(1000,20,4); > y <- rnorm(1000,10,2); Then I'd like to have x and y ploted on the same hist (I can do that already doing w <- c(x,y) then hist(w)) but the bars

The last two times I have originated message threads on R or Bioconductor I have received the message included below from someone named Patrick Connolly. Both times I was the originator of the message thread and used what I thought was a unique subject line that explained as best I could what my question was. Patrick seems to be implying that I am abusing the R and BioC help newsgroups in this

Een ingesloten tekst met niet-gespecificeerde tekenset is van het bericht gescrubt ... Naam: niet beschikbaar Url: https://stat.ethz.ch/pipermail/r-help/attachments/20080129/23f6cbc7/attachment.asc

I could use some advice regarding xyplot. I've got 2 time series. Both cover approximately the same period of time (ie, 1940 to 2009). But one series has annual data and the other has monthly data. One refers to university enrollment; the other to unemployment rates. Both are currently in the same data frame. I'd like to use the monthly times series as a light grayscale background for a

I could use some advice. I've got 2 time series. Both cover approximately the same period of time (ie, 1940 to 2009). But one series has annual data and the other has monthly data. One refers to university enrollment; the other to unemployment rates. Both are currently in the same data frame. I'd like to use the monthly times series as a light grayscale background for a plot of the

In the hypothetical example below, how do I add two segments() into the two panels, respectively? Say segments(x0=5, y0=10, x1=5, y1=20) on the left and segments(x0=15, y0=-10, x1=15, y1=-2) on the right? Many thanks in advance, Yuelin Li. ps. part of the code came from a solution given by Deepayan Sarkar. ------------------- library(lattice) set.seed(12345) x <- 0:20 y.male.obs <- - 1.2

Hi, I'm trying to recreate a sensitivity-specificity graph using the histbackback function. The only problem is that these graphs are typically drawn with vertical rather than horizontal bar plots (and the histbackback function only seems to work with horiz=TRUE argument, using "horiz=FALSE" doesn't work). Does anyone know if: 1) there's a different graphing function that

I would want to know if it's possible to plot two histograms on the same graphic in order to compare the bins one to one, like a Excel graphic. If it is possible can you help me to do it, because I don't know the exact R-commands. Thanks for your help Damiano > > >Martin Maechler <maechler at stat.math.ethz.ch> http://stat.ethz.ch/~maechler/ >Seminar fuer Statistik,

How do you label Xtick and Ytick marks with Currency symbols: $2000 instead of 2000? I would like to add dollar symbols to tickmarks on boxplots, histograms and back-to-back histograms. My Examples (requiring the lattice and Hmisc packages): data(case0102, package="Sleuth2") str(case0102) boxplot(Salary~Sex, case0102) histogram(~ Salary | Sex, data=case0102) require(Hmisc) #

I have some data which I am trying to fit with a negative binomial distribution. I have found the glm.nb function from MASS. I have reason to believe that the mean parameter mu depends on certain factors, and that the shape parameter theta depends on others. If, say, the factors are P and Q, it might be that mu ~ P:Q and theta ~ P (where mu ~ P:Q means that mu is a function of the pair (P,Q))

I am making frequency histograms using the histbackback function on my 2 datasets. However when I try to use the brks function: foo<-histbackback(log(fie11), log(fie86),ylim=c(0,9),probability=FALSE,axes=TRUE,ylab=("log10 Parcel Size"),brks=16) The graphic results in a 'NA' label for the y axis (no intervals are returned) Also when I use 'summary(foo)' the

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