情节传说中的R平方和n

时间:2016-06-18 08:37:29

标签: r plot legend

我希望绘图图例既包含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 可以通过以下方式实现(在此处{@ 3}}归功于@DirkEddelbuettel)

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)

这似乎应该是直截了当的,但我没有快乐。非常感谢任何帮助,谢谢。

1 个答案:

答案 0 :(得分:0)

这是一种方法,使用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的另一个解决方案,您可能会对堆栈中的其他两个帖子及其中的所有链接感兴趣:

Expression and new line in plot labels

Line break in expression()?