代码老是出错，求大神指点迷津

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

以下为源代码：

# 读取taxonomy，并添加各列名称
taxonomy = read.delim("rep_seqs_tax.txt", row.names= 1,header=F, sep="\t")
colnames(taxonomy) = c("kingdom","phylum","class","order","family","genus","species","evalue")
# 标记差异OTU类型
x$level = as.factor(ifelse(x$sig==1, "enriched",ifelse(x$sig==-1, "depleted","nosig")))
x$otu = rownames(x)
# 转换Pvalue为负对数
x$neglogp = -log(x$PValue)...

完整版如下： 完整版R代码.txt (1.8 KB)

2018-6-23 16:40:03 上传

R代码

上面是网络上的源代码，一直想用自己数据试试，但是运行到第二步老是出错。

是没有定义x么？如何解决？求教各位大神，刚刚接触，想学做个曼哈顿图进行生信分析，瞎忙了两天没有点进展根据网络教程所述最后运行成功的图如下：

附上源文件：
rep_seqs_tax.txt (90.73 KB)

2018-6-23 16:31:41 上传

源文件

真心求教！！！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：R语言

按照你的代码，变量“x”是忽然出现的，没有该变量存在。

good observation