Hi,
How can I overlap two density plots?
A <- c(8,10,10,11,11,11,12,12,12,12,11,11,11,10,10,7)
B <- c(11,13,13,14,14,14,15,15,15,15,14,14,14,13,13,10)
plot(density(A))
plot(density(B))
Regards,
Emre
P.S. There's a overlap.density function in package DAAG
Is there another way to do this?
--
---
Emre ÜNAL
http://www.geocities.com/dusemre
[[alternative HTML version deleted]]
"Emre Unal" <emre.unal at gmail.com> wrote in news:10a006e80804250057y15067d0m49d9cca11519bb4f at mail.gmail.com:> Hi, > > How can I overlap two density plots? > > > A <- c(8,10,10,11,11,11,12,12,12,12,11,11,11,10,10,7) > B <- c(11,13,13,14,14,14,15,15,15,15,14,14,14,13,13,10) > > plot(density(A))Instead: lines(density(B)) -- David Winsemius> plot(density(B)) > > Regards, > Emre > > > P.S. There's a overlap.density function in package DAAG > Is there another way to do this? > > > >
On Fri, Apr 25, 2008 at 10:57:27AM +0300, Emre Unal wrote:> > How can I overlap two density plots? > > A <- c(8,10,10,11,11,11,12,12,12,12,11,11,11,10,10,7) > B <- c(11,13,13,14,14,14,15,15,15,15,14,14,14,13,13,10)plot(density(A)) lines(density(B)) cu Philipp -- Dr. Philipp Pagel Lehrstuhl f?r Genomorientierte Bioinformatik Technische Universit?t M?nchen Wissenschaftszentrum Weihenstephan 85350 Freising, Germany and Institut f?r Bioinformatik und Systembiologie / MIPS Helmholtz Zentrum M?nchen - Deutsches Forschungszentrum f?r Gesundheit und Umwelt Ingolst?dter Landstrasse 1 85764 Neuherberg, Germany http://mips.gsf.de/staff/pagel
Very n?ively, you could do something like this,
plot(density(A))
lines(density(B),col=2)
, and tinker your xlim and ylim as suitable. The Cairo library gives
a pretty example,
data(iris)
attach(iris)
plot(Petal.Length, rep(-0.03,length(Species)), xlim=c(1,7),
ylim=c(0,1.7), xlab="Petal.Length",
ylab="Density",
pch=21, cex=1.5, col="#00000001", main = "Iris (yet
again)",
bg=c("#ff000020","#00ff0020","#0000ff20")[unclass(Species)])
for (i in 1:3)
polygon(density(Petal.Length[unclass(Species)==i],bw=0.2),col=c("#ff000040","#00ff0040","#0000ff40")[i])