首先,既然你要画残差(residual)的qqplot 那一定要先做回归。
你没说明回归自变量和因变量都是什么,假设第二列是因变量pressure[2],第三列是自变量pressure[3]
1.回归
>lm.fit<-lm(pressure[2]~1+pressure[3])
2.做qq图
这里不用qqplot这个指令,用qqnorm和qqline就好
>qqnorm(lm.fit$res)
>qqliine(lm.fit$res)
提取?就直接write.table()
比如你要第二列和第三列
>write.table(pressure[c(2,3)],file="c:\test.txt",quote=F,row.names=F,sep="\t")
第一个参数pressure[c(2,3)]是你要提取的变量
第二个参数file= 是你要存放的地址和文件名
第三个参数quote=F 是你打出来的数字都没有引号(默认是字符型)
第四个参数row.names=F 是不要行名称(为什么不要?自己试试就知道了)
第五个参数sep="\t" 每列数据之间空一个tab的距离,也就是8个字节,清晰明了。
追问回归什么的没学到。。
书就给个公式,residuals<- with(pressure,pressure-(0.168+0.007*temperature)^(20/3))
然后要用这些residuals建一个正常的QQplot,然后看分布。
那应该怎样写?
追答>qqnorm(residuals)
>qqline(residuals)
就好了,两幅图会一起显示,这就是qq图
本回答被提问者采纳