R 做差异分析,离散值计算总是出错

R做差异表达分析问题

library(limma)
library(edgeR)
setwd(“G:Data treating/R/Difference analyse”)#将目标文件所在目录通过setwd命令设置为当前工作目录???
getwd()#查看目录是否设置成功
library(“readxl”)
rawdata1<-read_excel(“Practical theoretical difference.xlsx”,sheet = “Eukaryota”)
head(rawdata1)
rawdata<-as.matrix(t(rawdata1),colnames=NULL)#转置
write.table(rawdata, file = “Eukaryota.txt”,col.names = FALSE,row.names = TRUE, quote = F, sep="\t")
rawdata2<-read.delim(“Eukaryota.txt”,header = TRUE)
colnames(rawdata2)
y<-DGEList(counts=rawdata2[,2:477],genes=rawdata2[,1])
y<-DGEList(counts=y c o u n t s , g e n e s = y counts,genes=y counts,genes=ygenes)
y_norm<-calcNormFactors(y)#默认为TMM标准化
#lcpm <- cpm(y)#每百万次计算计数(CPM)或每千字节读取(RPKM)
y<-plotMDS(y_norm,col=c(rep(“black”,238), rep(“red”,238)), dim = c(1, 2))#绘制二维散点图
group<-factor(c(rep(‘practical’,238),rep(‘theoretical’,238)))
design <- model.matrix(~group)
y<-DGEList(counts=rawdata2[,2:477],genes=rawdata2[,1])
dge <- estimateDisp(y, design = design, robust = TRUE) #估算离散值

结果出现错误:Error in solve.default(designunique, t(beta)) : no right-hand side in 'b’
真心不明白是什么意思。还希望各位指教一下。

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值