美文网首页生信星球培训第六期
学习小组Day5笔记--徐沫沫

学习小组Day5笔记--徐沫沫

作者: 徐沫沫 | 来源:发表于2018-09-21 23:17 被阅读14次

    理解数据类型:向量和数据框

    image.png

    part1: 向量

    1. 向量

    • 区分标量和向量
      “元素” 指的是数字或者字符串(用chr表示)等,根据元素个数区分:
      标量:一个元素组成的变量,可以是一个数字或字符串,字符串使用时要带双引号;
      向量:多个元素组成的变量,多个元素的有序排列,例如:c(1,2,3);
      给变量 赋值(定义):
      x<-c(1,2,3) # 将 x 定义为由元素1,2,3组成的向量;
      x<-(1:10) # x 定义为1-10之间所有的整数组成的向量;
      x<-seq(1,10,by = 0.5) # x 定义为1-10之间每隔0.5取一个数组成的向量;
      x<-rep(1:3,times=2) # x 定义为1-3 之间的整数再重复2次,组成的向量;
      image.png

    遇到了问题,我的变量不能查看,,,,哭哭哭~~~~

    2. 从向量中提取元素

    (1)根据元素位置,方括号内是元素的位置,- 表示除了:

    x[4] # x第4个元素
    x[-4] # 排除法,除了第4个元素之外剩余的元素
    x[2:4] # 第2到4个元素
    x[-(2:4)] # 除了……
    x[c(1,5)] # 第1个和第5个元素
    

    (2)根据值的大小

    x[x==10] # 等于10的元素
    x[x %in% c(1,2,5)] # 存在于向量c(1,2,5)中的元素
    

    Part2:数据框

    注意 将示例数据放在自己的 RData 文件夹下!!!

    获得示例数据:
    (1)新建txt文件,例如doudou.txt,输入内容,NA表示空值;
    用以下命令即可获得示例数据框:
    X<-read.csv('doudou.txt')

    (1)读取本地数据

    read.table(file = "huahua.txt",sep = "\t",header =T) #读取文件,分隔符为Tab,有表头,文件内容显示在控制台;
    a<-read.table(file = "huahua.txt",sep = "\t",header =T) #把这个数据框赋给一个变量a;,输入a,查看a的内容;

    (2)设置行名和列名

    X<-read.csv('doudou.txt') #变量X是一个数据框
    colnames(X) #查看列名
    rownames(X) #查看行名,默认值的行名就是行号,1.2.3.4...
    colnames(X)[1]<-"species" #有的公司返回数据,左上角第一格为空,R会自动补为x,用这个命令来修改
    X<-read.csv(file = "huahua.txt",sep = "\t",header =T,row.names=1) # 最后row.names=1 的意思是修改第一列为行名

    (3)数据框的导出

    write.table(X,file = "yu.txt",sep = ",",quote=F) # 将X内容保存名为yu.txt,将分隔符改为逗号,字符串不加双引号(默认格式带由双引号)

    (4)变量的保存与重新加载

    保存的格式是RData。

    save.image(file="practice.RData") # 保存当前所有变量,名为practice.RData
    save(frame1,file="frame1.RData") # 保存其中一个变量frame1,名为frame1.RData
    load("frame1.RData") # 再次使用RData时,加载frame1.RData。

    (5)提取元素

    X[x,y] # 提取变量X的第x行第y列
    X[x,] #第x行 注意行列的逗号的位置
    X[,y] #第y列 ,展示成一行
    X[y] #第y列 ,展示成一列
    X[a:b] #第a列到第b列
    X[c(a,b)] #第a列和第b列
    X$列名 #根据列名提取列

    (6)使用数据框中的变量

    提取某两列(例如:case列和values列)作散点图:

    c <-data.frame(case=paste("S",c(1:50)),values=runif(50))
    plot(c$case,c$values)
    

    --这种方法有个缺点:数据框名a在代码中重复出现,好像还能忍~~~~

    更简单的方法:

    方法1:attach

    将数据框名添加到搜索环境中:attach(c),作图时就只需输入列名,不用$。

    attach(c)
    plot(case,values)
    

    做完后将a删除出搜索环境 detach(c)
    局限性:两个以上数据框的列名有冲突时,attach会报错。

    方法2:with
    with(c,{
     plot(case,values)
     x<<-summary(values)   # 赋值给x,<<的意思是作为全局变量,也就是出了with循环仍有效。
    })
    x # 运行完后输出x
    

    相关文章

      网友评论

        本文标题:学习小组Day5笔记--徐沫沫

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