1、R-note一、基本函数1. 函数 c()向量,length() 长度,mode()众数,rbind() 组合,cbind()转置，mode() 属性（数值、字符等）2. 函数 mean( )-中位数, sum( )-求和, min( )-最小值, max( )-最大值, var( )-方差, sd( )-标准差, prod( ) 连乘3. 函数 help（）-帮助4. 正态分布函数 rnorm( ) 、泊松分布函数 rpois( ) 、指数分布函数 rexp( ) 、Gamma 分布函数 rgamma( ) 、均匀分布函数 runif( ) 、二项分布函数 rbinom( ) 、几何分布函数
2、rgeom( ) （一） 基本函数1 2:60*2+115 7 9 11。 。 。 （共 60 个数）2. a5:a 数列第 5 个数，a-5: 删除 a 数列第 5 位数a-(1:5): 删除 a 数列第 1-5 位数ac(2,4,7):a 数列第 2，4,7 位数aa which(a5)1 4 awhich(a5)1 77.rev（）函数-反转举例： a=1:10 rev(a)1 10 9 8 7 6 5 4 3 2 18.sort（）函数- 升序排列举例： a=c(1,4,2,5,3,5,4,7,4) sort(a)1 1 2 3 4 4 4 5 5 79.matrix（）函数 -将数据
3、转变成按行列排布举例： a=1:20 matrix(a,nrow=5,ncol=4),1 ,2 ,3 ,41, 1 6 11 162, 2 7 12 173, 3 8 13 184, 4 9 14 195, 5 10 15 2010.函数 t（） -矩阵行列反置举例： matrix(a,nrow=5,ncol=4),1 ,2 ,3 ,41, 1 6 11 162, 2 7 12 173, 3 8 13 184, 4 9 14 195, 5 10 15 20 t( matrix(a,nrow=5,ncol=4),1 ,2 ,3 ,4 ,51, 1 2 3 4 52, 6 7 8 9 103, 1
4、1 12 13 14 154, 16 17 18 19 2011. diag（）：矩阵对角元素向量或生成对角矩阵举例： a=matrix(1:16,nrow=4,ncol=4) a,1 ,2 ,3 ,41, 1 5 9 132, 2 6 10 143, 3 7 11 154, 4 8 12 16 diag(a)1 1 6 11 16 diag(diag(a),1 ,2 ,3 ,41, 1 0 0 02, 0 6 0 03, 0 0 11 04, 0 0 0 1612.统计分布每一种分布有四个函数：ddensity（密度函数） ，p分布函数，q分位数函数，r随机数函数。比如，正态分布的这四个函数
5、为 dnorm，pnorm，qnorm ，rnorm 。下面我们列出各分布后缀，前面加前缀 d、p、q 或 r 就构成函数名：norm：正态，t ：t 分布，f ：F 分布，chisq：卡方（包括非中心）unif：均匀，exp：指数，weibull：威布尔，gamma ：伽玛， beta：贝塔lnorm：对数正态，logis：逻辑分布，cauchy：柯西，binom：二项分布，geom ：几何分布， hyper：超几何，nbinom ：负二项，pois ：泊松signrank：符号秩，wilcox：秩和，tukey ：学生化极差 13.solve（）：矩阵求逆或解线性方程14. eigen（）
6、：矩阵的特征值分解15.data.frame():生成 excel 类似的数组举例： a=c(1,3,5,6,7) b=c(2,4,6,6,8) x=data.frame(a,b) xa b1 1 22 3 43 5 64 6 65 7 8 (x=data.frame(重量=a,运费 =b)重量 运费1 1 22 3 43 5 64 6 65 7 816.画散点图 plot（）（二）初级函数1.txt 文件读取：举例：(x=read.table(“D:/百度云/abc.txt“)V1 V2 V31 1 2 32 4 5 62.读取 excel 文件1）先将 excel 文件保存为 prn 文件
7、 y yage high weight1 18 150 502 17 160 602）安装 RODBC 安装包，导入 excle 文件举例： local(pkg library(RODBC) z (w for(i in 1:20)ai=i*2+3 a1 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 39 41 434.循环语句：while举例： a1=5 i=1 while(ai a1 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 3919 41 43 45 47 49 51 53 55 57
8、59 61 63 65 67 69 71 73 7537 77 79 81 83 85 87 89 91 93 95 97 99 101 103 105 107 109 11155 113 115 117 119 1215.source（）：数据来源，可将 txt 编辑好的程序直接读取，并运行。举例： source(“D:百度云a.txt“)1 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 3919 41 43 45 47 49 51 53 55 57 59 61 63 65 67 69 71 73 7537 77 79 81 83 85 87
9、 89 91 93 95 97 99 101 103 105 107 109 11155 113 115 117 119 1216.综合举例题目：模拟产生统计专业同学的名单（学号区分） ，记录数学分析，线性代数，概率统计三科成绩，然后进行一些统计分析 结果： num=seq(1008001,1008100) num1 1008001 1008002 1008003 1008004 1008005 1008006 1008007 1008008 100800910 1008010 1008011 1008012 1008013 1008014 1008015 1008016 1008017 10
10、0801819 1008019 1008020 1008021 1008022 1008023 1008024 1008025 1008026 100802728 1008028 1008029 1008030 1008031 1008032 1008033 1008034 1008035 100803637 1008037 1008038 1008039 1008040 1008041 1008042 1008043 1008044 100804546 1008046 1008047 1008048 1008049 1008050 1008051 1008052 1008053 100805
11、455 1008055 1008056 1008057 1008058 1008059 1008060 1008061 1008062 100806364 1008064 1008065 1008066 1008067 1008068 1008069 1008070 1008071 100807273 1008073 1008074 1008075 1008076 1008077 1008078 1008079 1008080 100808182 1008082 1008083 1008084 1008085 1008086 1008087 1008088 1008089 100809091
12、1008091 1008092 1008093 1008094 1008095 1008096 1008097 1008098 1008099100 1008100 x1=round(runif(100,min=80,max=100) x11 83 97 100 90 95 91 93 91 97 87 89 86 94 99 94 97 89 8319 82 98 98 84 100 87 89 97 99 94 88 93 89 90 89 83 93 9537 90 81 94 95 82 94 94 96 82 94 90 91 83 96 97 93 98 8155 87 80 94
13、 84 94 80 87 84 83 85 96 81 97 87 93 97 98 8273 85 90 89 90 88 89 98 88 81 81 96 85 98 87 96 91 86 9891 100 89 87 85 90 90 82 80 81 95 x2=round(rnorm(100,mean=80,sd=7) x21 73 79 88 88 92 79 77 90 87 83 93 85 80 80 84 87 63 72 78 77 85 81 80 7825 79 73 73 90 82 83 80 86 79 74 83 88 90 83 79 83 68 69
14、71 85 83 86 82 6749 84 72 66 91 74 83 73 78 86 77 81 81 89 84 83 82 91 81 80 69 92 85 82 7973 88 69 80 71 89 79 82 71 81 67 82 72 87 78 87 84 81 94 75 80 71 86 78 7597 75 72 95 75 x3=round(rnorm(100,mean=83,sd=18) x31 71 64 84 97 56 71 72 72 96 85 93 80 80 110 78 118 55 8519 99 83 92 113 65 91 76 12
15、0 63 78 72 70 94 101 69 84 100 8737 103 114 77 70 73 54 89 78 66 83 46 113 68 94 89 115 80 10555 67 100 73 89 80 64 115 74 81 99 82 90 92 90 73 101 88 9673 101 80 100 64 74 121 81 94 72 55 78 88 60 103 70 63 96 7291 104 93 92 110 74 82 80 111 64 92 x3which(x3100)=100 x31 71 64 84 97 56 71 72 72 96 8
16、5 93 80 80 100 78 100 55 8519 99 83 92 100 65 91 76 100 63 78 72 70 94 100 69 84 100 8737 100 100 77 70 73 54 89 78 66 83 46 100 68 94 89 100 80 10055 67 100 73 89 80 64 100 74 81 99 82 90 92 90 73 100 88 9673 100 80 100 64 74 100 81 94 72 55 78 88 60 100 70 63 96 7291 100 93 92 100 74 82 80 100 64
17、92 x=data.frame(num,x1,x2,x3) 数据合并 xnum x1 x2 x31 1008001 83 73 712 1008002 97 79 643 1008003 100 88 844 1008004 90 88 975 1008005 95 92 566 1008006 91 79 717 1008007 93 77 728 1008008 91 90 729 1008009 97 87 9610 1008010 87 83 8511 1008011 89 93 9312 1008012 86 85 8013 1008013 94 80 8014 1008014 99
18、 80 10015 1008015 94 84 7816 1008016 97 87 10017 1008017 89 63 5518 1008018 83 72 8519 1008019 82 78 9920 1008020 98 77 8321 1008021 98 85 9222 1008022 84 81 10023 1008023 100 80 6524 1008024 87 78 9125 1008025 89 79 7626 1008026 97 73 10027 1008027 99 73 6328 1008028 94 90 7829 1008029 88 82 7230 1
19、008030 93 83 7031 1008031 89 80 9432 1008032 90 86 10033 1008033 89 79 6934 1008034 83 74 8435 1008035 93 83 10036 1008036 95 88 8737 1008037 90 90 10038 1008038 81 83 10039 1008039 94 79 7740 1008040 95 83 7041 1008041 82 68 7342 1008042 94 69 5443 1008043 94 71 8944 1008044 96 85 7845 1008045 82 8
20、3 6646 1008046 94 86 8347 1008047 90 82 4648 1008048 91 67 10049 1008049 83 84 6850 1008050 96 72 9451 1008051 97 66 8952 1008052 93 91 10053 1008053 98 74 8054 1008054 81 83 10055 1008055 87 73 6756 1008056 80 78 10057 1008057 94 86 7358 1008058 84 77 8959 1008059 94 81 8060 1008060 80 81 6461 1008
21、061 87 89 10062 1008062 84 84 7463 1008063 83 83 8164 1008064 85 82 9965 1008065 96 91 8266 1008066 81 81 9067 1008067 97 80 9268 1008068 87 69 9069 1008069 93 92 7370 1008070 97 85 10071 1008071 98 82 8872 1008072 82 79 9673 1008073 85 88 10074 1008074 90 69 8075 1008075 89 80 10076 1008076 90 71 6
22、477 1008077 88 89 7478 1008078 89 79 10079 1008079 98 82 8180 1008080 88 71 9481 1008081 81 81 7282 1008082 81 67 5583 1008083 96 82 7884 1008084 85 72 8885 1008085 98 87 6086 1008086 87 78 10087 1008087 96 87 7088 1008088 91 84 6389 1008089 86 81 9690 1008090 98 94 7291 1008091 100 75 10092 1008092
23、 89 80 9393 1008093 87 71 9294 1008094 85 86 10095 1008095 90 78 7496 1008096 90 75 8297 1008097 82 75 8098 1008098 80 72 10099 1008099 81 95 64100 1008100 95 75 92 write.table(x,file=“D:百度云mark.txt“,col.names=F,row.names=F,sep=“) -保存数据7.data.frame（）：数据合并-合成举例：x=data.frame(num,x1,x2,x3)8.write.table
24、（）：将数据写入（保存）txt 文档举例：write.table(x,file=“D:百度云mark.txt“,col.names=F,row.names=F,sep=“ “) 9. mean():均值 举例：10.colMeans():列均值举例： colMeans(x)num x1 x2 x3 1008050.50 90.18 80.42 82.91 colMeans(x)c(“x1“,“x2“,“x3“)x1 x2 x3 90.18 80.42 82.9111.apply()：举例：apply(x,2,mean)num x1 x2 x3 1008050.50 90.18 80.42 82
25、.91 apply(x,2,max)num x1 x2 x3 1008100 100 95 100 apply(x,2,min)num x1 x2 x3 1008001 80 63 46 apply(xc(“x1“,“x2“,“x3“),1,sum)求总分1 227 240 272 275 243 241 242 253 280 255 275 251 254 279 256 284 207 240 259 25821 275 265 245 256 244 270 235 262 242 246 263 276 237 241 276 270 280 264 250 24841 223 21
26、7 254 259 231 263 218 258 235 262 252 284 252 264 227 258 253 250 255 22561 276 242 247 266 269 252 269 246 258 282 268 257 273 239 269 225 251 268 261 25381 234 203 256 245 245 265 253 238 263 264 275 262 250 271 242 247 237 252 240 262 which.max(apply(xc(“x1“,“x2“),1,sum)-求总分最大的人是第几个1 46 x$numwhic
27、h.max(apply(xc(“x1“,“x2“),1,sum)-求总分最大的人 num1 1008046第二讲12.hist（）：绘制直方图接上例：hist(x$x1)13.plot（）：绘制散点图接上例： plot(x1,x2) plot(x$x1,x$x2)14.table（）：列联函数 table(x$x1)-各数出现频数80 81 82 83 85 86 87 88 89 6 6 4 6 8 3 5 4 1 90 91 92 93 94 95 96 97 98 3 5 5 7 10 4 2 5 7 99 100 4 515.barplot（）：柱状图接上例： barplot(tabl
28、e(x$x1)16.pie（）：饼图接上例： pie(table(x$x1)17.boxplot():箱尾图接上例： boxplot(x$x1,x$x2,x$x3)1）箱子的上下横线为样本的 25%和 75%分位数 2）箱子中间的横线为样本的中位数 3）上下延伸的直线称为尾线，尾线的尽头为最高值和最低值 4）异常值 18.boxplot（）：箱线图接上例： boxplot(x$x1,x$x2,x$x3,col=c(“red“,“green“,“blue“),notch=T)或者： boxplot(x2:4,col=c(“red“,“green“,“blue“),notch=T)19.boxpl
29、ot（）：水平放置箱尾图接上例：boxplot(x$x1,x$x2,x$x3,horizontal=T)20.star（）：星象图1）每个观测单位的数值表示为一个图形 2）每个图的每个角表示一个变量，字符串类型会标注在图的下方 3）角线的长度表达值的大小 A：接上例： stars(xc(“x1“,“x2“,“x3“)B: stars(xc(“x1“,“x2“,“x3“),full=T,draw.segment=T)饼状分割C: stars(xc(“x1“,“x2“,“x3“),full=F,draw.segment=T)20.aplpack():脸谱图1)用五官的宽度和高度来描绘数值 2）人对
30、脸谱高度敏感和强记忆 3） 适合较少样本的情况 接上例： library(tcltk2) library(aplpack) faces(xc(“x1“,“x2“,“x3“)effect of variables:modified item Var “height of face “ “x1“width of face “ “x2“structure of face“ “x3“height of mouth “ “x1“width of mouth “ “x2“smiling “ “x3“height of eyes “ “x1“width of eyes “ “x2“height of hai
31、r “ “x3“width of hair “ “x1“style of hair “ “x2“height of nose “ “x3“width of nose “ “x1“width of ear “ “x2“height of ear “ “x3“21.其他脸谱： library(TeachingDemos) faces2(x)22 Stem（）：茎叶图 stem(x$x1)The decimal point is at the |80 | 00000000000082 | 000000000084 | 0000000086 | 0000000088 | 0000090 | 00000
32、00092 | 00000000000094 | 0000000000000096 | 000000098 | 00000000000100 | 0000023.QQ 图： qqnorm(),qqline(),1)可用于判断是否正态分布 2)直线的斜率是标准差，截距是均值 3)点的散布越接近直线，则越接近正态分布 接上例： qqnorm(x1) qqline(x1)24.散点图加深：plot （） plot(x$x1,x$x2,main=“数学分析与线性代数成绩的关系 “,xlab=“数学“,ylab=“ 线性“,xlim=c(0,100),ylim=c(0,100),xaxs=“i“,yax
33、s=“i“,col=“red“,pch=19)25.折线图：plot（）26.density（）：密度图 plot(density(rnorm(1000)27.内置数集：mtcars，iris,示例：iris sunflowerplot(iris,3:4,col=“red“,seg.col=“gold“)28.散点图集：pairs（）遍历样本中全部的变量配对，画出二元图 ，直观地了解所有变量之间的关系 示例： pairs(iris,1:4)2） plot(iris,1:4,main=“因素关系“,pch=19,col=“red“,cex=0.9)2）3） par 函数 par(mfrow=c(3,1) plot(x1,x2);plot(x2,x3);plot(x3,x1)29.绘图颜色：colors（）colors()30.绘图设备：dev.cur():31:绘图参数：mai():32.三维散点图：scatterplot3d() library(scatterplot3d) scatterplot3d(x2:4)33.地图：map （） map(“state“,interior=FALSE) map(“world“,fill=TRUE,col=heat.colors(10)
