r语言的下移命令

Posted 2023-02-24

参考技术A 1.getwd（）#查看当前工作目录

2.setwd（'E:graduation-design/OC/sctransform'）#设置工作路径，引号内的内容是一个例子，注意路径最好不要有中文和特殊符号。

3.install.packages.('包')#下载包

4.library（包）#加载包

5.update.packages(包)# 可以更新已经安装的包。

6.help(package = "包") # 获取该R包的帮助文件

7.head(object) #查看对象的前6行

8.tail()#查看对象的后6行

9.[,1]#访问数据的第一列

10.rm(object) #删除对象

11.save.image("myfile") #保存工作空间到文件myfile中（默认值为.RData）

12.example("foo")#函数foo的使用示例（引号可以省略）

13.class(object) #显示某个对象的类或类型

14.c()#定义向量

15.print() #展示数据

16.mode() # 查看数据类型

17.na.rm = TRUE# 去除缺失值

18.is.na() #检查是否有缺失值

19.na.omit() # 删除包含缺失值的行

20.nchar() #统计字符串长度

21.substr(x = ,start = 1, stop = 3) # 提取字符串（首字母大写）

22.toupper()#将字符串都转化成大写

23.tolower() #将字符串都转化成小写

24.gsub("^(\\w)","\\U\\1",tolower(temp),perl = T) #首字母大写

25.gsub("^(\\w)","\\L\\1",upper(temp),perl = T) #首字母小写

26.Sys.Date() # 显示目前系统的时间

27.read.table("c:/路径",seq = "") # 读取excel

28.read.table(网址) #读取网站内容

29.read.table("clipboard")#读取剪切板中的文件

30.read.table(gzfile("文件")) #读取压缩文件

31.write(x,file = "x,txt") #输出文本

32.write.table(x,gzfile("文件名"))# 输出压缩文件

R语言barplot柱状图+散点图+误差棒

今日绘图

今天做这个图

涉及到的很多细节问题：

1. 纵坐标标签水平显示

2. 横坐标下移

3. barplot叠加散点图

4. 柱状图加误差棒

5. excel计算mean(平均值)

6. excel计算se(标准误)

7. R语言创建矩阵和表格
   

所有的细节

都包含在这篇文章里

详情请看代码~

都有相应的解释说明~

所用到的数据

都在下方代码里~

可以复制粘贴到

R软件或RStudio中

直接出图~~~

代码如下

rm(list=ls()) #清除环境

op=par(mar=c(5,5,1,1))#设置图离边缘的距离

#做一个画柱状图的矩阵
data <- c(138302457,33047658) #需要用到的所有数据，用excel的AVERAGE函数算出的平均值（mean）
rnames <- c("PLC/PRF/5","HepG2") #行名
cnames <- c("") #列名
mymatrix <- matrix(data,nrow=2,ncol=1,byrow=F,#2排1列，按列填充
                   dimnames=list(rnames,cnames))
mymatrix #查看已做好的这个矩阵

#做一个画散点图的表格
value <- c(18167561,42749923,38225489,163789467,127506152,123611752) #需要用到的数据
group <- c("b","b","b","a","a","a") #数据所对应的类别，将按字母顺序出图，a对应的数据将出现在第一个
data2 <- data.frame(value,group) #将制作的表格命名为data2，第一列为value，第二列为group

#用矩阵画柱状图
barplot(mymatrix,col=c("white","white"),border = c("black","black"),
        #柱子白色填充，边框都设置为黑色，方便后期用图片处理软件（如：AI、PS）选择边框并将边框加粗
        xlab="",ylab="SOAT1 protein abundance",beside=T, #两个柱子平行
        horiz = F, #柱子垂直于x轴
        cex.lab=1.5, #坐标轴标题的大小
        font.lab=2, #加粗
        space=0.2, #两个柱子之间间隔0.2
        axes=F, #不画坐标轴
        ylim=c(-1*10^7,1.7*10^8), #y轴范围多设置一个刻度，这样才能使横坐标不与柱子底部重叠
        xlim=c(0.1,2.4)) #调整横坐标范围使柱子左侧不与纵坐标间隔太远,调整了起始点,就要根据柱宽/间隔调整终点

#用表格画散点图
stripchart(value~group,data=data2,method="jitter", #方法设置为jitter（随机分离所有点）
           #method默认为“overplot”（值相同的数据点将重合）；还可设置为“stack"（层叠堆放值相同的数据点）
           vertical = T,add=T, jitter=0.1, #垂直，叠加在已有的图上，散点分布的范围为0.1
           col=c("#f65c2a","#33a6e5"),pch=17,cex=2,at=c(0.7,1.9)) #颜色，散点形状，散点大小，散点出现的位置

#用画箭头的函数画误差棒
arrows(x0=0.7, y0=138302457, x1=0.7, y1=151095455, #用excel的STDEV(n1,n2,n3...)/SQRT(n)算出se(标准误)，y1=mean+se
       angle=90,col="#f65c2a", #箭头设置为90°就变成了“T”字形，设置颜色
       lwd = 3, length = 0.2) #线条粗细，横线长度
arrows(x0=1.9, y0=33047658, x1=1.9, y1=40601478, angle=90,col="#33a6e5",
       lwd = 3, length = 0.15)

#设置坐标轴
axis(2,at=c(0,5*10^7,1*10^8,1.5*10^8,1.7*10^8), #2是左侧纵坐标，at设置刻度位置
     label=c("0","5","1","1.5","1.7"), #太长的标签后期用图片处理软件（如：AI、PS）加上去，现在随便设置就行
     lwd=3,lwd.ticks = 3, #刻度标签，线粗细，ticks（小短线）粗细
     cex.axis=1, #标签大小
     las=1, #标签设置为水平
     font = 2) #标签加粗
axis(1,at=c(0.2,2.4),label=c("",""),lwd.ticks = NA,lwd=3) #1是下方横坐标，去掉ticks和标签，只是做一条齐平柱子的横线
axis(1,at=c(0.7,1.9),label=c("PLC/PRF/5","HepG2"),lwd=3,lwd.ticks = 3,cex.axis=1.3,font = 2) #设置横坐标

往期文章

END

排版 | 金忆秋

文字 | 金忆秋（侵删）

图片 | 金忆秋（不可商用）