Hi, all!

 

I'm pretty sure I'm missing something about this.

Is there a smart way of typping hat(R)^2 and it's value from a linear
regression?

 

I've just found this tricky one:

 

# Sample data

x <- sample(1:100,10)

y <- 2+3*x+rnorm(10)

 

# Run the regression

lm1 <- lm(y~x)

 

# Plotting

plot(x,y, main="Linear Regression", col="red")

abline(lm1, col="blue")

 

placex <- par("usr")[1]+.1*(par("usr")[2]-par("usr")[1])

placey1 <- par("usr")[3]+.9*(par("usr")[4]-par("usr")[3])

placey2 <- par("usr")[3]+.8*(par("usr")[4]-par("usr")[3])

 

# HERE: Is this the right way?

text(x=placex, y=placey1,

     bquote(R^2 == .(summary(lm1)$r.squared)), adj=c(0,0))

text(x=placex, y=placey2,

     bquote(R^2 == .(summary(lm1)$adj.r.squared)), adj=c(0,0))

text(x=placex,y=placey2,

     expression(hat(R)), adj=c(0,0))

 

In addition, when I save the plot as PDF, the expression hat(R) seems to be
somewhat displaced.

 

Any advice?

 

Thanks in advance,

 

Vicente

 

Aviso de confidencialidad

Este correo electrónico y, en su caso, cualquier fichero anexo, contiene
información confidencial exclusivamente dirigida a su(s) destinatario(s).
Toda copia o divulgación deberá ser autorizada por la Subdirección de
Análisis y Control Interno. Si ha recibido este mensaje por error, le
rogamos que nos lo comunique inmediatamente por esta misma vía y proceda a
su eliminación.
Antes de imprimir este mensaje, asegúrese de que es necesario. El medio
ambiente está en nuestras manos.

 


	[[alternative HTML version deleted]]

