向量:
a=c(1,2,3)
a=c(1,'a',2) # 最终都变为字符型
a=1:10 #创建一个向量,里面包括10个元素
给向量加一个维度,就会变成一个矩阵。
> a=1:10
> a
[1] 1 2 3 4 5 6 7 8 9 10
> dim(a) = c(5,2)
> a
[,1] [,2]
[1,] 1 6
[2,] 2 7
[3,] 3 8
[4,] 4 9
[5,] 5 10
class查看变量类型,str查看变量结构
$用来取元素的列
有选择性地选择索引
> a = 1:10
> dim(a)=c(2,5)
> a
[,1] [,2] [,3] [,4] [,5]
[1,] 1 3 5 7 9
[2,] 2 4 6 8 10
> #有选择性的选择索引
> a[,c(T,F,T,F,T)]
[,1] [,2] [,3]
[1,] 1 5 9
[2,] 2 6 10
grep用于数据搜索
例1:
#使用grep搜索数据框a中,Assay_Type里面行名为RNA-Seq的数据
grep('RNA-Seq',a$Assay_Type) #取出来的是下标
grepl('RNA-Seq',a$Assay_Type) #结果是一堆TRUE或者FALSE
例2:
#单细胞测序中的使用场景,用于将数据集分离
control <- colnames(seurat_obj1)[grep(".1",colnames(seurat_obj1))]
control_data <- seurat_obj1[,colnames(seurat_obj1) %in% control]