R help - Apr 2011 - Insert regression coefficients in a graphic.

Hello.

I did a graphic with the next code:

T=c(10,15,20,25,30)
p=c(1003,1005,1010,1011,1014)
plot(T, p, xlab=" ",, ylab=" ", pch=1, cex=2, lwd=2,
xlim=c(5,30), ylim=c(1000,1015))
abline(lm(p ~ T), lwd=2)
lm(p ~ T)

I want to see the regression coefficient in the graphic. ?How can I do that?

Thanks for your time, 

JJCV


--
View this message in context:
http://r.789695.n4.nabble.com/Insert-regression-coefficients-in-a-graphic-tp3468487p3468487.html
Sent from the R help mailing list archive at Nabble.com.

On Apr 22, 2011, at 2:16 PM, JJCV wrote:
> Hello.
>
> I did a graphic with the next code:
>
> T=c(10,15,20,25,30)
> p=c(1003,1005,1010,1011,1014)
> plot(T, p, xlab=" ",, ylab=" ", pch=1, cex=2, lwd=2,
> xlim=c(5,30), ylim=c(1000,1015))
> abline(lm(p ~ T), lwd=2)
> lm(p ~ T)
>
> I want to see the regression coefficient in the graphic. ?How can I  
> do that?
 > coef(lm(p ~ T)
+ )
(Intercept)           T
      997.40        0.56
 > text(10,1010, labels=paste("Slope =", coef(lm(p ~ T)  )[2]) )

I suppose you could format it with formatC() or sprintf() if you  
wanted it to fewer than 10 decimal places.


-- 
David Winsemius, MD
West Hartford, CT