美文网首页生信星球培训第四十一期
学习小组Day5笔记——山川石

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

作者: 山川石 | 来源:发表于2020-03-07 09:56 被阅读0次

    数据类型

    数据类型

    向量(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笔记——山川石

        本文链接:https://www.haomeiwen.com/subject/zmofdhtx.html