summary(exercise2_1)行业性别满意度电信业:38男:58不满意:75航空业:19女:62满意:45金融业:26旅游业:37二维列联表:>mytable1>addmargins(mytabl" />
习题
2.1
(1)简单频数分布表:
> load("D:\\工作总结\\人大\\R语言\\《统计学—基于R》(第3版)—例题和习题数据(公开资源)\\exercis e\\ch2\\exercise2_1.RData")
> summary(exercise2_1)
行业性别满意度
电信业:38 男:58 不满意:75
航空业:19 女:62 满意 :45
金融业:26
旅游业:37
二维列联表:
> mytable1<-table(exercise2_1$行业,exercise2_1$满意度)
> addmargins(mytable1) # 增加边界和
不满意满意 Sum
电信业 25 13 38
航空业 12 7 19
金融业 11 15 26
旅游业 27 10 37
Sum 75 45 120
三维列联表:
> mytable1<-ftable(exercise2_1, row.vars = c("性别","满意度"), col.var="行业");mytable1 行业电信业航空业金融业旅游业
性别满意度
男不满意 11 7 7 11
满意 6 3 7 6
女不满意 14 5 4 16
满意 7 4 8 4
(2)
条形图:
> count1<-table(exercise2_1$行业)
> count2<-table(exercise2_1$性别)
> count3<-table(exercise2_1$满意度)
> par(mfrow=c(1,3),mai=c(0.7,0.7,0.6,0.1),cex=0.7,cex.main=0.8)
> barplot(count1,xlab="行业",ylab="频数")
> barplot(count2,xlab="性别",ylab="频数")
> barplot(count3,xlab="满意度",ylab="频数")
帕累托图:
> count1<-table(exercise2_1$行业)
> par(mai=c(0.7,0.7,0.1,0.8),cex=0.8)
> x<-sort(count1,decreasing = T)
> bar<-barplot(x,xlab="行业",ylab="频数",ylim=c(0,1.2*max(count1)),col=2:5) > text(bar,x,labels = x,pos=3) # 条形图增加数值
> y<-cumsum(x)/sum(x) # cumsum累计求和
> par(new=T)
> plot(y,type="b",lwd=1.5,pch=15,axes=F)
> axis(4) # 右Y轴
> mtext("累积频率",side=4,line=3)
> mtext("累积分布曲线",line=-2.5,cex=0.8,adj=0.75)
复式条形图:
> mytable1<-table(exercise2_1$满意度,exercise2_1$行业)
> barplot(mytable1,xlab="行业",ylab="频数",legend=rownames(mytable1),args.legend=list(x= 13), beside = T)
脊形图:
> library(vcd)
> spine(行业~满意度,data=exercise2_1,xlab="满意度", ylab="行业",margins=c(4,3.5,1,2.5))
马赛克图:
> mosaicplot(~性别+行业+满意度,data=exercise2_1,col=2:3)