Entering edit mode
                    8.6 years ago
        miravet65
        
    
        •
    
    0
    I am having a problem in getting gene symbol after analysis GEO data.
I can get logFC, P.value but in gene.symbol only show <NA>
gr<-as.factor(gr)
gset$description <- gr
design <- model.matrix(~ description + 0, gset)
colnames(design) <- levels(gr)
fit <- lmFit(gset, design)
cont.matrix <- makeContrasts("cp-np", levels=design)
fit2 <- contrasts.fit(fit, cont.matrix)
fit2 <- eBayes(fit2, 0.01)
tT <- topTable(fit2, adjust="fdr", sort.by="B", number=Inf)
tT <- subset(tT, select=c("adj.P.Val","P.Value","logFC","Gene.symbol"))
write.table(tT, file=stdout(), row.names=F, sep="\t")
Output:
 Gene.symbol
 <NA>
  <NA>
<NA>
 FBXL21
      <NA>
   <NA>
                    
                
                
I think if you set gene symbols as data frames' rownames, It is not a problem.
This question lacks the required information to answer this. You should at least add some commands you ran, how can we guess what's wrong without?