Hi
this might be a little bit off topic, but here it goes: lets assume I have
the following:
set.seed(13)
dat1 <- rnorm(2000, mean=10, sd=10)
dat2 <- rnorm(100, mean=10, sd=20)
d.all <- density(dat, n=1024)
d.co <- density(x[[v]], , from=min(d.all$x), to=max(d.all$x),
n=1024)
d.diff <- list(
x = d.all$x,
y = d.all$y - d.co$y
)
ylim <- range(c(d.all$y, d.co$y, d.diff$y))
plot(
d.all,
ylim = ylim
)
abline(h=0)
lines(d.co, col="red")
lines(d.diff$x, d.diff$y, col="blue")
Now I would like to identify the areas where d.all is significantly larger
then d.co and where it is significantly smaller.
What is the easiest approach to do this? At the moment I am not doing any
tests, but I am sure there is a way to determine the ranges statistically?
Thanks,
Rainer
--
Rainer M. Krug, PhD (Conservation Ecology, SUN), MSc (Conservation Biology,
UCT), Dipl. Phys. (Germany)
Centre of Excellence for Invasion Biology
Stellenbosch University
South Africa
Tel : +33 - (0)9 53 10 27 44
Cell: +33 - (0)6 85 62 59 98
Fax (F): +33 - (0)9 58 10 27 44
Fax (D): +49 - (0)3 21 21 25 22 44
email: rai...@krugs.de
Skype: RMkrug
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
