Rcmdr 交互演示视频

视频地址http://v.youku.com/v_show/id_XODY0NzUwNzg4.html

HMDL文本和注释解读:

library(foreign, pos=15)
> abuse colnames(abuse) oldpar plot(fit1) #基本诊断图
> par(oldpar)
> qqPlot(fit1, simulate=TRUE, id.method="y", id.n=2) #残差qqplot检验正态性
> crPlots(fit1, span=0.5)
> library(zoo, pos=16)
> library(lmtest, pos=16)
> dwtest(relapses ~ depress + subst + stress + amt_ss + per_ss,
+ alternative="two.sided", data=abuse) #d-w test 检验 误差独立性
> outlierTest(fit1) #检验离群值
> anova(fit1, fit2) # 剔除stress变量,建立fit2回归模型。用anova()比较剔除stress后对回归有无影响。
> library(MASS, pos=18)
> stepwise(fit1, direction='backward/forward', criterion='BIC') #逐步回归和全子集回归(只介绍,不推荐)
> library(leaps, pos=19)
> plot(regsubsets(relapses ~ depress + subst + stress + amt_ss + per_ss,
+ data=abuse, nbest=1, nvmax=6), scale='bic')
> summary(fit2) # 选用fit2模型,解释回归分析结果
> summary(fit2) #对fit2中的变量标准化,比较各变量的标准回归系数

plotEx.data.frame 编写与示范~~



R作图(一) HDI Trend from 1990~2013

老师,我之前传上来的两份作图作业没有联接到Rweb上。。。所以重新提交一下。
数据来源:http://www.undp.org/content/undp/en/home/librarypage/hdr/2014-human-development-report/



R作图(二) CO2 Emission~GNP for 52 countries

这次作图有几点感受:
1. 这次数据太庞大了,虽然下载了excel格式的数据,但因为它是罗列了200多个国家的各种指标,按照首字母排序的,而且没有指明它们属于哪个region,所以一个一个分类太痛苦了。。。非常抱歉最后只呈现了一部分国家的数据。以后在数据整理上还要多学习。
2.不懂为什么Rweb上呈现的效果和我在Rstudio 做出来的不太一样。数据的颜色和位置都不太对啊,我用rainbow()设置了颜色,但是Rweb上看不出来差别。还有就是图例的排版也不太对。整体效果不如Rstudio和R上的图。

数据来源:http://www.undp.org/content/undp/en/home/librarypage/hdr/2014-human-development-report/