2021-10-11 统计学-基于R（第四版）第二章课后习题2.2、2.3记录及总结

本文链接：https://blog.csdn.net/m0_52461904/article/details/120708947

2.2 题目如下

（1）相关理解见注释，完整代码为：

> attach(faithful)
> par(mai=c(0.7,0.7,0.5,0.5),cex=0.7,font.main=1)
> hist(eruptions,prob=TRUE,col="pink",xlab="eruptions",ylab="密度",main="直方图")
#画直方图
> rug(jitter(eruptions))
#添加扰动点
> lines(density(eruptions),col="blue3",lwd=1)
#添加核密度曲线

画出来的效果图为：

（2）第二题又狠狠坑了我一把，答案的横坐标名称写错了还有课本上的奇怪格式也一直让我看不明白，心累，总之代码很简单：

> library(reshape2)
> library(DescTools)
> PlotMultiDens(faithful,xlab="观测数据",main="核密度比较曲线")

画出来的效果图为：

（3）相关理解见注释，完整代码为：

> par(mai=c(0.6,0.6,0.3,0.3),cex=0.7,font.main=1)
> boxplot(faithful,xlab="var",ylab="value",main="箱线图",col=palette)
> library(vioplot)
> vioplot(faithful,xlab="var",ylab="value",main="小提琴图",col=palette)

画出来的效果图为：

（4）画茎叶图完整代码为：

> library(aplpack)
> stem.leaf(faithful$eruptions)
1 | 2: represents 1.2
 leaf unit: 0.1
            n: 272
   12     s | 667777777777
   51    1. | 888888888888888888888888888899999999999
   71    2* | 00000000000011111111
   87     t | 2222222222333333
   92     f | 44444
   94     s | 66
   97    2. | 889
   98    3* | 0
  102     t | 3333
  108     f | 445555
  118     s | 6666677777
  134    3. | 8888888889999999
  (31)   4* | 0000000000000000111111111111111
  107     t | 22222222222233333333333333333
   78     f | 44444444444445555555555555555555555
   43     s | 6666666666677777777777
   21    4. | 88888888888899999
    4    5* | 0001

> dotchart(faithful$eruptions,xlab="eruptions",main="点图")

画出来的点图：

2.3 题目如下

（1）完整代码如下：

> plot(mtcars$mpg,mtcars$wt,pch=19,col="green4",xlab="mpg",ylab="wt",main="散点图")

画出的效果图为：

（2）完整代码如下：

> plot(mtcars,col="green4",main="散点图矩阵",gap=0.5)

画出的效果图为：

（3）完整代码如下：

> library(corrgram)
> corrgram(mtcars,main="相关系数矩阵",order=TRUE,lower.panel = panel.shade,upper.panel = panel.conf)
#order表示按相关系数排列变量，lower.panel表示对角线下方绘制阴影线，upper.panel用于对角线上方绘制饼图

画出的效果图为：

（4）完整代码如下：

> library(lattice)
> p1<-cloud(mpg~hp+wt,data = mtcars,par.setting=list(par.xlab.text=list(cex=0.6),par.ylab.text=list(cex=0.6),par.zlab.text=list(cex=0.6)),main="3D散点图")
#设置x、y、z轴标签字体大小和设置标题
> plot(p1)  #画3D散点图
> attach(mtcars)
> symbols(mpg,hp,wt,inches = 0.15,fg="black",bg="pink",main="气泡图")  #画气泡图
> mtext("气泡大小=wt",line=-2,cex=0.8,adj=0.5)  #添加注释

画出的效果图为：