Hi:

Try this:

co <- coef(llm)
bord <- pmin(y, co[1] + co[2] * x)
plot(x,y)
polygon(c(x, x[length(x)]), c(bord, bord[1]), col = 'red')
abline(co)

HTH,
Dennis

On Mon, Mar 14, 2011 at 4:02 PM, Marlin Keith Cox <marlink...@gmail.com>wrote:

> Hi,
> I would like to add a color under a lm line and not the plotted line.  Is
> this possible?  In the example, I do not want the area under the curve red,
> but rather under the llm line.
>
> x=seq(0,5,len=100)
> y=-(x-5)^2
> llm<-lm(y~x)
> plot(x,y)
>  polygon(c(x,x[length(x)]), c(y, y[1]), col='red')
> abline(llm)
> keith
>
> --
> M. Keith Cox, Ph.D.
> Alaska NOAA Fisheries, National Marine Fisheries Service
> Auke Bay Laboratories
> 17109 Pt. Lena Loop Rd.
> Juneau, AK 99801
> keith....@noaa.gov
> marlink...@gmail.com
> U.S. (907) 789-6603
>
>        [[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.
>

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

Reply via email to