R语言学习 day_8

• R语言实战书中有些代码比较旧，包已经更新但是代码没更。今天花了好多时间搞明白，qqPlot交互图。会使用help（)函数太重要了。

• 回归诊断

  回归诊断主要是对OLS回归的统计假设

  正态性、独立性、线性、同方差性

  1. 标准方法

     在R中可以使用plot()函数和par()函数，生成评价模型拟合情况的四幅图形

     #回归诊断，标准方法
     fit<-lm(weight~height,data=women)
     par(mfrow=c(2,2))
     plot(fit)
     
     

     

      

  2. 改进的方法

     car包

     正态性

     #改进的方法
     library(car)
     states<-as.data.frame(state.x77[,c("Murder","Population","Illiteracy","Income","Frost")])
     fit<-lm(Murder~.,data=states)
     qqPlot(fit,main="Q-Q Plot",id=list(method="identify",labels=row.names(states)),simulate=TRUE,pch=16)
     #car包已经更新原来的id.method="identify"已经不可用了，变成了id=list()
     #这里出现的“警告: 没有点在0.25英尺内”是Rstudio的问题，放入R本来的环境运行不存在该问题。建议进入原环境运行