Question: gene symbol in r
0
gravatar for miravet65
2.9 years ago by
miravet650
miravet650 wrote:

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>
R gene • 1.4k views
ADD COMMENTlink modified 17 months ago by zx87549.0k • written 2.9 years ago by miravet650

I think if you set gene symbols as data frames' rownames, It is not a problem.

ADD REPLYlink written 2.9 years ago by zjhzwang180

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?

ADD REPLYlink written 2.9 years ago by WouterDeCoster43k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1128 users visited in the last hour