在R语言绘图中标注带参数的公式

  • A+
所属分类:Script

在以往R语言的教材、帮助示例和各类读书笔记中,有关图形中标注公式的介绍只有指定符号和数字的介绍。如何像SPSS那样在计算结束后将带参数的模型计算公式一并绘制在图形中呢?其实提取模型参数很简单,关键是如何将模型变量名和公式提取出来并与参数重组公式,下面我给出一个简单的小例子说明这个过程。

   这是一个简单的回归过程:

## 输入数据
变量<-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])))

R语言线性回归

avatar

发表评论

:?: :razz: :sad: :evil: :!: :smile: :oops: :grin: :eek: :shock: :???: :cool: :lol: :mad: :twisted: :roll: :wink: :idea: :arrow: :neutral: :cry: :mrgreen: