生信星球培训第108期

生信星球小组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变量的位置,检查原因

上一篇 下一篇

猜你喜欢

热点阅读