This is a test version of Biostars. For the public version, visit https://www.biostars.org.
gene names for ballgown

Hi recently I performed an analysis on my data using stringtie and ballgown but i cant find my gene names, these are my commands:

    setwd("E:/share2/ball2")
    setRepositories
    library(ballgown)
    library(genefilter)
    library(dbplyr)
    library(magrittr)
    library(dplyr)
    library(readxl)
    remotes::install_github("alyssafrazee/RSkittleBrewer",force = TRUE)
pheno_data = read_excel("E:/my meta data.xlsx",col_names = T)
pheno_data<-as.data.frame(pheno_data)
bg_chrX = ballgown(dataDir = "ballgown", pData = pheno_data,samples = dir)
bg_chrX_filt = subset(bg_chrX,"rowVars(texpr(bg_chrX)) >1",genomesubset=TRUE)
results_transcripts <-stattest(bg_chrX_filt, feature='transcript', covariate='sex' 
          , getFC=TRUE, meas='FPKM')
results_transcripts <- data.frame(geneNames = ballgown::geneNames(bg_chrX_filt),
                                  geneIDs = ballgown::geneIDs(bg_chrX_filt),
                                  results_transcripts)

as you see in the result there is no even one gene name, look this:

    geneNames  geneIDs    feature  id        fc      pval      qval
9           .  MSTRG.2 transcript   9 0.8954503 0.5568966 0.9975262
40          .  MSTRG.2 transcript  40 1.0819969 0.7945403 0.9986724
50          . MSTRG.38 transcript  50 0.9637638 0.8832754 0.9986724
100         . MSTRG.50 transcript 100 1.1624027 0.4955879 0.9975262
101         . MSTRG.51 transcript 101 1.1214982 0.7252486 0.9985385
102         . MSTRG.52 transcript 102 0.7734894 0.2886172 0.9975262

is there any problem in commands?? how can I find my gene names directly??

rna-seq deg ballgown stringtie

Hello mo.ghasemi717!

Questions similar to yours can already be found at:

We have closed your question to allow us to keep similar content in the same thread.

If you disagree with this please tell us why in a reply below. We'll be happy to talk about it.

Cheers!

0 answers

No answers yet.

Log in to answer this question.