R语言:microeco:一个用于微生物群落生态学数据挖掘的R包,第六:trans_env class

#环境变量在分析微生物群落结构和组装机制方面是非常有用的。我们首先展示RDA分析(db-RDA和RDA)。

复制代码
> t1 <- trans_env$new(dataset = dataset, add_data = env_data_16S)
Env data is stored in object$data_env ...
> t1$cal_ordination(method = "RDA", use_measure = "bray")
No taxa_level provided; use Genus level automatically !
The original ordination result is stored in object$res_ordination ...
The R2 is stored in object$res_ordination_R2 ...
> t1$cal_ordination_envsquare()
Error: attempt to apply non-function
> t1$cal_ordination_envfit()
Result is stored in object$res_ordination_envfit ...
> t1$trans_ordination(adjust_arrow_length = TRUE, max_perc_env = 10)
The result list is stored in object$res_ordination_trans ...
> t1$plot_ordination(plot_color = "Group")

#Mantel检验可以检验环境变量与距离矩阵之间是否存在显著的相关关系。

复制代码
> t1$cal_mantel(use_measure = "bray")
The result is stored in object$res_mantel ...
> t1$res_mantel
   by_group     Variables mantel type Correlation method Correlation coefficient p.value  p.adjusted Significance
1       All      Latitude mantel test            pearson              0.52010091   0.001 0.001375000           **
2       All     Longitude mantel test            pearson              0.37732601   0.001 0.001375000           **
3       All      Altitude mantel test            pearson              0.22102713   0.001 0.001375000           **
4       All   Temperature mantel test            pearson              0.45198101   0.001 0.001375000           **
5       All Precipitation mantel test            pearson              0.27905384   0.001 0.001375000           **
6       All           TOC mantel test            pearson              0.13000213   0.001 0.001375000           **
7       All           NH4 mantel test            pearson             -0.05538846   0.918 0.918000000             
8       All           NO3 mantel test            pearson              0.06758353   0.044 0.048400000            *
9       All            pH mantel test            pearson              0.40853579   0.001 0.001375000           **
10      All  Conductivity mantel test            pearson              0.26425041   0.001 0.001375000           **
11      All            TN mantel test            pearson              0.13205237   
0.002 0.002444444           **

环境变量与分类群之间的相关性对分析和推断群落结构的影响因素具有重要意义。

在本例中,我们首先进行了差异丰度检验和随机森林分析,以获得重要属。然后利用这些分类群进行相关性分析。

> t2 <- trans_diff$new(dataset = dataset, method = "rf", group = "Group", rf_taxa_level = "Genus")

> t1 <- trans_env$new(dataset = dataset, add_data = env_data_16S[, 4:11])

> t1cal_cor(use_data = "Genus", p_adjust_method = "fdr", other_taxa = t2res_rf$Taxa[1:5])

> t1$plot_cor(pheatmap = FALSE)

图画的有些难看,这是临时跑的,大家有需求的话可以用自己的数据跑一下。没需求就跑着玩,熟悉一下就好。

相关推荐
超龄超能程序猿6 小时前
使用 Python 对本地图片进行图像分类
开发语言·人工智能·python·机器学习·分类·数据挖掘·scipy
LabEx1 天前
科研数据可视化核心技术:基于 AI 与 R 语言的热图、火山图及网络图绘制实践指南
人工智能·信息可视化·r语言·r语言绘图·乐备实·labex·科研数据绘图
Jet45051 天前
第100+43步 ChatGPT学习:R语言实现特征选择曲线图
学习·chatgpt·r语言
一条破秋裤1 天前
一份多光谱数据分析
笔记·数据挖掘·数据分析
cal_1 天前
数据分析中的拉链表解析
大数据·数据挖掘·数据分析
九章云极AladdinEdu1 天前
华为昇腾NPU与NVIDIA CUDA生态兼容层开发实录:手写算子自动转换工具链(AST级代码迁移方案)
人工智能·深度学习·opencv·机器学习·华为·数据挖掘·gpu算力
Chef_Chen1 天前
从0开始学习R语言--Day40--Kruskal-Wallis检验
开发语言·学习·r语言
王小王-1231 天前
基于Python的程序员数据分析与可视化系统的设计与实现
python·数据挖掘·数据分析·招聘数据分析·程序员数据分析·招聘薪资数据分析·智联招聘可视化
Deng9452013141 天前
基于数据挖掘的课程推荐系统研究
人工智能·数据挖掘·数据预处理·基于用户的协同过滤·文本特征提取
quant_19861 天前
R语言如何接入实时行情接口
开发语言·经验分享·笔记·python·websocket·金融·r语言