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.