R语言第三章 统计绘图表示第一节

 

R语言

 

一、综合例子

 

> num=seq(10378001,10378100)

> num

  [1]10378001 10378002 10378003 10378004 10378005 10378006 10378007 10378008

  [9]10378009 10378010 10378011 10378012 10378013 10378014 10378015 10378016

 [17]10378017 10378018 10378019 10378020 10378021 10378022 10378023 10378024

 [25]10378025 10378026 10378027 10378028 10378029 10378030 10378031 10378032

 [33]10378033 10378034 10378035 10378036 10378037 10378038 10378039 10378040

 [41]10378041 10378042 10378043 10378044 10378045 10378046 10378047 10378048

 [49]10378049 10378050 10378051 10378052 10378053 10378054 10378055 10378056

 [57]10378057 10378058 10378059 10378060 10378061 10378062 10378063 10378064

 [65]10378065 10378066 10378067 10378068 10378069 10378070 10378071 10378072

 [73]10378073 10378074 10378075 10378076 10378077 10378078 10378079 10378080

 [81]10378081 10378082 10378083 10378084 10378085 10378086 10378087 10378088

 [89]10378089 10378090 10378091 10378092 10378093 10378094 10378095 10378096

 [97]10378097 10378098 10378099 10378100

> x1=round(runif(100,min=80,max=100))

> x2=round(rnorm(100,mean=80,sd=7))

> x3=round(rnorm(100,mean=83,sd=18))

> x3[which(x3>100)]=100

> x=data.frame(num,x1,x2,x3)

>which.max(apply(x[c("x1","x2","x3")],1,sum))           #取出下标为67的数值(它是最高分)

[1] 67

>x$num[which.max(apply(x[c("x1","x2","x3")],1,sum))]    #取出最高分的学号

[1] 10378067

 

二、对x1进行直方图分析

 

绘制直方函数hist()

hist(x$x1)

 

 

 

三、散点图绘制函数plot()

 

plot(x1,x2);

plot(x$x1,x$x2);

 

 

 

四、列联表分析

 

列联函数table(),柱状图绘制函数barplot()

> table(x$x1);                            #列联函数统计分析各个分数线有多少人

81 82  83  84 85  86  87 88  89  90 91  92  93 94  95  96 97  98  99 100

 1  4   6   5  4   4   3  5   5   3  9   8   5  2   6   5  4   6  13   2

barplot(table(x$x1));  #柱状图显示

五、饼图绘制函数pie()

 

pie(table(x$x1));   #饼图函数pie使用

 

 

 

六、箱尾图

 

1、 箱子的上下横线为样本的25%和75%分位数

2、 箱子中间的横线为样本的中位数

3、 上下延伸的直线称为尾线，尾线的尽头为最高值和最低值

4、 异常值

 

Boxplot() 函数使用

boxplot(x$x1,x$x2,x$x3); #将x1、x2、x3的数据框使用箱尾图表示

 

七、水平放置箱尾图

boxplot(x$x1,x$x2,x$x3,horizontal=T);   #水平放置使用horizontal=T函数

 

 

八、星相图

         1、每个观测单位的数值表示为一个图形

         2、每个图的每个角表示一个变量，字符串会标注在图的下方

         3、角线的长度表达值的大小

 

                  stars(x[c("x1","x2","x3")]);   #使用stars函数算出星相图