数据类型
向量(vector)
标量与向量的区分
- 元素:数字或者字符串。
- 标量:单个元素组成的变量。
- 向量:多个元素组成的变量。
标量和向量都属于变量。
赋值 <-
x <- c(1,2,3) #向量赋值
x <- 1:10 #1-10的所有整数
x <- seq(1,10,by = 0.5) #1-10之间每隔0.5取一个数
x <- rep(1:3,times = 2) #1-3重复两次
引自CSDNR语言学习笔记(1)——rep函数(Kalinda_yu)
rep函数的用法:
- rep(x, …):将变量x的值循环n遍
- times:整个数组循环几遍
- each:每个元素循环几遍
- length.out 输出长度为多少,长了会被截掉,短了会根据前面的规则补上
向量中提取元素
- 根据元素位置
x[-4]#排除法,除了第4个元素之外剩余的元素
x[2:4]#第2到4个元素
x[-(2:4)]#除了第2-4个元素
x[c(1,5)] #第1个和第5个元素
- 根据值
x[x==10]#等于10的元素
x[x<0]
x[x %in% c(1,2,5)]#存在于向量c(1,2,5)中的元素
数据框(data frame)
读取本地数据
注意:一定要将数据放在工作目录里,否则报错
read.table(file = "huahua.txt",seq = "\t",header = T)
#header是逻辑值,指示文件的第一行是否包含变量的名称header = T代表第一行用于列名称,header=F则为第一行即为具体数据;seq是字段分隔符,代表每行数据由该字符分隔
a <- read.table(file = "huahua.txt",seq = "\t",header = T)
#将读取的数据赋给变量a
设置行名和列名
X <- read.csv('doudou.txt') #变量X是一个数据框
colnames(X) #查看列名
rownames(X) #查看行名,默认值的行名就是行号,1.2.3.4...
colnames(X)[1] <- "bioplanet"#修改第一列行名
X <- read.csv(file = "huahua.txt",sep = " ",header =T,row.names=1)#row.names的意思是修改第一列为行名
数据框的导出
write.table(X,file = "yu.txt",sep = ",",quote=F)#file = "导出文件名及类型"
变量的保存与重新加载
数据保存的格式是RData。
save.image(file="bioinfoplanet.RData")#保存当前所有变量
save(X,file="test.RData")#保存其中一个变量
load("test.RData")#再次使用RData时的加载命令
提取元素
X[x,y]#第x行第y列
X[x,]#第x行
X[,y]#第y列
X[y] #也是第y列
X[a:b]#第a列到第b列
X[c(a,b)]#第a列和第b列
X$列名#提取指定列
直接使用数据框中的变量
- attach
a <-data.frame(case=paste("S",1:50)),values=runif(50))
attach(a)
plot(case,values)
attch意为锁定数据框;detach则是解锁数据框。
使用两个以上时,数据框中不要出现相同的的列名。
- with
with(a,{
plot(case,values)
x<<-summary(values) #求和并赋值给x,<<的意思是作为全局变量,也就是出了大括号仍有效。
})
x
with就是把操作限制在数据框上,“,”前为数据框名称,“{}”内为具体操作。
问题
save(X,file="test.RData")这句代码如果报错X not found,是为什么,应该怎么解决?
字母的大小写是严格区分的,可能之前变量名为x,检查environment里面是否有X。