QQ图实例(R软件).doc

下载文档 降价啦

14
0
约1.49千字
约 10页
2017-04-24 发布于重庆
举报
版权申诉
保障服务

QQ图实例(R软件).doc

1、本文档共10页，可阅读全部内容。
2、有哪些信誉好的足球投注网站（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

QQ图实例(R软件)

rm(list=ls()) DA=read.csv(file=file.choose(),head=T) lm.reg=lm(DA$Y~X1+X2+X3+X4+X5+X6+X7+X8+X9+X10+X11+X12+X13+X14+X15+X16+X17+X18+X19+X20+X21+X22+X23+X24+X25+X26+X27+X28+X29+X30+X31+X32+X33+X34+X35+X36+X37+X38+X39+X40+X41+X42+X43+X44+X45+X46+X47+X48+X49+X50,data=DA) summary(lm.reg) plot(lm.reg,2) XX-cor(DA[1:50]) kappa(XX,exact=TRUE) eigen(XX) p=1 n=nrow(DA) d1=dffits(lm.reg) Y=1:n Y[abs(d1)2*sqrt((p+1)/n)] infl=lm.influence(lm.reg,do.coef=FALSE) D=cooks.distance(lm.reg,infl=lm.influence(lm.reg,do.coef=FALSE), + rs=weighted.residuals(lm.reg), + sd=sqrt(deviance(lm.reg)/df.residual(lm.reg)), + hat=infl$hat) sort(D) D=abs(1-covratio(lm.reg, infl=lm.influence(lm.reg, do.coef = FALSE),res = + weighted.residuals(lm.reg))) sort(D) Reg_Diag-function(fm){ + n-nrow(fm$model); df-fm$df.residual + p-n-df-1; s-rep( , n); + res-residuals(fm); s1-s; s1[abs(res)==max(abs(res))]-* + sta-rstandard(fm); s2-s; s2[abs(sta)2]-* + stu-rstudent(fm); s3-s; s3[abs(sta)2]-* + h-hatvalues(fm); s4-s; s4[h2*(p+1)/n]-* + d-dffits(fm); s5-s; s5[abs(d)2*sqrt((p+1)/n)]-* + c-cooks.distance(fm); s6-s; s6[c==max(c)]-* + co-covratio(fm); abs_co-abs(co-1) + s7-s; s7[abs_co==max(abs_co)]-* + data.frame(residual=res, s1, standard=sta, s2, + student=stu, s3, hat_matrix=h, s4, + DFFITS=d, s5,cooks_distance=c, s6, + COVRATIO=co, s7)} Reg_Diag(lm.reg)