赞
踩
点击蓝字↑↑↑“微生态”,轻松关注不迷路
本文由阿童木根据实践经验而整理,希望对大家有帮助。
原创微文,欢迎转发转载。
导读R语言lm函数可对两组数据进行回归分析。geom_point函数可以将数据绘制成散点图,geom_smooth函数可以继续在图中添加回归直线。下面展示在一个图中绘制散点和多条回归直线的方法。
一、模拟输入数据
set.seed(1995) # 随机种子data=data.frame(matrix(abs(round(rnorm(40, mean=20, sd=5))), 10, 4)) # 随机正整数,5行,4列colnames(data)=c("y", "x1", "x2", "x3") # 列名data # 查看数据,如下
二、数据编辑
scale(data)# scale函数标准化数据
data.frame(scale(data))
library(reshape)data2=melt(data.frame(scale(data)), id="y")# 使用melt函数将数据框里的数据以“y”变量为ID进行分组
三、ggplot画图
library(ggplot2)# 加载画图包pic labs(x=" ", y="y") +# xy轴标签geom_point() +# 添加散点geom_smooth(method="lm") +# 添加线性回归直线theme(panel.grid=element_blank(), panel.background=element_rect(fill='transparent', color='black')) +# 去点网格、去掉背景、添加边框guides(color=guide_legend(title="x"))# 修改legend标题ggsave(pic, file="pic.png", height=7, width=10)# 打开结果png文件,如下:
感谢阅读~
你可能还喜欢
1 技术贴 | 16S专题 | 简单介绍如何用自己的笔记本处理高通量16S数据
2 技术贴 | 宏基因组专题 | 组装工具盘点和比较
3 技术贴 | R语言菌群Alpha多样性分析和绘图
4 技术贴 | 宏转录组专题 | DDBJ数据库:宏转录组测序数据下载
5 技术贴 | R语言pheatmap聚类分析和热图
微生态科研学术群期待与您交流更多微生态科研问题
(联系微生态老师即可申请入群)
了解更多菌群知识,请关注“微生态”。
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。