向量和数据框
1、向量
- 标量和向量的区分
元素:数字或者字符串(用chr表示)等
标量:一个元素组成的变量
向量:多个元素组成的变量,并且是一排有序排列的元素,可以把一个向量作为数据框中的一列
赋值:给变量定义,意思是赋予这个变量一个数值(不只是数值,包括字符串/数据框等等),第一次赋值被第二次的覆盖了,以此类推,最后一次赋值为最终结果。
> x<- c(1,2,3) #将x定义为由元素1,2,3组成的向量。
> x
[1] 1 2 3
> x<- 1:10 #将1-10之间所有的整数赋值给x
> x
[1] 1 2 3 4 5 6 7 8 9 10
> x<- seq(1,10,by = 1) #1-10之间每隔1取一个数
> x
[1] 1 2 3 4 5 6 7 8 9 10
> x<- rep(1:3,times=2) #1-3 重复2次
> x
[1] 1 2 3 1 2 3
- 从向量中提取元素
(1)根据元素位置
> x
[1] 1 2 3 1 2 3
> x[4] # x第4个元素
[1] 1
> x[-4] #排除法,除了第4个元素之外剩余的元素
[1] 1 2 3 2 3
> x[2:4]#第2到4个元素
[1] 2 3 1
> x[-(2:4)]#除了第2-4个元素
[1] 1 2 3
> x[c(1,5)] #第1个和第5个元素
[1] 1 2
(2)根据值
> x<- seq(-2,10,by = 1) #-2到10之间每隔1取一个数
> x[x==10]#等于10的元素
[1] 10
> x[x<0]
[1] -2 -1
> x[x %in% c(1,2,5)]#存在于向量c(1,2,5)中的元素
[1] 1 2 5
2、数据框
- 读取本地数据
setwd("C:/Users/TL/Desktop/新建文件夹")
a <- read.table("huahua.txt",sep = "\t",header = T)#sep数据按照什么符号隔开,header设置表头(T)与否(F)
- 设置行名和列名
X <- read.csv('doudou.txt') #变量X是一个数据框
> View(X)
> colnames(X) #查看列名
[1] "X1" "X2"
> rownames(X) #查看行名,默认值的行名就是行号,1.2.3.4...
[1] "1" "2" "3" "4" "5"
> colnames(X)[1]<-"bioplanet"#有的数据,左上角第一格为空,R会自动补为x,这个命令来修改
> View(X)
> X<-read.csv("huahua.txt",sep = " ",header =T,row.names=1)#row.names将第一列设为行名
- 数据框的导出
write.table(X,"lixiaoli.txt",sep = ",",quote=F)#sep = ","分隔符改为逗号,quote=F字符串不加双引号
- 变量的保存与重新加载
这次没有处理完的数据下次想接着用,要学会保存和重新加载,保存的格式是RData。
save.image(file="bioinfoplanet.RData")#保存当前所有变量
> save(X,file="test.RData")#保存其中一个变量
> rm(list=ls())#清空环境变量
> options(stringsAsFactors = F)
> load("bioinfoplanet.RData")#再次使用RData时的加载命令
a
X1 X2
1 A 1
2 B NA
3 C NA
4 D 3
5 E NA
> X
X2
A 1
B NA
C NA
D 3
E NA
- 提取元素
> X <- read.csv('doudou.txt') #变量X是一个数据框
> X[1,2] #第1行第2列
[1] 1
> X[1,]#第1行
X1 X2
1 A 1
> X[,2]#第2列
[1] 1 NA NA 3 NA
> X[1] #也是第y列
X1
1 A
2 B
3 C
4 D
5 E
> X[1:2]#第1列到第2列
X1 X2
1 A 1
2 B NA
3 C NA
4 D 3
5 E NA
> X[c(1,2)]#第a列和第b列
X1 X2
1 A 1
2 B NA
3 C NA
4 D 3
5 E NA
> X$x1#也可以提取列(这个命令不用写括号,并且支持Tab自动补全,不过只能提取一列)
NULL
> X$X1#也可以提取列(这个命令不用写括号,并且支持Tab自动补全,不过只能提取一列)
[1] "A" "B" "C" "D" "E"
- 直接使用数据框中的变量
attach
提取某两列case和values作散点图:
将数据框名添加到搜索环境中:attach(a),作图时就只需输入列名(连$都不用了)。
做完后将a删除出搜索环境 detach(a)
options(stringsAsFactors = T)
> a <-data.frame(case=paste0("S",1:9),values=runif(9))
> attach(a)
> plot(case,values)
> detach(a)
Question:save(X,file="test.RData")这句代码如果报错object X not found,是为什么,应该怎么解决?
找不到对象X,说明变量X没有被赋值或者不存在X这个变量,看是否存在名为X的变量,如果没有就重新对X进行赋值;还需要看看大小写有没有输错。