学习小组Day5笔记——山川石

数据类型

数据类型

向量(vector)

标量与向量的区分

  • 元素:数字或者字符串。
  • 标量:单个元素组成的变量。
  • 向量:多个元素组成的变量。
    标量和向量都属于变量。

赋值 <-

  1. x <- c(1,2,3) #向量赋值
  2. x <- 1:10 #1-10的所有整数
  3. x <- seq(1,10,by = 0.5) #1-10之间每隔0.5取一个数
  4. x <- rep(1:3,times = 2) #1-3重复两次

引自CSDNR语言学习笔记(1)——rep函数(Kalinda_yu)

rep函数的用法:

  • rep(x, …):将变量x的值循环n遍
  • times:整个数组循环几遍
  • each:每个元素循环几遍
  • length.out 输出长度为多少,长了会被截掉,短了会根据前面的规则补上

向量中提取元素

  1. 根据元素位置
x[-4]#排除法,除了第4个元素之外剩余的元素
x[2:4]#第2到4个元素
x[-(2:4)]#除了第2-4个元素
x[c(1,5)] #第1个和第5个元素
  1. 根据值
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。

你可能感兴趣的:(学习小组Day5笔记——山川石)