2
我想在使用ggplot2創建的圖表上打印線性迴歸的F統計量和df。我已經非常接近了,但我似乎無法在自由度(1和8)之間出現逗號。在R中的等式中包含逗號
#generate some data
x <- 1:10
y <- 1:10 + rnorm(10, 0.1, 0.2)
#run linear regression
mod1 <- lm(y~x)
#get r2 and df and round for better printing
Fstat <- round(summary(mod1)$fstatistic, 3)[1]
df <- round(summary(mod1)$fstatistic, 1)[2:3]
#write out the equation to be printed to plot
eq <- paste("F[",df[1],"]","[,]","[",df[2],"]==",Fstat)
#plot the data and label with the equation
qplot(x, y) + annotate("text", x=2, y=8, label=eq, parse=TRUE)