生信星球培训第十六期

学习小组day5 ying

2019-06-10  本文已影响0人  Ying_a20b

显示工作路径 getwd()

设置工作路径 setwd("C:/Users/PoPo/Desktop/Rdata")

函数或者命令不会用时,除了百度/谷歌搜索以外,用这个命令查看帮助:?read.table,调出对应的帮助文档,翻到example部分研究一下。

1. 向量和数据框是两种重要的数据类型

向量是由多个元素组成的变量,标量是一个元素组成的变量,元素可以是数字或者字符串,字符串必须加引号

给变量赋值

x<-c(1,2,3) #常用的向量定义写法,意为将x定义为由元素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 重复2次

2. 从向量中提取元素

根据元素的位置

x[4]#x第4个元素

x[-4]#排除法,除了第4个元素之外剩余的元素

x[2:4]#第2到4个元素

x[-(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)中的元素


数据框

在Rdata新建doudou.txt

X1,X2

A,1

B,

C,

D,3

E,

小抄截图中显示的NA表示空值,所以新建的时候像我一样空着就好。

用以下命令即可获得示例数据框:

X<-read.csv('doudou.txt')



读取本地数据

read.table(file ="huahua.txt",sep ="\t",header =T)#读取文件,分隔符为Tab,有表头 

a<-read.table(file ="huahua.txt",sep ="\t",header =T)#把这个数据框赋给一个变量a 



设置行名和列名

X<-read.csv('doudou.txt') #在示例数据里有doudou.txt 注意这里的变量X是一个数据框

colnames(X)#查看列名

rownames(X)#查看行名,默认值的行名就是行号,1.2.3.4...

colnames(X)[1]<-"bioplanet"#有的公司返回数据,左上角第一格为空,R会自动补为x,用这个命令来修改

X<-read.csv(file ="huahua.txt",sep ="\t",header =T,row.names=1)#最后row.names的意思是修改第一列为行名


数据框的导出

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



变量的保存与重新加载

#这次没有处理完的数据下次想接着用怎么办?--学会保存和重新加载。保存的格式是RData。

save.image(file="bioinfoplanet.RData")#保存当前所有变量

save(frame1,file="frame1.RData")#保存其中一个变量

load("frame1.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$列名#也可以提取列(优秀写法,而且这个命令还优秀到不用写括号的地步,并且支持Tab自动补全哦,不过只能提取一列)



直接 使用数据框中的变量

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

> plot(c$case,c$values)

上一篇下一篇

猜你喜欢

热点阅读