Well you need to recalculate the x values and need to interpolate for
the position where you lines cross the m lines ....
Uwe Ligges
On 16.06.2011 23:35, Muhammad Rahiz wrote:
Hi all,
I have the following script which fills the values which are less than
the mean of a given timeseries.
If you look closely, the colored regions are "out of line".
Any suggestions how I can rectify this?
Thanks
Muhammad
# -----
#rm(list=ls())
x <- abs(rnorm(100))
tt <- 1:100
m <- mean(x)
w <- which(x>=m)
x1 <- x ; x2 <- x ; x3 <- x
x1[w] <- m
x2[1:length(x)] <- m
tx <- c(tt,rev(tt)) ; ty <- c(x1,rev(x2))
par(mfrow=c(2,1))
yy <- c(0,3) # y-limit
plot(tt,x,type="l",ylim=yy)
abline(h=m)
plot(tx,ty,type="n",ylim=yy)
polygon(tx,ty,col="red")
lines(x)
# -----
______________________________________________
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.
______________________________________________
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.
