统计学基于R语言第二版CH2例题代码.docx
《统计学基于R语言第二版CH2例题代码.docx》由会员分享，可在线阅读，更多相关《统计学基于R语言第二版CH2例题代码.docx（14页珍藏版）》请在冰点文库上搜索。
统计学基于R语言第二版CH2例题代码
统计学基于R语言第二版
Ch2例题代码
第2章
设置路径
setwd（"example//ch2//第2章例题CSV格式数据"）
文本框2--1
example2_1<-read.csv（'example2_1.csv'）
save（example2_1,file='example2_1.RData'）
load（'example2_1.RData'）
example2_1
summary（example2_1）
count1<-table（example2_1$社区）
count1
prop.table（count1）*100
count2<-table（example2_1$性别）
count2
prop.table（count2）*100
count3<-table（example2_1$态度）
prop.table（count3）*100
文本框2--2
mytable1<-table（example2_1$社区,example2_1$性别）
mytable1
addmargins（mytable1）
addmargins（prop.table（mytable1））*100
mytable2<-table（example2_1$社区,example2_1$态度）
mytable2
addmargins（mytable2）
addmargins（prop.table（mytable2））*100
mytable3<-table（example2_1$性别,example2_1$态度）
mytable3
addmargins（mytable3）
addmargins（prop.table（mytable3））*100
文本框2--3
#install.packages（'gmodels'）
library（gmodels）
CrossTable（example2_1$性别,example2_1$态度）
文本框2--4
mytable1<-ftable（example2_1）
mytable1
mytable2<-ftable（example2_1,row.vars=c（'性别','态度'）,col.vars=c（'社区'））
mytable2
addma<-addmargins（table（example2_1））
ftable（addma）
tab1<-table（example2_1）
prop1<-prop.table（tab1）*100
addma2<-addmargins（prop1）
ftable（addma2,row.vars=c（'性别','态度'）,col.vars=c（'社区'））
文本框2--5
example2_2<-read.csv（'example2_2.csv'）
save（example2_2,file='example2_2.RData'）
load（'example2_2.RData'）
example2_2
vector2_2<-as.vector（example2_2$销售额）
vector2_2
library（plyr）
count<-table（round_any（vector2_2,10,floor））
count<-as.numeric（count）
pcount<-prop.table（count）*100
coumsum<-cumsum（pcount）
name<-paste（seq（160,270,by=10）,'-',seq（170,280,by=10）,seq=''）
gt<-data.frame（count,pcount,coumsum）
colnames（gt）<-c（'频数','百分比','累计百分比'）
rownames（gt）<-name
gt<-round（gt,4）
gt
d<-table（cut（vector2_2,breaks=10*（16:
28）,right=F））
percent<-prop.table（d）
cpercent<-cumsum（percent）*100
mytable<-cbind（d,percent,cpercent）
mytable
文本框2--6
load（'example2_1.RData'）
count1<-table（example2_1$社区）
count2<-table（example2_1$性别）
count3<-table（example2_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='社区',horiz=T,main='水平条形图',col=2:
5）
barplot（count2,xlab='性别',ylab='频数',main='垂直条形图'）
barplot（count3,xlab='态度',ylab='频数',main='垂直条形图',col=2:
3）
文本框2--7
load（'example2_1.RData'）
count1<-table（example2_1$社区）
par（mai=c（0.7,0.7,0.1,0.8）,cex=0.8）
x<-sort（count1,decreasing=T）
barplot（x,xlab='社区',ylab='频数',col=2:
5）
y<-cumsum（x）/sum（x）
par（new=T）
plot（y,type='b',lwd=1.5,pch=15,axes=F,xlab=,main=''）
axis（4）
par（las=0）
mtext（'累计频率',side=4,line=3）
mtext（'累计分布曲线',line=-2.5,cex=0.8,adj=0.75）
文本框2--8
load（'example2_1.RData'）
mytable<-table（example2_1$态度,example2_1$社区）
par（mfrow=c（2,2）,cex=0.6）
barplot（mytable,xlab='社区',ylab='频数',ylim=c（0,16）,col=c（'green','blue'）,legend=rownames（mytable）,args.legend=list（x=12）,beside=T,main='社区并列条形图'）
barplot（mytable,xlab='社区',ylab='频数',ylim=c（0,30）,col=c（'green','blue'）,legend=rownames（mytable）,args.legend=list（x=4.8）,beside=F,main='社区并列条形图'）
mytable<-table（example2_1$态度,example2_1$性别）
barplot（mytable,xlab='性别',ylab='频数',ylim=c（0,30）,col=c（'green','blue'）,legend=rownames（mytable）,args.legend=list（x=4.5）,beside=T,main='性别并列条形图'）
barplot（mytable,xlab='性别',ylab='频数',ylim=c（0,60）,col=c（'green','blue'）,legend=rownames（mytable）,args.legend=list（x=1）,beside=F,main='性别并列条形图'）
文本框2--9
load（'example2_1.RData'）
library（vcd）
spine（社区~性别,data=example2_1,xlab='性别',ylab='社区',margins=c（4,3.5,1,2.5））
spine（态度~社区,data=example2_1,xlab='态度',ylab='社区',margins=c（4,3.5,1,2.5））
文本框2--10
load（'example2_1.RData'）
mosaicplot（~性别+社区+态度,data=example2_1,col=2:
3,main=''）
文本框2--11
load（'example2_1.RData'）
count1<-table（example2_1$社区）
name<-names（count1）
percent<-prop.table（count1）*100
label<-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=label,init.angle=90）
library（plotrix）
pie3D（count1,labels=label,explode=0.1,labelcex=0.7）
文本框2--12
library（plotrix）
fan.plot（count1,labels=label,ticks=200,col=c（'gray30','gray70','gray50','gray90'））
文本框2--13
load（'example2_2.RData'）
par（mfrow=c（2,2）,mai=c（0.6,0.6,0.4,0.1）,cex=0.7）
hist（example2_2$销售额,xlab='销售额',ylab='频数',main='普通'）
hist（example2_2$销售额,xlab='销售额',ylab='频数',breaks=20,col='blue',main='分20组'）
hist（example2_2$销售额,xlab='销售额',ylab='密度',breaks=20,col='lightblue',prob=T,main='增加轴须线和核密度线'）
rug（example2_2$销售额）
lines（density（example2_2$销售额）,col='red'）
hist（example2_2$销售额,xlab='销售额',ylab='密度',breaks=20,col='lightblue',prob=T,main='增加正态密度线'）
curve（dnorm（x,mean（example2_2$销售额）,sd（example2_2$销售额））,add=T,col='red'）
rug（jitter（example2_2$销售额））
文本框2--14
hist（faithful$eruptions,breaks=20,prob=T,xlab='喷发持续时间',col='lightblue',main=''）
rug（faithful$eruptions）
lines（density（faithful$eruptions,bw=0.1）,type='l',lwd=2,col='red'）
curve（dnorm（x,mean=mean（faithful$eruptions）,sd=sd（faithful$eruptions））,add=T,col='blue',lwd=2,lty=6）
文本框2--15
load（'example2_2.RData'）
stem（example2_2$销售额）
library（aplpack）
stem.leaf（example2_2$销售额,1,2）
stem.leaf.backback（example2_2$销售额[1:
60],example2_2$销售额[61:
120]）
文本框2--16
par（mfrow=c（3,1）,mai=c（0.4,0.2,0.3,0.2））
x<-rnorm（1000,50,5）
boxplot（x,range=1.5,col='red',horizontal=T,main='相邻值与箱子连线的箱线图（range=1.5）'）
boxplot（x,range=3,col='green',horizontal=T,main='相邻值与箱子连线的箱线图（range=3）'）
boxplot（x,range=0,varwidth=T,col='pink',horizontal=T,main='极值与箱子连线的箱线图'）
文本框2--17
load（'example2_2.RData'）
x<-example2_2$销售额
layout（matrix（c（1,2）,ncol=1）,heights=c（2,1））
par（mai=c（0.15,0.4,0.2,0.2）,cex=0.8）
hist（x,freq=F,col='lightblue',breaks=15,xlab='',ylab='',main=''）
rug（x,col='blue4'）
abline（v=quantile（x）,col='blue4',lwd=2,lty=6）
points（quantile（x）,c（0,0,0,0,0）,lwd=5,col='red2'）
lines（density（x）,col='red',lwd=2）
par（mai=c（0.35,0.42,0.2,0.43）,cex=0.8）
boxplot（x,col='pink',lwd=2,horizontal=T）
rug（x,ticksize=0.1,col='blue4'）
abline（v=quantile（x）,col='blue4',lwd=2,lty=6）
文本框2--18
example2_3<-read.csv（'example2_3.csv'）
save（example2_3,file='example2_3.RData'）
load（'example2_3.RData'）
example2_3
class（example2_3）
boxplot（example2_3,col='lightblue',xlab='运动员',ylab='射击环数',cex.lab=0.8,cex.axis=0.6）
文本框2--19
load（'example2_3.RData'）
library（vioplot）
vioplot（example2_3,col='lightblue'）
文本框2--20
load（'example2_3.RData'）
example2_3<-cbind（example2_3,id=factor（1:
20））
example2_3
library（reshape）
example2_3_1<-melt（example2_3,id.vars=c（'id'）,variable_name='运动员'）
example2_3_1<-rename（example2_3_1,c（value='射击环数'））
save（example2_3_1,file='example2_3_1.RData'）
head（example2_3_1）
tail（example2_3_1）
load（'example2_3_1.RData'）
example2_3_1
dotchart（example2_3_1$射击环数,example2_3_1$运动员,xlab='射击环数',pch=20）
library（lattice）
dotplot（射击环数~运动员,data=example2_3_1,col='blue',pch=19）
文本框2--21
load（'example2_3_1.RData'）
library（lattice）
par（cex=0.7,mai=c（0.6,0.6,0.1,0.1））
densityplot（~射击环数|运动员,data=example2_3_1,col='blue',cex=0.5）
attach（example2_3_1）
library（sm）
pare（射击环数,运动员,lty=1:
6,col=1:
6,lwd=1.5）
legend（'topleft',legend=levels（运动员）,lty=1:
6,col=1:
6）
文本框2--22
x<-seq（0,25,length=100）
y<-4+0.5*x+rnorm（100,0,2）
d<-data.frame（x,y）
plot（x,y）
polygon（d[chull（d）,],col='pink',lty=3,lwd=2）
points（d）
abline（lm（y~x,lwd=2,col=4））
abline（v=mean（x）,h=mean（y）,lty=2,col='gray70'）
文本框2--23
example2_4<-read.csv（'example2_4.csv'）
save（example2_4,file='example2_4.RData'）
load（'example2_4.RData'）
par（mfcol=c（1,2）,mai=c（0.7,0.7,0.3,0.1）,cex=0.8,cex.main=0.8）
plot（广告费用,销售收入,main='普通带网格线',type='n'）
grid（）
points（广告费用,销售收入,main='普通带网格线'）
rug（广告费用,side=1,col=4）
rug（销售收入,side=2,col=4）
plot（广告费用,销售收入,main='带有拟合直线'）
abline（lm（销售收入~广告费用,data=example2_4）,col='red'）
rug（广告费用,side=1,col=4）
rug（销售收入,side=2,col=4）
文本框2--24
par（fig=c（0,0.8,0,0.8）,mai=c（0.9,0.9,0.1,0.1））
plot（广告费用,销售收入,xlab='广告费用',ylab='销售收入',cex.lab=0.7,cex.axis=0.7）
abline（lm（销售收入~广告费用,data=example2_4）,col='blue'）
par（fig=c（0,0.8,0.52,1）,new=T）
boxplot（广告费用,horizontal=T,axes=F）
par（fig=c（0.52,1,0,0.9）,new=T）
boxplot（销售收入,axes=F）
文本框2--25
plot（广告费用,销售收入,xlab='',ylab='销售收入'）
abline（lm（销售收入~广告费用,data=example2_4））
points（销售网点数,销售收入,pch=2,col='blue'）
abline（lm（销售收入~销售网点数,data=example2_4）,col='blue'）
points（销售人员数,销售收入,pch=3,col='blue'）
abline（lm（销售收入~销售人员数,data=example2_4）,col='blue'）
legend（'bottomright',legend=c（'广告费用','销售网点数','销售人员数'）,pch=1:
3,col=c（'black','blue','red'））
文本框2--26
library（scatterplot3d）
attach（example2_4）
s3d<-scatterplot3d（销售人员数,销售网点数,销售收入,pch=16,highlight.3d=T,type='h',cex.lab=0.7）
fig<-lm（销售收入~销售人员数+销售网点数）
s3d$plane3d（fig,col='blue'）
文本框2--27
plot（example2_4,cex=0.8,gap=0.5）
library（car）
attach（example2_4）
scatterplotMatrix（~销售收入+销售网点数+销售人员数+广告费用,diagonal='histogram',gap=0.5）
文本框2--28
n<-30
x<-rnorm（n）
y<-rnorm（n）
z<-abs（rnorm（n））+5:
1
plot（x,y,cex=z,col='pink',pch=19）
points（x,y,cex=z,lwd=1.5）
文本框2--29
load（'example2_4.RData'）
attach（example2_4）
r<-sqrt（销售收入/pi）
symbols（广告费用,销售网点数,circle=r,inches=0.3,fg='white',bg='lightblue',ylab='销售网点数',xlab='广告费用'）
text（广告费用,销售网点数,rownames（example2_4）,cex=0.6）
文本框2--30
example2_5<-read.csv（'example2_5.csv'）
save（example2_5,file='example2_5.RData'）
load（'example2_5.RData'）
par（mai=c（0.7,0.7,0.1,0.1））
matplot（t（example2_5[2:
9]）,type='b',lty=1:
7,xlab='消费项目',ylab='支出金额',pch=1,xaxt='n'）
cn<-as.character（colnames（example2_5）[2:
9]）
axis（side=1,at=1:
8,labels=cn,cex.axis=0.6）
文本框2--31
load（'example2_5.RData'）
#install.packages（'fmsb'）
library（fmsb）
example2_5
radarchart（example2_5[,2:
9],axistype=0,seg=4,maxmin=F,vlabels=names（example2_5[,2:
9]）,pcol=1:
7,plwd=1.5）
rn<-as.character（example2_5[,1]）
legend（x='topleft',legend=rn,col=1:
7,lwd=1,text.width=0.5,cex=0.7）
文本框2--32
load（'example2_5.RData'）
matrix2_5<-as.matrix（example2_5[,2:
9]）
rownames（matrix2_5）<-example2_5[,1]
save（matrix2_5,file='matrix2_5.RData'）
stars（matrix2_5,key.loc=c
升级会员 