R语言讲义(包括各种回归)说课材料.ppt

R 语言讲义;免费(没有权力和铜臭) 资源公开, 可改变代码(不是黑盒子,也不是吝啬鬼, 透明是防止“腐败”的最好方式) 容易学习。可编程以实行复杂的课题 可扩展: 通过数千个网上提供的适用于不同领域、不同目的、不同方法的软件包来实现你的目标。也可以把你的方法贡献出来 功能强大(绘图功能, 优秀的内在帮助系统, R社区的支持，不断更新，不断修正) 没有任何一个商业软件有如此多和如此新的算法;点击CRAN得到一批镜像网站;点击镜像网站比如Berkeley;;;;Packages （每个都有大量数据和可以读写修改的函数/程序）;Packages （继续）;Packages （网上）;所有这些Packages可以自由下载;安装Packages;;;关机时是否保存?;几个有用的函数;赋值和运算 ;赋值和运算 ;从高到低的运算次序;一些基本运算例子;一些基本运算例子;一些基本运算例子;一些基本运算例子;一些基本运算例子;一些基本运算例子;一些基本运算例子;一些基本运算例子;一些基本运算例子;分布和产生随机数 ;a=factor(letters[1:10]) a[3]=w#不行 a=as.character(a) a[3]=w a=factor(a) a;输入输出数据;序列和向量;向量?矩阵;矩阵;矩阵/高维数组;矩阵/高维数组/scale;Data.frame;缺失值问题等;List;Categorical data A survey asks people if they smoke or not. The data is Yes, No, No, Yes, Yes x=c(Yes,No,No,Yes,Yes) table(x);x factor(x);Barplot:Suppose, a group of 25 people are surveyed as to their beer-drinking preference. The categories were (1) Domestic can, (2) Domestic bottle, (3) Microbrew and (4) import. The raw data is 3 4 1 1 3 4 3 3 1 3 2 1 2 1 2 3 2 3 1 1 1 1 4 3 1 beer = scan() 3 4 1 1 3 4 3 3 1 3 2 1 2 1 2 3 2 3 1 1 1 1 4 3 1 barplot(beer) # this isnt correct barplot(table(beer)) # Yes, call with summarized data barplot(table(beer)/length(beer)) # divide by n for proportiontable(beer)/length(beer);Table/categorical data;smokes = c(Y,N,N,Y,N,Y,Y,Y,N,Y) amount = c(1,2,2,3,3,1,2,1,3,2) (tmp=table(smokes,amount)) # store the table options(digits=3) # only print 3 decimal places prop.table(tmp,1) # the rows sum to 1 now prop.table(tmp,2) # the columns sum to 1 now #上两行等价于下面两行 sweep(tmp, 1, margin.table(tmp, 1), /) sweep(tmp, 2, margin.table(tmp, 2), /) prop.table(tmp)#amount # all the numbers sum to 1 options(digits=7) # restore the number of digits;array/matrix??table??data.frame;写函数;#几个图一起: par(mfrow=c(2,4))#par(mfcol=c(2,4)) layout(matrix(c(1,1,1,2,3,4,2,3,4),nr=3,byrow=T)) hist(rnorm(100),col=Red,10) hist(rnorm(100),col=Blue,8) hist(rnorm(100),col=Green) hist(rnorm(100),col=Brown) #par(mar = c(bottom, left, top, right))设置边缘 #缺省值c(5, 4, 4, 2) + 0.1 (英寸) spring