R语言操作n

1.加载安装vegan包

2.查看data(varechem)和data(varespec),探索其维度和结构

3.基于varespec构建物种互作网络,输出gml文件并采用gephi可视化为图片,输出pdf,阈值为r>0.6,p<0.05

4.基于varespec和varechem构建物种-环境互作网络,输出gml文件并采用gephi可视化为图片,输出pdf,阔值为r>0.6,p<0.05

r 复制代码
# 1
install.packages("vegan")
library(vegan)

# 2
data(varechem)
data(varespec)
dim(varechem)
dim(varespec)
str(varechem)
str(varespec)

# 3
install.packages(c("Hmisc", "igraph"))
library(Hmisc)
library(igraph)
cor_spec <- rcorr(as.matrix(varespec), type = "pearson")
species <- colnames(varespec)
edges <- data.frame()
for (i in 1:(ncol(varespec)-1)) {
  for (j in (i+1):ncol(varespec)) {
    r_val <- cor_spec$r[i, j]
    p_val <- cor_spec$P[i, j]
    if (abs(r_val) > 0.6 && p_val < 0.05) {
      edges <- rbind(edges, data.frame(from = species[i], to = species[j]))
    }
  }
}
g <- graph_from_data_frame(edges, directed = FALSE)
write_graph(g, "species_network.gml", format = "gml")

# 4
library(psych)
cor_env_spec <- corr.test(varespec, varechem, method = "pearson", adjust = "none")
edges_env <- data.frame()
env_vars <- colnames(varechem)
for (i in 1:ncol(varespec)) {
  for (j in 1:ncol(varechem)) {
    r_val <- cor_env_spec$r[i, j]
    p_val <- cor_env_spec$p[i, j]
    if (abs(r_val) > 0.6 && p_val < 0.05) {
      edges_env <- rbind(edges_env, data.frame(from = species[i], to = env_vars[j]))
    }
  }
}
nodes_env <- data.frame(name = c(species, env_vars),
                        type = c(rep("species", length(species)),
                                 rep("environment", length(env_vars))))
g_env <- graph_from_data_frame(edges_env, directed = FALSE, vertices = nodes_env)
write_graph(g_env, "species_env_network.gml", format = "gml")

gephi安装网址

相关推荐
KY_chenzhao3 分钟前
基于R语言的空间异质性数据分析技术
数据分析·r语言·地理加权回归·判别分析·主成份分析
南瓜胖胖2 小时前
R语言开始绘图--柱状图
开发语言·r语言
纪伊路上盛名在1 天前
python、R、shell兼容1
开发语言·python·jupyter·r语言·shell·生物信息·效率
十三画者2 天前
【工具】grcMalaria用于处理和分析“斑点疟疾基因报告卡”的R软件包
数据挖掘·数据分析·r语言·数据可视化
Teacher.chenchong2 天前
R语言空间分析实战:地理加权回归联合主成份与判别分析破解空间异质性难题
开发语言·回归·r语言
KY_chenzhao2 天前
AI赋能R-Meta分析核心技术:从热点挖掘到高级模型
人工智能·r语言·论文·生态·meta分析·热点
KY_chenzhao2 天前
基于R语言的贝叶斯网络模型实践技术应用:开启科研新视角
r语言·贝叶斯·贝叶斯网络·因果推断·复杂系统
Yolo566Q2 天前
基于R语言地理加权回归、主成份分析、判别分析等空间异质性数据分析实践技术应用
经验分享·数据分析·回归·r语言
还没想好1235 天前
R包gbmt轨迹模型指标计算的函数
r语言
Takoony5 天前
正则表达式r前缀使用指南
开发语言·正则表达式·r语言