NMDS作图报错

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

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="")
  )
前面运行过程都未报错，但是在运行绘图时报错（如图）请问各位大神这种情况处理？是数据整理出错了吗？（示例数据上传）还是作图代码写错了？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：NMDS MDS Packages install package

dfGroup = read.csv("nmdsgroup.csv",row.names = 1)#这里用了read.csv
data$group<- dfGroup$Group# 这里只添加Group这一列，而不是整个data.frame
就OK了
另外数据名字尽量不要和参数名字一致。比如ggplot里有group这个参数，尽量不要在数据里也用同样的名字命名。

handdyyy 发表于 2022-12-8 20:04
dfGroup = read.csv("nmdsgroup.csv",row.names = 1)#这里用了read.csv
data$group

非常感谢，本来就是copy别人的代码，但是没看清楚。感谢大神的提醒。