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.
Insert regression coefficients in a graphic.
3 messages · JJCV, David Winsemius, csrabak
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
Em 22/4/2011 21:48, David Winsemius escreveu:
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.
Wouldn't round() work as well?