• R语言
  • 线性模型在R中的实现(1)

觉得自己在这个论坛得到不少帮助,自己似乎没有啥贡献,有些惭愧。



最近自己在学习R,主要是线性模型方面的,学了点皮毛,在这里把自己学到的一点点东西在这里向大家报告一下,或许对象我这样的初学者有些帮助,不对的地方也好请高手们指正:



//*从excel表中读取数据,先将数据文件存成csv文件,再用read.csv去读取*/

data<-read.csv(file="e:\\work\\2009下半年\\modern statistics-2009new\\ding learning note\\linear model\\例2.2.csv")



/*将data的第二列数值赋给变量x,第三列赋给变量y*/

x<-data[,2]

y<-data[,3]



/*做散点图*/

plot(x,y)



/*线性回归用的是函数lm()*/

z<-lm(y~x)

summary(z)



/*以下代码是没有截距的回归*/

z<-lm(y~x-1)

z<-lm(y~0+x)





library(MASS)

z<-lm(y~x)

summary(z)

plot(x,y)

plot(fitted(z),studres(z))/*残差分析,拟合值与学生化残差图*/

qqnorm(studres(z))/*对残差进行正态性检验*/

qqline(studres(z))



这里调用library(mass)我还是不太明白,估计是将回归的一些结果拿出来用吧,因为之前我没有这个的时候,做不了plot(fitted(z),studres(z))



这是一元线性回归的一点点内容,多元的我正在学,有机会在向大家报告
library 是把某个包载入,以便使用其中提供的函数或数据。
你不是用require(MASS)的时候,你试试?studres

R默认加载的包里是没有这个命令的