remove NAs values from gset
1
0
Entering edit mode
@mortezahadizade-16792
Last seen 29 days ago
Iran

Hi everyone

Sometimes in microarray analysis, “ex” matrix is contained any NA values, I remove NAs with “na.omit” command, after removing NAs from “ex” matrix, do I need to remove NAs from “gset”???

I would be very pleased if you could guide me.

 

gset <- getGEO(dataset, GSEMatrix =TRUE, AnnotGPL=TRUE, destdir = "C:/R Database /Data/")

if (length(gset) > 1) idx <- grep(platform, attr(gset, "names")) else idx <- 1

gset <- gset[[idx]]

ex <- exprs(gset)

> max(ex)

[1] NA

> dim(ex)

[1] 18999    24

> ex<- na.omit(ex)

> max(ex)

[1] 7.683

> dim(ex)

[1] 16764    24

> gset<- na.omit(gset)?????

microarray R • 943 views
ADD COMMENT
1
Entering edit mode
@james-w-macdonald-5106
Last seen 7 hours ago
United States

I don't believe there is a na.omit function that will dispatch correctly on an ExpressionSet, so you will need to do things slightly differently. But do note that an ExpressionSet is intended to act just like a matrix, so subsetting with the [ operator works just like you would think.

naind <- apply(exprs(gset), 1, function(x) any(is.na(x)))

gset <- gset[!naind,]
ADD COMMENT
0
Entering edit mode

Or probably faster

naind <- is.na(rowSums(exprs(gset)))

gset <- gset[!naind,]
ADD REPLY
0
Entering edit mode

Thank you so much for your valuable guidance, your answer was reasonable and accurate.

ADD REPLY

Login before adding your answer.

Traffic: 492 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