gephi插件不可用_利用Gephi软件绘制网络图

作者:中国科学院天津工业生物技术研究所 王敬敬 博士

利用Gephi软件绘制网络图

1. 生成物种相关性矩阵

此步骤需要在R语言环境下运行,依赖psych包,输入文件为典型的OTU表或属水平丰度矩阵,示例如下。

be22f5efe8e67b4e14a2ef2253f31629.png

输入文件FH_CK.txt文件格式, CK为空白对照的5个重复;FH为处理组的5个重复;G1-G148为丰度大于0.2%的属。

# 安装需要的包
install.packages("psych")

# 加载包
library(psych)

# 读取otu-sample矩阵,行为sample,列为otu
otu=read.table("FH_CK.txt", head=T, row.names=1)

# 计算OTU间两两相关系数矩阵
# 数据量小时可以用psych包corr.test求相关性矩阵,数据量大时,可应用WGCNA中corAndPvalue, 但p值需要借助其他函数矫正
occor = corr.test(otu,use="pairwise",method="spearman",adjust="fdr",alpha=0.05)
occor.r = occor$r # 取相关性矩阵R值
occor.p = occor$p # 取相关性矩阵p值

# 确定物种间存在相互作用关系的阈值,将相关性R矩阵内不符合的数据转换为0
occor.r[occor.p>0.05|abs(occor.r)<0.6] = 0

# 将occor.r保存为csv文件
write.csv(occor.r,file="FH_CK_0.05_occor.csv")

2. Gephi生成点、边文件

从官网https:/

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值