R 进行生存分析时出现log2问题
我计划使用TCGAbiolinks通过使用以下代码执行生存分析:R 进行生存分析时出现log2问题,r,bioconductor,R,Bioconductor,我计划使用TCGAbiolinks通过使用以下代码执行生存分析: library(TCGAbiolinks) library(SummarizedExperiment) query_RNAseq <- GDCquery(project = "TCGA-PAAD", data.category = "Gene expression", data.type = "Gene expression
library(TCGAbiolinks)
library(SummarizedExperiment)
query_RNAseq <- GDCquery(project = "TCGA-PAAD",
data.category = "Gene expression",
data.type = "Gene expression quantification",
platform = "Illumina HiSeq",
file.type = "results",
experimental.strategy = "RNA-Seq",
legacy = TRUE)
GDCdownload(query_RNAseq, method = "api",chunks.per.download = 10)
PAADRnaseqSE <- GDCprepare(query_RNAseq)
# Survival Analysis SA
clinical_patient_Cancer <- GDCquery_clinic("TCGA-PAAD","clinical")
dataPAADcomplete <- log2(PAADRnaseqSE)
库(TCGAbiolinks)
图书馆(总结实验)
查询\u RNAseq代码来自哪里?你期望最后一行做什么log2()
对数字有效,而paddrnase
是一个复杂的数据结构…只是猜测,但是log2(as.numeric(PAADRnaseqSE))
有效吗?@Ben Bolkeryes,谢谢,我使用SummaredExperiment包中的分析功能将PadRNaseqse
转换为矩阵后解决了这个问题:PadMatrix@Mist,它不起作用,但我已经解决了这个问题。谢谢你的提示
> dataPAADcomplete <- log2(PAADRnaseqSE)
Error in log(<S4 object of class "RangedSummarizedExperiment">, 2) :
unused argument (2)