Entering edit mode
3.4 years ago
eleramos10
•
0
Hello! I hope you can help me!
Im trying to work on a miRNA dataset with colorectal cancer samples. I found many useful options of datasets of this paper https://europepmc.org/article/pmc/5990389 (table 1). So for example, to load the data I write:
gset <- getGEO("GSE120300", GSEMatrix=TRUE, AnnotGPL=FALSE)
if(length(gset)>1) idx <- grep("GPL16016",attr(gset,"-
names")) else idx <- 1
gset <- gset[[idx]]
dim(gset)
However, when I try to remove the NA values from the exprs(gset) I get 0 features! I do not know what to do!
This is the code Im using to remove the NA values:
RMV = which(apply(exprs(gset), 1, function(x) any ( is.na(x))))
gset = gset[-RMV,]
dim(gset)
Can you help me? Thanks in advance!