生信星球小组Day3
2021-05-16 本文已影响0人
生信小白
向量
1标量和向量
x <- c(1,2,3) # 将x定义为元素1,2,3组成的向量
x <- 1:10 # 将x定义为1-10之间所有的整数
x <- seq(1,10, by = 0.5) # 1-10之间每隔0.5取一个数
x <- rep(1:3, times = 2) # 1-3,并且重复两次
2从向量中提取元素
## 根据位置提取
x[4] #提取第四个
x[-4] #除了第四个都提取
x[2:4] #提取2-4元素
x[c(1,5)] #提取第一个和第五个元素
## 根据值提取
x[x==10] #提取等于10的元素
x[x<0] #提取小于0的元素
x[x %in% c(1,2,5)] #提取存在于1,2,5向量中的元素
数据框
# 01_设置工作路径
getwd()
setwd("/Users/apple/Desktop/")
# 读取doudou文件
x <- read.csv("doudou.txt")
# 读取huahua文件
y <- read.csv("huahua.txt",sep = "",header = T,row.names = 1)
# 02_查看行列名
colnames(x)
rownames(x)
# 补充左上1
colnames(x)[1] <- "bioplanet"
# 03_数据框导出
write.table(x, # 变量x
file = "yu.txt", #输出文件名
sep = ",", # 分隔符,
quote = F) # 默认带双引号,选择FALSE
# 04_保存变量
save.image(file = "bioinfoplanet.Rdata") #保存所有变量
save(x,y,file = "bioinfoplanet2.Rdata") #保存x,y两个变量
load("bioinfoplanet2.Rdata")
# 05_提取元素
x[x,y] #提取x行y列
x[x,] #提取x行
x[,y] #提取y列
x[y] #默认列
x[a:b] #提取a-b列
x[c(a,b)] #提取a,b两列
x$X1 #提取x的X1列,指定列名提取
报错save(x,file="test.Rdata")
说明x这个变量不存在,找到定义x变量的位置,检查原因