【生信学习笔记】差异分析如何解决DESeq中

先来观察一下完整流程

library(DESeq2)
DESeq2group<-c(data[,3])
#DESeq2group<-c(rep("Control",15),rep("AD",10))

#或者
colData<-data.frame(row.names =names(testing),
                    condition=factor(DESeq2group,levels = c("Control","AD")))
dds<-DESeqDataSetFromMatrix(countData =sapply(test,as.integer),colData = colData,design=~condition)

dds<-dds[rowSums(counts(dds))>1,]

dds<-DESeq(dds)

res2<-results(dds)

res2<-as.data.frame(res2)

#test=as.data.frame(testing)

#先as.data.frame, 直接numeric会将数据降维

#dda<-DESeqDataSetFromMatrix(countData =sapply(test,as.integer),colData = colData,design=~condition)

#sapply和as.integer用来解决数据类型问题

如有错误请大佬多多指正

 

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值