搜索
人大经济论坛 附件下载

附件下载

所在主题:
文件名:  nmdsdata.xlsx
资料下载链接地址: https://bbs.pinggu.org/a-3813941.html
附件大小:
10.88 KB   举报本内容
NMDS作图
#install.packages("ggplot2")
#install.packages("permute")
#install.packages("lattice")
#install.packages("vegan")
library(ggplot2)
library(permute)
library(lattice)
library(vegan)

df = read.csv("C:/Users/tangyj/Desktop/NMDS/nmdsdata.csv",header = T, row.names = 1)
#df=t(df) # 对数据进行转置
dfGroup = read.delim("C:/Users/tangyj/Desktop/NMDS/nmdsgroup.csv",row.names = 1)
dfNmds<-metaMDS(df,distance="bray",k = 2)
# 绘图前的数据整理
data = data.frame(dfNmds$points)
data$group = dfGroup$Group

# 绘图
#使用ggplot2画图,逐步添加
ggplot(data,aes(x = MDS1,
y = MDS2,
color = group,
group = group,
fill = group)
)+
geom_point(size=2)+
theme_classic()+
stat_ellipse( # 添加置信区间
geom = "polygon",
level = 0.95,
alpha=0.3)+
geom_text( # 添加文本标签
aes(label=rownames(data)),
vjust=1.5,
size=2,
color = "black"
)+
labs( # 在副标题处添加stress
subtitle = paste("stress=",round(dfNmds$stress,3),sep="")
)
前面运行过程都未报错,但是在运行绘图时报错(如图)请问各位大神这种情况处理?是数据整理出错了吗?(示例数据上传)还是作图代码写错了?



    熟悉论坛请点击新手指南
下载说明
1、论坛支持迅雷和网际快车等p2p多线程软件下载,请在上面选择下载通道单击右健下载即可。
2、论坛会定期自动批量更新下载地址,所以请不要浪费时间盗链论坛资源,盗链地址会很快失效。
3、本站为非盈利性质的学术交流网站,鼓励和保护原创作品,拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知,将积极的采取必要措施;同时,本站也将在技术手段和能力范围内,履行版权保护的注意义务。
(如有侵权,欢迎举报)
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

GMT+8, 2026-1-15 21:58