1

我有一个情节，想在情节区域添加一些回归统计数据（例如 F、R2、p）。我熟悉text()，但一直无法找到全面的信息来源，其中包含有关如何使用数学符号、子和上标等构建文本字符串的示例。非常感谢任何具有详细示例的来源。

例如，我有一个简单的线性回归，我想从中提取统计数据并将它们添加到我的绘图中。例如

reg1 <- lm(WW1 ~ PC1, data = WW_Data)

我想要像F1,69 = 14.38, p = < 0.001, R2adj = 0.16“1,69”和“adj”是下标，而“p”是斜体的东西。

编辑：

感谢您@Backlin对我的问题的出色回答。我已经扩展了一点，如果你得到一个非常重要的 p 值，代码会替换为“< 0.001”并将所有统计数据四舍五入到小数点后 2 位，除了我四舍五入到 3 的 p 值。

WW_Data <- data.frame(WW1=rnorm(10), PC1=1:10)
reg1 <- lm(WW1~PC1, WW_Data)
sreg1 <- summary(reg1)
plot(0, 0)
text(0, .2, eval(substitute(
    expression(list(F[list(fn,fd)]==fv,italic(p)==pv,R[adj]^2==R2adj)),
        list(fv = round(sreg1$fstatistic[1],2), fn = sreg1$fstatistic[2],
             fd = sreg1$fstatistic[3], pv = ifelse(sreg1$coefficients["PC1",4] < 0.001, "< 0.001",round(sreg1$coefficients["PC1",4],3)),
             R2adj = round(sreg1$adj.r.squared,2)))))
4

1 回答 1

3

我自己也为此苦苦挣扎，但实际上都是?plotmath. 你的表达方式如下，

# Fixed expression
text(x, y, expression(list(F[list(1,69)]==14.38,italic(p)<0.001,R[adj]^2==0.16)))

# Using the values of your lm
sreg1 <- summary(reg1)
text(x, y, eval(substitute(
    expression(list(F[list(fn,fd)]==fv,italic(p)==pv,R[adj]^2==R2adj)),
        list(fv = sreg1$fstatistic[1], fn = sreg1$fstatistic[2],
             fd = sreg1$fstatistic[3], pv = sreg1$coefficients["PC1",4],
             R2adj = sreg1$adj.r.squared))))

这是一个看起来像的虚拟示例。

WW_Data <- data.frame(WW1=rnorm(10), PC1=1:10)
reg1 <- lm(WW1~PC1, WW_Data)
sreg1 <- summary(reg1)
plot(0, 0)
text(0, .2, eval(...)) # The expression above

在此处输入图像描述

于 2012-08-09T08:34:23.983 回答