生信学习Day5笔记(数据结构)--立小里

向量和数据框

向量和数据框

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进行赋值;还需要看看大小写有没有输错。

©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容