统计学基于R语言第二版CH1例题代码.docx
《统计学基于R语言第二版CH1例题代码.docx》由会员分享，可在线阅读，更多相关《统计学基于R语言第二版CH1例题代码.docx（7页珍藏版）》请在冰豆网上搜索。
统计学基于R语言第二版CH1例题代码
统计学基于R语言第二版
Ch1例题代码
第一章
设置路径
setwd（"C:
//Users//ch1//第1章例题CSV格式数据"）
读入数据
example_1<-read.csv（"example1_1.csv"）
文本框1--1
x<-8
y<-c（80,87,98,73,100）
z<-example_1
文本框1--2
sum（y）
mean（y）
var（y）
barplot（y,xlab="类别",ylab="计数",col=c（1:
5））
文本框1--3
help（mean）
help（package="stats"）
文本框1--4
install.packages（"car"）
install.packages（"vcd"）
library（car）
文本框1--5
install.packages（c（"agricolae","aplpack","BSDA","car","corrgram"））
install.packages（c（"DescTools","e1071","fmsb","forecast","gmodels"））
install.packages（c（"gplots","HH","Hmisc","lm.beta","lsr"））
install.packages（c（"pastecs","plotrix","plyr","psych","reshape"））
install.packages（c（"scatterplot3d","sm","TeachingDemos","vcd","vioplot"））
文本框1--6
names<-c（'张青松','王宇翔','田思宇','徐丽娜','张志杰','赵颖颖','王志强','宋丽媛','袁芳芳','张建国'）
stas<-c（68,85,74,88,63,78,90,80,58,63）
math<-c（85,91,74,100,82,84,78,100,51,70）
mark<-c（84,63,61,49,89,51,59,53,79,91）
mana<-c（89,76,80,71,78,60,72,73,91,85）
acco<-c（86,66,69,66,80,60,66,70,85,82）
class<-c（'统计学','数学','营销学','管理学','会计学'）
table1_1<-data.frame（names,stas,math,mark,mana,acco）
colnames（table1_1）<-c（"学生姓名",class）
#table1_1<-data.frame（学生姓名=names,统计学=stas,数学=math,营销学=mark,管理学=mana,会计学=acco）
table1_1
matrix1_1<-matrix（cbind（stas,math,mark,mana,acco）,ncol=5）
dimnames（matrix1_1）<-list（names,class）
save（matrix1_1,file="matrix1_1.RData"）
matrix1_1
文本框1--7
load（"mmatrix1_1.RData"）
文本框1--8
read.csv（'example1_1.csv'）
read.csv（'example1_1.csv',header=F）
文本框1--9
#install.packages（'excel'）
#library（xlsx）
table1_1<-read.csv（'table1_1.csv'）
table1_1
文本框1--10
save（table1_1,file='table1_1.RData'）
write.csv（table1_1,file='table1_1.csv'）
文本框1--11
load（'table1_1.RData'）
table1_1
head（table1_1,3）
tail（table1_1）
t（table1_1）
t
文本框1--12
newdata<-table1_1[order（table1_1$学生姓名）,]
newdata
newdata1<-table1_1[order（table1_1$统计学,decreasing=TRUE）,]
newdata1
文本框1--13
load（'table1_1.RData'）
mean（table1_1$统计学）
mean（table1_1[,2]）
文本框1--14
load（'matrix1_1.RData'）
rowsum1<-colSums（matrix1_1[,1:
5]）
rbind（matrix1_1,totals=rowsum1）
colsum1<-rowSums（matrix1_1）
colsum1
cbind（matrix1_1,totals=colsum1）
matrix1_1
mean（matrix1_1）
apply（matrix1_1,1,sum）
apply（matrix1_1,2,sum）
apply（matrix1_1,1,mean）
apply（matrix1_1,2,mean）
apply（matrix1_1,2,sd）
table1_1
apply（table1_1[,2:
6],2,sum）
apply（table1_1[,2:
6],2,mean）
apply（table1_1[,2:
6],2,sd）
apply（table1_1[,-1,1:
10],1,mean）
文本框1--15
load（"table1_1.RData"）
library（reshape）
table1_1
rename（table1_1,c（学生姓名='姓名',统计学='统计'））
文本框1--16
load（"table1_1.RData"）
table1_1_1<-edit（table1_1）
save（table1_1_1,file='table1_1_1.RData'）
load（"table1_1_1.RData"）
table1_1_1
文本框1--17
newdata<-data.frame（性别=character（0）,年龄=numeric（0））
newdata<-edit（newdata）
save（newdata,file='newdata.RData'）
load（'newdata.RData'）
newdata
文本框1--18
x<-c（2,4,6,NA）
is.na（x）
文本框1--19
x<-c（2,4,6,NA）
sum（x）
sum（x,na.rm=T）
文本框1--20
table1_1_2<-read.csv（'table1_1_2.csv'）
save（table1_1_2,file='table1_1_2.RData'）
load（'table1_1_2.RData'）
table1_1_2
table1_1_3<-na.omit（table1_1_2）
table1_1_3
save（table1_1_3,file='table1_1_3.RData'）
文本框1--21
load（'table1_1.RData'）
table1_1
matrix1_1<-as.matrix（table1_1[,2:
6]）
matrix1_1
rownames（matrix1_1）<-table1_1[,1]
matrix1_1
table1<-as.data.frame（matrix1_1）
table1
文本框1--22
table1_1
as.vector（table1_1$统计学）
as.vector（as.matrix（table1_1[,2:
6]））
文本框1--23
rnorm（10）
set.seed（15）
rnorm（10,50,5）
runif（10,0,100）
rexp（10）
rchisq（10,15）
set.seed（1000）
rnorm（10）
文本框1--24
example_1<-read.csv（'example1_1.csv'）
save（example_1,file='example1_1.RData'）
load（'example1_1.RData'）
example_1
sample（example_1$姓名,10,replace=F）
sample（example_1$姓名,100,replace=T）
sample（example_1$姓名[example_1$考试分数<60]）
sample（example_1$考试分数[example_1$考试分数>90]）
文本框1--25
文本框1--26
load（'example1_1.RData'）
example_1
x<-example_1[,4]
x
myfun<-function（x）{
n<-length（x）
mean<-sum（x）/n
median<-median（x）
r<-max（x）-min（x）
s<-sd（x）
summ<-data.frame（c（mean,median,r,s））
rownames（summ）=c（'平均数','中位数','极差','标准差'）
names（summ）<-'值'
return（summ）
}
myfun（x）
文本框1--27
layout（matrix（c（1,2,3,3）,nrow=2,ncol=2）,widths=c（2,1））
par（mai=c（0.6,0.6,0.1,0.1）,cex=0.7）
x<-rnorm（5000）
y<-rchisq（5000,10）
hist（x,prob=T,col='lightblue',xlab='x',ylab='Density',ylim=c（0,0.4）,main=''）
hist（y,freq=F,col='pink',xlab='y',ylab='Density',main=''）
boxplot（x,col='red',lwd=1）
文本框1--28
layout（matrix（c（1,2,3,3）,nrow=2,ncol=2,byrow=T）,heights=c（2,1））
par（mai=c（0.6,0.6,0.1,0.1）,cex=0.7）
x<-rnorm（1000）
y<-rchisq（1000,10）
hist（x,prob=T,col='lightblue',xlab='x',ylab='Density',ylim=c（0,0.4）,main=''）
plot（x,y,xlab='x',ylab='y'）
boxplot（y,col='red',lwd=1,horizontal=T,varwidth=T）
升级会员 