统计计算R操作.pdf

《统计计算R操作.pdf》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、R操作李丽贵州大学理学院2016年3月3日星期四R简介•R是一种用于统计计算和统计图形的语言及环境•免费、源代码开放的软件，软件包一直在更新。•R的思想是：它可以提供一些集成的统计工具，但更大量的是它提供各种数学计算、统计计算的函数，从而使使用者能灵活机动的进行数据分析，甚至创造出符合需要的新的统计计算方法。2016年3月3日星期四如何使用R软件？1、help("**")?sample帮助-R函数帮助（文本）-**2、查询函数属于哪个软件包??lm.glm2016年3月3日星期四demo•手动执行##Displayademo,pa

2、usingbetweenpagesdemo(lm.glm,package="stats",ask=TRUE)•自动执行##Displayitwithoutpausingdemo(lm.glm,package="stats",ask=FALSE)2016年3月3日星期四2.2R向量•x1<-c(1,2)•1:5•1:nn:1n=5•x2<-c(1:3)•x3<-c(1:3,11:13)•x<-c(x1,x2)•length(x)2016年3月3日星期四向量运算•c(5,3,2)-c(2,1,0)•c(5,NA,2)-c(NA,1,0

3、)•c(6,18,9)/3•a=c(6,18,9,16)/c(3,2,3,4)•options(digits=3)#保留3位有效数字，默认7•sum(a)•sd(a)#标准差•var(a)#样本方差•range(a)#最大值和最小值2016年3月3日星期四•sort(a)#从小到大排列•sort(a,decreasing=T)•order(a)#从小到大排列的元素下标•a[order(a)]•a[order(a)]==sort(a)2016年3月3日星期四查看变量•ls()objects()•ls.str()#moreinform

4、ation•rm(c1)•rm(list=ls())其他-删除所有对象2016年3月3日星期四数列•a<-numeric(n)•seq(10)•seq(-2,3)•seq(0,1,.1)•seq(0,by=0.1,1)•seq(by=0.1,0,1)•rep(5,4)•rep(c(2,6,8),4)•rep(c(10,28),c(2,3))2016年3月3日星期四2016年3月3日星期四字符型向量•c1=c("x","tan(x)","年龄")•c1=c("Ilovehomeland")•paste("I","love","hom

5、eland")•paste(c("X","Y"),"=",1:4)•paste("x",1:4,sep="")2016年3月3日星期四向量下标•x=seq(1,8,by=2)•x[2]•seq(1,8,by=2)[2]•(c(2,4,6)+(-2))[2]a=c(2,4,6)+(-2)a[2]2016年3月3日星期四•也可单独改变某个元素的值•>x[3]=-1•访问向量的一部分的几种方法：•(1)正整数下标向量•x[c(1,3)]•x[1:2]•(2)负整数下标向量•>x[-1]＃表示删除第一个元素•(3)逻辑下标向量•>x[x<

6、10]＃表示选取x<10的元素•>x[x<1]2016年3月3日星期四(4)字符型下标向量定义向量时可以给元素加上名字>ages=c(Li=23,Zhang=33,Wang=45)>agesLiZhangWang233345访问时可以用通常的方法，还可以用元素名访问>ages[“Zhang”]Zhang332016年3月3日星期四画余弦函数•x1<-0:100•x2<-x1*2*pi/100•y=cos(x2)•plot(x2,y,type="l")2016年3月3日星期四1.00.5y0.0-0.5-1.00123456x220

7、16年3月3日星期四下面定义一个向量的分段函数1xx,0y1xx,0程序：x<-c(-2.3,4,-5,7)y=numeric(length(x))y[x<0]=1-x[x<0]y[x>=0]=1+x[x>=0]y2016年3月3日星期四2.3矩阵和多维数组•1、矩阵函数格式•matrix(data=NA,nrow=1,ncol=1,byrow=FALSE,dimnames=NULL)•data为数组的数据向量，nrow为函数，ncol为列数，byrow=TRUE表示按行写入数据，否则按列。•R中默认是按列填数据。

8、•A=matrix(1:12,ncol=4,byrow=TRUE)•c(A)＃显示A的所有向量，是按列拉直的向量。2016年3月3日星期四•cbind()把向量一列一列拼成一个大矩阵•rbind()把向量一行一行拼成一个大矩阵。•B=rbind(c