EnhancedVolcanop plot using GSE data
0
0
Entering edit mode
3 months ago
mnx0723 • 0

I am very lacking in skills. Please advise on the code. volcanoplot is not drawn properly.

test <- getGEO("GSE58144", GSEMatrix = TRUE)
test <- test[[1]]
exprs_data <- exprs(test)
exprs_data[exprs_data < 0] <- 0
exprs_data <- round(exprs_data)
sample_info <- pData(test)
sample_info <- dplyr::select(sample_info,source_name_ch1, characteristics_ch1)


colData <- data.frame(condition = test$`condition:ch1`,
stringsAsFactors = FALSE)


dds <- DESeqDataSetFromMatrix(countData = exprs_data, colData = colData, design = ~condition)
dds <-DESeq(dds)
res <- results(dds)

EnhancedVolcano(
 res,
 lab = rownames(res),
 x = 'log2FoldChange',
 y = 'padj',
 xlim = c(-2, 2),
 title = 'Volcano Plot',
 pCutoff = 0.05,
 FCcutoff = 0.5,
 pointSize = 2.5)

enter image description here

microarray GSE EnhancedVolcano • 231 views
ADD COMMENT
0
Entering edit mode

Hey, you should not perform DE analysis using DESeq2 for these data. They do not have raw count values.

ADD REPLY

Login before adding your answer.

Traffic: 2732 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6