使用R,ggplot2绘制NMDS图
- 输入文件准备
mapping.txt :样品名与分组信息
otu.txt : otu表
文件下载 : https://github.com/tangshanshaolu/R_test.git
- 代码部分
library(vegan)
library(ggplot2)
pdf("nmds.pdf")
otu<-read.table("otu.txt",header=T,sep="\t")
rownames(otu)<-otu[,1]
otu<-otu[,-1]
otu<-data.frame(t(otu))
map<-read.table("mapping.txt",header=F,sep="\t")
sites<-as.vector(map[,1])
amt<-as.vector(map[,2

最低0.47元/天 解锁文章
1193

被折叠的 条评论
为什么被折叠?



