在以往R语言的教材、帮助示例和各类读书笔记中,有关图形中标注公式的介绍只有指定符号和数字的介绍。如何像SPSS那样在计算结束后将带参数的模型计算公式一并绘制在图形中呢?其实提取模型参数很简单,关键是如何将模型变量名和公式提取出来并与参数重组公式,下面我给出一个简单的小例子说明这个过程。
[code lang="R"]
## 输入数据
变量<-c(0.1,0.11,0.12,0.13,0.14,0.15,0.16,0.17,0.18,0.20,0.21,0.23)
GDP<-c(42,43.5,45,45.5,45,47.5,49,53,50,55,55,60)
## 回归分析
lm.sol<-lm(GDP ~ 变量)
summary(lm.sol)
## 提取参数
c <- round(coef(lm.sol),3)
# 提取公式变量名
fm <- strsplit(deparse(formula(lm.sol))," " )
## 绘图
plot(lm.sol$model ,main = "回归结果")
lines(lm.sol$fitted.values,lm.sol$model[,2])
# 在图上标注回归公式
text(51 ,0.13,as.expression_r(paste(fm[[1]][1], "=", c[1],"+",c[2],"*",fm[[1]][3])))
[/code]