楼主: weishanyun4
2135 2

[问答] 请帮我看看这个代码醋错误在哪个地方 [推广有奖]

  • 0关注
  • 0粉丝

高中生

22%

还不是VIP/贵宾

-

威望
0
论坛币
2929 个
通用积分
0.0000
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
202 点
帖子
5
精华
0
在线时间
46 小时
注册时间
2017-3-27
最后登录
2025-8-13

楼主
weishanyun4 发表于 2020-8-28 10:22:55 |AI写论文
20论坛币
library(ggpubr)
setwd("C:\\Users\\ALIENWARE\\")            
rt=read.table("risk.txt",sep="\t",header=T,row.names=1,check.names=F)        
rt=rt[,3:(ncol(rt)-2)]
Type=read.table("clinical.txt",sep="\t",check.names=F,header=T)
row.names(Type)=Type[,1]
sameSample=intersect(row.names(Type),row.names(rt))
rt=rt[sameSample,]
Type=Type[sameSample,]
#准备
data=data.frame()
for(i in colnames(rt)){
  data=rbind(data,cbind(expression=rt[,i],gene=i,Type=as.vector(Type[,2])))
}
write.table(data,file="data.txt",sep="\t",row.names=F,quote=F)
#作图
data=read.table("data.txt",sep="\t",header=T,check.names=F)      
p=ggboxplot(data, x="gene", y="expression", color = "Type",
     ylab="lncRNA expression",
     xlab="",
     palette = rainbow(length(levels(data$Type))) )
p=p+rotate_x_text(45)
pdf(file="boxplot.pdf",width=7,height=5)                          
p+stat_compare_means(aes(group=Type),symnum.args=list(cutpoints = c(0, 0.001, 0.01, 0.05, 1), symbols = c("***", "**", "*", "ns")),label = "p.signif")
dev.off()

关键词:Expression intersect cutpoints Clinical express

沙发
llb_321 在职认证  发表于 2020-8-29 19:17:52
ggboxplot之后的代码应该没有问题,除非palette语句可能报错。
之前的数据准备和处理,由于没有你的原始数据,无法判断。
而且,你也没有说明到底是哪种报错信息,所以,你给大家出了个难题欧
已有 1 人评分论坛币 收起 理由
cheetahfly + 10 热心帮助其他会员

总评分: 论坛币 + 10   查看全部评分

藤椅
llb_321 在职认证  发表于 2020-8-29 19:21:59
另外,data=read.table("data.txt",sep="\t",header=T,check.names=F)    这行代码没有必要,因为data本来就在的,你重新读入数据,反倒可能会导致数据类型的改变。

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-30 03:20