Thank you, but I need only the area under the llm colored. What if there were two lm lines? Is it possible to color only between them. I have search and found shaded polygons, but not shades that incorporate lm lines.
Here is the updated data example. The color on this follows both the lm and the plotted lines x=seq(0,5,len=100) y=-(x-5)^2 llm<-lm(y~x) 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) On Mon, Mar 14, 2011 at 3:24 PM, Dennis Murphy <[email protected]> wrote: > 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 > <[email protected]>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 >> [email protected] >> [email protected] >> U.S. (907) 789-6603 >> >> [[alternative HTML version deleted]] >> >> ______________________________________________ >> [email protected] mailing list >> https://stat.ethz.ch/mailman/listinfo/r-help >> PLEASE do read the posting guide >> http://www.R-project.org/posting-guide.html<http://www.r-project.org/posting-guide.html> >> and provide commented, minimal, self-contained, reproducible code. >> > > -- M. Keith Cox, Ph.D. Alaska NOAA Fisheries, National Marine Fisheries Service Auke Bay Laboratories 17109 Pt. Lena Loop Rd. Juneau, AK 99801 [email protected] [email protected] U.S. (907) 789-6603 [[alternative HTML version deleted]] ______________________________________________ [email protected] 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.
