仅供参考。建议读一读《Applied.Econometrics.with.R》,链接是:
https://bbs.pinggu.org/thread-402964-1-1.html
祝你好运了。
- 最小二乘法
- 有三种方式可以实现最小二乘法的简单线性回归,假设数据byu
- (1)lm (byu$salary ~ byu$age + byu $ exper)
- (2)lm (salary ~ age + exper, data= byu)
- (3)attach(byu)
- lm(salary ~ age+exper)
- lm(salary~.,byu) #利用全部自变量做线性回归
- lm()只能得出回归系数,要想得到更为详尽的回归信息,应该将结果作为数据保存或者使用“拟合模型”(fitted model)
- result<-lm(salary~age+ exper + age*exper, data=byu)
- summary(result)
- myresid<-result$resid #获得残差
- vcov(result) #针对于拟合后的模型计算方差-协方差矩阵
- shapiro.test(b) #做残差的正太性检验
- qqnorm(b$res); qqline(b$res) #做残差的QQ图
- 检验回归图形:
- par(mfrow=c(2,2))
- plot(result)
- car包回归诊断函数
- 1.正态性
- qqplot
- 2.误差的独立性
- durbinWatsonTest( ) #p值不显著,说明无自相关,误差项之间独立
- 3.线性
- crplot( )
复制代码