统计学--基于R（第3版）（基于R应用的统计学丛书）作者：贾俊平 习题答案 第三章

3.1

#3_1
load('C:/exercise/ch3/exercise3_1.RData')
exercise3_1
#(1)
#计算平均数
mean(exercise3_1$网购金额)
#计算标准差
sd(exercise3_1$网购金额)
#计算极差
range<-max(exercise3_1$网购金额)-min(exercise3_1$网购金额)
range
#计算四分位差
IQR(exercise3_1$网购金额,type=6)

#(2)
#计算10%、25%、50%、75%、90%的分位数
quantile(exercise3_1$网购金额,probs=c(0.1,0.25,0.5,0.75,0.9),type=6)

#(3)
#计算标准分数，检测数据的离群点
as.vector(round(scale(exercise3_1$网购金额),4))
#通常将3个标准差之外的数据确定为离群点，由表中数据可见只有3.2108一个，即第14个数据(2270)为离群点

#(4)计算偏度系数和峰度系数，分析网购金额分布特点
#library(agricolae)#若该package无法下载，可用package(e1071)

library(e1071)
#偏度系数
skewness(exercise3_1$网购金额)#中等偏斜分布
#峰度系数
kurtosis(exercise3_1$网购金额)#尖峰分布

3.2

#3_2
#选择适当的图形和统计量比较3种方法组装产品数量的分布特点
load('C:/exercise/ch3/exercise3_2.RData')
exercise3_2
#平均数
mean(exercise3_2$方法A)#最大
mean(exercise3_2$方法B)
mean(exercise3_2$方法C)
#偏度系数
library(e1071)
skewness(exercise3_2$方法A)#轻微偏斜（右偏）
skewness(exercise3_2$方法B)#轻微偏斜（左偏）
skewness(exercise3_2$方法C)#严重偏斜（左偏）
#峰度系数
kurtosis(exercise3_2$方法A)#扁平分布
kurtosis(exercise3_2$方法B)#扁平分布
kurtosis(exercise3_2$方法C)#尖峰分布
#箱线图
boxplot(exercise3_2,col="lightblue",xlab="方法",ylab="产品数量")
#直方图
par(mfrow=c(2,2),mai=c(0.6,0.6,0.4,0.1),cex=0.7)
hist(exercise3_2$方法A,breaks="Sturges",xlab="产品数量",ylab="频数",col="lightblue",main="方法A")
hist(exercise3_2$方法B,breaks="Sturges",xlab="产品数量",ylab="频数",col="green",main="方法B")
hist(exercise3_2$方法C,breaks="Sturges",xlab="产品数量",ylab="频数",col="red",main="方法C")
#茎叶图
library(aplpack)
stem.leaf(exercise3_2$方法A)
stem.leaf(exercise3_2$方法B)
stem.leaf(exercise3_2$方法C)

3.3

#3_3
#绘制背靠背茎叶图，计算有关的描述统计量，比较分析男女学生身高的特点和差异
load('C:/exercise/ch3/exercise3_3.RData')
exercise3_3
library(aplpack)
#背靠背茎叶图
stem.leaf.backback(exercise3_3$男生,exercise3_3$女生)

#平均身高
mean(exercise3_3$男生)
mean(exercise3_3$女生)

#众数
mode<-function(x){
sx<-sort(unique(x)) #列出单一的数据并排序，再存入sx
tab<-tabulate(match(x,sx)) #比较x和sx中相同的数值，列出它们在sx中的位置，再计算每个位置的频数放入对象tab中
sx[tab==max(tab)] #找出sx对象中频数最多的元素
}
mode(exercise3_3$男生)
mode(exercise3_3$女生)

#或利用which.max函数计算众数
which.max(table(exercise3_3$男生))
which.max(table(exercise3_3$女生))

#中位数
median(exercise3_3$男生)
median(exercise3_3$女生)