R语言第十周作业.pdf

R语言第十周作业

Preparedon22November2020

实验十数据可视化（绘图二）

实验目的：

1.熟悉R条形图的制作

2.熟悉R直方图的制作

3.熟悉用Q-Q图来检验分布的特征

4.了解相关概率计算和随机抽样的R实现

实验要求：完成每个实验内容的同时，在每题后面附上程序代码。要求独立完成。

实验内容：

一．

1.对于iris数据，用每类花（iris$Speciees）的样本数作为高度，制作条形图。

heights-table(iris$Species)

barplot(heights)

2.用每类花的、、、的平均值分别制作条形图，四图同显。

par(mfrow=c(2,2))

attach(iris)

SLength-tapply,Species,mean)SLength-tapply,Species,mean)

SWidth-tapply,Species,mean)SWidth-tapply,Species,mean)

PLength-tapply,Species,mean)PLength-tapply,Species,mean)

PWidth-tapply,Species,mean)PWidth-tapply,Species,mean)

3.分别制作、、、的直方图（用密度值做代表，设置prob=T）,添加拟合的密度曲线，四

图同显。

par(mfrow=c(2,2))

attach(iris)

hist,prob=T,col=5)

lines(density)

hist,prob=T,col=2)

lines(density)

hist,prob=T,col=3)

lines(density)

hist,prob=T,col=4)

lines(density)

4.分别制作、、、的正态Q-Q图，四图同显，并说明他们的正态性。

par(mfrow=c(2,2))

attach(iris)attach(iris)

qqline,col=5)qqline,col=5)

qqline,col=6)qqline,col=6)

qqline,col=9)qqline,col=9)

qqline,col=8)

二．用随机模拟的方法计算下面积分的数值，并和实际值作比较。

ᵰ

∫sinᵆᵅᵆ

0

m-100000

sam(x=runif(m,0,pi),y=runif(m,0,1))

head(sam$y=sin(sam$x))

plot(sam)plot(sam)plot(sam)

a-4*sum(sam$y=sin(sam$x))/m

a

[1]

实际值

b--(cos(pi)-cos(0))

b

[1]2

三．模拟来自指数分布随机变量的样本均值的极限分布。

x-rexp(1000)

n-c(2,5,10,60)

m-10000

xbar()

for(iin1:length(n)){

for(jin1:m)

xbar[j,i]-mean(sample(x,n[i],rep=T))

}

par(mfrow=c(2,2))

for(iin1:4)