rug(jitter(eruptions))#添加扰动点> lines(density(eruptio.._统计学导论基于r课后题答案">
赞
踩
(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="点图")
画出来的点图:
(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) #添加注释
画出的效果图为:
先到这!
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。