Off topic:gene names for ballgown
0
0
Entering edit mode
3.7 years ago

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 • 526 views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 2813 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