情节传说中的R平方和n

Question

我想一个情节图例包含R平方值，以及另一个值（数据点的数量， n ，在这种情况下，但它并不特别重要）。 我希望R格式为R ^ 2 =，其中R是斜体，2是上标。

在以下使用bquote示例中，可以通过这种方式自行添加R平方

DF <- data.frame(VAR1=rnorm(100), VAR2=rnorm(100))
with(DF, plot(VAR1, VAR2))
fit <- lm(VAR2 ~ VAR1, data=DF)
r2 <- summary(fit)$adj.r.squared
mylabel = bquote(italic(R)^2 == .(format(r2, digits = 2)))
legend("bottomright", legend = mylabel ) # display legend

添加R平方值和n可以通过以下方式实现（在此处获得@DirkEddelbuettel的信用我如何使用R在我的散点图上绘制我的R平方值？ ）

legend("topleft", legend = c( paste("R2 =", format(r2, digits = 2)),
                              paste("n = ", nrow(DF)) ) )

但是，我无法找到一种方法来修改它以便以不粘贴调用的方式格式化R ^ 2。 例如，

legend("topleft", legend = c( expression(paste(italic(R)^2, " = ", format(r2, digits = 2))),
                          paste("n = ", nrow(DF)) ) )

将成功格式化R ^ 2，还会粘贴format(r2, digits = 2)

这似乎应该是直截了当的，但我没有快乐。 任何帮助将不胜感激，谢谢。

Answer 1

这是一种方法，使用atop在几行内分割。

# your code
DF <- data.frame(VAR1=rnorm(100), VAR2=rnorm(100))
with(DF, plot(VAR1, VAR2))
fit <- lm(VAR2 ~ VAR1, data=DF)
r2 <- summary(fit)$adj.r.squared
mylabel = bquote(italic(R)^2 == .(format(r2, digits = 2)))

# the legend
legend("topright", legend=bquote(atop(italic(R)^2 == .(format(r2, digits = 2)),"n" == .(nrow(DF)))),bty="n")

请注意，我必须删除图例周围的框（ bty=n" ），因为它管理得不好。对于涉及mtext另一个解决方案，您可能会对堆叠中的其他两个帖子以及所有链接感兴趣：

图标签中的表达式和新行

表达式（）中的换行符？