29

我试图将回归的结果(即 R2)放在图表中,但似乎无法弄清楚如何从表达式中调用变量(它粘贴变量名称)。

这是我的代码。

R2Val <- signif(summary(sMod_pre90)$r.squared[1], 2)
text(92, 4, expression(paste(R^2, " = ", R2Val)), adj = 0, cex = 0.85)
4

3 回答 3

32

使用bquote(). 这是一个带有虚拟数据的示例:

set.seed(1)
DF <- data.frame(A = rnorm(100), B = rnorm(100))
mod <- lm(B ~ A, data = DF)
R2Val<-signif(summary(mod)$r.squared[1], 2)

被包裹的表达式部分.()在环境中被求值,即 的值R2Val被替换。

plot(B ~ A, data = DF)
text(1.5, 2, labels = bquote(R^2 == .(R2Val)), adj = 0, cex = 0.85)

另一个潜在的解决方案是substitute,它的工作原理类似:

plot(B ~ A, data = DF)
text(1.5, 2, labels = substitute(R^2 == A, list(A = R2Val)), adj = 0, cex = 0.85)
于 2010-12-12T18:25:31.770 回答
1

要结合substitute()and paste(),我使用以下代码;

waic <-1;
chisquare <-2;
plot(x=1:2,y=1:2, 
     main =  substitute(paste(
    "chi^2 goodness of fit ", 
       chi^2*(D/theta) , "=",  chisquare ,
       ".  WAIC =",waic),


list(chisquare=chisquare,waic=waic)  
)
)

结果如下;

在此处输入图像描述

于 2019-06-22T13:08:30.420 回答
0

我已经设法使用替代功能将它组合在一起。

R2Val <- signif(summary(sMod_pre90)$r.squared[1], 2)
text(92, 4, substitute(R^2 ~ "=" ~ R2Val), adj = 0, cex = 0.85)

一切都很好。

于 2010-12-12T18:22:13.680 回答