这是一个简单的回归过程:
## 输入数据
变量<-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(paste(fm[[1]][1], "=", c[1],"+",c[2],"*",fm[[1]][3])))
结果: