[R] color under lm line

[Eik Vettorazzi]

Hi Keith,
how about this:

fit<-function(x,coefs) x*coefs[2]+coefs[1]
plot(1:10,1:10,type="n")
x<-par("usr")[1:2]
cf1<-c(a=3,b=.5) #eg coefs from a lm-object
cf2<-c(a=1,b=.9) #another line, coloring only intersection
polygon(c(x,rev(x)),c(fit(x,cf1),fit(rev(x),cf2)),col="red")

setting one of the coef vectors to c(0,0) colors the area under the
other line.

hth.

Am 15.03.2011 00:02, schrieb Marlin Keith Cox:
> 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
> 

-- 
Eik Vettorazzi
Institut für Medizinische Biometrie und Epidemiologie
Universitätsklinikum Hamburg-Eppendorf

Martinistr. 52
20246 Hamburg

T ++49/40/7410-58243
F ++49/40/7410-57790