图鬼狐:跟着生信星球和生信技能树的教程与视频学了一周后,小白菜对收割生信的信心又加大了。今日开始,小白菜记录每天所学,学他一个月,那时应该快到丰收的季节了。(星星眼)
(R和Rstudio的下载不讲了)↓
R常用指令:这是个外部网站,对于使用较频繁的指令整理的比较好,但是要开VPN才能点开链接了。
R的基本变量类型
# 1.向量
a=c(1,2,3)#c <- creat
a=c(1,"a",2)
a=1:10
a=seq(1,10)
a=LETTERS[1:7]
#2.数组
a=1:10
dim(a)=c(2,5)#两行五列,改变a的维度
pheatmap::pheatmap(a)#画热图
class(a)#a的类型
str(a)#查看a的结构
is.matrix(a)#判断a是否是矩阵
as.character(a)#改变a的结构
#3.数据框
b=as.data.frame(a)#a和b数据结构不一样了
str(b)
#4.列表
b$V5#9 10
d=options()#option是一个内置列表
length(d)#88个元素
e=d$repos#CRAN "https://mirrors.tuna.tsinghua.edu.cn/CRAN/"
mode(e)#e是"character"
as.numeric(lapply(d,length))#查看d里面每个元素的长度</pre>
读取文本数据
#取元素
b[,3]
b[c(T,F),3]
#读取文本数据
a=read.table("GSE28735_series_matrix.txt.gz",comment.char='!',header = T,sep = '\t')
#comment.char='!'是注释符号的意思,表明!为注释内容,不读
?read.table #不会就用它查看帮助
row.names(a)=a$ID_REF##去掉行名
a=a[,-1]#去掉第一列
a=log2(a)
pheatmap::pheatmap(a[1:10],)#画热图</pre>