[R] Polygon question

[Uwe Ligges]

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 at 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.