统计学-基于R第3版习题答案(第二章)
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
习题
2.1
(1)简单频数分布表:
工作总结人大\\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)
(3)
饼图:
> count1<-table(exercise2_1$行业)
> name<-names(count1)
> percent<-prop.table(count1)*100
> label1<-paste(name," ",percent,"%",sep="")
> par(pin=c(3,3),mai=c(0.1,0.4,0.1,0.4),cex=0.8) # 圆的大小> pie(count1,labels=label1,init.angle = 90)
扇形图:
> count1<-table(exercise2_1$行业)
> name<-names(count1)
> percent<-count1/sum(count1)*100
> labs<-paste(name," ",percent,"%",sep="")
> library(plotrix)
> fan.plot(count1,labels=labs,ticks=200)
2.2
(1)分10组,绘制频数分布表
工作总结人大\\R语言《统计学—基于R》(第3版)—例题和习题数据(公开资源)\\exercis e\\ch2\\exercise2_2.RData")
> library(actuar)
> v<-as.vector(exercise2_2$灯泡寿命)
> gd1<-grouped.data(v, breaks = 10, right = FALSE)
> table1<-data.frame(gd1);table1
Var.1 v
1 [2600, 2800) 1
2 [2800, 3000) 4
3 [3000, 3200) 12
4 [3200, 3400) 13
5 [3400, 3600) 27
6 [3600, 3800) 20
7 [3800, 4000) 19
8 [4000, 4200) 4
(2)直方图
> d<-exercise2_2$灯泡寿命
> hist(d,breaks=10,xlab="寿命",ylab="频数")