This is a test version of Biostars. For the public version, visit https://www.biostars.org.
RNA sequencing analysis

Hi,

I'm trying to perform RNA sequencing analysis using data download from NCBI: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE184199

I'm following this workflow: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4934518/

These are the commands I've done so far:

FileURL <- paste(           "https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE184199",
            "format=file",
            "file=GSE184199_ALL_exp_FPKM.txt.gz",
            sep="&")
          download.file(FileURL, "GSE184199_ALL_exp_FPKM.txt.gz")

However, when I try to run the next command line I get an error message:

GenewiseCounts <- read.delim("GSE184199_ALL_exp_FPKM.txt.gz", row.names="EntrezGeneID")
Error in data[[rowvar]] : 
  attempt to select less than one element in get1index

Does anyone know how I might correct this?

Thanks!

sequencing rna
• 1,396 views • link
updated by Ram Ram • written by bioinformatics • bioinformatics

1 answer

Your tutorial does not say to use FPKM. EdgeR takes raw counts, not FPKM. Even if you could get it to work...don't.

If you must use this data, you can't use this tutorial. You need to look up a limma or limma-zoom tutorial.

• 0 views • link
written by swbarnes2 swbarnes2

Ok thanks for your advise. I will try to follow a limma-zoom tutorial.

• 0 views • link
written by bioinformatics • bioinformatics

Log in to answer this question.

More posts like this